Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 631 records · Page 35Linked to original sources

Probability distributions of bed load particle velocities, accelerations, hop distances, and travel times informed by Jaynes's principle of maximum entropy

We describe the most likely forms of the probability distributions of bed load particle velocities, accelerations, hop distances, and travel times, in a manner that formally appeals to inferential statistics while honoring mechanical and kinematic constraints imposed by equilibrium transport conditions. The analysis is based on E. Jaynes's elaboration of the implications of the similarity between the Gibbs entropy in statistical mechanics and the Shannon entropy in information theory. By maximizing the information entropy of a distribution subject to known constraints on its moments, our choice of the form of the distribution is unbiased. The analysis suggests that particle velocities and travel times are exponentially distributed and that particle accelerations follow a Laplace distribution with zero mean. Particle hop distances, viewed alone, ought to be distributed exponentially. However, the covariance between hop distances and travel times precludes this result. Instead, the covariance structure suggests that hop distances follow a Weibull distribution. These distributions are consistent with high-resolution measurements obtained from high-speed imaging of bed load particle motions. The analysis brings us closer to choosing distributions based on our mechanical insight.

Journal of Geophysical Research F: Earth Surface↗

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface↗

Constraining the relative importance of raindrop- and flow-driven sediment transport mechanisms in postwildfire environments and implications for recovery time scales

Mountain watersheds recently burned by wildfire often experience greater amounts of runoff and increased rates of sediment transport relative to similar unburned areas. Given the sedimentation and debris flow threats caused by increases in erosion, more work is needed to better understand the physical mechanisms responsible for the observed increase in sediment transport in burned environments and the time scale over which a heightened geomorphic response can be expected. In this study, we quantified the relative importance of different hillslope erosion mechanisms during two postwildfire rainstorms at a drainage basin in Southern California by combining terrestrial laser scanner-derived maps of topographic change, field measurements, and numerical modeling of overland flow and sediment transport. Numerous debris flows were initiated by runoff at our study area during a long-duration storm of relatively modest intensity. Despite the presence of a well-developed rill network, numerical model results suggest that the majority of eroded hillslope sediment during this long-duration rainstorm was transported by raindrop-induced sediment transport processes, highlighting the importance of raindrop-driven processes in supplying channels with potential debris flow material. We also used the numerical model to explore relationships between postwildfire storm characteristics, vegetation cover, soil infiltration capacity, and the total volume of eroded sediment from a synthetic hillslope for different end-member erosion regimes. This study adds to our understanding of sediment transport in steep, postwildfire landscapes and shows how data from field monitoring can be combined with numerical modeling of sediment transport to isolate the processes leading to increased erosion in burned areas.

California↗

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

Lateral flow separation occurs in rivers where banks exhibit strong curvature. In canyon-boundrivers, lateral recirculation zones are the principal storage of fine-sediment deposits. A parallelized,three-dimensional, turbulence-resolving model was developed to study the flow structures along lateralseparation zones located in two pools along the Colorado River in Marble Canyon. The model employs thedetached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacingin the interior of the flow. The DES-3D model is validated using Acoustic Doppler Current Profiler flowmeasurements taken during the 2008 controlled flood release from Glen Canyon Dam. A point-to-pointvalidation using a number of skill metrics, often employed in hydrological research, is proposed here forfluvial modeling. The validation results show predictive capabilities of the DES model. The model reproducesthe pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position ofthe primary and secondary eddy cells, and return current. The lateral recirculation zone is open, havingcontinuous import of fluid upstream of the point of reattachment and export by the recirculation returncurrent downstream of the point of separation. Differences in magnitude and direction of near-bed andnear-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between therecirculation return current and the main flow is dynamic, with large temporal changes in flow direction andmagnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shearlayer becoming three-dimensional without preferred orientation downstream.

Arizona↗

Debris flow runup on vertical barriers and adverse slopes

Runup of debris flows against obstacles in their paths is a complex process that involves profound flow deceleration and redirection. We investigate the dynamics and predictability of runup by comparing results from large-scale laboratory experiments, four simple analytical models, and a depth-integrated numerical model (D-Claw). The experiments and numerical simulations reveal the important influence of unsteady, multidimensional flow on runup, and the analytical models highlight key aspects of the underlying physics. Runup against a vertical barrier normal to the flow path is dominated by rapid development of a shock, or jump in flow height, associated with abrupt deceleration of the flow front. By contrast, runup on sloping obstacles is initially dominated by a smooth flux of mass and momentum from the flow body to the flow front, which precedes shock development and commonly increases the runup height. D-Claw simulations that account for the emergence of shocks show that predicted runup heights vary systematically with the adverse slope angle and also with the Froude number and degree of liquefaction (or effective basal friction) of incoming flows. They additionally clarify the strengths and limitations of simplified analytical models. Numerical simulations based on a priori knowledge of the evolving dynamics of incoming flows yield quite accurate runup predictions. Less predictive accuracy is attained in ab initio simulations that compute runup based solely on knowledge of static debris properties in a distant debris flow source area. Nevertheless, the paucity of inputs required in ab initio simulations enhances their prospective value in runup forecasting.

Journal of Geophysical Research F: Earth Surface↗

Numerical simulation of large-scale bed load particle tracer advection-dispersion in rivers with free bars

Asymptotic characteristics of the transport of bed load tracer particles in rivers have been described by advection-dispersion equations. Here we perform numerical simulations designed to study the role of free bars, and more specifically single-row alternate bars, on streamwise tracer particle dispersion. In treating the conservation of tracer particle mass, we use two alternative formulations for the Exner equation of sediment mass conservation: the flux-based formulation, in which bed elevation varies with the divergence of the bed load transport rate, and the entrainment-based formulation, in which bed elevation changes with the net deposition rate. Under the condition of no net bed aggradation/degradation, a 1-D flux-based deterministic model that does not describe free bars yields no streamwise dispersion. The entrainment-based 1-D formulation, on the other hand, models stochasticity via the probability density function (PDF) of particle step length, and as a result does show tracer dispersion. When the formulation is generalized to 2-D to include free alternate bars, however, both models yield almost identical asymptotic advection-dispersion characteristics, in which streamwise dispersion is dominated by randomness inherent in free bar morphodynamics. This randomness can result in a heavy-tailed PDF of waiting time. In addition, migrating bars may constrain the travel distance through temporary burial, causing a thin-tailed PDF of travel distance. The superdiffusive character of streamwise particle dispersion predicted by the model is attributable to the interaction of these two effects.

Journal of Geophysical Research F: Earth Surface↗

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗

A model integrating longshore and cross-shore processes for predicting long-term shoreline response to climate change

We present a shoreline change model for coastal hazard assessment and management planning. The model, CoSMoS-COAST (Coastal One-line Assimilated Simulation Tool), is a transect-based, one-line model that predicts short-term and long-term shoreline response to climate change in the 21 st century. The proposed model represents a novel, modular synthesis of process-based models of coastline evolution due to longshore and cross-shore transport by waves and sea-level rise. Additionally, the model uses an extended Kalman filter for data assimilation of historical shoreline positions to improve estimates of model parameters and thereby improve confidence in long-term predictions. We apply CoSMoS-COAST to simulate sandy shoreline evolution along 500 km of coastline in Southern California, which hosts complex mixtures of beach settings variably backed by dunes, bluffs, cliffs, estuaries, river mouths, and urban infrastructure, providing applicability of the model to virtually any coastal setting. Aided by data assimilation, the model is able to reproduce the observed signal of seasonal shoreline change for the hindcast period of 1995-2010, showing excellent agreement between modeled and observed beach states. The skill of the model during the hindcast period improves confidence in the model's predictive capability when applied to the forecast period (2010-2100) driven by GCM-projected wave and sea-level conditions. Predictions of shoreline change with limited human intervention indicate that 31% to 67% of Southern California beaches may become completely eroded by 2100 under sea-level rise scenarios of 0.93 to 2.0 m.

California↗

Links between atmosphere, ocean, and cryosphere from two decades of microseism observations on the Antarctic Peninsula

The lack of landmasses, climatological low pressure, and strong circumpolar westerly winds between the latitudes of 50°S to 65°S produce exceptional storm‐driven wave conditions in the Southern Ocean. This combination makes the Antarctic Peninsula one of Earth's most notable regions of high‐amplitude wave activity and thus, ocean‐swell‐driven microseism noise in both the primary (direct wave‐coastal region interactions) and secondary (direct ocean floor forcing due to interacting wave trains) period bands. Microseism observations are examined across 23 years (1993–2015) from Palmer Station (PMSA), on the west coast of the Antarctic Peninsula, and from East Falkland Island (EFI). These records provide a spatially integrative measure of both Southern Ocean wave amplitudes and the interactions between ocean waves and the solid Earth in the presence of sea ice, which can reduce wave coupling with the continental shelf. We utilize a spatiotemporal correlation‐based approach to illuminate how the distribution of sea ice influences seasonal microseism power. We characterize primary and secondary microseism power due to variations in sea ice and find that primary microseism energy is both more sensitive to sea ice and more capable of propagating across ocean basins than secondary microseism energy. During positive phases of the Southern Annular Mode, sea ice is reduced in the Bellingshausen Sea and overall storm activity in the Drake Passage increases, thus strongly increasing microseism power levels.

Antarctic Penninsula↗

Terrestrial ecosystem model performance in simulating productivity and its vulnerability to climate change in the northern permafrost region

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m −2 yr −1 ), most models produced higher NPP (309 ± 12 g C m −2 yr −1 ) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m −2 yr −1 ), which mainly resulted from differences in simulated maximum monthly GPP (GPP max ). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C ( V c max_25 ), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO 2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPP max as well as their sensitivity to climate change.

Journal of Geophysical Research: Biogeosciences↗

DOM composition and transformation in boreal forest soils: The effects of temperature and organic-horizon decomposition state

The boreal region stores large amounts of organic carbon (C) in organic-soil horizons, which are vulnerable to destabilization via warming and disturbance. Decomposition of soil organic matter (SOM) contributes to the production and turnover of dissolved organic matter (DOM). While temperature is a primary control on rates of SOM and DOM cycling, little is known about temperature effects on DOM composition in soil leachate. Here we conducted a 30 day incubation to examine the effects of temperature (20 versus 5°C) and SOM decomposition state (moss versus fibric versus amorphous horizons) on DOM composition in organic soils of interior Alaska. We characterized DOM using bulk dissolved organic C (DOC) concentration, chemical fractionation, optical properties, and ultrahigh-resolution mass spectrometry. We observed an increase in DOC concentration and DOM aromaticity in the 20°C treatment compared to the 5°C treatment. Leachate from fibric horizons had higher DOC concentration than shallow moss or deep amorphous horizons. We also observed chemical shifts in DOM leachate over time, including increases in hydrophobic organic acids, polyphenols, and condensed aromatics and decreases in low-molecular weight hydrophilic compounds and aliphatics. We compared ultrahigh-resolution mass spectrometry and optical data and observed strong correlations between polyphenols, condensed aromatics, SUVA 254 , and humic-like fluorescence intensities. These findings suggest that biolabile DOM was preferentially mineralized, and the magnitude of this transformation was determined by kinetics (i.e., temperature) and substrate quality (i.e., soil horizon). With future warming, our findings indicate that organic soils may release higher concentrations of aromatic DOM to aquatic ecosystems.

Journal of Geophysical Research: Biogeosciences↗

Nitrate reduction mechanisms and rates in an unconfined eogenetic karst aquifer in two sites with different redox potential

This study integrates push-pull tracer tests (PPTT) with microbial characterization of extracted water via quantitative polymerase chain reaction (qPCR) and reverse transcriptase qPCR (RT-qPCR) of selected functional N transformation genes to quantify nitrate reduction mechanisms and rates in sites with different redox potential in a karst aquifer. PPTT treatments with nitrate (AN) and nitrate-fumarate (ANC) were executed in two wells representing anoxic and oxic geochemical end-members. Oxic aquifer zero-order nitrate loss rates (mmol L −1 h −1 ) were similar for AN and ANC treatment, ranging from 0.03 ± 0.01 to 0.05 ± 0.01. Anoxic aquifer zero-order nitrate loss rates ranged from 0.03 ± 0.02 (AN) to 0.13 ± 0.02 (ANC). Microbial characterization indicates mechanisms influencing these rates were dissimilatory nitrate reduction to ammonium (DNRA) at the anoxic site with AN treatment, assimilatory reduction of nitrate to ammonium (ANRA) with ANC treatment in the water column at both sites, and additional documented nitrate reduction that occurred in unsampled biofilms. With carbon treatment, total numbers of microbes (16S rRNA genes) significantly increased (fourteenfold to thirtyfold), supporting stimulated growth with resulting ANRA. Decreased DNRA gene concentrations ( nrfA DNA) and increased DNRA activity ratio (nrfA-cDNA/DNA) supported the assertion that DNRA occurred in the anoxic zone with AN and ANC treatment. Furthermore, decreased DNRA gene copy numbers at the anoxic site with ANC treatment suggests that DNRA microbes in the anoxic site are chemolithoautotrophic. Increased RT-qPCR denitrification gene expression ( nirK and nirS ) was not observed in water samples, supporting that any observed NO 3 -N loss due to denitrification may be occurring in unsampled microbial biofilms.

Florida↗

Interactions among vegetation, climate, and herbivory control greenhouse gas fluxes in a subarctic coastal wetland

High-latitude ecosystems are experiencing the most rapid climate changes globally, and in many areas these changes are concurrent with shifts in patterns of herbivory. Individually, climate and herbivory are known to influence biosphere-atmosphere greenhouse gas (GHG) exchange; however, the interactive effects of climate and herbivory in driving GHG fluxes have been poorly quantified, especially in coastal systems that support large populations of migratory waterfowl. We investigated the magnitude and the climatic and physical controls of GHG exchange within the Yukon-Kuskokwim Delta in western Alaska across four distinct vegetation communities formed by herbivory and local microtopography. Net CO 2 flux was greatest in the ungrazed Carex meadow community (3.97 ± 0.58 [SE] µmol CO 2 m −2 s −1 ), but CH 4 flux was greatest in the grazed community (14.00 ± 6.56 nmol CH 4 m −2 s −1 ). The grazed community is also the only vegetation type where CH 4 was a larger contributor than CO 2 to overall GHG forcing. We found that vegetation community was an important predictor of CO 2 and CH 4 exchange, demonstrating that variation in regional gas exchange is best explained when the effect of grazing, determined by the difference between grazed and ungrazed communities, is included. Further, we identified an interaction between temperature and vegetation community, indicating that grazed regions could experience the greatest increases in CH 4 emissions with warming. These results suggest that future GHG fluxes could be influenced by both climate and by changes in herbivore population dynamics that expand or contract the vegetation community most responsive to future temperature change.

Alaska↗

Estimating belowground carbon stocks in isolated wetlands of the Northern Everglades Watershed, central Florida, using ground penetrating radar (GPR) and aerial imagery

Peat soils store a large fraction of the global soil carbon (C) pool and comprise 95% of wetland C stocks. While isolated freshwater wetlands in temperate and tropical biomes account for more than 20% of the global peatland C stock, most studies of wetland soil C have occurred in expansive peatlands in northern boreal and subarctic biomes. Furthermore, the contribution of small depressional wetlands in comparison to larger wetland systems in these environments is very uncertain. Given the fact that these wetlands are numerous and variable in terms of their internal geometry, innovative methods are needed for properly estimating belowground C stocks and their overall C contribution to the landscape. In this study, we use a combination of ground penetrating radar (GPR), aerial imagery, and direct measurements (coring) in conjunction with C core analysis to develop a relation between C stock and surface area, and estimate the contribution of subtropical depressional wetlands to the total C stock of pine flatwoods at the Disney Wilderness Preserve (DWP), Florida. Additionally, GPR surveys were able to image collapse structures underneath the peat basin of depressional wetlands, depicting lithological controls on the formation of depressional wetlands at the DWP. Results indicate the importance of depressional wetlands as critical contributors to the landscape C budget at the DWP and the potential of GPR-based approaches for (1) rapidly and noninvasively estimating the contribution of depressional wetlands to regional C stocks and (2) evaluating the formational processes of depressional wetlands.

Florida↗

Field-scale observations of a transient geobattery resulting from natural attenuation of a crude oil spill

We present evidence of a geobattery associated with microbial degradation of a mature crude oil spill. Self-potential measurements were collected using a vertical array of nonpolarizing electrodes, starting at the land surface and passing through the smear zone where seasonal water table fluctuations have resulted in the coating of hydrocarbons on the aquifer solids. These passive electrical potential measurements exhibit a dipolar pattern associated with a current source. The anodic and cathodic reactions of this natural battery occur below and above the smear zone, respectively. The smear zone is characterized by high magnetic susceptibility values associated with the precipitation of semiconductive magnetic iron phase minerals as a by-product of biodegradation, facilitating electron transfer between the anode and the cathode. This geobattery response appears to have a transient nature, changing on a monthly scale, probably resulting from chemical and physical changes in subsurface conditions such as water table fluctuations.

Journal of Geophysical Research G: Biogeosciences↗

Influence of atmospheric rivers on vegetation productivity and fire patterns in the southwestern U.S.

In the southwestern U.S., the meteorological phenomenon known as atmospheric rivers (ARs) has gained increasing attention due to its strong connections to floods, snowpacks, and water supplies in the West Coast states. Relatively less is known about the ecological implications of ARs, particularly in the interior Southwest, where AR storms are less common. To address this gap, we compared a chronology of AR landfalls on the west coast between 1989 and 2011 and between 25°N and 42.5°N to annual metrics of the normalized difference vegetation index (NDVI; an indicator of vegetation productivity) and daily resolution precipitation data to assess influences of AR-fed winter precipitation on vegetation productivity across the southwestern U.S. We mapped correlations between winter AR precipitation during landfalling ARs and (1) annual maximum NDVI and (2) area burned by large wildfires summarized by ecoregion during the same year as the landfalls and during the following year. Interannual variations of AR precipitation strongly influenced both NDVI and area burned by wildfire in some dryland ecoregions. The influence of ARs on dryland vegetation varied significantly depending on the latitude of landfall, with those ARs making landfall below 35°N latitude more strongly influencing these systems, and with effects observed as far as 1300 km from the landfall location. As climatologists' understanding of the synoptic patterns associated with the occurrence of ARs continues to evolve, an increased understanding of how AR landfalls, in aggregate, influence vegetation productivity and associated wildfire activity in dryland ecosystems may provide opportunities to better predict ecological responses to climate and climate change.

Arizona, California, Nevada, Utah↗

Spatial heterogeneity of within-stream methane concentrations

Streams, rivers, and other freshwater features may be significant sources of CH 4 to the atmosphere. However, high spatial and temporal variabilities hinder our ability to understand the underlying processes of CH 4 production and delivery to streams and also challenge the use of scaling approaches across large areas. We studied a stream having high geomorphic variability to assess the underlying scale of CH 4 spatial variability and to examine whether the physical structure of a stream can explain the variation in surface CH 4 . A combination of high-resolution CH 4 mapping, a survey of groundwater CH 4 concentrations, quantitative analysis of methanogen DNA, and sediment CH 4 production potentials illustrates the spatial and geomorphic controls on CH 4 emissions to the atmosphere. We observed significant spatial clustering with high CH 4 concentrations in organic-rich stream reaches and lake transitions. These sites were also enriched in the methane-producing mcrA gene and had highest CH 4 production rates in the laboratory. In contrast, mineral-rich reaches had significantly lower concentrations and had lesser abundances of mcrA . Strong relationships between CH 4 and the physical structure of this aquatic system, along with high spatial variability, suggest that future investigations will benefit from viewing streams as landscapes, as opposed to ecosystems simply embedded in larger terrestrial mosaics. In light of such high spatial variability, we recommend that future workers evaluate stream networks first by using similar spatial tools in order to build effective sampling programs.

Journal of Geophysical Research G: Biogeosciences↗

Deposition of mercury in forests across a montane elevation gradient: Elevational and seasonal patterns in methylmercury inputs and production

Global mercury contamination largely results from direct primary atmospheric and secondary legacy emissions, which can be deposited to ecosystems, converted to methylmercury, and bioaccumulated along food chains. We examined organic horizon soil samples collected across an elevational gradient on Whiteface Mountain in the Adirondack region of New York State, USA to determine spatial patterns in methylmercury concentrations across a forested montane landscape. We found that soil methylmercury concentrations were highest in the midelevation coniferous zone (0.39 ± 0.07 ng/g) compared to the higher elevation alpine zone (0.28 ± 0.04 ng/g) and particularly the lower elevation deciduous zone (0.17 ± 0.02 ng/g), while the percent of total mercury as methylmercury in soils decreased with elevation. We also found a seasonal pattern in soil methylmercury concentrations, with peak methylmercury values occurring in July. Given elevational patterns in temperature and bioavailable total mercury (derived from mineralization of soil organic matter), soil methylmercury concentrations appear to be driven by soil processing of ionic Hg, as opposed to atmospheric deposition of methylmercury. These methylmercury results are consistent with spatial patterns of mercury concentrations in songbird species observed from other studies, suggesting that future declines in mercury emissions could be important for reducing exposure of mercury to montane avian species.

New York↗