Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth Systems”

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 307 records · Page 17

Earth system dynamics: The interrelation of atmospheric, ocean and solid earth dynamics

The research work performed during the time period 16 Oct. 1992 through 31 Dec. 1993 is summarized. The overall research activity, including a list of the major findings of the EOS IDS research to date, is described, the publications and presentations are listed, and a budget request for the subsequent year is attached. Specifically, the report covers: EOS panel activities; major findings of research; team member contributions; new research directions; EOS restructuring effect; changes in requirements; plans for using existing data; collaborations with other EOS and non-EOS investigations; EOS instrument team interaction; instrument development verification and validation; interaction with EOSDIS and DAAC's; team coordination; overall management; summary of response to site review questions and comments; science computing facility; and additional new research activities.

Tapley, Byron D.↗

Orbital Noise in the Earth System is a Common Cause of Climate and Greenhouse-Gas Fluctuation

The mismatch between fossil isotopic data and climate models known as the cool-tropic paradox implies that either the data are flawed or we understand very little about the climate models of greenhouse warming. Here we question the validity of the climate models on the scientific background of orbital noise in the Earth system. Our study shows that the insolation pulsation induced by orbital noise is the common cause of climate change and atmospheric concentrations of carbon dioxide and methane. In addition, we find that the intensity of the insolation pulses is dependent on the latitude of the Earth. Thus, orbital noise is the key to understanding the troubling paradox in climate models.

Liu, H. S.↗

Analysis of secondary organic aerosol simulation bias in the Community Earth System Model (CESM2.1)

Organic aerosol (OA) has been considered as one of the most important uncertainties in climate modeling due to the complexity in presenting its chemical production and depletion mechanisms. To better understand the capability of climate models and probe into the associated uncertainties in simulating OA, we evaluate the Community Earth System Model version 2.1 (CESM2.1) configured with the Community Atmosphere Model version 6 (CAM6) with comprehensive tropospheric and stratospheric chemistry representation (CAM6-Chem) through a long-term simulation (1988–2019) with observations collected from multiple datasets in the United States. We find that CESM generally reproduces the interannual variation and seasonal cycle of OA mass concentration at surface layer with a correlation of 0.40 compared to ground observations and systematically overestimates (69%) in summer and underestimates (-19%) in winter. Through a series of sensitivity simulations, we reveal that modeling bias is primarily related to the dominant fraction of monoterpene-formed secondary organic aerosol (SOA), and a strong positive correlation of 0.67 is found between monoterpene emission and modeling bias in the eastern US during summer. In terms of vertical profile, the model prominently underestimates OA and monoterpene concentrations by 37%–99% and 82%–99%, respectively, in the upper air (>500 m) as validated against aircraft observations. Our study suggests that the current volatility basis set (VBS) scheme applied in CESM might be parameterized with monoterpene SOA yields that are too high, which subsequently results in strong SOA production near the emission source area. We also find that the model has difficulty in reproducing the decreasing trend of surface OA in the southeastern US probably because of employing pure gas VBS to represent isoprene SOA which is in reality mainly formed through multiphase chemistry; thus, the influence of aerosol acidity and sulfate particle change on isoprene SOA formation has not been fully considered in the model. This study reveals the urgent need to improve the SOA modeling in climate models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗

Performance and Accuracy Implications of Parallel Split Physics-Dynamics Coupling in the Energy Exascale Earth System Atmosphere Model

Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.

97 MATHEMATICS AND COMPUTING↗

The Tropospheric Emission Spectrometer: From Discovery Mission to Earth System Sounding

The Tropospheric Emission Spectrometer (TES) was launched aboard the Aura spacecraft in 2004 and decommissioned in February 2018. Over more than a decade of observations, TES has shown the potential of high-resolution Fourier Transform Spectrometer (FTS) measurements to alter our understanding of atmospheric composition and its relationship to the Earth System. We survey some of the major contributions TES has made to understanding air quality, chemistry-climate interactions, atmospheric hydrology, and carbon cycle science. We then show how innovations in retrieval methods pioneered by the TES team are now being applied to multi-spectral, multi-instrument records, establishing a new paradigm for composition measurements into the future.

Bowman, Kevin W.↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

Benthic Kinetic Energy in an Earth System Model and Observations

The benthic layer of the ocean is closely correlated with every upper layer above it. Close analysis of benthic layer stress can yield an accurate depiction of the state of the ocean, as well as the atmosphere above it. Through the use of the Department of Energy’s Energy Exascale Earth System Model (E3SM), we took a look at how these features changed over time in a simulation spanning the historical period from 1958 to 2018 in the Arctic Ocean north of 60° N. Making use of E3SM’s regional refinement capability to place higher resolution in regions of interest, the E3SM-Arctic configuration has ~10 km resolution in the ocean and sea ice components in this region, coarsening up to ~60 km outside of the Arctic (Veneziani et al. 2022). In particular, we narrowly focused on the benthic kinetic energy levels, and how trends changed over time. To validate the fidelity of the E3SM model, we compared model current speed outputs over a 369-day period in the Bering Strait to actual recorded observations during this time period. Our analysis of the benthic kinetic energy mean, standard deviation, and trend throughout this time period can be used to characterize the effects of environmental impacts in the Arctic Ocean.

54 ENVIRONMENTAL SCIENCES↗

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 1: Experimental protocols and surface changes

There is now substantial literature on climate model studies of equatorial or tropical stratospheric SO 2 injections that aim to counteract the surface warming produced by rising concentrations of greenhouse gases. Here we present the results from the first systematic intercomparison of climate responses in three Earth system models wherein the injection of SO 2 occurs at different latitudes in the lower stratosphere: CESM2-WACCM6, UKESM1.0 and GISS-E2.1-G. The first two use a modal aerosol microphysics scheme, while two versions of GISS-E2.1-G use a bulk aerosol (One-Moment Aerosol, OMA) and a two-moment (Multiconfiguration Aerosol TRacker of mIXing state, MATRIX) microphysics approach, respectively. Our aim in this work is to determine commonalities and differences between the climate model responses in terms of the distribution of the optically reflective sulfate aerosols produced from the oxidation of SO 2 and in terms of the surface response to the resulting reduction in solar radiation. A focus on understanding the contribution of characteristics of models transport alongside their microphysical and chemical schemes, and on evaluating the resulting stratospheric responses in different models, is given in the companion paper (Bednarz et al., 2023). The goal of this exercise is not to evaluate these single-point injection simulations as stand-alone proposed strategies to counteract global warming; instead we determine sources and areas of agreement and uncertainty in the simulated responses and, ultimately, the possibility of designing a comprehensive intervention strategy capable of managing multiple simultaneous climate goals through the combination of different injection locations. We find large disagreements between GISS-E2.1-G and the CESM2-WACCM6 and UKESM1.0 models regarding the magnitude of cooling per unit of aerosol optical depth (AOD) produced, which varies from 4.7 K per unit of AOD in CESM2-WACCM6 to 16.7 K in the GISS-E2.1-G version with two-moment aerosol microphysics. By normalizing the results with the global mean response in each of the models and thus assuming that the amount of SO 2 injected is a free parameter that can be managed independently, we highlight some commonalities in the overall distributions of the aerosols, in the inter-hemispheric surface temperature response and in shifts to the Intertropical Convergence Zone, as well as some areas of disagreement, such as the extent of the aerosol confinement in the equatorial region and the efficiency of the transport to polar latitudes. In conclusion, we demonstrate that it is possible to use these simulations to produce more comprehensive injection strategies in multiple climate models. However, large differences in the injection magnitudes can be expected, potentially increasing inter-model spreads in some stratospheric quantities (such as aerosol distribution) while reducing the spread in the surface response in terms of temperature and precipitation; furthermore, the selection of the injection locations may be dependent on the models' specific stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗

Automation is all you need: Faster Earth system models with AI/ML

Focal Area: Data acquisition and assimilation enabled by machine learning (ML), artificial intelligence (AI) and advanced methods. Science Challenge: Tropical cyclones can in- duce extreme water cycle events through dramatic precipitation and storm surge. More reliable models of intensity will translate into better prediction of the impact of extreme events in large scale Earth systems simulations. We demonstrate and describe AI/ML methodologies for rapid assimilation of new, in situ data products.

54 ENVIRONMENTAL SCIENCES↗

Diurnal rainfall response to the physiological and radiative effects of CO 2 in tropical forests in the Energy exascale earth system model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO$_2$ concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO$_2$. The plant physiological and radiative impacts of CO$_2$ on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO$_2$ levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO$_2$ radiative and CO$_2$ plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime CAPE (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

54 ENVIRONMENTAL SCIENCES↗

Bridging New Observational Capabilities and Process-Level Simulation: Insights into Aerosol Roles in the Earth System

The spatial distribution of ambient aerosol particles significantly impacts aerosol–radiation–cloud interactions, which contribute to the largest uncertainty in global anthropogenic radiative forcing estimations. However, the atmospheric boundary layer and lower free troposphere have not been adequately sampled in terms of spatiotemporal resolution, hindering a comprehensive characterization of various atmospheric processes and impeding our understanding of the Earth system. To address this research data gap, we have leveraged the development of uncrewed aerial systems (UAS) and advanced measurement techniques to obtain mesoscale spatial data on aerosol microphysical and optical properties around the U.S. Southern Great Plains (SGP) atmospheric observatory. Our study also benefits from state-of-the-art laboratory facilities that include three-dimensional molecular imaging techniques enabled by secondary ion mass spectrometry and nanogram-level chemical composition analysis via micronebulization aerosol mass spectrometry. Through our study, we have developed a framework for observation–modeling integration, enabling an examination of how various assumptions about the organic–inorganic components mixing state, inferred from chemical analysis, affect clouds and radiation in observation-constrained model simulations. By integrating observational constraints (derived from offline chemical analysis of the aerosol surface using collected samples) with in situ UAS observations, we have identified a prominent role of organic-enriched nanometer layers located at the surface of aerosol particles in determining profiles of aerosol optical and hygroscopic properties over the SGP observatory. Furthermore, we have improved the agreement between predicted clouds and ground-based cloud lidar measurements. This UAS–model–laboratory integration exemplifies how these new advanced capabilities can significantly enhance our understanding of aerosol–radiation–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗

On the Generation of Weddell Sea Polynyas in a High-Resolution Earth System Model

Larger Weddell Sea polynyas (WSPs), differentiated in this study from the smaller Maud Rise Polynyas (MRPs), forming to the east of the prime meridian in the proximity of the Maud Rise seamount, have last been observed in the 1970s. We investigate WSPs that grow realistically out of MRPs in a high-resolution (HR) preindustrial simulation with the Energy Exascale Earth System Model version 0.1. The formation of MRPs requires HR to simulate the detailed flow around Maud Rise, while the realistic formation of WSPs requires a model to produce MRPs. Furthermore, WSPs tend to follow periods of a prolonged build-up of a heat reservoir at depth and weakly negative wind-stress curl in association with the core of the southern hemisphere westerlies at an anomalously northern position. While this scenario also leads to drier conditions over the central Weddell Sea, which some literature claims to be a necessary condition for the formation of WSPs, our model results indicate that open-ocean polynyas do not occur during periods of weakly negative wind-stress curl despite drier atmospheric conditions. Our study supports the hypothesis noted in earlier studies that a shift from a weakly negative to a strongly negative wind-stress curl over the Weddell Sea is a prerequisite for WSPs to form, together with a large heat reservoir at depth. However, the ultimate trigger is a pronounced MRP; whose associated convection creates high surface salinity anomalies that propagate westward with the flow of the Weddell Gyre. If large enough, these anomalies trigger the formation of a WSP and a pulse of newly formed Antarctic Bottom Water.

54 ENVIRONMENTAL SCIENCES↗

Reduced ENSO Variability at the LGM Revealed by an Isotope-Enabled Earth System Model

Studying the El Nino Southern Oscillation (ENSO) in the past can help us better understand its dynamics and improve its future projections. However, both paleoclimate reconstructions and model simulations of ENSO strength at the Last Glacial Maximum (LGM; 21 ka B.P.) have led to contradicting results. Here we perform model simulations using the recently developed water isotope-enabled Community Earth System Model (iCESM). For the first time, model-simulated oxygen isotopes are directly compared with those from ENSO reconstructions using the individual foraminifera analysis (IFA). We find that the LGM ENSO is most likely weaker comparing with the preindustrial. The iCESM suggests that total variance of the IFA records may only reflect changes in the annual cycle instead of ENSO variability as previously assumed. Furthermore, the interpretation of subsurface IFA records can be substantially complicated by the habitat depth of thermocline-dwelling foraminifera and their vertical migration with a temporally varying thermocline.

El Niño–Southern Oscillation;isotope-enabled si↗

Future Changes in Active and Inactive Atlantic Hurricane Seasons in the Energy Exascale Earth System Model

Abstract North Atlantic tropical cyclones (TCs) have considerable interannual variability, with La Niña and the positive phase of the Atlantic Meridional Mode (AMM) tending to drive active hurricane seasons, and El Niño and the negative AMM often driving inactive seasons. Here, we analyze how active and inactive Atlantic hurricane seasons may change in the future using the high resolution Energy Exascale Earth System Model (E3SM). We performed atmosphere‐only simulations forced by sea‐surface temperature patterns characteristic of La Niña and the positive AMM jointly, and El Niño and the negative AMM jointly, in historical and future climates. Projected Atlantic TCs become more frequent in the future by approximately 34% during El Niño and negative AMM and by 66% during La Niña and positive AMM, with a significant increase in the portion of intense TCs. Warmer SSTs increase TC potential intensity, with reduced wind shear and increased mid‐tropospheric humidity further supporting TC activity.

54 ENVIRONMENTAL SCIENCES↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C. Andrews↗

Representing Global Soil Erosion and Sediment Flux in Earth System Models

Abstract Soil erosion produces enormous amounts of sediment, carbon, and nutrient fluxes from land to rivers, thus playing crucial roles in global biogeochemical cycles and food security. To predict soil erosion in the context of climate and land use changes, we explicitly parameterize cropland management actions (i.e., irrigation, conserved agriculture, and crop residue management) and geological factors (i.e., lithology and glacier) in the Energy Exascale Earth System Model (E3SM) soil erosion module. The erosion model is calibrated using a global‐scale regionalized parameter calibration method. The spatial variabilities of the modeled and the Revised Universal Soil Loss Equation (RUSLE) based benchmark soil erosion are consistent across vegetation, climate and soil properties. Compared with independent data, our model shows a bias reduction in 59% of the observations relative to the RUSLE‐based soil erosion, with 53% of the bias reduction exceeding 50%. This improvement is mainly due to a better representation of the topographic effect on soil erosion. Our results indicate that conserved agriculture practices have effectively reduced soil erosion in cropland by over 25% in the United States and Argentina. In contrast, irrigation has increased soil erosion in many Asian countries. For upland sediment flux, our model is consistent with the WBMsed benchmark data in inter‐basin variability but could be more skillful in simulating intra‐basin variability because it couples soil erosion and sediment flux explicitly. The developed model provides useful skills for more realistic predictions of soil erosion and river sediment dynamics under environmental changes.

58 GEOSCIENCES↗