SEARCH · USGS Science

Results for “Atmospheric Research”

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

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49Linked to original sources

USGS1999

Tritium and 14C concentrations in unsaturated-zone gases at test hole UZB-2, Amargosa Desert Research Site, 1994-98: A section in U.S. Geological Survey Toxic Substances Hydrology Program: Proceedings of the technical meeting, Charleston, South Carolina, March 8-12, 1999: Volume 3 (Part C) (WRI 99-4018C)>

Tritium concentrations have been determined yearly since April 1994 from water-vapor samples collected at test hole UZB-2. The hole was drilled about 100 m (meters) south of the southwest corner of a commercial burial site for low-level radioactive wastes in September 1993. UZB-2 is equipped with ten 2.5-cm (centimeters) diameter air ports permanently installed in the unsaturated zone between the depths of 5.5 and 108.8 m below land surface. Depth to ground water is about 110 m. Additional sampling ports were driven by hand to depths of 0.5, 1.0 and 1.5 m in May 1997. Initial samples of water vapor collected in April 1994 showed elevated tritium concentrations of more than 100 TU (tritium units) from all 10 air ports, with a maximum concentration of 762±10 TU from an air port at a depth of 24.1 m. Subsequent tritium concentrations increased in all air ports, although tritium concentrations at depths of less than 34.1 m have remained relatively constant since July 1995. The largest observed increase in tritium has been at a depth of 47.9 m. There, tritium concentration has increased from 198±5 TU in April 1994 to 2,570±30 TU in June 1998. Large increases also have been measured in samples collected from air ports at depths of 106.4 and 108.8 m, just above the water table. During September and October 1998, carbon dioxide samples were collected from all ten air ports in UZB-2 and at a depth of 1.5 m, and analyzed for radioactive carbon-14 ( 14 C). 14 C concentrations are highest in air ports at depths less than 6 m where they exceed 2,000 pmc (percent modern carbon). Concentrations decrease rapidly in air ports at depth and are about 20 pmc below 94.2 m. However, at 47.9 meters, the 14 C concentration is 205±1 pmc, which is 2 to 4 times higher than concentrations in air ports immediately above and below. This depth corresponds to the largest tritium increase in UZB-2. Concentrations of both tritium and 14 C are greater than what could be expected from atmospheric fallout. The distribution of tritium and 14 C likely represent a complex pattern of lateral and vertical transport through the unsaturated zone from buried wastes to UZB-2.

Nevada↗

USGS2021

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): A new observational database to support explosive eruptive column model validation and development

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Journal of Volcanology and Geothermal Research↗

USGS2010

Testing mixing models of old and young groundwater in a tropical lowland rain forest with environmental tracers

We tested three models of mixing between old interbasin groundwater flow (IGF) and young, locally derived groundwater in a lowland rain forest in Costa Rica using a large suite of environmental tracers. We focus on the young fraction of water using the transient tracers CFC‐11, CFC‐12, CFC‐113, SF 6 , 3 H, and bomb 14 C. We measured 3 He, but 3 H/ 3 He dating is generally problematic due to the presence of mantle 3 He. Because of their unique concentration histories in the atmosphere, combinations of transient tracers are sensitive not only to subsurface travel times but also to mixing between waters having different travel times. Samples fall into three distinct categories: (1) young waters that plot along a piston flow line, (2) old samples that have near‐zero concentrations of the transient tracers, and (3) mixtures of 1 and 2. We have modeled the concentrations of the transient tracers using (1) a binary mixing model (BMM) of old and young water with the young fraction transported via piston flow, (2) an exponential mixing model (EMM) with a distribution of groundwater travel times characterized by a mean value, and (3) an exponential mixing model for the young fraction followed by binary mixing with an old fraction (EMM/BMM). In spite of the mathematical differences in the mixing models, they all lead to a similar conceptual model of young (0 to 10 year) groundwater that is locally derived mixing with old (>1000 years) groundwater that is recharged beyond the surface water boundary of the system.

Water Resources Research↗

USGS2024

Responses of marginal and intrinsic water-use efficiency to changing aridity using FLUXNET observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

Journal of Geophysical Research Biogeosciences↗

USGS2012

Analysis of time series of glacier speed: Columbia Glacier, Alaska

During the summers of 1984 and 1985, laser measurements were made of the distance from a reference location to markers on the surface of the lower reach of Columbia Glacier, Alaska. The data have numerous gaps, mostly because of inclement weather. The laser measurements were corrected for variations in atmospheric temperature and pressure and for the trajectory of the marker. The marker speed was calculated by using a combination of cubic splines and digital filters to arrive at a data set with 1-hour time intervals. The speed varies from 7 to 15 m/d and has three noteworthy components: (1) a low-frequency perturbation in speed with a time scale of days related to increased precipitation, (2) semidiurnal and diurnal variations related to sea tides, and (3) diurnal variations related to glacier surface melt. The low frequency and tidal period variations are separated by low-pass filtering and subtracting the filtered data from the original data. The variations from melt runoff and tides have the same frequency range and are examined using harmonic analysis and multiple regression with the various forcing functions. The high-frequency portion of the ice speed signal is dominated by tidal influences and to a lesser extent by meltwater influences. The low-frequency portion of the signal is dominated by the effects of precipitation.

Journal of Geophysical Research Solid Earth↗

USGS2000

Importance of the Mississippi River Basin for investigating agricultural–chemical contamination of the hydrologic cycle

This special issue is devoted to recent and ongoing research relating to the fate and transport of agricultural chemicals in the Mississippi River Basin by the US Geological Survey Toxic Substances Hydrology (Toxics) Program. The Mississippi River Basin drains approximately 3 200 000 km 2 representing 41% of the United States. This is the largest river in the United States and the third largest in the world. The Mississippi River discharges an average of 19 920 m 3 /s of water into the Gulf of Mexico. The river is an extensively used resource, supplying drinking water to 70 cities in the United States. The Mississippi River Basin has undergone dramatic land use and cultural changes over the last 150 years. Approximately 70 million people now live within the basin, representing approximately 27% of the nation's population. This basin has also become one of the most productive agricultural regions in the world in terms of both crops and livestock grown. Approximately 65% of the nation's harvested cropland is grown in this basin, with more than 100 000 metric tons (t) of pesticides and approximately 6 500 000 t of commercial nitrogen fertilizers applied to cropland within the basin annually. The drainage of more than 20 000 000 ha within the basin has been enhanced by means of tile lines and ditches to lower the water table to make the cropland more productive. While removing the water from the soil as intended, this practice also leads to more rapid transport of contaminants to the river, and ultimately the Gulf of Mexico. Furthermore, the extensive chemical use in the Mississippi River Basin has led to the transport of pesticides and nitrate into the region’s streams, aquifers, and atmosphere. An estimated 1 000 000 t of nitrate-N is transported from the Mississippi River Basin into the Gulf of Mexico annually. The peak annual load of herbicides to the Gulf of Mexico has been documented at 1920 t. The fundamental goal of the papers presented in this volume is to provide a scientific basis for decisions necessary to promote sound and efficient agricultural practices and protect the quality of the nation's water resources.

Science of the Total Environment↗

USGS2020

Surface methane concentrations along the mid-Atlantic bight driven by aerobic subsurface production rather than seafloor gas seeps

Relatively minor amounts of methane, a potent greenhouse gas, are currently emitted from the oceans to the atmosphere, but such methane emissions have been hypothesized to increase as oceans warm. Here, we investigate the source, distribution, and fate of methane released from the upper continental slope of the U.S. Mid-Atlantic Bight, where hundreds of gas seeps have been discovered between the shelf-break and ~1600 m water depth. Using physical, chemical, and isotopic analyses, we identify two main sources of methane in the water column: seafloor gas seeps and in situ aerobic methanogenesis which primarily occurs at 100 – 200 m depth in the water column. Stable isotopic analyses reveal that water samples collected at all depths were significantly impacted by aerobic methane oxidation, the dominant methane sink in this region, with more than 50% of the methane being oxidized, on average. Due to methane oxidation in the deeper water column, below 200 m depth, surface concentrations of methane are influenced more by methane sources found near the surface (0 – 10 m depth) and in the subsurface (10 - 200 m depth), rather than seafloor emissions at greater depths.

Atlantic margin↗

USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2013

Mercury bioaccumulation in fishes from subalpine lakes of the Wallowa-Whitman National Forest, northeastern Oregon and western Idaho

Mercury (Hg) is a globally distributed pollutant that poses considerable risks to human and wildlife health. Over the past 150 years since the advent of the industrial revolution, approximately 80 percent of global emissions have come from anthropogenic sources, largely fossil fuel combustion. As a result, atmospheric deposition of Hg has increased by up to 4-fold above pre-industrial times. Because of their isolation, remote high-elevation lakes represent unique environments for evaluating the bioaccumulation of atmospherically deposited Hg through freshwater food webs, as well as for evaluating the relative importance of Hg loading versus landscape influences on Hg bioaccumulation. The increase in Hg deposition to these systems over the past century, coupled with their limited exposure to direct anthropogenic disturbance make them useful indicators for estimating how changes in Hg emissions may propagate to changes in Hg bioaccumulation and ecological risk. In this study, we evaluated Hg concentrations in fishes of high-elevation, sub-alpine lakes in the Wallowa-Whitman National Forest in northeastern Oregon and western Idaho. Our goals were to (1) assess the magnitude of Hg contamination in small-catchment lakes to evaluate the risk of atmospheric Hg to human and wildlife health, (2) quantify the spatial variability in fish Hg concentrations, and (3) determine the ecological, limnological, and landscape factors that are best correlated with fish total mercury (THg) concentrations in these systems. Across the 28 study lakes, mean THg concentrations of resident salmonid fishes varied as much as 18-fold among lakes. Importantly, our top statistical model explained 87 percent of the variability in fish THg concentrations among lakes with four key landscape and limnological variables— catchment conifer density (basal area of conifers within a lake’s catchment), lake surface area, aqueous dissolved sulfate, and dissolved organic carbon. The basal area of conifers within a lake’s catchment was by far the most important variable explaining fish THg concentrations, with an increase in THg concentrations of more than 400 percent across the forest density spectrum. Across all study lakes, fish THg concentrations ranged from 0.004 to 0.438 milligrams per kilogram wet weight (mg/kg ww). Only a single individual fish sample exceeded the U.S. Environmental Protection Agency’s (USEPA) human health tissue residue criteria of 0.3 mg/kg ww. However, 54 percent of fish (N=177) exceeded the more stringent tissue residue criteria (0.04 mg/kg ww) adopted by the Oregon Department of Environmental Quality to better protect subsistence fishers. Additionally, 2 and 10 percent of fish exceeded levels associated with reduced common loon reproduction and behavior, respectively. Whereas 25 and 68 percent of fish sampled exceeded concentrations deemed protective of mink and kingfisher, respectively. These results suggest that THg concentrations may be present in these lakes at levels associated with ecological risk. It is important to note however, that accurate inference on potential impairment cannot be made within the context of this study design and further research is needed to better quantify these risks.

Oregon;Idaho↗

USGS2025

The ghost plume phenomenon and its impact on zenith-facing remote sensing measurements of volcanic SO2 emission rates

A large source of error in SO 2 emission rates derived from mobile Differential Optical Absorption Spectroscopy (DOAS) of volcanic gas plumes is the uncertainty in atmospheric light paths between the sun and the instrument, particularly under non-ideal atmospheric conditions, such as the presence of low clouds. DOAS instruments measure the SO 2 column density along the effective light path, so changes to that pathway directly affect the measured SO 2 signal. Due to complex radiative transfer mechanisms when a cloud is between the DOAS viewing position and a volcanic plume, measured plumes can appear spatially offset from their true location, a phenomenon informally referred to as “ghost plumes.” In addition to the appearance of ghost plumes, DOAS measurements recorded in non-ideal conditions have poorly characterized errors and are often discarded, limiting the data available to characterize volcanic degassing. In this study we simulate the radiative transfer associated with zenith-facing mobile DOAS traverses using the McArtim radiative transfer model for scenarios when there is a cloud layer between the instrument and the volcanic plume. In total, 217 permutations of atmospheric optical conditions are considered with varying cloud opacities (AOD = 0, 1, 2, 4, 8, 20), plume opacities (AOD = 0, 1, 2, 4, 8), solar zenith angles (SZA = 1°, 30°, 60°), and cloud thicknesses (200, 400, 800 m). We first develop objective criteria for selecting SO 2 baseline absorption levels and plume spatial extents. The simulated plume traverses are then integrated to obtain the SO 2 cross-sectional burdens which, after multiplication with the wind speed, yield SO 2 emission rates. We find large modification in the shape of the modeled cross-sectional burdens even under translucent (low AOD) cloud conditions in our modeled scenarios. Despite modification of the plume shape, the presence of a low cloud layer is typically not a large source of error in the SO 2 cross-sectional burden or emission rate obtained from zenith-facing DOAS traverses. We find that all measured cross-sectional burdens simulated using an aerosol-free plume in the above conditions and SZA ≤ 30° are within ±25% of the true value.

Journal of Volcanology and Goethermal Research↗

USGS2009

Hydrologic support of carbon dioxide flux revealed by whole-lake carbon budgets

Freshwater lakes are an important component of the global carbon cycle through both organic carbon (OC) sequestration and carbon dioxide (CO 2 ) emission. Most lakes have a net annual loss of CO 2 to the atmosphere and substantial current evidence suggests that biologic mineralization of allochthonous OC maintains this flux. Because net CO 2 flux to the atmosphere implies net mineralization of OC within the lake ecosystem, it is also commonly assumed that net annual CO 2 emission indicates negative net ecosystem production (NEP). We explored the relationship between atmospheric CO 2 emission and NEP in two lakes known to have contrasting hydrologie characteristics and net CO 2 emission. We calculated NEP for calendar year 2004 using whole-lake OC and inorganic carbon (IC) budgets, NEP oc and NEP IC , respectively, and compared the resulting values to measured annual CO 2 flux from the lakes. In both lakes, NEP Ic and NEP Ic were positive, indicating net autotrophy. Therefore CO 2 emission from these lakes was apparently not supported by mineralization of allochthonous organic material. In both lakes, hydrologie CO 2 inputs, as well as CO 2 evolved from netcalcite precipitation, could account for the net CO 2 emission. NEP calculated from diel CO 2 measurements was also affected by hydrologie inputs of CO 2 . These results indicate that CO 2 emission and positive NEP may coincide in lakes, especially in carbonate terrain, and that all potential geologic, biogeochemical, and hydrologie sources of CO 2 need to be accounted for when using CO 2 concentrations to infer lake NEP. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

USGS2011

Secular trends in the geologic record and the supercontinent cycle

Geologic secular trends are used to refine the timetable of supercontinent assembly, tenure, and breakup. The analysis rests on what is meant by the term supercontinent, which here is defined broadly as a grouping of formerly dispersed continents. To avoid the artificial pitfall of an all-or-nothing definition, quantitative measures of “supercontinentality” are presented: the number of continents, and the area of the largest continent, which both can be gleaned from global paleogeographic maps for the Phanerozoic. For the secular trends approach to be viable in the deep past when the very existence of supercontinents is debatable and reconstructions are fraught with problems, it must first be calibrated in the Phanerozoic against the well-constrained Pangea supercontinent cycle. The most informative geologic variables covering both the Phanerozoic and Precambrian are the abundances of passive margins and of detrital zircons. Both fluctuated with size of the largest continent during the Pangea supercontinent cycle and can be quantified back to the Neoarchean. The tenure of Pangea was a time represented in the rock record by few zircons and few passive margins. Thus, previously documented minima in the abundance of detrital zircons (and orogenic granites) during the Precambrian (Condie et al., 2009a, Gondwana Research 15, 228–242) now can be more confidently interpreted as marking the tenures of supercontinents. The occurrences of carbonatites, granulites, eclogites, and greenstone-belt deformation events also appear to bear the imprint of Precambrian supercontinent cyclicity. Together, these secular records are consistent with the following scenario. The Neoarchean continental assemblies of Superia and Sclavia broke up at ca. 2300 and ca. 2090 Ma, respectively. Some of their fragments collided to form Nuna by about 1750 Ma; Nuna then grew by lateral accretion of juvenile arcs during the Mesoproterozoic, and was involved in a series of collisions at ca. 1000 Ma to form Rodinia. Rodinia broke up in stages from ca. 1000 to ca. 520 Ma. Before Rodinia had completely come apart, some of its pieces had already been reassembled in a new configuration, Gondwana, which was completed by 530 Ma. Gondwana later collided with Laurentia, Baltica, and Siberia to form Pangea by about 300 Ma. Breakup of Pangea began at about 180 Ma (Early Jurassic) and continues today. In the suggested scenario, no supercontinent cycle in Earth history corresponded to the ideal, in which all the continents were gathered together, then broke apart, then reassembled in a new configuration. Nuna and Gondwana ended their tenures not by breakup but by collision and name change; Rodinia's assembly overlapped in time with its disassembly; and Pangea spalled Tethyan microcontinents throughout much of its tenure. Many other secular trends show a weak or uneven imprint of the supercontinent cycle, no imprint at all. Instead, these secular trends together reveal aspects of the shifting background against which the supercontinents came and went, making each cycle unique. Global heat production declined; plate tectonics sped up through the Proterozoic and slowed down through the Phanerozoic; the atmosphere and oceans became oxidized; life emerged as a major geochemical agent; some rock types went extinct or nearly so (BIF, massif-type anorthosite, komatiite); and other rock types came into existence or became common (blueschists, bioclastic limestone, coal).

Earth-Science Reviews↗

USGS2001

Density of basalt core from Hilo drill hole, Hawaii

Density measurements of 1600 samples of core from 889 to 3097 m depth below sea level in the Hawaii Scientific Drilling Program hole near Hilo, Hawaii show marked differences between the basaltic rock types and help define stratigraphy in the hole. Water-saturated densities of subaerial lava flows (occurring above 1079 m depth) have the broadest range because of the large density variation within a single lava flow. Water-saturated densities commonly range from 2.0 to 3.0 with an average of 2.55 ?? 0.24 g/cc. Dikes and sills range from 2.8 to 3.1 g/cc). Densities of hyaloclastite commonly range from 2.3 to 2.7, with an overall average of about 2.5 g/cc. The low-density of most hyaloclastite is due primarily to palagonitization of abundant glass and presence of secondary minerals in the interstices between fragments. Four principal zones of pillow lava, separated by hyaloclastite, occur in the drill core. The shallowest (1983-2136 m) is paradoxically the densest, averaging 3.01 ?? 0.10 g/cc. The second (2234-2470 m) is decidedly the lightest, averaging 2.67 ?? 0.13 g/cc. The third (2640-2790 m) and fourth (2918-bottom at 3097 m) are high, averaging 2.89 ?? 0.17 and 2.97 ?? 0.08 g/cc, respectively. The first pillow zone includes degassed pillows i.e. lava erupted on land that flowed into the sea. These pillows are poor in vesicles, because the subaerial, one-atmosphere vesicles were compressed when the flow descended to deeper water and higher pressure. The second (low-density, non-degassed) pillow zone is the most vesicle-rich, apparently because it was erupted subaqueously at a shallow depth. The higher densities of the third and fourth zones result from a low vesicularity of only a few percent and an olivine content averaging more than 5% for the third zone and about 10% for the fourth zone. The uppermost hyaloclastite extending about 400 m below the bottom of the subaerial basalt is poorly cemented and absorbs up to 6 wt% of water when immersed. Progressing downward the hyaloclastite absorbs less water and becomes better cemented. This change is apparently due to palagonitization of glass and addition of secondary minerals in the deeper older hyaloclastite, a process favored by the increase of temperature with depth. The cementation is largely complete at 1800 m depth where the temperature attains about 20??C. The zone of freshest, uncemented hyaloclastite represents the weakest rock in the drill hole and is a likely level for tectonic or landslide disruption. ?? 2001 Published by Elsevier Science B.V.

Journal of Volcanology and Geothermal Research↗

USGS2018

Regional variability of nitrate fluxes in the unsaturated zone and groundwater, Wisconsin, USA

Process-based modeling of regional NO3− fluxes to groundwater is critical for understanding and managing water quality, but the complexity of NO3− reactive transport processes make implementation a challenge. This study introduces a regional vertical flux method (VFM) for efficient estimation of reactive transport of NO3− in the vadose zone and groundwater. The regional VFM was applied to 443 well samples in central-eastern Wisconsin. Chemical measurements included O2, NO3−, N2 from denitrification, and atmospheric tracers of groundwater age including carbon-14, chlorofluorocarbons, tritium, and tritiogenic helium. VFM results were consistent with observed chemistry, and calibrated parameters were in-line with estimates from previous studies. Results indicated that (1) unsaturated zone travel times were a substantial portion of the transit time to wells and streams (2) since 1945 fractions of applied N leached to groundwater have increased for manure-N, possibly due to increased injection of liquid manure, and decreased for fertilizer-N, and (3) under current practices and conditions, approximately 60% of the shallow aquifer will eventually be affected by downward migration of NO3−, with denitrification protecting the remaining 40%. Recharge variability strongly affected the unsaturated zone lag times and the eventual depth of the NO3− front. Principal components regression demonstrated that VFM parameters and predictions were significantly correlated with hydrogeochemical landscape features. The diverse and sometimes conflicting aspects of N management (e.g. limiting N volatilization versus limiting N losses to groundwater) warrant continued development of large-scale holistic strategies to manage water quality and quantity.

Wisconsin↗

USGS2016

Risk, liability, and economic issues with long-term CO 2 storage—A review

Given a scarcity of commercial-scale carbon capture and storage (CCS) projects, there is a great deal of uncertainty in the risks, liability, and their cost implications for geologic storage of carbon dioxide (CO 2 ). The probabilities of leakage and the risk of induced seismicity could be remote, but the volume of geologic CO 2 storage (GCS) projected to be necessary to have a significant impact on increasing CO 2 concentrations in the atmosphere is far greater than the volumes of CO 2 injected thus far. National-level estimates of the technically accessible CO 2 storage resource (TASR) onshore in the United States are on the order of thousands of gigatons of CO 2 storage capacity, but such estimates generally assume away any pressure management issues. Pressure buildup in the storage reservoir is expected to be a primary source of risk associated with CO 2 storage, and only a fraction of the theoretical TASR could be available unless the storage operator extracts the saltwater brines or other formation fluids that are already present in the geologic pore space targeted for CO 2 storage. Institutions, legislation, and processes to manage the risk, liability, and economic issues with CO 2 storage in the United States are beginning to emerge, but will need to progress further in order to allow a commercial-scale CO 2 storage industry to develop in the country. The combination of economic tradeoffs, property rights definitions, liability issues, and risk considerations suggests that CO 2 storage offshore of the United States may be more feasible than onshore, especially during the current (early) stages of industry development.

Natural Resources Research↗

USGS2008

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS2009

Aeolian dunes as ground truth for atmospheric modeling on Mars

Martian aeolian dunes preserve a record of atmosphere/surface interaction on a variety of scales, serving as ground truth for both Global Climate Models (GCMs) and mesoscale climate models, such as the Mars Regional Atmospheric Modeling System (MRAMS). We hypothesize that the location of dune fields, expressed globally by geographic distribution and locally by dune centroid azimuth (DCA), may record the long-term integration of atmospheric activity across a broad area, preserving GCM-scale atmospheric trends. In contrast, individual dune morphology, as expressed in slipface orientation (SF), may be more sensitive to localized variations in circulation, preserving topographically controlled mesoscale trends. We test this hypothesis by comparing the geographic distribution, DCA, and SF of dunes with output from the Ames Mars GCM and, at a local study site, with output from MRAMS. When compared to the GCM: 1) dunes generally lie adjacent to areas with strongest winds, 2) DCA agrees fairly well with GCM modeled wind directions in smooth-floored craters, and 3) SF does not agree well with GCM modeled wind directions. When compared to MRAMS modeled winds at our study site: 1) DCA generally coincides with the part of the crater where modeled mean winds are weak, and 2) SFs are consistent with some weak, topographically influenced modeled winds. We conclude that: 1) geographic distribution may be valuable as ground truth for GCMs, 2) DCA may be useful as ground truth for both GCM and mesoscale models, and 3) SF may be useful as ground truth for mesoscale models. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS2022

Ready for real time: Performance of Global Navigation Satellite System in 2019 Mw 7.1 Ridgecrest, California, rapid response products

Global Navigation Satellite Systems (GNSSs) have undergone notable advancement in the last few decades, leading to the availability of a dataset with capabilities well beyond its original intended purpose. The proliferation of high‐rate (1 Hz or greater) GNSS receivers in areas of seismological interest now allows for routine consideration of dynamic earthquake ground motions, with centimeter‐level displacement accuracy via precise point positioning methods. Real‐time (RT) GNSS observations, from stations that are both telemetered and processed to displacement with minimal latency, have lower accuracy compared to post‐processed (PP) GNSS displacements due to imprecise knowledge of atmospheric conditions, satellite clocks, and satellite orbits in RT. Whether the quality of RT high‐rate GNSS is sufficient for use in rapid response products remains to be thoroughly examined. Here, we highlight RT GNSS displacement time series processed during the 2019 M w "> M w Mw 7.1 Ridgecrest, California, earthquake in the context of common rapid‐response products, magnitude estimation, and kinematic fault‐slip models. We discuss how these data can be used to supplement RT seismic data for rapid characterization of significant earthquakes. We find that kinematic fault‐slip models using RT GNSS data retain the general spatiotemporal characteristics of those with PP data, with subtle differences in size and amplitude of modeled slip asperities. We demonstrate the effect of these rapid seismic source models using RT GNSS data on the U.S. Geological Survey product ShakeMap—a downstream ground‐motion prediction algorithm informed by the rupture dimensions estimated in the slip model. Discrepancies in the ShakeMap estimate are minor, within ±12% change, with the most severe variation at the fault edges. Our analysis suggests that, when used in conjunction with available seismic data sources, RT GNSS is sufficient and valuable for rapid earthquake characterization.

California↗