Search USGSSearch

SEARCH · Search USGS

Results for “Fast Times”

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 271 records · Page 15Linked to original sources

I. Thermal evolution of Ganymede and implications for surface features. II. Magnetohydrodynamic constraints on deep zonal flow in the giant planets. III. A fast finite-element algorithm for two-dimensional photoclinometry

The work is divided into three independent papers: PAPER I: Thermal evolution models are presented for Ganymede, assuming a mostly differentiated initial state of a water ocean overlying a rock layer. The only heat sources are assumed to be primordial heat (provided by accretion) and the long-lived radiogenic heat sources in the rock component. As Ganymede cools, the ocean thins, and two ice layers develop, one above composed of ice I, and the other below composed of high-pressure polymorphs of ice. Subsolidus convection proceeds separately in each ice layer, its transport of heat calculated using a simple parameterized convection scheme and the most recent data on ice rheology. The model requires that the average entropy of the deep ice layer exceed that of the ice I layer. If the residual ocean separating these layers becomes thin enough, then a Rayleigh-Taylor-like ("diapiric") instability may ensue, driven by the greater entropy of the deeper ice and merging the two ice mantles into a single convective layer. This instability is not predicted by linear analysis but occurs for plausible finite amplitude perturbations associated with large Rayleigh number convection. The resulting warm ice diapirs may lead to a dramatic "heat pulse" at the surface and to fracturing of the lithosphere, and may be directly or indirectly responsible for resurfacing and grooved terrain formation on Ganymede. The timing of this event depends rather sensitively on poorly known rheological parameters but could be consistent with chronologies deduced from estimated cratering rates. Irrespective of the occurrence or importance of the heat pulse, we find that lithospheric fracturing requires rapid stress loading (on a timescale ≾ 10 4 ) years). Such a timescale can be realized by warm ice diapirism, but not directly by gradual global expansion. In the absence of any quantitative and self-consistent model for the resurfacing of Ganymede by liquid water, we favor resurfacing by warm ice flows,which we demonstrate to be physically possible, a plausible consequence of our models, compatible with existing observations, and a hypothesis testable by Galileo. We discuss core formation as an alternative driver for resurfacing, and conclude that it is less attractive. We also consider anew the puzzle of why Callisto differs so greatly from Ganymede, offering several possible explanations. The models presented do not provide a compelling explanation for all aspects of Ganymedean geological evolution, since we have identified several potential problems, most notably the apparently extended period of grooved terrain formation (several hundred million years), which is difficult to reconcile with the heat pulse phenomenon. PAPER II: The observed zonal flows of the giant planets will, if they penetrate below the visible atmosphere, interact significantly with the planetary magnetic field outside the metalized core. The appropriate measure of this interaction is the Chandrasekhar number Q = ( H 2 )/(4πρνα 2 λ) (where H = radial component of the magnetic field, ν = eddy viscosity, λ = magnetic diffusivity, α -1 = lengthscale on which λ varies); at depths where Q ≳ 1 the velocity will be forced to oscillate on a small lengthscale or decay to zero. We estimate the conductivity due to semiconduction in H 2 (Jupiter, Saturn) and ionization in H 2 O (Uranus, Neptune) as a function of depth; the value λ ≃ 10 10 cm 2 s -1 needed for Q = 1 is readily obtained well outside the metallic core (where λ ≃ 10 2 cm 2 s -1 ). These assertions are quantified by a simple model of the equatorial zonal jet in which the flow is assumed uniform on cylinders concentric with the spin axis, and the viscous and magnetic torques on each cylinder are balanced. We solve this "Taylor constraint" simultaneously with the dynamo equation to obtain the velocity and magnetic field in the equatorial plane. With this model we reproduce the widely differing jet widths of Jupiter and Saturn (though not the flow at very high or low latitudes) using ν = 2500 cm 2 s -1 , consistent with the requirement that viscous dissipation not exceed the specific luminosity. A model Uranian jet consistent with the limited Voyager data can also be constructed, with appropriately smaller ν, but only if one assumes a two-layer interior. We tentatively predict a wide Neptunian jet. For Saturn (but not Jupiter or Uranus) the model has a large magnetic Reynolds number where Q = 1 and hence exhibits substantial axisymmetrization of the field in the equatorial plane . This effect may or may not persist at higher latitudes. The one-dimensional model presented is only a first step. Variation of the velocity and magnetic field parallel to the spin axis must be modeled in order to answer several important questions, including: 1) What is the behavior of flows at high latitudes, whose Taylor cylinders are interrupted by the region with Q ≳ 1? 2) To what extent is differential rotation in the envelope responsible for the spin-axisymmetry of Saturn's magnetic field? PAPER III: It is shown that the problem of two-dimensional photoclinometry (PC) -- the reconstruction of a surface z ( x , y ) from a brightness image B( x , y ) -- may be formulated in a natural way in terms of finite elements. The resulting system of equations is underdetermined as a consequence of the lack of boundary conditions for z , but a unique solution may be chosen by minimizing a function S expressing the "roughness" of the surface. An efficient PC algorithm based on this formulation is presented, requiring ~ 10.66 (four-byte) memory locations and ~10 4 floating multiplications/additions per pixel, and incorporating: 1) Minimization of the roughness by the penalty method, which yields the smallest set of equations. 2) Iterative solution of the nonlinear equations by Newton's method. 3) Solution of the linearized equations by an inner iterative cycle of successive over-relaxation, which takes advantage of the extreme sparseness of the system. 4) Multigridding, in which the solutions to the smaller problems obtained by reducing the resolution are used recursively to greatly speed convergence at the higher resolutions, and 5) A rapid noniterative initial estimate of z obtained by exploiting the special symmetry of the equations obtained in the first linearization. The algorithm is extensively demonstrated on 200 by 200 pixel synthetic "images" generated from digital topographic data for northern Utah over a range of phase angles. Rms error in the solution is ~ 22 m, out of ~ 660 m total relief. The error is dominated by "stripes" with the same azimuth as the light source, resulting from use of the roughness criterion in lieu of boundary conditions; the rms error along profiles parallel to the stripes is only ~ 2-8 m, depending on the phase angle. Satisfactory solutions are obtained even in the presence of quantization error, noise, and moderate blur in the image. Applications of the PC algorithm to both remote sensing and photomicrography are sketched; a photoclinometric map of a low-relief Precambrian era fossil is presented as an example of the latter. Prospects for dealing with photometrically inhomogeneous surfaces, and an extension of the method to the analysis of side-looking radar data ("radarclinometry") are also discussed.

Thesis

Composition and variation of noise recorded at the Yellowknife Seismic Array, 1991-2007

We analyze seismic noise recorded on the 18 short-period, vertical component seismometers of the Yellowknife Seismic Array (YKA). YKA has an aperture of 23 km and is sited on cratonic lithosphere in an area with low cultural noise. These properties make it ideal for studying natural seismic noise at periods of 1-3 s. We calculated frequency-wave number spectra in this band for over 6,000 time windows that were extracted once per day for 17 years (1991-2007). Slowness analysis reveals a rich variety of seismic phases originating from distinct source regions: R g waves from the Great Slave Lake; L g waves from the Atlantic, Pacific, and Arctic Oceans; and teleseismic P waves from the north Pacific and equatorial mid-Atlantic regions. The surface wave energy is generated along coastlines, while the body wave energy is generated at least in part in deep-water, pelagic regions. Surface waves tend to dominate at the longer periods and, just as in earthquake seismograms, L g is the most prominent arrival. Although the periods we study are slightly shorter than the classic double-frequency microseismic band of 4-10 s, the noise at YKA has clear seasonal behavior that is consistent with the ocean wave climate in the Northern Hemisphere. The temporal variation of most of the noise sources can be well fit using just two Fourier components: yearly and biyearly terms that combine to give a fast rise in microseismic power from mid-June through mid-October, followed by a gradual decline. The exception is the R g energy from the Great Slave Lake, which shows a sharp drop in noise power over a 2-week period in November as the lake freezes. The L g noise from the east has a small but statistically significant positive slope, perhaps implying increased ocean wave activity in the North Atlantic over the last 17 years. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Lower crustal deformation beneath the central Transverse Ranges, southern California: Results from the Los Angeles Region Seismic Experiment

We present a P wave velocity model derived from active source seismic data collected during the 1994 Los Angeles Region Seismic Experiment. Our model extends previously published upper crustal velocity models to mantle depths. Our model was developed by both ray tracing through a layered model and calculating travel times through a gridded model. It includes an 8-km-thick crustal root centered beneath the surface trace of the San Andreas fault, north of the highest topography in the San Gabriel Mountains. A simple mass balance calculation suggests that ∼36 km of north-south shortening across the San Andreas fault in the central Transverse Ranges could have formed this root. If north-south compression began when the “Big Bend” in the San Andreas fault formed at ∼5 Ma, 36 km of shortening implies a north-south contraction rate of ∼7.1 mm/yr across the central Transverse Ranges. If, instead, north-south compression began when the Transverse Ranges formed at 3.4–3.9 Ma, 36 km of shortening implies a contraction rate of 9.2–10.6 mm/yr. North of the San Andreas fault, the Mojave Desert crust has a low-velocity (6.3 km/s) mid and lower crust and a 28-km-deep Moho. South of the San Andreas fault, beneath the Los Angeles and San Gabriel Valley basins, there is a fast (6.6–6.8 km/s), thick (10–12 km) lower crust with a 27-km-deep Moho. Farther south still, the lower crust of the Continental Borderland is fast (6.6–6.8 km/s) and thin (5 km) with a shallow (22 km deep) Moho.

California

Drainage reversals in Mono Basin during the late pliocene and Pleistocene

Mono Basin, on the eastern flank of the central Sierra Nevada, is the highest of the large hydrographically closed basins in the Basin and Range province. We use geomorphic features, shoreline deposits, and basalt-filled paleochannels to reconstruct an early to middle Pleistocene record of shorelines and changing spillways of Lake Russell in Mono Basin. During this period of time, Lake Russell repeatedly attained altitudes between 2205 and 2280 m-levels far above the present surface of Mono Lake (~1950 m) and above its last overflow level (2188 m). The spill point of Lake Russell shifted through time owing to late Tertiary and Quaternary faulting and volcanism. During the early Pleistocene, the lake periodically discharged through the Mount Hicks spillway on the northeastern rim of Mono Basin and flowed northward into the Walker Lake drainage basin via the East Walker River. Paleochannels recording such discharge were incised prior to 1.6 Ma, possibly between 1.6 and 1.3 Ma, and again after 1.3 Ma (ages of basaltic flows that plugged the paleochannels). Faulting in the Adobe Hills on the southeastern margin of the basin eventually lowered the rim in this area to below the altitude of the Mount Hicks spillway. Twice after 0.76 Ma, and possibly as late as after 0.1 Ma, Lake Russell discharged southward through the Adobe Hills spillway into the Owens-Death Valley system of lakes. This study supports a pre-Pleistocene aquatic connection through Mono Basin between the hydrologically distinct Lahontan and Owens-Death Valley systems, as long postulated by biologists, and also confirms a probable link during the Pleistocene for species adapted to travel upstream in fast-flowing water.

Geological Society of America Bulletin

Temporal influences on selenium partitioning, trophic transfer, and exposure in a major U.S. river

Hydrologic and irrigation regimes mediate the timing of selenium (Se) mobilization to rivers, but the extent to which patterns in Se uptake and trophic transfer through recipient food webs reflect the temporal variation in Se delivery is unknown. We investigated Se mobilization, partitioning, and trophic transfer along approximately 60 river miles of the selenium-impaired segment of the Lower Gunnison River (Colorado, USA) during six sampling trips between June 2015 and October 2016. We found temporal patterns in Se partitioning and trophic transfer to be independent of those in dissolved Se concentrations and that the recipient food web sustained elevated Se concentrations from earlier periods of high Se mobilization. Using an ecosystem-scale Se accumulation model tailored to the Lower Gunnison River, we predicted that the endangered Razorback Sucker ( Xyrauchen texanus ) and Colorado Pikeminnow ( Ptychocheilus lucius ) achieve whole-body Se concentrations exceeding aquatic life protection criteria during periods of high runoff and irrigation activity (April–August) that coincide with susceptible phases of reproduction and early-life development. The results of this study challenge assumptions about Se trophodynamics in fast-flowing waters and introduce important considerations for the management of Se risks for biota in river ecosystems.

Colorado

Nitrogen excess in North American ecosystems: Predisposing factors, ecosystem responses, and management strategies

Most forests in North America remain nitrogen limited, although recent studies have identified forested areas that exhibit symptoms of N excess, analogous to overfertilization of arable land. Nitrogen excess in watersheds is detrimental because of disruptions in plant/soil nutrient relations, increased soil acidification and aluminum mobility, increased emissions of nitrogenous greenhouse gases from soil, reduced methane consumption in soil, decreased water quality, toxic effects on freshwater biota, and eutrophication of coastal marine waters. Elevated nitrate ( NO 3 − ) loss to groundwater or surface waters is the primary symptom of N excess. Additional symptoms include increasing N concentrations and higher N:nutrient ratios in foliage (i.e., N:Mg, N:P), foliar accumulation of amino acids or NO 3 − , and low soil C:N ratios. Recent nitrogen-fertilization studies in New England and Europe provide preliminary evidence that some forests receiving chronic N inputs may decline in productivity and experience greater mortality. Long-term fertilization at Mount Ascutney, Vermont, suggests that declining and slow N-cycling coniferous stands may be replaced by fast-growing and fast N-cycling deciduous forests. Symptoms of N saturation are particularly severe in high-elevation, nonaggrading spruce–fir ecosystems in the Appalachian Mountains and in eastern hardwood watersheds at the Fernow Experimental Forest near Parsons, West Virginia. In the Los Angeles Air Basin, mixed conifer forests and chaparral watersheds with high smog exposure are N saturated and exhibit the highest streamwater NO 3 − concentrations for wildlands in North America. High-elevation alpine watersheds in the Colorado Front Range and a deciduous forest in Ontario, Canada, are N saturated, although N deposition is moderate (∼8 kg·ha −1 ·yr −1 ). In contrast, the Harvard Forest hardwood stand in Massachusetts has absorbed >900 kg N/ha during 8 yr of N amendment studies without significant NO 3 − leaching, illustrating that ecosystems vary widely in the capacity to retain N inputs. Overly mature forests with high N deposition, high soil N stores, and low soil C:N ratios are prone to N saturation and NO 3 − leaching. Additional characteristics favoring low N retention capacity include a short growing season (reduced plant N demand) and reduced contact time between drainage water and soil (i.e., porous coarse-textured soils, exposed bedrock or talus). Temporal patterns of hydrologic fluxes interact with biotic uptake and internal cycling patterns in determining ecosystem N retention. Soils are the largest storage pool for N inputs, although vegetation uptake is also important. Recent studies indicate that nitrification may be widespread in undisturbed ecosystems, and that microbial assimilation of NO 3 − may be a significant N retention mechanism, contrary to previous assumptions. Further studies are needed to elucidate the sites, forms, and mechanisms of N retention and incorporation into soil organic matter, and to test potential management options for mitigating N losses from forests. Implementation of intensive management practices in N-saturated ecosystems may only be feasible in high-priority areas and on a limited scale. Reduction of N emissions would be a preferable solution, although major reductions in the near future are unlikely in many areas due to economic, energy-use, policy, and demographic considerations.

Ecological Applications

Prey patchiness and larval fish growth and survival: inferences from an individual-based model

We used an individual-based simulation model to evaluate how prey patchiness and fish swimming behavior affect larval fish survival and mortality source (predation or starvation). Simulations revealed that cohort survival increased linearly with greater average patch residence times and that patch residence times for individual fish with different fates (survived, starved, eaten) diverged substantially during the first few days of feeding. Further, by examining the interaction of patch spatial distribution (uniform random versus clumped) with three possible swimming behaviors, we found that swimming behavior, via its effect on prey encounter and feeding rates, affected both cohort survival rates and whether fish died from predation or starvation, but that the spatial distribution of patches (fine-scale (100s m)) only influenced whether fish died from predation or starvation. Within a particular patch spatial arrangement, however, patch intensity (division of food between patches and non-patches) had a major effect on survival. Except at high food levels, fish did not survive when there were no patches, suggesting that average, well-mixed prey densities will not support sufficiently rapid growth for survival. As patch intensity increased, survival increased to a maximum and then declined, with peak survival occurring at higher patch intensity as average food concentration declined. Finally, the degree of patchiness also determined the intensity of selection on growth rates. In patchier environments, there was stronger selection for fast growth rates leading to eight-fold differences in average cohort growth rates after only three days of growth. In general, survival was not directly related to fast cohort growth the best survival occurred with high average prey densities and weak patchiness. Prey patchiness, by influencing the average as well as the variance in individual growth rates, can have a substantial impact on survival rates of larval fish cohorts.

Ecological Modelling

Migration trends for king and common eiders and yellow-billed loons past Point Barrow in a rapidly changing environment

Most of the king ( Somateria spectabilis ) and common eiders ( S. mollissima v-nigra ) nesting in northern Alaska and northwestern Canada migrate past Point Barrow, Alaska, during the spring and fall migration. Yellow-billed loons ( Gavia adamsii ) also migrate past Point Barrow and are a species of international conservation concern. Spring migration counts of eiders have been conducted approximately every ten years at Point Barrow since 1976, and indicated that both eider species experienced population declines of approximately 50% between 1976 and 1996, and that the declines had stabilized by 2004. Population estimates derived from migration counts have not been previously estimated for yellow-billed loons. We conducted spring counts of eiders and loons in 2015 and 2016 to obtain population estimates to compare with those from 1994, 1995, 2003, and 2004 in order to evaluate long-term and current trends. We estimated (95% confidence intervals) that 796,419 (304,011) king and 96,775 (39,913) common eiders passed Point Barrow in 2015, and 322,381 (145,833) king and 130,390 (34,548) common eiders passed Point Barrow in 2016. Both and king and common eider population estimates increased from 1994 through 2016, however, the increase over time was not significant (F < 5.07, P > 0.087, df = 1). Our population estimates for king eiders were very different between the two years of this study, possibly due to a very short and intense migration peak in 2016, resulting in a population count that was biased low because sampling periods did not adequately capture the peak of migration. The numbers of common eiders were similar between the two years, as well as for the 12 years since the previous count. Photo analysis of flocks indicated that observer counts were on average 4% lower than photo counts (paired t-test; |t| = 3.26, df = 297, P < 0.001) for flocks less than 1400 individuals (observer count). Estimates of yellow-billed loon populations were very variable and are biased low as numbers of loons passing Pt. Barrow were still high when our counts ended in late May. It is important that counts continue to be conducted for these species of conservation and subsistence importance, but that techniques be refined to reduce bias and variability, and to find solutions to the increasing difficulty of conducting a count from the shore-fast ice in spring.

Alaska

Concentrations of environmental DNA (eDNA) reflect spawning salmon abundance at fine spatial and temporal scales

Developing fast, cost-effective assessments of wild animal abundance is an important goal for many researchers, and environmental DNA (eDNA) holds much promise for this purpose. However, the quantitative relationship between species abundance and the amount of DNA present in the environment is likely to vary substantially among taxa and with ecological context. Here, we report a strong quantitative relationship between eDNA concentration and the abundance of spawning sockeye salmon in a small stream in Alaska, USA, where we took temporally- and spatially-replicated samples during the spawning period. This high-resolution dataset suggests that (1) eDNA concentrations vary significantly day-to-day, and likely within hours, in the context of the dynamic biological event of a salmon spawning season; (2) eDNA, as detected by species-specific quantitative PCR probes, seems to be conserved over short distances (tens of meters) in running water, but degrade quickly over larger scales (ca. 1.5 km); and (3) factors other than the mere presence of live, individual fish — such as location within the stream, live/dead ratio, and water temperature — can affect the eDNA-biomass correlation in space or time. A multivariate model incorporating both biotic and abiotic variables accounted for over 75% of the eDNA variance observed, suggesting that where a system is well-characterized, it may be possible to predict species' abundance from eDNA surveys, although we underscore that species- and system-specific variables are likely to limit the generality of any given quantitative model. Nevertheless, these findings provide an important step toward quantitative applications of eDNA in conservation and management.

Alaska

Modeling the complex impacts of timber harvests to find optimal management regimes for Amazon tidal floodplain forests

At the Amazon estuary, the oldest logging frontier in the Amazon, no studies have comprehensively explored the potential long-term population and yield consequences of multiple timber harvests over time. Matrix population modeling is one way to simulate long-term impacts of tree harvests, but this approach has often ignored common impacts of tree harvests including incidental damage, changes in post-harvest demography, shifts in the distribution of merchantable trees, and shifts in stand composition. We designed a matrix-based forest management model that incorporates these harvest-related impacts so resulting simulations reflect forest stand dynamics under repeated timber harvests as well as the realities of local smallholder timber management systems. Using a wide range of values for management criteria (e.g., length of cutting cycle, minimum cut diameter), we projected the long-term population dynamics and yields of hundreds of timber management regimes in the Amazon estuary, where small-scale, unmechanized logging is an important economic activity. These results were then compared to find optimal stand-level and species-specific sustainable timber management (STM) regimes using a set of timber yield and population growth indicators. Prospects for STM in Amazonian tidal floodplain forests are better than for many other tropical forests. However, generally high stock recovery rates between harvests are due to the comparatively high projected mean annualized yields from fast-growing species that effectively counterbalance the projected yield declines from other species. For Amazonian tidal floodplain forests, national management guidelines provide neither the highest yields nor the highest sustained population growth for species under management. Our research shows that management guidelines specific to a region&rsquo;s ecological settings can be further refined to consider differences in species demographic responses to repeated harvests. In principle, such fine-tuned management guidelines could make management more attractive, thus bridging the currently prevalent gap between tropical timber management practice and regulation.

Amazon River, Mazagão watershed

Frac sand in the United States: a geological and industry overview

A new mineral rush is underway in the upper Midwest of the United States, especially in Wisconsin and Minnesota, for deposits of high-quality frac sand that the mining industry calls &ldquo;Northern White&rdquo; sand or &ldquo;Ottawa&rdquo; sand. Frac sand is a specialized type of sand that is added to fracking fluids that are injected into unconventional oil and gas wells during hydraulic fracturing (fracking or hydrofracking), a process that enhances petroleum extraction from tight (low permeability) reservoirs. Frac sand consists of natural sand grains with strict mineralogical and textural specifications that act as a proppant (keeping induced fractures open), extending the time of release and the flow rate of hydrocarbons from fractured rock surfaces in contact with the wellbore. The current sand mining surge has been driven by the boom in unconventional oil and gas production that has been largely spurred by advancements in technology promoting the expansion of hydraulic fracturing and horizontal drilling over the past decade. Because of its superior quality, the sand of the upper Midwest not only supports the majority of domestic oil and gas production, but it also supplies frac sand to some of Canada&rsquo;s western shale basins. The principal sources of &ldquo;Northern White&rdquo; or &ldquo;Ottawa&rdquo; sand in the upper Midwest are the Middle and Upper Ordovician St. Peter Sandstone and the Lower Ordovician and Upper Cambrian Jordan Formation, with the Upper Cambrian Wonewoc and Mount Simon Formations gaining in importance. Additional frac sand sources to the south include the Upper Cambrian Hickory Sandstone Member of the Riley Formation in Texas, which is referred to informally as &ldquo;Brown&rdquo; or &ldquo;Brady&rdquo; sand, and the Middle Ordovician Oil Creek Formation in Oklahoma. More than 40 United States industry operators are involved in the mining, processing, transportation, and distribution of frac sand to a robust market that is fast-growing in the United States and throughout the world. In addition to the abrupt rise in frac sand mining and distribution, a new industry has emerged from the production of alternative proppants, such as coated sand and synthetic beads. Alternative proppants, developed through new technologies, are competing with supplies of natural frac sand. In the long term, the vitality of both industries will be tied to the future of hydraulic fracturing of tight oil and gas reservoirs, which will be driven by the anticipated increases in global energy consumption.

Open-File Report

Atmospheric Deposition and Surface-Water Chemistry in Mount Rainier and North Cascades National Parks, U.S.A., Water Years 2000 and 2005-2006

High-elevation aquatic ecosystems in Mount Rainier and North Cascades National Parks are highly sensitive to atmospheric deposition of nitrogen and sulfur. Thin, rocky soils promote fast hydrologic flushing rates during snowmelt and rain events, limiting the ability of basins to neutralize acidity and assimilate nitrogen deposited from the atmosphere. Potential effects of nitrogen and sulfur deposition include episodic or chronic acidification of terrestrial and aquatic ecosystems. In addition, nitrogen deposition can cause eutrophication of water bodies and changes in species composition in lakes and streams. This report documents results of a study performed by the U.S. Geological Survey, in cooperation with the National Park Service, of the effects of atmospheric deposition of nitrogen and sulfur on surface-water chemistry in Mount Rainier and North Cascades National Parks. Inorganic nitrogen in wet deposition was highest in the vicinity of North Cascades National Park, perhaps due to emissions from human sources and activities in the Puget Sound area. Sulfur in wet deposition was highest near the Pacific coast, reflecting the influence of marine aerosols. Dry deposition generally accounted for less than 30 percent of wet plus dry inorganic nitrogen and sulfur deposition, but occult deposition (primarily fog) represents a potentially substantial unmeasured component of total deposition. Trend analyses indicate inorganic nitrogen in wet deposition was relatively stable during 1986-2005, but sulfur in wet deposition declined substantially during that time, particularly after 2001, when emissions controls were added to a large powerplant in western Washington. Surface-water sulfate concentrations at the study site nearest the powerplant showed a statistically significant decrease between 2000 and 2005-06, but there was no statistically significant change in alkalinity, indicating a delayed response in surface-water alkalinity. Seasonal patterns in surface-water chemistry and streamflow are strongly influenced by melting of seasonal snowpacks, which release large amounts of dilute, slightly acidic water to terrestrial and aquatic ecosystems during spring snowmelt. Concentrations of sulfate, alkalinity, and base cations in surface water declined rapidly during snowmelt, then gradually recovered during summer and fall. Preferential elution of acidic solutes from the snowpack at the beginning of snowmelt may cause episodic acidification in small alpine streams; evidence is provided by a stream sample collected at one of the sites during spring 2006 that was acidic (pH = 4.8, alkalinity = -18 microequivalents per liter) and had high concentrations of nitrate and sulfate and low concentrations of weathering products. Rain-on-snow events caused sharp declines in specific conductance, which was measured continuously using an in-stream sensor. A strong correlation was observed between measured specific conductance and measured alkalinity (r2 = 0.76), permitting estimation of alkalinity from specific-conductance data using a regression equation. Estimated alkalinity declined by an order of magnitude during the rain-on-snow events, in one case to 8 microequivalents per liter. Actual declines in alkalinity might be greater because the regression equation accounts only for dilution effects; at low concentrations, the relation between specific conductance and alkalinity is likely to be nonlinear and have a negative intercept (negative alkalinity). Thus, episodic acidification is possible during rain-on-snow events. The scale of episodic acidification is unknown, but if it occurs, it could have detrimental effects on aquatic life and amphibians. Historical lake-survey data indicate that most lakes are oligotrophic and have low nitrogen and phosphorus concentrations. Nitrogen limitation is more common in lakes in Mount Rainier National Park than in North Cascades National Park due to higher nitrate concentrations at North Cascades. T

Scientific Investigations Report

Molluscan isotope sclerochronology in marine palaeoclimatology: Taxa, technique and timespan issues

Study of the accretionary biomineralised hardparts of organisms (sclerochronology) can make a useful contribution to palaeoclimatology. Ontogenetic sequences of isotopic data (δ 18 O and Δ 47 values) from the shells of marine molluscs are a source of information on seasonal sea-surface temperatures that can be used for detailed and rigorous evaluation of the outputs of numerical climate models. In situations where there is significant seasonality, and where shell preservation is adequate, accurate information about winter and summer surface temperature can be obtained from shallow-water benthic forms (bivalves and gastropods), in particular the early ontogeny of fast-growing species. Accurate information about winter surface temperature can also be obtained from individuals that lived at mid-shelf depths (20–40 m), but summer seafloor values from these need upward adjustment to derive a plausible surface temperature if thermal stratification of the water column occurs in this season. Ontogenetic δ 18 O profiles from planktonic pteropod gastropods are a potential source of insight into seasonal surface temperatures in the ocean basins; these organisms merit investigation for provision of information to complement shelf data. Temperature profiles constructed from shell δ 18 O require an estimate of the δ 18 O value of ambient seawater, which can be derived by back-calculation from the Δ 47 -temperature supplied by the same shell material. Alternatively, through appropriate sampling and data processing, seasonal temperatures can be obtained directly from Δ 47 profiles. Climate parameters are defined in terms of the mean state over a period of 30 consecutive years, a statistic (e.g., for seasonal temperatures) which can be derived from the long isotopic temperature records obtainable from bivalve species that live for many tens or hundreds of years. Efforts should be made to acquire such records, especially averaged data from crossdated shells, to specify climate parameters for precise times in the past. Information for precise times would be of particular value for icehouse intervals like the late Cenozoic, characterised by high frequency (high amplitude) climate fluctuation. Short records from non-crossdated shells can nevertheless provide useful insights into climate, particularly if a large dataset is obtained, supplying a reliable picture of the mean state and range of variation in climate parameters over the interval represented by the shells.

Quaternary Science Reviews

Slow recovery of headwater-stream fishes following a catastrophic poisoning event

Accidental spills of chemicals and other pollutants can decimate populations of stream-dwelling species. Recovery from such accidents can be relatively fast and complete when the affected stream reaches can be recolonized from upstream and downstream sources. However, faunal recoveries from accidental spills that extirpate populations from entire headwater streams have not been extensively documented, and understanding resilience of headwater-stream biota is relevant for assessing threats to at-risk species. We assessed recovery of fish populations in a 5.7-km long headwater stream in the southeastern United States following a complete, or nearly complete, fish-kill caused by a chemical spill near the source of the stream. We sampled for fishes at five stream locations, two downstream and three upstream from a perched, culverted road-crossing located 2.4 km upstream from the stream mouth, over a period of 18.5 months following the poisoning event. We observed 11 fish species, representing <65% of the fish species expected based on occurrences in nearby tributary streams. In post-poisoning sampling, only three of these taxa were observed upstream of the culvert; all 11 species, including the federally threatened Cherokee Darter Etheostoma scotti, were found downstream of the culvert but were mostly represented by a few, large individuals. In contrast, dead individuals of at least eight taxa including the Cherokee Darter were observed upstream of the culvert at the time of the fish-kill. These observations provide evidence of slow recovery of a headwater fish fauna, and especially upstream of a barrier to fish movement, where the recolonization sources are primarily downstream. Additional case studies may reveal whether this result applies generally to headwater streams. Slow recovery could make species that primarily inhabit or maintain greatest abundances in headwaters, including multiple at-risk fishes, particularly vulnerable to the threat of accidental spills that result in local population extirpation.

Journal of Fish and Wildlife Management

NBII-SAIN Data Management Toolkit

The Strategic Plan for the U.S. Geological Survey Biological Informatics Program (2005-2009) recognizes the need for effective data management: Though the Federal government invests more than $600 million per year in biological data collection, it is difficult to address these issues because of limited accessibility and lack of standards for data and information...variable quality, sources, methods, and formats (for example observations in the field, museum specimens, and satellite images) present additional challenges. This is further complicated by the fast-moving target of emerging and changing technologies such as GPS and GIS. Even though these technologies offer new solutions, they also create new informatics challenges (Ruggiero and others, 2005). The USGS National Biological Information Infrastructure program, hereafter referred to as NBII, is charged with the mission to improve the way data and information are gathered, documented, stored, and accessed. The central objective of this project is a direct reflection of the purpose of NBII as described by John Mosesso, Program Manager of the U.S. Geological Survey-Biological Informatics Program-GAP Analysis: At the outset, the reason for bringing about NBII was that there were significant amounts of data and information scattered all over the U.S., not accessible, in incompatible formats, and that NBII was tasked with addressing this problem...NBII's focus is to pull data together that truly matters to someone or communities. Essentially, the core questions are: 1) what are the issues, 2) where is the data, and 3) how can we make it usable and accessible (John Mosesso, U.S. Geological Survey, oral commun., 2006). Redundancy in data collection can be a major issue when multiple stakeholders are involved with a common effort. In 2001 the U.S. General Accounting Office (USGAO) estimated that about 50 percent of the Federal government's geospatial data at the time was redundant. In addition, approximately 80 percent of the cost of a spatial information system is associated with spatial data collection and management (U.S. General Accounting Office, 2003). These figures indicate that the resources (time, personnel, money) of many agencies and organizations could be used more efficiently and effectively. Dedicated and conscientious data management coordination and documentation is critical for reducing such redundancy. Substantial cost savings and increased efficiency are direct results of a pro-active data management approach. In addition, details of projects as well as data and information are frequently lost as a result of real-world occurrences such as the passing of time, job turnover, and equipment changes and failure. A standardized, well documented database allows resource managers to identify issues, analyze options, and ultimately make better decisions in the context of adaptive management (National Land and Water Resources Audit and the Australia New Zealand Land Information Council on behalf of the Australian National Government, 2003). Many environmentally focused, scientific, or natural resource management organizations collect and create both spatial and non-spatial data in some form. Data management appropriate for those data will be contingent upon the project goal(s) and objectives and thus will vary on a case-by-case basis. This project and the resulting Data Management Toolkit, hereafter referred to as the Toolkit, is therefore not intended to be comprehensive in terms of addressing all of the data management needs of all projects that contain biological, geospatial, and other types of data. The Toolkit emphasizes the idea of connecting a project's data and the related management needs to the defined project goals and objectives from the outset. In that context, the Toolkit presents and describes the fundamental components of sound data and information management that are common to projects involving biological, geospatial, and other related data

Open-File Report

Characterization of the porosity distribution in the upper part of the karst Biscayne aquifer using common offset ground penetrating radar, Everglades National Park, Florida

The karst Biscayne aquifer is characterized by a heterogeneous spatial arrangement of porosity and hydraulic conductivity, making conceptualization difficult. The Biscayne aquifer is the primary source of drinking water for millions of people in south Florida; thus, information concerning the distribution of karst features that concentrate the groundwater flow and affect contaminant transport is critical. The principal purpose of the study was to investigate the ability of two-dimensional ground penetrating radar (GPR) to rapidly characterize porosity variability in the karst Biscayne aquifer in south Florida. An 800-m-long GPR transect of a previously investigated area at the Long Pine Key Nature Trail in Everglades National Park, collected in fast acquisition common offset mode, shows hundreds of diffraction hyperbolae. The distribution of diffraction hyperbolae was used to estimate electromagnetic (EM) wave velocity at each diffraction location and to assess both horizontal and vertical changes in velocity within the transect. A petrophysical model (complex refractive index model or CRIM) was used to estimate total bulk porosity. A set of common midpoint surveys at selected locations distributed along the common-offset transect also were collected for comparison with the common offsets and were used to constrain one-dimensional (1-D) distributions of porosity with depth. Porosity values for the saturated Miami Limestone ranged between 25% and 41% for common offset GPR surveys, and between 23% and 39% for common midpoint GPR surveys. Laboratory measurements of porosity in five whole-core samples from the saturated part of the aquifer in the study area ranged between 7.1% and 41.8%. GPR estimates of porosity were found to be valid only under saturated conditions; other limitations are related to the vertical resolution of the GPR signal and the volume of the material considered by the measurement methodology. Overall, good correspondence between GPR estimates and the direct porosity values from the whole-core samples confirms the ability of GPR common offset surveys to provide rapid characterization of porosity variability in the Biscayne aquifer. The common offset survey method has several advantages: (1) improved time efficiency in comparison to other GPR acquisition modes such as common midpoints; and (2) enhanced lateral continuity of porosity estimates, particularly when compared to porosity measurements on 1-D samples such as rock cores. The results also support the presence of areas of low EM wave velocity or high porosity under saturated conditions, causing velocity pull-down areas and apparent sag features in the reflection record. This study shows that GPR can be a useful tool for improving understanding of the petrophysical properties of highly heterogeneous systems such as karst aquifers, and thus may assist with the development of more accurate groundwater flow models, such as those used for restoration efforts in the Everglades.

Florida

Uranium-series disequilibrium in tuffs from Yucca Mountain, Nevada, as evidence of pore-fluid flow over the last million years

Samples of tuff from boreholes drilled into fault zones in the Exploratory Studies Facility (ESF) and relatively unfractured rock of the Cross Drift tunnels, at Yucca Mountain, Nevada, have been analysed by U-series methods. This work is part of a project to verify the finding of fast flow-paths through the tuff to ESF level, indicated by the presence of 'bomb' 36Cl in pore fluids. Secular radioactive equilibrium in the U decay series, (i.e. when the radioactivity ratios 234U/238U, 230Th/234U and 226Ra/230Th all equal 1.00) might be expected if the tuff samples have not experienced radionuclide loss due to rock-water interaction occurring within the last million years. However, most fractured and unfractured samples were found to have a small deficiency of 234U (weighted mean 234U/238U=0.95??0.01) and a small excess of 230Th (weighted mean 230Th/234U 1.10??0.02). The 226Ra/230Th ratios are close to secular equilibrium (weighted mean = 0.94??0.07). These data indicate that 234U has been removed from the rock samples in the last ???350 ka, probably by pore fluids. Within the precision of the measurement, it would appear that 226Ra has not been mobilized and removed from the tuff, although there may be some localised 226Ra redistribution as suggested by a few ratio values that are significantly different from 1.0. Because both fractured and unfractured tuffs show approximately the same deficiency of 234U, this indicates that pore fluids are moving equally through fractured and unfractured rock, More importantly, fractured rock appears not to be a dominant pathway for groundwater flow (otherwise the ratio would be more strongly affected and the Th and Ra isotopic ratios would likely also show disequilibrium). Application of a simple mass-balance model suggests that surface infiltration rate is over an order of magnitude greater than the rate indicated by other infiltration models and that residence time of pore fluids at ESF level is about 400 a. Processes of U sorption, precipitation and re-solution are believed to be occurring and would account for these anomalous results but have not been included in the model. Despite the difficulties, the U-series data suggest that fractured rock, specifically the Sundance and Drill Hole Wash faults, are not preferred flow paths for groundwater flowing through the Topopah Spring tuff and, by implication, rapid-flow, within 50 a, from the surface to the level of the ESF is improbable. ?? 2002 Elsevier Science Ltd. All rights reserved.

Applied Geochemistry

Quick and dirty (and accurate) 3-D paleoseismic trench models using coded scale bars

Structure‐from‐motion (SfM) modeling has dramatically increased the speed of generating geometrically accurate orthophoto mosaics of paleoseismic trenches, but some aspects of this technique remain time and labor intensive. Model accuracy relies on control points to establish scale, reduce distortion, and orient 3D models. Traditional SfM methods use total station or Global Navigation Satellite System (GNSS) surveys to constrain models, but collecting control points along a vertical trench wall is often inhibited by poor line of sight to the survey sensor or limited sky view and requires many hours in the field and office. We used physical scale bars printed with coded targets to constrain SfM models of a dusty, 46‐m‐long trench excavation across the Teton fault (Wyoming, U.S.A.). We present a workflow for generating quick and accurate 3D SfM models and orthophoto mosaics and compare the effectiveness of using scale bar, GNSS, and total‐station control in the models. Our results show that the scale bar model deviates from total station survey points by an average of 3.1 cm (maximum of 5.3 cm). In addition, the scale‐bar model only deviates an average of 1.7 cm (maximum 3.5 cm) when compared to the best model alternative, the SfM model controlled by the total station survey. Scale bars eliminate several hours needed to collect and incorporate control points from total station or GNSS surveys and significantly simplify the workflow, at the cost of slightly increased 3D model and orthophoto mosaic error. Our results further suggest that trench models can be constrained with at least four physical scale bars, but using five to six physical scale bars provides redundant control for field deployment and model optimization. The scale bar method for paleoseismic trenches proves to be portable and fast, minimizes the need for specialized survey equipment, and maintains model accuracy needed for mapping trench walls.

Seismological Research Letters