Search USGS⌕ Search

SEARCH · Search USGS

Results for “Complexity”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16Linked to original sources
USGS2003

Measurements of thermal updraft intensity over complex terrain using American white pelicans and a simple boundary-layer forecast model

An examination of boundary-layer meteorological and avian aerodynamic theories suggests that soaring birds can be used to measure the magnitude of vertical air motions within the boundary layer. These theories are applied to obtain mixed-layer normalized thermal updraft intensity over both flat and complex terrain from the climb rates of soaring American white pelicans and from diagnostic boundary-layer model-produced estimates of the boundary-layer depth z i and the convective velocity scale w * . Comparison of the flatland data with the profiles of normalized updraft velocity obtained from previous studies reveals that the pelican-derived measurements of thermal updraft intensity are in close agreement with those obtained using traditional research aircraft and large eddy simulation (LES) in the height range of 0.2 to 0.8 z i . Given the success of this method, the profiles of thermal vertical velocity over the flatland and the nearby mountains are compared. This comparison shows that these profiles are statistically indistinguishable over this height range, indicating that the profile for thermal updraft intensity varies little over this sample of complex terrain. These observations support the findings of a recent LES study that explored the turbulent structure of the boundary layer using a range of terrain specifications. For terrain similar in scale to that encountered in this study, results of the LES suggest that the terrain caused less than an 11% variation in the standard deviation of vertical velocity.

Boundary-Layer Meteorology↗
USGS1998

The Fremont complex: A behavioral perspective

The Fremont complex is composed of farmers and foragers who occupied the Colorado Plateau and Great Basin region of western North America from about 2100 to 500 years ago. These people included both immigrants and indigenes who shared some material culture and symbolic attributes, but also varied in ways not captured by definitions of the Fremont as a shared cultural tradition. The complex reflects a mosaic of behaviors including full-time farmers, full-time foragers, part-time farmer/foragers who seasonally switched modes of production, farmers who switched to full-time foraging, and foragers who switched to full-time farming. Farming defines the Fremont, but only in the sense that it altered the matrix in which both farmers and foragers lived, a matrix which provided a variety of behavioral options to people pursuing an array of adaptive strategies. The mix of symbiotic and competitive relationships among farmers and between farmers and foragers presents challenges to detection in the archaeological record. Greater clarity results from use of a behavioral model which recognizes differing contexts of selection favoring one adaptive strategy over another. The Fremont is a case where the transition from foraging to farming is followed by a millennium of adaptive diversity and terminates with the abandonment of farming. As such, it serves as a potential comparison to other cases in the world during the early phases of the food producing transition.

Journal of World Prehistory↗
USGS1999

Reactive solute transport in streams: A surface complexation approach for trace metal sorption

A model for trace metals that considers in-stream transport, metal oxide precipitation-dissolution, and pH-dependent sorption is presented. Linkage between a surface complexation submodel and the stream transport equations provides a framework for modeling sorption onto static and/or dynamic surfaces. A static surface (e.g., an iron- oxide-coated streambed) is defined as a surface with a temporally constant solid concentration. Limited contact between solutes in the water column and the static surface is considered using a pseudokinetic approach. A dynamic surface (e.g., freshly precipitated metal oxides) has a temporally variable solid concentration and is in equilibrium with the water column. Transport and deposition of solute mass sorbed to the dynamic surface is represented in the stream transport equations that include precipitate settling. The model is applied to a pH-modification experiment in an acid mine drainage stream. Dissolved copper concentrations were depressed for a 3 hour period in response to the experimentally elevated pH. After passage of the pH front, copper was desorbed, and dissolved concentrations returned to ambient levels. Copper sorption is modeled by considering sorption to aged hydrous ferric oxide (HFO) on the streambed (static surface) and freshly precipitated HFO in the water column (dynamic surface). Comparison of parameter estimates with reported values suggests that naturally formed iron oxides may be more effective in removing trace metals than synthetic oxides used in laboratory studies. The model's ability to simulate pH, metal oxide precipitation-dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between trace metal chemistry and hydrologic transport at the field scale.

Water Resources Research↗
USGS2005

Modeling and measuring the nocturnal drainage flow in a high-elevation, subalpine forest with complex terrain

The nocturnal drainage flow of air causes significant uncertainty in ecosystem CO 2 , H 2 O, and energy budgets determined with the eddy covariance measurement approach. In this study, we examined the magnitude, nature, and dynamics of the nocturnal drainage flow in a subalpine forest ecosystem with complex terrain. We used an experimental approach involving four towers, each with vertical profiling of wind speed to measure the magnitude of drainage flows and dynamics in their occurrence. We developed an analytical drainage flow model, constrained with measurements of canopy structure and SF 6 diffusion, to help us interpret the tower profile results. Model predictions were in good agreement with observed profiles of wind speed, leaf area density, and wind drag coefficient. Using theory, we showed that this one‐dimensional model is reduced to the widely used exponential wind profile model under conditions where vertical leaf area density and drag coefficient are uniformly distributed. We used the model for stability analysis, which predicted the presence of a very stable layer near the height of maximum leaf area density. This stable layer acts as a flow impediment, minimizing vertical dispersion between the subcanopy air space and the atmosphere above the canopy. The prediction is consistent with the results of SF 6 diffusion observations that showed minimal vertical dispersion of nighttime, subcanopy drainage flows. The stable within‐canopy air layer coincided with the height of maximum wake‐to‐shear production ratio. We concluded that nighttime drainage flows are restricted to a relatively shallow layer of air beneath the canopy, with little vertical mixing across a relatively long horizontal fetch. Insight into the horizontal and vertical structure of the drainage flow is crucial for understanding the magnitude and dynamics of the mean advective CO 2 flux that becomes significant during stable nighttime conditions and are typically missed during measurement of the turbulent CO 2 flux. The model and interpretation provided in this study should lead to research strategies for the measurement of these advective fluxes and their inclusion in the overall mass balance for CO 2 at this site with complex terrain.

Journal of Geophysical Research D: Atmospheres↗
USGS2007

Triggering mechanism and tsunamogenic potential of the Cape Fear Slide complex, U.S. Atlantic margin

Analysis of new multibeam bathymetry data and seismic Chirp data acquired over the Cape Fear Slide complex on the U.S. Atlantic margin suggests that at least 5 major submarine slides have likely occurred there within the past 30,000 years, indicating that repetitive, large-scale mass wasting and associated tsunamis may be more common in this area than previously believed. Gas hydrate deposits and associated free gas as well as salt tectonics have been implicated in previous studies as triggers for the major Cape Fear slide events. Analysis of the interaction of the gas hydrate phase boundary and the various generations of slides indicates that only the most landward slide likely intersected the phase boundary and inferred high gas pressures below it. For much of the region, we believe that displacement along a newly recognized normal fault led to upward migration of salt, oversteepening of slopes, and repeated slope failures. Using new constraints on slide morphology, we develop the first tsunami model for the Cape Fear Slide complex. Our results indicate that if the most seaward Cape Fear slide event occurred today, it could produce waves in excess of 2 m at the present-day 100 m bathymetric contour.

Cape Fear Slide↗
USGS2011

Hydrodynamics of a bathymetrically complex fringing coral reef embayment: Wave climate, in situ observations, and wave prediction

This paper examines the relationship between offshore wave climate and nearshore waves and currents at Hanalei Bay, Hawaii, an exposed bay fringed with coral reefs. Analysis of both offshore in situ data and numerical hindcasts identify the predominance of two wave conditions: a mode associated with local trade winds and an episodic pattern associated with distant source long-period swells. Analysis of 10 months of in situ data within the bay show that current velocities are up to an order of magnitude greater during long-period swell episodes than during trade wind conditions; overall circulation patterns are also fundamentally different. The current velocities are highly correlated with incident wave heights during the swell episodes, while they are not during the modal trade wind conditions. A phase-averaged wave model was implemented with the dual purpose of evaluating application to bathymetrically complex fringing reefs and to examine the propagation of waves into the nearshore in an effort to better explain the large difference in observed circulation during the two offshore wave conditions. The prediction quality of this model was poorer for the episodic condition than for the lower-energy mode, however, it illustrated how longer-period swells are preferentially refracted into the bay and make available far more nearshore wave energy to drive currents compared to waves during modal conditions. The highly episodic circulation, the nature of which is dependent on complex refraction patterns of episodic, long-period swell has implications for flushing and sediment dynamics for incised fringing reef-lined bays that characterize many high islands at low latitudes around the world.

Hawai'i↗
USGS2022

Observation-constrained multicycle dynamic models of the southern San Andreas and the northern San Jacinto Faults: Addressing complexity in paleoearthquake extent and recurrence with realistic 2D fault geometry

Understanding mechanical conditions that lead to complexity in earthquakes is important to seismic hazard analysis. In this study, we simulate physics-based multicycle dynamic models of the San Andreas fault (Carrizo through San Bernardino sections) and the San Jacinto fault (Claremont and Clark strands). We focus on a complex fault geometry based on the Southern California Earthquake Center Community Fault Model and its effect over multiple earthquake cycles. Using geodetically derived strain rates, we validate the models against geologic slip rates and recurrence intervals at various paleoseismic sites. We find that the interactions among fault geometry, dynamic rupture and interseismic stress accumulation produce stress heterogeneities, leading to rupture segmentation and variability in earthquake recurrence. Our models produce earthquakes with rupture extents similar to a recent comprehensive paleoseismic catalog. The “earthquake gates” of the Big Bend and the Cajon Pass occasionally impede dynamic ruptures. The angle of compression, which is the subtraction of the maximum shear strain rate direction from the local fault strike, can better determine the likelihood of the impedance of restraining bends to dynamic ruptures. Because the Big Bend has an angle of compression of ∼20°, ruptures that traverse the Big Bend, like the 1857 Fort Tejon earthquake, are more frequent than expected based on empirical relations which predict the ∼40° restraining bend to terminate most ruptures. Our models indicate that large ruptures tend to initiate north of the Big Bend and propagate southwards, similar to the 1857 earthquake, providing critical information for ground shaking assessment in the region.

California↗
USGS2010

Kinematic evidence for extensional unroofing of the Franciscan Complex along the Coast Range Fault, Northern Diablo Range, California

Franciscan metagraywacke immediately below the Del Puerto ophiolite, an outlier of the Coast Range ophiolite in the northern Diablo Range, was sheared during top‐to‐the‐east displacement on the Coast Range fault. This represents normal faulting and extensional offset. It was accompanied by attenuation of the Coast Range ophiolite and Great Valley sequence in the hanging wall along layer‐parallel normal faults that sole into the Coast Range fault. Extension occurred as the Franciscan Complex moved relatively west, out from under North American lithosphere and across the subducting ocean plate below. This effected a lengthening and thinning in the wedge of material above the down‐going plate, presumably in response to instability brought about by subduction shallowing (Krueger and Jones, 1989) and accretion of the Franciscan Central belt in the latest Cretaceous to early Paleocene. As a result, blueschist facies terranes of the uppermost part of the Franciscan Complex are now juxtaposed directly against hanging wall units that bear only low‐grade metamorphism.

Tectonics↗
USGS2012

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth↗
USGS2010

Age and origin of Cretaceous planktonic foraminifers from limestone of the Franciscan Complex near Laytonville, California

Detailed planktonic foraminiferal biostratigraphy from eight measured sections of Cretaceous limestone near Laytonville, California, indicates a composite sequence that extends in age from late Albian to early Turonian. The sequence contains seven biozones and two subzones based on the first and last appearance datums of planktonic foraminifers examined in thin section. Unequivocal biostratigraphic facing directions show four sections are right side up and four are reversed, and confirm the stratigraphic polarity employed in the paleomagnetic studies of Alvarez et al. (1980) and Tarduno et al. (1986). Temporal changes in the microfauna deviate from the global trend. Early Albian through early Cenomanian planktonic foraminiferal assemblages dominated by opportunistic, eurytopic species of Hedbergella and Globigerinelloides give way in the late Cenomanian to assemblages containing more abundant larger, heavily ornamented stenotopic species of Rotalipora and Praeglobotruncana typical of stratified Tethyan oceans. Radiolarian replacement chert increases slightly in the late Albian to early Cenomanian part of the Laytonville sequence and then decreases in the late Cenomanian and early Turonian. Benthic foraminifers show a similar decrease in abundance and diversity. We propose that this complex of biogenic and lithogenic patterns records transit via oceanic plate motion from a depositional site in the southern part of the paleoequatorial zone of high productivity, characterized by pronounced upwelling and habitat destabilization, perhaps augmented by topographic upwelling, to the central part of the equatorial zone dominated by biogenic calcite deposition. Our model of northward transit from below the equator supports the paleomagnetic determinations of Alvarez et al. (1980) and Tarduno et al. (1986). Arrival at the paleoequator coincided with the onset of stratification in the world ocean during the middle Cenomanian, indicated by the diversification of the more complex rotaliporids, and the subsequent intensification of oxygen-depleted intermediate waters in the latest Cenomanian to early Turonian. Two pulses in heterohelicid abundance signify the onset or intensification of an oxygen-minimum layer if Cretaceous heterohelicids occupied habitats analagous to Tertiary biserial heterohelicids. The first pulse beginning in the middle Cenomanian Rotalipora reicheli Zone may indicate a minor expansion of the oxygen-minimum zone or regional upwelling due to transit of the site beneath the equatorial divergence. The second pulse in the late Cenomanian Dicarinella algeriana Subzone signals the onset of upwelling of deeper oceanic water masses that characterized the succeeding Whiteinella archaeocretacea Zone. Organic-rich black shales typical of the Whiteinella archaeocretacea Zone are missing at Laytonville, although samples do contain the low-diversity, partially dissolved, planktonic assemblages that characterize this zone. The lack of black shales suggests that deposition occurred at depths greater than the oxygen-depleted intermediate water depths or alternately indicates paleoceanographic conditions unique to the Pacific.

California↗
USGS2010

Remagnetization and northward coastwise transport of Franciscan Complex rocks, northern California: A reinterpretation of the paleomagnetic data

Previous paleomagnetic studies of accreted oceanic rocks within the Franciscan Complex of northern California have concluded that these rocks originated far to the south of their present positions with respect to the North American continent. Based on positive “fold” tests, the characteristic remanent magnetizations were inferred to predate accretion-related deformation and metamorphism, and to have been acquired during or soon after deposition of these rocks. Thus, the paleomagnetic data were thought to provide direct information on ancient oceanic plate motions. However, the plate motions implied by some of these paleomagnetic data are problematic (e.g., exceptional plate velocities), and uniform-polarity magnetizations in almost all of these rocks indicate the possibility of remagnetization. Recent work on oceanic rocks in similar subduction complexes of Japan and Mexico have shown that they were most likely chemically remagnetized during accretion prior to disruption of the original stratigraphic sequences. Modern analogs indicate that the oceanic rocks in Mexico were probably remagnetized while still part of a shallow-dipping subducting slab (<10°) at the base of an accretionary prism. Assuming these rocks were near horizontal at the time of remagnetization, paleolatitudes at which these rocks were subducted and subsequent arc-parallel displacements along the western margin of North America can be inferred. In this paper, Franciscan rocks in northern California are reinterpreted as also having been remagnetized prior to accretion-related deformation. This scenario satisfies both geologic and paleomagnetic constraints for these rocks, and resolves conflicts between data indicating both remagnetization and tectonic displacement. Transport of the Laytonville Limestone from the southern hemisphere is not required. Paleolatitudes of subduction and remagnetization in the northern hemisphere (12° to 33°) appear to be inversely proportional to age of accretion (middle Cretaceous to Oligocene) for the Franciscan rocks. Subsequent northward diplacements (800 to 3700 km) and clockwise rotations (56° and 154°) of these rocks inferred from the paleomagnetic data are consistent with potential displacements along the western margin of North America during late Mesozoic and Cenozoic time calculated using examples of modern subduction zones and current plate reconstruction models.

California↗
USGS1983

Sm-Nd age and isotopic systematics of the bimodal suite, ancient gneiss complex, Swaziland

Studies of the development and stabilization of the Archaean crust often focus on the relative temporal relationships between the metamorphosed basaltic to ultramafic volcanic units (greenstone belts) and the sialic gneiss terrains that make up the oldest sections of the terrestrial crust. At the heart of this interest are the questions of the processes responsible for crust formation in the Archaean and whether or not the various units of an Archaean crustal section represent new additions to the crust from the mantle or are products of the reprocessing of even older crustal materials. One area where this controversy has been particularly pronounced is the Archaean crustal section of south-west Africa1-6. The oldest rocks in the Kaapvaal craton consist of the Onverwacht Group of mafic to ultramafic metavolcanics of the Barberton greenstone belt and a grey-gneiss complex termed the ancient gneiss complex (AGC) of Swaziland. We report here the results of a whole-rock Sm-Nd isotopic study of the AGC and the implications these data may have for crustal evolution in the Kaapvaal craton. ?? 1983 Nature Publishing Group.

Nature↗
USGS1985

Franciscan complex calera limestones: Accreted remnants of farallon plate oceanic plateaus

The Calera Limestone, part of the Franciscan Complex of northern California, may have formed in a palaeoenvironment similar to Hess and Shatsky Rises of the present north-west Pacific1. We report here new palaeomagnetic results, palaeontological data and recent plate-motion models that reinforce this assertion. The Calera Limestone may have formed on Farallon Plate plateaus, north of the Pacific-Farallon spreading centre as a counterpart to Hess or Shatsky Rises. In one model2, the plateaus were formed by hotspots close to the Farallon_Pacific ridge axis. On accretion to North America, plateau dissection in the late Cretaceous to Eocene (50-70 Myr) could explain the occurrence of large volumes of pillow basalt and exotic blocks of limestone in the Franciscan Complex. Partial subduction of the plateaus could have contributed to Laramide (70-40 Myr) compressional events3. ?? 1985 Nature Publishing Group.

Nature↗
USGS1987

Increased solubility of quartz in water due to complexing by organic compounds

Quartz is the most stable natural solid phase of silica. It weathers extremely slowly at the Earth's surface 1 , and often resists weathering even after all other silicate minerals have been degraded. However, there is ample evidence from both ancient and modern environments indicating enhanced dissolution and mobility of silica under conditions that cannot easily be explained by the inorganic controls of quartz solubility 2 . Increased solubility of quartz has been observed particularly in soils rich in organic material; however, no direct link between dissolved organic carbon and dissolved silica has been identified 3 . Here we present evidence for an increase in the solubility of quartz in a natural water brought about by dissolved organic compounds. These compounds were produced by the biodegradation of petroleum, and consist largely of a complex mixture of organic acids. We propose that silica is being complexed and mobilized by these organic acids in waters having close to neutral p H.

Nature↗
USGS2018

Strike-slip 23 January 2018 MW 7.9 Gulf of Alaska rare intraplate earthquake: Complex rupture of a fracture zone system

Large intraplate earthquakes in oceanic lithosphere are rare and usually related to regions of diffuse deformation within the oceanic plate. The 23 January 2018 M W 7.9 strike-slip Gulf of Alaska earthquake ruptured an oceanic fracture zone system offshore Kodiak Island. Bathymetric compilations show a muted topographic expression of the fracture zone due to the thick sediment that covers oceanic basement but the fracture zone system can be identified by offset N-S magnetic anomalies and E-W linear zones in the vertical gravity gradient. Back-projection from global seismic stations reveals that the initial rupture at first propagated from the epicenter to the north, likely rupturing along a weak zone parallel to the ocean crustal fabric. The rupture then changed direction to eastward directed with most energy emitted on Aka fracture zone resulting in an unusual multi-fault earthquake. Similarly, the aftershocks show complex behavior and are related to two different tectonic structures: (1) events along N-S trending oceanic fabric, which ruptured mainly strike-slip and additionally, in normal and oblique slip mechanisms and (2) strike-slip events along E-W oriented fracture zones. To explain the complex faulting behavior we adopt the classical stress and strain partitioning concept and propose a generalized model for large intra-oceanic strike-slip earthquakes of trench-oblique oriented fracture zones/ocean plate fabric near subduction zones. Taking the Kodiak asperity position of 1964 maximum afterslip and outer-rise Coulomb stress distribution into account, we propose that the unusual 2018 Gulf of Alaska moment release was stress transferred to the incoming oceanic plate from co- and post-processes of the nearby great 1964 M W 9.2 megathrust earthquake.

Alaska↗
USGS2019

Modeling spatially and temporally complex range dynamics when detection is imperfect

Species distributions are determined by the interaction of multiple biotic and abiotic factors, which produces complex spatial and temporal patterns of occurrence. As habitats and climate change due to anthropogenic activities, there is a need to develop species distribution models that can quantify these complex range dynamics. In this paper, we develop a dynamic occupancy model that uses a spatial generalized additive model to estimate non-linear spatial variation in occupancy not accounted for by environmental covariates. The model is flexible and can accommodate data from a range of sampling designs that provide information about both occupancy and detection probability. Output from the model can be used to create distribution maps and to estimate indices of temporal range dynamics. We demonstrate the utility of this approach by modeling long-term range dynamics of 10 eastern North American birds using data from the North American Breeding Bird Survey. We anticipate this framework will be particularly useful for modeling species’ distributions over large spatial scales and for quantifying range dynamics over long temporal scales.

Scientific Reports↗
USGS2001

Evaluation of flash-flood discharge forecasts in complex terrain using precipitation

Operational prediction of flash floods produced by thunderstorm (convective) precipitation in mountainous areas requires accurate estimates or predictions of the precipitation distribution in space and time. The details of the spatial distribution are especially critical in complex terrain because the watersheds are generally small in size, and small position errors in the forecast or observed placement of the precipitation can distribute the rain over the wrong watershed. In addition to the need for good precipitation estimates and predictions, accurate flood prediction requires a surface-hydrologic model that is capable of predicting stream or river discharge based on the precipitation-rate input data. Different techniques for the estimation and prediction of convective precipitation will be applied to the Buffalo Creek, Colorado flash flood of July 1996, where over 75 mm of rain from a thunderstorm fell on the watershed in less than 1 h. The hydrologic impact of the precipitation was exacerbated by the fact that a significant fraction of the watershed experienced a wildfire approximately two months prior to the rain event. Precipitation estimates from the National Weather Service's operational Weather Surveillance Radar-Doppler 1988 and the National Center for Atmospheric Research S-band, research, dual-polarization radar, colocated to the east of Denver, are compared. In addition, very short range forecasts from a convection-resolving dynamic model, which is initialized variationally using the radar reflectivity and Doppler winds, are compared with forecasts from an automated-algorithmic forecast system that also employs the radar data. The radar estimates of rain rate, and the two forecasting systems that employ the radar data, have degraded accuracy by virtue of the fact that they are applied in complex terrain. Nevertheless, the radar data and forecasts from the dynamic model and the automated algorithm could be operationally useful for input to surface-hydrologic models employed for flood warning. Precipitation data provided by these various techniques at short time scales and at fine spatial resolutions are employed as detailed input to a distributed-parameter hydrologic model for flash-flood prediction and analysis. With the radar-based precipitation estimates employed as input, the simulated flood discharge was similar to that observed. The dynamic-model precipitation forecast showed the most promise in providing a significant discharge-forecast lead time. The algorithmic system's precipitation forecast did not demonstrate as much skill, but the associated discharge forecast would still have been sufficient to have provided an alert of impending flood danger.

Journal of Hydrologic Engineering↗
USGS2009

Quantifying ecological, morphological, and genetic variation to delimit species in the coast horned lizard species complex (Phrynosoma)

Lineage separation and divergence form a temporally extended process whereby populations may diverge genetically, morphologically, or ecologically, and these contingent properties of species provide the operational criteria necessary for species delimitation.We inferred the historical process of lineage formation in the coast horned lizard (Phrynosoma coronatum) species complex by evaluating a diversity of operational species criteria, including divergence in mtDNA (98 specimens; 2,781 bp) and nuclear loci (RAG-1, 1,054 bp; BDNF 529 bp), ecological niches (11 bioclimatic variables; 285 unique localities), and cranial horn shapes (493 specimens; 16 landmarks). A phylogenetic analysis of mtDNA recovers 5 phylogeographic groups arranged latitudinally along the Baja California Peninsula and in California. The 2 southern phylogeographic groups exhibit concordance between genetic, morphological, and ecological divergence; however, differentiation is weak or absent at more recent levels defined by phylogeographic breaks in California. Interpreting these operational species criteria together suggests that there are 3 ecologically divergent and morphologically diagnosable species within the P. coronatum complex. Our 3-species taxonomic hypothesis invokes a deep coalescence event when fitting the mtDNA genealogy into the species tree, which is not unexpected for populations that have diverged recently. Although the hypothesis that the 3 phylogeographic groups distributed across California each represent distinctive species is not supported by all of the operational species criteria evaluated in this study, the conservation status of the imperiled populations represented by these genealogical units remains critical.

Proceedings of the National Academy of Sciences of↗