USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied Sciences”

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 865 records · Page 48Linked to original sources

Prediction, time variance, and classification of hydraulic response to recharge in two karst aquifers

Many karst aquifers are rapidly filled and depleted and therefore are likely to be susceptible to changes in short-term climate variability. Here we explore methods that could be applied to model site-specific hydraulic responses, with the intent of simulating these responses to different climate scenarios from high-resolution climate models. We compare hydraulic responses (spring flow, groundwater level, stream base flow, and cave drip) at several sites in two karst aquifers: the Edwards aquifer (Texas, USA) and the Madison aquifer (South Dakota, USA). A lumped-parameter model simulates nonlinear soil moisture changes for estimation of recharge, and a time-variant convolution model simulates the aquifer response to this recharge. Model fit to data is 2.4% better for calibration periods than for validation periods according to the Nash–Sutcliffe coefficient of efficiency, which ranges from 0.53 to 0.94 for validation periods. We use metrics that describe the shapes of the impulse-response functions (IRFs) obtained from convolution modeling to make comparisons in the distribution of response times among sites and between aquifers. Time-variant IRFs were applied to 62% of the sites. Principal component analysis (PCA) of metrics describing the shapes of the IRFs indicates three principal components that together account for 84% of the variability in IRF shape: the first is related to IRF skewness and temporal spread and accounts for 51% of the variability; the second and third largely are related to time-variant properties and together account for 33% of the variability. Sites with IRFs that dominantly comprise exponential curves are separated geographically from those dominantly comprising lognormal curves in both aquifers as a result of spatial heterogeneity. The use of multiple IRF metrics in PCA is a novel method to characterize, compare, and classify the way in which different sites and aquifers respond to recharge. As convolution models are developed for additional aquifers, they could contribute to an IRF database and a general classification system for karst aquifers.

South Dakota, Texas↗

Informing amphibian conservation efforts with abundance-based metapopulation models

Science-based management strategies are needed to halt or reverse the global decline of amphibians. In many cases, sound management requires reliable models built using monitoring data. Historically, monitoring and statistical modeling efforts have focused on estimating occupancy using detection–nondetection data. Spatial occupancy models are useful for studying colonization–extinction dynamics, but richer insights can be gained from estimating abundance and density-dependent demographic rates. We developed an integrated abundance-based metapopulation model of the processes contributing to spatiotemporal variation in patch population density. We fit our model to a combination of detection–nondetection and count data from a 14-yr study of a reintroduced metapopulation of federally threatened Chiricahua Leopard Frogs (Lithobates chiricahuensis). Pond-specific population growth rate was influenced by pond hydroperiod and frog density, such that permanent and semipermanent ponds with low densities of adult frogs experienced the highest annual population growth rates. Immigration rate declined as the distance among ponds increased. After reintroduction in 2003, metapopulation-level abundance increased and appeared to stabilize around 1300 adult frogs (95% CI = 1192–1471) by year 2015. Further, changes in metapopulation abundance were driven mostly by changes in abundance at a few ponds. These high-density populations, which would not have been identifiable with traditional occupancy-based metapopulation models, are likely especially important for species recovery in the area. Abundance-based metapopulation models can be widely applied to inform conservation efforts, by providing higher quality information needed to prioritize habitat patches for management and can be used to make more accurate predictions of metapopulation extinction risk.

Herpetologica↗

Limits of mountain and continental glaciations east of the Continental Divide in northern Montana and north-western North Dakota, U.S.A.

This chapter provides an overview of the limits of glaciations and glacial history in, and east and south-east of, Glacier National Park, Montana, and on the Northern Plains further east in Montana and north-western North Dakota. The term “Laurentide glacier” was applied to a continental ice sheet east of the Rocky Mountains in North America. It describes Laurentide Ice Sheet as any Quaternary continental ice sheet east of the Rocky Mountains in the United States and Canada. Laurentide till refers to till deposited by a Laurentide Ice Sheet. A Laurentide continental ice sheet is distinguished from a Cordilleran continental ice sheet in the Cordilleran region in parts of Washington, Idaho, and Montana in the United States and in adjacent Canada. Clague indicated that Cordilleran Ice Sheets formed several times during the Pleistocene. The chapter also reviews that the base for the digital map is simplified. Selected hydrographic features, selected towns and cities, selected physiographic features, and a grid of 1° × 2° topographic quadrangles are included to aid the reader in location of the glacial limits and other features depicted here on other maps at different scales.

Montana, North Dakota↗

Use and environmental occurrence of pharmaceuticals in freestall dairy farms with manured forage fields

Environmental releases of antibiotics from concentrated animal feeding operations (CAFOs) are of increasing regulatory concern. This study investigates the use and occurrence of antibiotics in dairy CAFOs and their potential transport into first-encountered groundwater. On two dairies we conducted four seasonal sampling campaigns, each across 13 animal production and waste management systems and associated environmental pathways: application to animals, excretion to surfaces, manure collection systems, soils, and shallow groundwater. Concentrations of antibiotics were determined using on line solid phase extraction (OLSPE) and liquid chromatography-tandem mass spectrometry (LC/MS/MS) with electrospray ionization (ESI) for water samples, and accelerated solvent extraction (ASE) LC/MS/MS with ESI for solid samples. A variety of antibiotics were applied at both farms leading to antibiotics excretion of several hundred grams per farm per day. Sulfonamides, tetracyclines, and their epimers/isomers, and lincomycin were most frequently detected. Yet, despite decades of use, antibiotic occurrence appeared constrained to within farm boundaries. The most frequent antibiotic detections were associated with lagoons, hospital pens, and calf hutches. When detected below ground, tetracyclines were mainly found in soils, whereas sulfonamides were found in shallow groundwater reflecting key differences in their physicochemical properties. In manure lagoons, 10 compounds were detected including tetracyclines and trimethoprim. Of these 10, sulfadimethoxine, sulfamethazine, and lincomycin were found in shallow groundwater directly downgradient from the lagoons. Antibiotics were sporadically detected in field surface samples on fields with manure applications, but not in underlying sandy soils. Sulfadimethoxine and sulfamethazine were detected in shallow groundwater near field flood irrigation gates, but at highly attenuated levels.

Environmental Science & Technology↗

Design of exploration and minerals-data-collection programs in developing areas

This paper considers the practical problem of applying economic analysis to designing minerals exploration and data collection strategies for developing countries. Formal decision rules for the design of government exploration and minerals-data-collection programs are derived by using a minerals-industry planning model that has been extended to include an exploration function. Rules derived are applicable to centrally planned minerals industries as well as market-oriented minerals sectors. They pertain to the spatial allocation of exploration effort and to the allocation of activities between government and private concerns for market-oriented economies. Programs characterized by uniform expenditures, uniform information coverage across regions, or uniform-density grid drilling progrmas are shown to be inferior to the strategy derived. Moreover, for market-oriented economies, the economically optimal mix in exploration activities between private and government data collection would require that only private firms assess local sites and that government agencies carry out regional surveys.

Socio-Economic Planning Sciences↗

Nitrapyrin in streams: The first study documenting off-field transport of a nitrogen stabilizer compound

Nitrapyrin is a bactericide that is co-applied with fertilizer to prevent nitrification and enhance corn yields. While there have been studies of the environmental fate of nitrapyrin, there is no documentation of its off-field transport to streams. In 2016, 59 water samples from 11 streams across Iowa were analyzed for nitrapyrin and its degradate, 6-chloropicolinic acid (6-CPA), along with three widely used herbicides, acetochlor, atrazine, and metolachlor. Nitrapyrin was detected in seven streams (39% of water samples) with concentrations ranging from 12 to 240 ng/L; 6-CPA was never detected. The herbicides were ubiquitously detected (100% of samples, 28–16000 ng/L). Higher nitrapyrin concentrations in streams were associated with rainfall events following spring fertilizer applications. Nitrapyrin persisted in streams for up to 5 weeks. These results highlight the need for more research focused on the environmental fate and transport of nitrapyrin and the potential toxicity this compound could have on nontarget organisms.

Environmental Science & Technology Letters↗

Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes

Reliable estimation of sensible heat flux ( H ) is important in energy balance models for quantifying evapotranspiration (ET). This study was conducted to evaluate the value of adding the Priestley-Taylor (PT) equation to the METRIC (Mapping Evapotranspiration at high Resolution with Internalized Calibration) model. METRIC was used to estimate energy fluxes for 10 Landsat images from the 2005, 2006 and 2007 crop growing seasons in south-central Nebraska, USA, where each image owing to recent rainfall exhibited high residual moisture content even at the hot pixel. The METRIC model performed satisfactorily for net radiation ( R n ) and soil heat flux ( G ) estimation with a root mean square error (RMSE) of 52 and 24 W m -2 , respectively. A RMSE of 122 W m -2 for H indicated the limitation of the METRIC model in estimating H for high residual moisture content of the hot pixel (Alfalfa reference ET fraction, ET r F > 0.15). The modified METRIC model (wet METRIC or wMETRIC) incorporating the PT equation was applied to calculate H at the anchor pixels (hot and cold) for high residual moisture content of the hot pixel. The α coefficient of the PT equation was locally calibrated using hourly meteorological data from an automatic weather station and R n and G data from a Bowen ratio flux tower. The mean α coefficient value was 1.14. The wMETRIC model reduced the RMSE of H from 122 to 64 W m -2 and that of latent heat flux, LE, from 163 to 106 W m -2 . The RMSE of daily ET decreased from 1.7 to 1.1 mm d -1 with wMETRIC. The results indicate that treatment of anchor pixels for high residual moisture content with the PT approach gives improved estimation of H , LE and daily ET. It is recommended that the wMETRIC model be used for estimating ET if the hot pixel has high residual moisture (i.e. reference ET fraction > 0.15).

Nebraska↗

Behavioral responses of sea lamprey (Petromyzon marinus) and white sucker (Catostomus commersonii) to turbulent flow during fishway passage attempts

An understanding of how undesirable and desirable fish species respond behaviorally to turbulent flow in fishways would guide development of selective fish passage techniques. We applied high-resolution computational fluid dynamics modeling and competing risks analysis towards the development of predictive selective passage models. Sea lamprey ( Petromyzon marinus ; an invasive fish in the Great Lakes Basin, North America) upstream passage probability declined from 0.73 to 0.03 as flow conditions became increasingly turbulent, while declines in white sucker ( Catostomus commersonii , a native fish in the region) upstream passage probability were less substantial (0.53 to 0.44). Deploying a sea lamprey trap in the fishway did not effectively reduce sea lamprey upstream passage probability, though capture rate increased during trials with cooler water temperature and low total kinetic energy. Bifurcated fishways that maintain low turbulent flow in the entrapment route and high turbulent flow in the upstream passage route could increase the effectiveness of trapping sea lamprey in fishways as a means to advance selective passage goals.

Wisconsin↗

Food webs of the Delta, Suisun Bay and Suisun Marsh: an update on current understanding and possibilities for management

This paper reviews and highlights recent research findings on foodweb processes since an earlier review by Kimmerer et al. (2008). We conduct this review within a conceptual framework of the Delta-Suisun food web, which includes both temporal and spatial components. The temporal component is based on knowledge that the landscape has changed markedly from historical conditions. The spatial component of our framework acknowledges that the food web is not spatially static; it varies regionally and across habitat types within regions. The review highlights the idea of a changing baseline with respect to foodweb function. New research also indicates that interactions between habitat-specific food webs vary across the current landscape. For example, based on early work in the South Delta, the food web associated with submerged aquatic vegetation was thought to provide little support to species of concern; however, data from other regions of the estuary suggest that this conceptual model may not apply across the entire region. Habitat restoration has been proposed as a method of re-establishing historic foodweb processes to support species of concern. Benefits are likely for species that directly access such restored habitats, but are less clear for pelagic species. Several topics require attention to further improve the knowledge of food webs needed to support effective management, including: 1) synthesis of factors responsible for low pelagic biomass; 2) monitoring and research on effects of harmful algal blooms; 3) broadening the scope of long-term monitoring; 4) determining benefits of tidal wetland restoration to species of concern, including evaluations of interactions of habitat-specific food webs; and 5) interdisciplinary analysis and synthesis. The only certainty is that food webs will continue to change in response to the changes in the physical environment and new species invasions.

Californina↗

A framework tool that applies weight-of-evidence integration to the analysis of existing datasets to guide freshwater conservation

The overarching issue we address here is how to extract clear and actionable ecological and management insights from real-world field data that often do not satisfy traditional statistical assumptions. Toward this goal, we developed a general 12+6 step adaptive management framework tool. We applied this framework tool to existing biodiversity monitoring data to create a proof-of-concept result that addresses the overarching question of “ why might a specific native stream fish taxon be present or absent at specific locations? ” Our multi-step framework tool links established steps and steps that are unique to our framework through weight-of evidence (WOE) integration, an approach that combines quantitative results from multiple visualization and statistical procedures. The systematic use of all steps in our framework can provide improved conservation outcomes compared to a single analysis. Advantages accrue from our approach because our framework tool refines the overarching goal into related sub-questions, applies a specific quantitative procedure to each sub-step, combines results from all sub-questions using a WOE integration, identifies testable questions that elucidate ambiguities and gaps revealed through WOE integration, and proposes practical field methods for obtaining this clarifying information through future research and data collections. The process of considering multiple visualizations and analyses as individual pieces of a shared puzzle offers a new way to approach the use of existing data. Our team-based approach transforms the collection and analysis of existing data into a series of field tests that can guide future actions (e.g., data collection-analysis events, restoration initiatives, research). Habitat and impact regressors will vary with taxa and system, but our structured process tool has broad generality for a range of conservation issues in which freshwater systems are threatened by human impacts.

Kansas↗

A comprehensive freshwater mussel database for the Duck River Drainage, Tennessee

We have developed a comprehensive database for freshwater mussels for the Duck River drainage in Tennessee, including its largest tributary, the Buffalo River. This database is intended to serve as an expandable template that could be applied statewide. The Duck River is one of the most biologically diverse rivers in the world, with historically over 70 mussel species, and it has been selected as a priority watershed by multiple management and conservation entities. The database for this system compiles over 7,000 mussel records, spanning 200 years, from multiple Federal, State, academic, and private entities, representing 77 native species. The database is spatially explicit and includes temporal and methodological data for each record, and notes of negative survey data were made when possible. The database can facilitate the creation of distribution maps for each species and temporal maps of species richness to show watershed-wide trends. This project addresses the present lack of a centralized mussel database in Tennessee for a critical system. It will be available to facilitate species status assessments, inform conservation planning, and serve as a model for similar databases for other Tennessee watersheds.

Tennessee↗

An index of floodplain surface complexity

Floodplain surface topography is an important component of floodplain ecosystems. It is the primary physical template upon which ecosystem processes are acted out, and complexity in this template can contribute to the high biodiversity and productivity of floodplain ecosystems. There has been a limited appreciation of floodplain surface complexity because of the traditional focus on temporal variability in floodplains as well as limitations to quantifying spatial complexity. An index of floodplain surface complexity (FSC) is developed in this paper and applied to eight floodplains from different geographic settings. The index is based on two key indicators of complexity, variability in surface geometry (VSG) and the spatial organisation of surface conditions (SPO), and was determined at three sampling scales. FSC, VSG, and SPO varied between the eight floodplains and these differences depended upon sampling scale. Relationships between these measures of spatial complexity and seven geomorphological and hydrological drivers were investigated. There was a significant decline in all complexity measures with increasing floodplain width, which was explained by either a power, logarithmic, or exponential function. There was an initial rapid decline in surface complexity as floodplain width increased from 1.5 to 5 km, followed by little change in floodplains wider than 10 km. VSG also increased significantly with increasing sediment yield. No significant relationships were determined between any of the four hydrological variables and floodplain surface complexity.

Hydrology and Earth System Sciences↗

Investigating the chemical space coverage of multiple chromatographic and ionization methods using non-targeted analysis on surface and drinking water collected using passive sampling

Multiple non-targeted analysis tools were used to look for a broad range of possible chemical contaminants present in surface and drinking water using liquid chromatography separation and high-resolution mass spectrometry detection, including both quadrupole time of flight (Q-ToF) and Orbitrap instruments. Two chromatographic techniques were evaluated on an LC-Q-ToF with electrospray ionization in both positive and negative modes: (1) the traditionally used reverse phase C18 and (2) the hydrophilic interaction liquid chromatography (HILIC) aimed to capture more polar contaminants that may be present in water. Multiple ionization modes were evaluated with an LC-Orbitrap, including electrospray (ESI) and atmospheric pressure chemical ionization (APCI), also in both positive and negative modes. A suspect screening library of over 1300 possible environmental contaminants, including pesticides, pharmaceuticals, personal care products, illicit drugs/drugs of abuse, and various anthropogenic markers was made with experimentally collected data with the LC-Q-ToF with both column types, with 227 chemicals being retained by the HILIC column. The non-targeted methods using multiple chromatographic and ionization modes were applied to environmental water samples collected with polar organic chemical integrative samplers (POCIS), including surface water upstream and downstream from wastewater effluent discharge, and the downstream drinking water intake and treated drinking water for three distinct sampling events. For the LC-Q-ToF, 442 chemical features were detected on the C18 column and 91 with the HILIC column in the POCIS extracts, while 556 features were found on the Orbitrap workflow by ESI and 131 features detected by APCI. Over 100 chemicals were tentatively identified by suspect screening and database searching. The comprehensive and systematic evaluation of these methods serve as a step in characterizing the chemical space covered when utilizing different chromatography and ionization methods, or different instrument workflows on complex environmental mixtures.

Science of the Total Environment↗

Natural background concentrations of nutrients in streams and rivers of the conterminous United States

Determining natural background concentrations of nutrients in watersheds in the developed world has been hampered by a lack of pristine sampling sites covering a range of climatic conditions and basin sizes. Using data from 63 minimally impacted U.S. Geological Survey reference basins, we developed empirical models of the background yield of total nitrogen (TN) and total phosphorus (TP) from small watersheds as functions of annual runoff, basin size, atmospheric nitrogen deposition rate, and region-specific factors. We applied previously estimated in-stream loss rates to yields from the small watershed models to obtain estimates of background TN and TP yield and concentration throughout the stream/river network in 14 ecoregions of the conterminous United States. Background TN concentration varies from less than 0.02 mg L-1 in the xeric west to more than 0.5 mg L-1 along the southeastern coastal plain. Background TP concentration varies from less than 0.006 mg L-1 in the xeric west to more than 0.08 mg L-1 in the great plains. TN concentrations in U.S. streams and rivers currently exceed natural background levels by a much larger factor (6.4) than do TP concentrations (2.0). Because of local variation in runoff and other factors, the range of background nutrient concentrations is very large within some nutrient ecoregions. It is likely that background concentrations in some streams in these regions exceed proposed nutrient criteria.

Environmental Science & Technology↗

Quantifying groundwater’s role in delaying improvements to Chesapeake Bay water quality

A study has been undertaken to determine the time required for the effects of nitrogen-reducing best management practices (BMPs) implemented at the land surface to reach the Chesapeake Bay via groundwater transport to streams. To accomplish this, a nitrogen mass-balance regression (NMBR) model was developed and applied to seven watersheds on the Delmarva Peninsula. The model included the distribution of groundwater return times obtained from a regional groundwater-flow (GWF) model, the history of nitrogen application at the land surface over the last century, and parameters that account for denitrification. The model was (1) able to reproduce nitrate concentrations in streams and wells over time, including a recent decline in the rate at which concentrations have been increasing, and (2) used to forecast future nitrogen delivery from the Delmarva Peninsula to the Bay given different scenarios of nitrogen load reduction to the water table. The relatively deep porous aquifers of the Delmarva yield longer groundwater return times than those reported earlier for western parts of the Bay watershed. Accordingly, several decades will be required to see the full effects of current and future BMPs. The magnitude of this time lag is critical information for Chesapeake Bay watershed managers and stakeholders.

Delaware, Maryland↗

Habitat use of anadromous and amphidromous sturgeons in North America: A systematic review

Sturgeons are among the most endangered fishes in the world. Identifying habitat use characteristics to inform restoration projects is crucial for recovery. However, small sample sizes, inadequate replication of studies, and limited spatial extents complicate our ability to effectively apply the findings of single studies to endangered species conservation across the larger riverscape. We synthesized information from amphidromous and anadromous sturgeons in North America to identify species-specific knowledge gaps and conduct a quantitative comparison of species–habitat relationships. We provided a qualitative summary of substrate use and synthesized estimates of depth and velocity during spawning and non-spawning activity. Generalized patterns among species were identified, such as spawning in fast water on hard substrate and then using slow water with soft substrate areas when not spawning. We noted species-specific variability during spawning that may be attributed to historical maximum length, egg characteristics, and watershed features. This study provides some of the first estimates of habitat use that can be adapted for many populations. Results can contribute to empirically grounded decision-support tools used to prioritize information needs for recovery.

Canadian Journal of Fisheries and Aquatic Sciences↗

Book review: Bird census techniques, Second edition

Conservation concerns, federal mandates to monitor birds, and citizen science programs have spawned a variety of surveys that collect information on bird populations. Unfortunately, all too frequently these surveys are poorly designed and use inappropriate counting methods. Some of the flawed approaches reflect a lack of understanding of statistical design; many ornithologists simply are not aware that many of our most entrenched counting methods (such as point counts) cannot appropriately be used in studies that compare densities of birds over space and time. It is likely that most of the readers of The Condor have participated in a bird population survey that has been criticized for poor sampling methods. For example, North American readers may be surprised to read in Bird Census Techniques that the North American Breeding Bird Survey 'is seriously flawed in its design,' and that 'Analysis of trends is impossible from points that are positioned along roads' (p. 109). Our conservation efforts are at risk if we do not acknowledge these concerns and improve our survey designs. Other surveys suffer from a lack of focus. In Bird Census Techniques, the authors emphasize that all surveys require clear statements of objectives and an understanding of appropriate survey designs to meet their objectives. Too often, we view survey design as the realm of ornithologists who know the life histories and logistical issues relevant to counting birds. This view reflects pure hubris: survey design is a collaboration between ornithologists, statisticians, and managers, in which goals based on management needs are met by applying statistical principles for design to the biological context of the species of interest. Poor survey design is often due to exclusion of some of these partners from survey development. Because ornithologists are too frequently unaware of these issues, books such as Bird Census Techniques take on added importance as manuals for educating ornithologists about the relevance of survey design and methods and the often subtle interdisciplinary nature of surveys. Review info: Bird Census Techniques, Second Edition. By Colin J. Bibby, Neil D. Burgess, David A. Hill, and Simon H. Mustoe. 2000. Academic Press, London, UK. xvii 1 302 pp. ISBN 0- 12-095831-7.

Condor↗

Lead isotopes in North Pacific deep water - Implications for past changes in input sources and circulation patterns

The sources of non-anthropogenic Pb in seawater have been the subject of debate. Here we present Pb isotope time-series that indicate that the non-anthropogenic Pb budget of the northernmost Pacific Ocean has been governed by ocean circulation and riverine inputs, which in turn have ultimately been controlled by tectonic processes. Despite the fact that the investigated locations are situated within the Asian dust plume, and proximal to extensive arc volcanism, eolian contributions have had little impact. We have obtained the first high-resolution and high-precision Pb isotope time-series of North Pacific deep water from two ferromanganese crusts from the Gulf of Alaska in the NE Pacific Ocean, and from the Detroit Seamount in the NW Pacific Ocean. Both crusts were dated applying 10 Be/9Be ratios and yield continuous time-series for the past 13.5 and 9.6 Myr, respectively. Lead isotopes show a monotonic evolution in 206Pb/204Pb from low values in the Miocene (??? 18.57) to high values at present day (??? 18.84) in both crusts, even though they are separated by more than 3000 km along the Aleutian Arc. The variation exceeds the amplitude found in Equatorial Pacific deep water records by about three-fold. There also is a striking similarity in 207Pb/204Pb and 208Pb/ 204Pb ratios of the two crusts, indicating the existence of a local circulation cell in the sub-polar North Pacific, where efficient lateral mixing has taken place but only limited exchange (in terms of Pb) with deep water from the Equatorial Pacific has occurred. Both crusts display well-defined trends with age in Pb-Pb isotope mixing plots, which require the involvement of at least four distinct Pb sources for North Pacific deep water. The Pb isotope time-series reveal that eolian supplies (volcanic ash and continent-derived loess) have only been of minor importance for the dissolved Pb budget of marginal sites in the deep North Pacific over the past 6 Myr. The two predominant sources have been young volcanic arcs, one located in the northeastern part and one located in the northwestern part of the Pacific margin, from where material has been eroded and delivered to the ocean, most likely via riverine pathways. ?? 2003 Elsevier Science B.V. All rights reserved.

Earth and Planetary Science Letters↗