Search NASA⌕ Search

SEARCH · Search NASA

Results for “Primary Production”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 667 records · Page 37

Climate, Hydrology, and Nutrients Control the Seasonality of Si Concentrations in Rivers

Abstract The seasonal behavior of fluvial dissolved silica (DSi) concentrations, termedDSi regime, mediates the timing of DSi delivery to downstream waters and thus governs river biogeochemical function and aquatic community condition. Previous work identified five distinct DSi regimes across rivers spanning the Northern Hemisphere, with many rivers exhibiting multiple DSi regimes over time. Several potential drivers of DSi regime behavior have been identified at small scales, including climate, land cover, and lithology, and yet the large‐scale spatiotemporal controls on DSi regimes have not been identified. We evaluate the role of environmental variables on the behavior of DSi regimes in nearly 200 rivers across the Northern Hemisphere using random forest models. Our models aim to elucidate the controls that give rise to (a) average DSi regime behavior, (b) interannual variability in DSi regime behavior (i.e., Annual DSi regime), and (c) controls on DSi regime shape (i.e., minimum and maximum DSi concentrations). Average DSi regime behavior across the period of record was classified accurately 59% of the time, whereas Annual DSi regime behavior was classified accurately 80% of the time. Climate and primary productivity variables were important in predicting Average DSi regime behavior, whereas climate and hydrologic variables were important in predicting Annual DSi regime behavior. Median nitrogen and phosphorus concentrations were important drivers of minimum and maximum DSi concentrations, indicating that these macronutrients may be important for seasonal DSi drawdown and rebound. Our findings demonstrate that fluctuations in climate, hydrology, and nutrient availability of rivers shape the temporal availability of fluvial DSi.

Environmental Sciences & Ecology↗

Insights Into Nature-Based Climate Solutions: Managing Forests for Climate Resilience and Carbon Stability

Successful implementation of forest management as a nature-based climate solution is dependent on the durability of management-induced changes in forest carbon storage and sequestration. As forests face unprecedented stability risks in the face of ongoing climate change, much remains unknown regarding how management will impact forest stability, or how interactions with climate might shift the response of forests to management across spatiotemporal scales. Here, we used a process-based model to simulate multidecadal projections of forest dynamics in response to changes in management and climate. Simulations were conducted across gradients in forest type, edaphic factors, and management intensity under two alternate radiative forcing scenarios (RCP4.5 and RCP8.5). This allowed for the quantification of forest stability shifts in response to climate change, and the role of management in modulating that response, where ecosystem stability is characterized as the resilience and temporal stability of net primary production, aboveground biomass, and soil carbon. Our results indicate that forest structure is primarily shaped by management, but the same management strategy often produced divergent structures over time, due to interactions with regional climate change. We found that management can be used to increase stability and minimize the release of stored carbon by reducing mortality, but also highlight the regional dependency of management-induced changes in resilience to climate change.

54 ENVIRONMENTAL SCIENCES↗

Coupling Remote Sensing With a Process Model for the Simulation of Rangeland Carbon Dynamics

Rangelands provide significant environmental benefits through many ecosystem services, which may include soil organic carbon (SOC) sequestration. However, quantifying SOC stocks and monitoring carbon (C) fluxes in rangelands are challenging due to the considerable spatial and temporal variability tied to rangeland C dynamics as well as limited data availability. We developed the Rangeland Carbon Tracking and Management (RCTM) system to track long-term changes in SOC and ecosystem C fluxes by leveraging remote sensing inputs and environmental variable data sets with algorithms representing terrestrial C-cycle processes. Bayesian calibration was conducted using quality-controlled C flux data sets obtained from 61 Ameriflux and NEON flux tower sites from Western and Midwestern US rangelands to parameterize the model according to dominant vegetation classes (perennial and/or annual grass, grass-shrub mixture, and grass-tree mixture). The resulting RCTM system produced higher model accuracy for estimating annual cumulative gross primary productivity (GPP) (R 2 > 0.6, RMSE <390 g C m -2 ) relative to net ecosystem exchange of CO 2 (NEE) (R 2 > 0.4, RMSE <180 g C m -2 ). Model performance in estimating rangeland C fluxes varied by season and vegetation type. The RCTM captured the spatial variability of SOC stocks with R 2 = 0.6 when validated against SOC measurements across 13 NEON sites. Model simulations indicated slightly enhanced SOC stocks for the flux tower sites during the past decade, which is mainly driven by an increase in precipitation. Future efforts to refine the RCTM system will benefit from long-term network-based monitoring of vegetation biomass, C fluxes, and SOC stocks.

54 ENVIRONMENTAL SCIENCES↗

Incorporating the Acclimation of Photosynthesis and Leaf Respiration in the Noah‐MP Land Surface Model: Model Development and Evaluation

Abstract Realistic simulation of leaf photosynthetic and respiratory processes is needed for accurate prediction of the global carbon cycle. These two processes systematically acclimate to long‐term environmental changes by adjusting photosynthetic and respiratory traits (e.g., the maximum photosynthetic capacity at 25°C ( V cmax,25 ) and the leaf respiration rate at 25°C ( R 25 )) following increasingly well‐understood principles. While some land surface models (LSMs) now account for thermal acclimation, they do so by assigning empirical parameterizations for individual plant functional types (PFTs). Here, we have implemented an Eco‐Evolutionary Optimality (EEO)‐based scheme to represent the universal acclimation of photosynthesis and leaf respiration to multiple environmental effects, and that therefore requires no PFT‐specific parameterizations, in a standard version of the widely used LSM, Noah MP. We evaluated model performance with plant trait data from a 5‐year experiment and extensive global field measurements, and carbon flux measurements from FLUXNET2015. We show that observed R 25 and V cmax,25 vary substantially both temporally and spatially within the same PFT ( C.V. >20%). Our EEO‐based scheme captures 62% of the temporal and 70% of the spatial variations in V cmax,25 (73% and 54% of the variations in R 25 ). The standard scheme underestimates gross primary production by 10% versus 2% for the EEO‐based scheme and generates a larger spread in r (correlation coefficient) across flux sites (0.79 ± 0.16 vs. 0.84 ± 0.1, mean ± S.D.). The standard scheme greatly overestimates canopy respiration (bias: ∼200% vs. 8% for the EEO scheme), resulting in less CO 2 uptake by terrestrial ecosystems. Our approach thus simulates climate‐carbon coupling more realistically, with fewer parameters.

Ren, Yanghang [Department of Earth System Science ↗

The Imprint of Southern Ocean Storms on Modeled Surface Chlorophyll, Their Drivers and Satellite Biases

Southern Ocean (SO) phytoplankton chlorophyll is highly variable on sub-seasonal time scales. Although the SO is the windiest ocean basin globally, it is not conclusively understood how storms impact SO phytoplankton dynamics. Much of our existing knowledge stems from satellites, but biases due to data gaps from cloud cover and low solar angles remain unquantified. Here, we use ocean–sea-ice simulations with the Community Earth System Model to quantify the climatological 1997–2018 imprint of storms on chlorophyll and phytoplankton dynamics in the ice-free SO. Additionally, by comparing the full-field model output to synthetic satellite observations, we quantify sampling biases in satellite-derived estimates. We find that both the sign and the magnitude of the average surface chlorophyll imprint vary substantially across storms but last for at least 4 days after the storm passing. Based on our analysis, more than one third of the storms explain the majority of local non-seasonal chlorophyll variability, but satellite-derived storm imprints are often too large in magnitude. On the day of the storm passing, changes in vertical mixing predominantly cause surface chlorophyll anomalies, and reduced light availability due to enhanced cloud cover outweighs the enhanced nutrient availability due to entrainment. Interestingly, storms imprint differently on total net primary production than on surface chlorophyll, demonstrating the difficulty to derive carbon-cycle impacts from a surface-chlorophyll assessment. With SO future storm activity projected to increase, complementing satellite observations with other observing technologies, for example, profiling floats, is necessary to better constrain how storms impact biological carbon cycling in the SO.

Nissen, Cara [University of Colorado, Boulder, CO ↗

Ubiquity and Causes of Soil Water Preferential Flow Across 17 Ecoregions

Abstract Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high‐frequency, multi‐depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation. Plain Language Summary Water can bypass part of the soil's matrix through a process called preferential flow (PF). This quick transport of water through the soil brings with it nutrients and contaminants and eventually makes it to groundwater and streams. To ensure ample amounts of good quality groundwater and surface water we need to understand when and where PF occurs. We inferred when PF occurred at 40 different sites across 17 ecoregions in the USA using soil moisture and rainfall data. We found that PF happened at all sites and in up to 60% of rainfall events ≥2 mm. Preferential flow was most likely at sites with high rainfall intensities, high clay content in soils, low variability in soil moisture, and high vegetation productivity. As rainfall intensities are predicted to increase due to climate change and vegetation becomes more productive, PF occurrence becomes more important for predicting groundwater recharge, water quality, and streamflow generation. Key Points Preferential flow (PF) is ubiquitous across the USA and occurs in up to 60% of all rainfall events ≥2 mm Rainfall intensity, soil texture, and antecedent soil moisture emerge as critical in generating PF across diverse ecoregions Two different PF detection approaches show similar relationships between key drivers and occurrence of PF

Li, Bonan↗

Geochemical Phosphorus Sequestration in Tundra Soils Impedes Delivery of Bioavailable Phosphorus to the Kuparuk River, Alaska, USA: Implications for the Broader Arctic Region

Long-term river monitoring of the Kuparuk River (North Slope, Alaska, USA) confirms significant increases in solutes that are indicative of active layer thickening due to thawing permafrost. However, there is no evidence of an increase in total dissolved phosphorus (TDP) or soluble reactive phosphorus (SRP), the nutrient that limits primary production in this and similar rivers in the region. Here, we show that Mehlich-3 extractable iron (Fe) and aluminum (Al) in active layer soils impart high P geochemical sorption capacities across a range of landscape features that we would expect to promote lateral movement of water and solutes to headwater streams in our study watershed. Reanalysis of a recently published pan-Arctic soils database that includes active layer and permafrost soil samples suggests that this high P sorption capacity could be common in other parts of the Arctic region. We conclude that soil minerals enhance P retention on hillslopes and propose pedogenic secondary Fe and Al minerals may continue to retain P in these soils and limit biological productivity in the adjacent river even as active layer thickening increases potential P mobility in the watershed. We suggest that similar interactions may occur in other areas of the Arctic where comparable geochemical conditions prevail.

Sutor, Frederick W. [Univ. of Vermont, Burlington,↗

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape

Arctic warming is altering vegetation and carbon dynamics with global implications, yet Earth System Model (ESM) predictions in the Arctic remain highly uncertain, in part due to historically limited data for model parameterization and validation. As such, ESMs typically represent Arctic ecosystems in an oversimplified manner. Recently, nine plant functional types (PFTs) designed to realistically represent tundra vegetation were integrated into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) and parameterized using plot-scale observations from a single site. Additional evaluation was needed to determine their transferability across the Arctic. Here, in this study, we evaluated whether refined representation of tundra vegetation improved model accuracy by conducting spatially explicit 100 × 100 m resolution ELM simulations on Alaska's Seward Peninsula. Simulations with the default two-PFT configuration and with the nine Arctic-specific PFTs were benchmarked against observations of net ecosystem exchange, gross primary production, and aboveground biomass from multiple data streams including an eddy covariance flux tower, flux chambers, and aircraft and unoccupied aerial system hyperspectral remote sensing. Evaluation revealed that Arctic-specific PFT simulations produced more realistic landscape-level carbon exchanges, and better captured observed heterogeneity in biomass and productivity, explaining 60%–70% of spatial variance (R 2 = 0.6–0.7) compared to just 12%–18% (R 2 = 0.12–0.18) with the default configuration. However, the refined model failed to reproduce observed aboveground biomass for highly productive alder-willow communities, requiring further evaluation of carbon allocation parameterizations for tall shrubs that are increasingly expanding across tundra landscapes. Our results demonstrate that enhanced representation of vegetation heterogeneity boosts predictive understanding of tundra carbon dynamics, facilitating regional to pan-Arctic model and remote-sensing scaling.

Murphy, Bailey A. [Oak Ridge National Laboratory (↗

Investigating Coastal Vegetation Dynamics and Ecosystem Impacts Under Elevated CO 2 and Temperature: A Process‐Based Approach

Coastal forests are increasingly vulnerable to climate change and sea-level rise, with flooding and salinity driving transitions to marsh-dominated ecosystems. Using the coastal version of FATES-Hydro, we conducted 30-year simulations at two coastal forest sites—a broadleaf swamp white oak stand at Lake Erie and a conifer loblolly pine stand at Chesapeake Bay—under historical climate and elevated CO 2 (+100 ppm) and temperature (+1.5°C) scenarios. Elevated CO 2 increased net primary productivity at both sites, while warming alone intensified hydraulic stress and accelerated mortality, particularly in the conifer stand. Simulations show that elevated temperatures intensify vapor pressure deficit and hydraulic stress on trees already experiencing salinity- and submersion-driven water stress, increasing tree mortality beyond what would be expected in a non-water-limited environment. Marsh expansion partially compensated for tree loss at the Lake Erie site but reduced ecosystem productivity in the conifer forest at Chesapeake Bay. In conclusion, our results highlight how differences in stand structure, phenology, and local hydrology modulate ecosystem trajectories under climate change, emphasizing the importance of demographic and community-level processes for predicting the fate of coastal forests.

Ding, Junyan [Barcelona Supercomputing Center (BSC↗

Contrasting Parametric Sensitivities in Two Global Vegetation Models Using Parameter Perturbation Ensembles

Uncertainty in land model projections remains high and the roles of parametric and structural uncertainty are difficult to disentangle. To compare parametric sensitivity across model structures we present two parameter perturbation ensembles using the Community Land Model (CLM) operating in satellite phenology mode. The ensembles contrast two vegetation modules: (a) the default CLM vegetation module and (b) the Functionally Assembled Terrestrial Ecosystem Simulator (CLM-FATES). We perturbed over 300 parameters and quantified their effects on biophysical fluxes globally and across biomes. Most parameters have minimal impact on biophysical fluxes, with only a few substantially influencing results. While both models exhibit similar parameter sensitivity for some fluxes, CLM-FATES shows larger spread in gross primary productivity (GPP), driven by strong sensitivity to carboxylation rate. CLM-FATES also shows a weaker GPP response to soil hydrology parameters and exhibits higher water use efficiency (WUE). Cross-model comparisons reveal similar sensitivities for some parameters (e.g., leaf dimension) but divergent responses to others (e.g., stomatal intercept), highlighting underlying structural differences. Differences in WUE and sensitivity to hydrology and stomatal conductance parameters underscore how model structure fundamentally alters parametric sensitivity. The data sets generated from these ensembles can be used to identify influential parameters and guide future calibration efforts.

Foster, A. C. [NSF National Center for Atmospheric↗

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗

The impact of nitrogen oxides on electrochemical carbon dioxide reduction

The electroreduction of carbon dioxide offers a promising avenue to produce valuable fuels and chemicals using greenhouse gas carbon dioxide as the carbon feedstock. Because industrial carbon dioxide point sources often contain numerous contaminants, such as nitrogen oxides, understanding the potential impact of contaminants on carbon dioxide electrolysis is crucial for practical applications. Herein, we investigate the impact of various nitrogen oxides, including nitric oxide, nitrogen dioxide, and nitrous oxide, on carbon dioxide electroreduction on three model electrocatalysts (i.e., copper, silver, and tin). We demonstrate that the presence of nitrogen oxides (up to 0.83%) in the carbon dioxide feed leads to a considerable Faradaic efficiency loss in carbon dioxide electroreduction, which is caused by the preferential electroreduction of nitrogen oxides over carbon dioxide. The primary products of nitrogen oxides electroreduction include nitrous oxide, nitrogen, hydroxylamine, and ammonia. Despite the loss in Faradaic efficiency, the electrocatalysts exhibit similar carbon dioxide reduction performances once a pure carbon dioxide feed is restored, indicating a negligible long-term impact of nitrogen oxides on the catalytic properties of the model catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adaptive responses of marine diatoms to zinc scarcity and ecological implications

Scarce dissolved surface ocean concentrations of the essential algal micronutrient zinc suggest that Zn may influence the growth of phytoplankton such as diatoms, which are major contributors to marine primary productivity. However, the specific mechanisms by which diatoms acclimate to Zn deficiency are poorly understood. Using global proteomic analysis, we identified two proteins (ZCRP-A/B, Zn/Co Responsive Protein A/B) among four diatom species that became abundant under Zn/Co limitation. Characterization using reverse genetic techniques and homology data suggests putative Zn/Co chaperone and membrane-bound transport complex component roles for ZCRP-A (a COG0523 domain protein) and ZCRP-B, respectively. Metaproteomic detection of ZCRPs along a Pacific Ocean transect revealed increased abundances at the surface (<200 m) where dZn and dCo were scarcest, implying Zn nutritional stress in marine algae is more prevalent than previously recognized. These results demonstrate multiple adaptive responses to Zn scarcity in marine diatoms that are deployed in low Zn regions of the Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗

High impact of bacterial predation on cyanobacteria in soil biocrusts

Diverse bacteria lead a life as pathogens or predators of other bacteria in many environments. However, their impact on emerging ecological processes in natural settings remains to be assessed. Here we describe a novel type of obligate, intracellular predatory bacterium of widespread distribution that preys on soil cyanobacteria in biocrusts. The predator, Candidatus Cyanoraptor togatus, causes localized, cm-sized epidemics that are visible to the naked eye, obliterates cyanobacterial net primary productivity, and severely impacts crucial biocrust properties like nitrogen cycling, dust trapping and moisture retention. The combined effects of high localized morbidity and areal incidence result in decreases approaching 10% of biocrust productivity at the ecosystem scale. Our findings show that bacterial predation can be an important loss factor shaping not only the structure but also the function of microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗

Increased photosynthesis during spring drought in energy-limited ecosystems

Abstract Drought is often thought to reduce ecosystem photosynthesis. However, theory suggests there is potential for increased photosynthesis during meteorological drought, especially in energy-limited ecosystems. Here, we examine the response of photosynthesis (gross primary productivity, GPP) to meteorological drought across the water-energy limitation spectrum. We find a consistent increase in eddy covariance GPP during spring drought in energy-limited ecosystems (83% of the energy-limited sites). Half of spring GPP sensitivity to precipitation was predicted solely from the wetness index (R 2 = 0.47, p < 0.001), with weaker relationships in summer and fall. Our results suggest GPP increases during spring drought for 55% of vegetated Northern Hemisphere lands ( >30° N). We then compare these results to terrestrial biosphere model outputs and remote sensing products. In contrast to trends detected in eddy covariance data, model mean GPP always declined under spring precipitation deficits after controlling for air temperature and light availability. While remote sensing products captured the observed negative spring GPP sensitivity in energy-limited ecosystems, terrestrial biosphere models proved insufficiently sensitive to spring precipitation deficits.

54 ENVIRONMENTAL SCIENCES↗

Catalytic carbon–carbon bond cleavage in lignin via manganese–zirconium-mediated autoxidation

Abstract Efforts to produce aromatic monomers through catalytic lignin depolymerization have historically focused on aryl–ether bond cleavage. A large fraction of aromatic monomers in lignin, however, are linked by various carbon–carbon (C–C) bonds that are more challenging to cleave and limit the yields of aromatic monomers from lignin depolymerization. Here, we report a catalytic autoxidation method to cleave C–C bonds in lignin-derived dimers and oligomers from pine and poplar. The method uses manganese and zirconium salts as catalysts in acetic acid and produces aromatic carboxylic acids as primary products. The mixtures of the oxygenated monomers are efficiently converted to cis,cis -muconic acid in an engineered strain of Pseudomonas putida KT2440 that conducts aromatic O -demethylation reactions at the 4-position. This work demonstrates that autoxidation of lignin with Mn and Zr offers a catalytic strategy to increase the yield of valuable aromatic monomers from lignin.

09 BIOMASS FUELS↗

Implicating the red body of Nannochloropsis in forming the recalcitrant cell wall polymer algaenan

Stramenopile algae contribute significantly to global primary productivity, and one class, Eustigmatophyceae, is increasingly studied for applications in high-value lipid production. Yet much about their basic biology remains unknown, including the nature of an enigmatic, pigmented globule found in vegetative cells. Here, we present an in-depth examination of this “red body,” focusing on Nannochloropsis oceanica. During the cell cycle, the red body forms adjacent to the plastid, but unexpectedly it is secreted and released with the autosporangial wall following cell division. Shed red bodies contain antioxidant ketocarotenoids, and overexpression of a beta-carotene ketolase results in enlarged red bodies. Infrared spectroscopy indicates long-chain, aliphatic lipids in shed red bodies and cell walls, and UHPLC-HRMS detects a C32 alkyl diol, a potential precursor of algaenan, a recalcitrant cell wall polymer. We propose that the red body transports algaenan precursors from plastid to apoplast to be incorporated into daughter cell walls.

59 BASIC BIOLOGICAL SCIENCES↗

Reduced CO 2 uptake and growing nutrient sequestration from slowing overturning circulation

Current Earth system models (ESMs) project dramatic slowing (28–42% by 2100) of Atlantic Meridional Overturning Circulation and Southern Meridional Overturning Circulation (SMOC) across a range of climate scenarios, with a complete shutdown of SMOC possible by year 2300. Slowing meridional overturning circulation (MOC) differentially impacts the ocean biological and solubility carbon pumps, leaving the net impact on ocean carbon uptake uncertain. Here using a suite of ESMs, we show that slowing MOC reduces anthropogenic carbon uptake by the solubility pump but increases deep-ocean storage of carbon and nutrients by the biological pump. The net effect reduces ocean uptake of anthropogenic CO 2 . The deep-ocean nutrient sequestration will increasingly depress global-scale, marine net primary production over time. As a result, MOC slowdown represents a positive feedback that could extend or intensify peak-warmth climate conditions on multi-century timescales.

58 GEOSCIENCES↗