Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vegetation dynamics modeling”

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 253 records · Page 14

Global rules for translating land-use change (LUH2) to land-cover change for CMIP6 using GLM2

Anthropogenic land-use and land-cover change activities play a critical role in Earth system dynamics through significant alterations to biogeophysical and biogeochemical properties at local to global scales. To quantify the magnitude of these impacts, climate models need consistent land-cover change time series at a global scale, based on land-use information from observations or dedicated land-use change models. However, a specific land-use change cannot be unambiguously mapped to a specific land-cover change. Here, nine translation rules are evaluated based on assumptions about the way land-use change could potentially impact land cover. Utilizing the Global Land-use Model 2 (GLM2), the model underlying the latest Land-Use Harmonization dataset (LUH2), the land-cover dynamics resulting from land-use change were simulated based on multiple alternative translation rules from 850 to 2015 globally. For each rule, the resulting forest cover, carbon density and carbon emissions were compared with independent estimates from remote sensing observations, U.N. Food and Agricultural Organization reports, and other studies. The translation rule previously suggested by the authors of the HYDE 3.2 dataset, that underlies LUH2, is consistent with the results of our examinations at global, country and grid scales. This rule recommends that for CMIP6 simulations, models should (1) completely clear vegetation in land-use changes from primary and secondary land (including both forested and non-forested) to cropland, urban land and managed pasture; (2) completely clear vegetation in land-use changes from primary forest and/or secondary forest to rangeland; (3) keep vegetation inland-use changes from primary non-forest and/or secondary non-forest to rangeland. Our analysis shows that this rule is one of three (out of nine)rules that produce comparable estimates of forest cover, vegetation carbon and emissions to independent estimates and also mitigate the anomalously high carbon emissions from land-use change observed in previous studies in the 1950s. According to the three translation rules, contemporary global forest area is estimated to be 37.42×106 km2, within the range derived from remote sensing products. Likewise, the estimated carbon stock is in close agreement with reference biomass datasets, particularly over regions with more than 50 % forest cover

58 GEOSCIENCES↗

A deep learning hybrid predictive modeling (HPM) approach for estimating evapotranspiration and ecosystem respiration

Climate change is reshaping vulnerable ecosystems, leading to uncertain effects on ecosystem dynamics, including evapotranspiration (ET) and ecosystem respiration (R eco ). However, accurate estimation of ET and R eco still remains challenging at sparsely monitored watersheds, where data and field instrumentation are limited. In this study, we developed a hybrid predictive modeling approach (HPM) that integrates eddy covariance measurements, physically based model simulation results, meteorological forcings, and remote-sensing datasets to estimate ET and R eco in high space–time resolution. HPM relies on a deep learning algorithm and long short-term memory (LSTM) and requires only air temperature, precipitation, radiation, normalized difference vegetation index (NDVI), and soil temperature (when available) as input variables. We tested and validated HPM estimation results in different climate regions and developed four use cases to demonstrate the applicability and variability of HPM at various FLUXNET sites and Rocky Mountain SNOTEL sites in Western North America. To test the limitations and performance of the HPM approach in mountainous watersheds, an expanded use case focused on the East River Watershed, Colorado, USA. The results indicate HPM is capable of identifying complicated interactions among meteorological forcings, ET, and R eco variables, as well as providing reliable estimation of ET and R eco across relevant spatiotemporal scales, even in challenging mountainous systems. The study documents that HPM increases our capability to estimate ET and R eco and enhances process understanding at sparsely monitored watersheds.

54 ENVIRONMENTAL SCIENCES↗

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Fine-scale vegetation composition and structure shape spatiotemporal variation in surface albedo across a low Arctic tundra landscape

The unprecedented rate of warming in the Arctic is driving changes in the structure and composition of tundra vegetation. Increases in deciduous tall shrub cover, height, and density are of particular concern, as these changes alter local surface albedo in ways that could amplify effects on the regional surface energy budget (SEB). Despite this importance, significant uncertainties remain in understanding the interplay between fine-scale vegetation patterns and emergent albedo dynamics across space and time. Here, we address these uncertainties by (1) quantifying spatiotemporal variation in surface shortwave albedo and (2) determining the relative influence of fine-scale vegetation composition, structure, and environmental conditions on albedo across a representative low-Arctic tundra landscape on Alaska’s Seward Peninsula. To do this, we synthesized multi-scale, multi-platform remote sensing observations, including a novel Landsat-derived albedo time series, a fine-scale map of Arctic plant functional type (PFT) fractional cover, and airborne LiDAR estimates of canopy height and topography. We show that there are substantial reductions in winter albedo for pixels dominated by tall, woody PFTs (28.13%) relative to pixels dominated by non-woody vegetation, but almost no change in summer albedo (3% increase). Further, we identified a unimodal trend in the relationship between canopy height and the timing of the springtime transition from high (snowy) to low (leafy) albedo (peak at 5.5 m), possibly because of competing ‘snow-fence’ and ‘protrusion’ snow-shrub interactions. To explore the primary drivers of albedo, we constructed a random forest model and found that canopy height and the fractional cover of woody PFTs were as- or more important predictors of winter albedo than topographic features. These findings provide strong evidence for the impacts of local vegetation characteristics on regional surface albedo, highlighting the need for better quantification of snow-shrub interactions to accurately predict the Arctic’s SEB under future environmental change.

Arctic↗

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

Modeling Snow Dynamics and Stable Water Isotopes Across Mountain Landscapes

Abstract A coupled hydrologic and snowpack stable water isotope model assesses controls on isotopic inputs across a mountainous basin. Annually, the most depleted isotope conditions occur in the upper subalpine where snow accumulation is high, and rainfall is low. Snowmelt isotopic evolution over time indicates fractionation processes account for <25% snowmelt enrichment. Meltwater isotopic inputs are largely determined by controls on the amount, phase and isotopic mass of precipitation coincident with the ablation period. Effect of vapor loss from the snowpack on d‐excess in snowmelt is a balance between energy and snow‐availability. It is highest above treeline, and in the grass and aspen‐dominated portions of the upper montane where vegetation shading is low. Deep snowpack in conifer forests limit the influence of vapor loss in the subalpine. Wet years reduce the effects of vapor loss on snowmelt across the basin, except in the lower montane where added snowfall bolsters snow‐limited conditions.

58 GEOSCIENCES↗

Recent Widespread Deceleration of Global Surface Urban Heat Islands Unveiled by Satellites

Tracking the temporal dynamics of urban heat island (UHI) is critical for urban heat adaptation and mitigation strategies. However, whether UHI trends have shifted recently and their underlying drivers remain unknown. Here we investigate the variabilities in surface UHI trends and their associated determinants in 2,104 cities worldwide from 2000 to 2022. Our findings reveal that approximately half of the world's cities have experienced notable shifts in surface UHI trends, predominantly characterized by UHI deceleration. These shifts can be primarily attributed to alterations in vegetation trends during the day and to modifications in surface albedo and local warming trends at night. Our study challenges the conventional linear models commonly employed to estimate surface UHI trends, suggesting potential biases in such estimates. Our findings underscore the need for nuanced policies to curtail UHI growth by considering changes in urban underlying surfaces and background climate, particularly from a nonlinear perspective.

54 ENVIRONMENTAL SCIENCES↗

Comparing deuterium excess to large-scale precipitation recycling models in the tropics

Precipitation recycling is essential to sustaining regional ecosystems and water supplies, and it is impacted by land development and climate change. This is especially true in the tropics, where dense vegetation greatly influences recycling. Unfortunately, large-scale models of recycling often exhibit high uncertainty, complicating efforts to estimate recycling. Here, we examine how deuterium excess (d-excess), a stable-isotope quantity sensitive to recycling effects, acts as an observational proxy for recycling. While past studies have connected variability in d-excess to precipitation origins at local or regional scales, our study leverages >3000 precipitation isotope samples to quantitatively compare d-excess against three contemporary recycling models across the global tropics. Using rank-correlation, we find statistically significant agreement ( $\bar \tau = 0.52$ τ ¯ = 0.52 to $0.70$ 0.70 ) between tropical d-excess and recycling that is strongly mediated by seasonal precipitation, vegetation density, and scale mismatch. Our results detail the complex relationship between d-excess and precipitation recycling, suggesting avenues for further investigation.

54 ENVIRONMENTAL SCIENCES↗

Detecting forest response to droughts with global observations of vegetation water content

Droughts in a warming climate have become more common and more extreme, making understanding forest responses to water stress increasingly pressing. Analysis of water stress in trees has long focused on water potential in xylem and leaves, which influences stomatal closure and water flow through the soil-plant-atmosphere continuum. At the same time, changes of vegetation water content (VWC) are linked to a range of tree responses, including fluxes of water and carbon, mortality, flammability, and more. Unlike water potential, which requires demanding in situ measurements, VWC can be retrieved from remote sensing measurements, particularly at microwave frequencies using radar and radiometry. Here, we highlight key frontiers through which VWC has the potential to significantly increase our understanding of forest responses to water stress. To validate remote sensing observations of VWC at landscape scale and to better relate them to data assimilation model parameters, we introduce an ecosystem-scale analogue of the pressure-volume curve, the non-linear relationship between 44 average leaf or branch water potential and water content commonly used in plant hydraulics. The 45 sources of variability in these ecosystem-scale pressure-volume curves and their relationship to 46 forest response to water stress are discussed. We further show to what extent diel, seasonal, and 47 decadal dynamics of VWC reflect variations in different processes relating the tree response to 48 water stress. VWC can also be used for inferring belowground conditions – which are difficult to 49 impossible to observe directly. Lastly, we discuss how a dedicated geostationary spaceborne 50 observational system for VWC, when combined with existing datasets, can capture diel and 51 seasonal water dynamics to advance the science and applications of global forest vulnerability to 52 future droughts.

54 ENVIRONMENTAL SCIENCES↗

Convolutional Non-Homogeneous Poisson Process and its Application to Wildfire Ignition Risk Quantification for Power Delivery Networks

To quantify wildfire ignition risks on power delivery networks, the current practice predominantly relies on the empirically calculated fire danger indices, which may not well capture the effects of dynamically changing environmental factors. This article proposes a spatio-temporal point process model, known as the Convolutional Non-homogeneous Poisson Process (cNHPP), and applies the model to quantify wildfire ignition risks for power delivery networks. The proposed model captures both the current (i.e., instantaneous) and cumulative (i.e., historical) effects of key environmental processes (i.e., covariates) on wildfire risks, as well as the spatio-temporal dependency among different segments of the power delivery network. The computation and interpretation of the intensity function are thoroughly investigated. We apply the proposed approach to estimate wildfire ignition risks on major transmission lines in California, using historical fire data, meteorological and vegetation data obtained from the National Oceanic and Atmospheric Administration and National Aeronautics and Space Administration. Here, a comprehensive comparison study is performed to show the applicability and predictive capability of the proposed approach.

Non-homogeneous Poisson Process↗

Rising Water Levels and Vegetation Shifts Drive Substantial Reductions in Methane Emissions and Carbon Dioxide Uptake in a Great Lakes Coastal Freshwater Wetland

ABSTRACT Coastal freshwater wetlands are critical ecosystems for both local and global carbon cycles, sequestering substantial carbon while also emitting methane (CH 4 ) due to anoxic conditions. Estuarine freshwater wetlands face unique challenges from fluctuating water levels, which influence water quality, vegetation, and carbon cycling. However, the response of CH 4 fluxes and their drivers to altered hydrology and vegetation remains unclear, hindering mechanistic modeling. To address these knowledge gaps, we studied an estuarine freshwater wetland in the Great Lakes region, where rising water levels led to a vegetation shift from emergent Typha dominance in 2015–2016 to floating‐leaved species in 2020–2022. Using eddy covariance flux measurements during the peak growing season (June–September) of both periods, we observed a 60% decrease in CH 4 emissions, from 81 ± 4 g C m −2 in 2015–2016 to 31 ± 3 g C m −2 in 2020–2022. This decline was driven by two main factors: (1) higher water levels, which suppressed ebullitive fluxes via increased hydrostatic pressure and extended CH 4 residence time, enhancing oxidation potential in the water column; and (2) reduced CH 4 conductance through plants. Net carbon dioxide (CO 2 ) uptake decreased by 90%, from −267 ± 26 g C m −2 in 2015–2016 to −27 ± 49 g C m −2 in 2020–2022. Additionally, diel CH 4 flux patterns shifted, with a distinct morning peak observed in 2015–2016 but absent in 2020–2022, suggesting changes in plant‐mediated transport and a potential decoupling from photosynthesis. The dominant factors influencing CH 4 fluxes shifted from water temperature and gross primary productivity in 2015–2016 to atmospheric pressure in 2020–2022, suggesting an increased role of ebullition as a primary transport pathway. Our results demonstrate that changes in water levels and vegetation can substantially alter CH 4 and CO 2 fluxes in coastal freshwater wetlands, underscoring the critical role of hydrological shifts in driving carbon dynamics in these ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Declining groundwater storage expected to amplify mountain streamflow reductions in a warmer world

Abstract Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Carroll, Rosemary W. H. (ORCID:0000000293028074)↗

Packaging “vegetable oils”: Insights into plant lipid droplet proteins

Abstract Plant neutral lipids, also known as “vegetable oils”, are synthesized within the endoplasmic reticulum (ER) membrane and packaged into subcellular compartments called lipid droplets (LDs) for stable storage in the cytoplasm. The biogenesis, modulation, and degradation of cytoplasmic LDs in plant cells are orchestrated by a variety of proteins localized to the ER, LDs, and peroxisomes. Recent studies of these LD-related proteins have greatly advanced our understanding of LDs not only as steady oil depots in seeds but also as dynamic cell organelles involved in numerous physiological processes in different tissues and developmental stages of plants. In the past 2 decades, technology advances in proteomics, transcriptomics, genome sequencing, cellular imaging and protein structural modeling have markedly expanded the inventory of LD-related proteins, provided unprecedented structural and functional insights into the protein machinery modulating LDs in plant cells, and shed new light on the functions of LDs in nonseed plant tissues as well as in unicellular algae. Here, we review critical advances in revealing new LD proteins in various plant tissues, point out structural and mechanistic insights into key proteins in LD biogenesis and dynamic modulation, and discuss future perspectives on bridging our knowledge gaps in plant LD biology.

Cai, Yingqi (ORCID:0000000203575809)↗

Whole-Ecosystem Warming Increases Plant-Available Nitrogen and Phosphorus in an Ombrotrophic Bog

Warming is expected to increase the net release of carbon from peatland soils, contributing to future warming. This positive feedback may be moderated by the response of peatland vegetation to rising atmospheric [CO 2 ] or to increased soil nutrient availability. We asked whether a gradient of whole-ecosystem warming (from + 0 °C to + 9 °C) would increase plant-available nitrogen and phosphorus in an ombrotrophic bog in northern Minnesota, USA, and whether elevated [CO 2 ] would modify the nutrient response. We tracked changes in plant-available nutrients across space and through time and in comparison with other nutrient pools, and assessed whether nutrient warming responses were captured by a point version of the land-surface model, ELM-SPRUCE. We found that warming exponentially increased plant-available ammonium and phosphate, but that nutrient dynamics were unaffected by elevated [CO 2 ]. The warming response increased by an order of magnitude between the first and fourth year of the experimental manipulation, perhaps because of dramatic mortality of Sphagnum mosses in the surface peat of the warmest treatments. However, neither the magnitude nor the temporal dynamics of the responses were captured by ELM-SPRUCE. Here, relative increases in plant-available ammonium and phosphate with warming were similar, but the response varied across raised hummocks and depressed hollows and with peat depth. Plant-available nutrient dynamics were only loosely correlated with inorganic and organic porewater nutrients, likely representing different processes. Future predictions of peatland nutrient availability under climate change scenarios must account for dynamic changes in nutrient acquisition by plants and microbes, as well as microtopography and peat depth.

54 ENVIRONMENTAL SCIENCES↗

Boreal–Arctic wetland methane emissions modulated by warming and vegetation activity

Abstract Wetland methane (CH 4 ) emissions over the Boreal–Arctic region are vulnerable to climate change and linked to climate feedbacks, yet understanding of their long-term dynamics remains uncertain. Here, we upscaled and analysed two decades (2002–2021) of Boreal–Arctic wetland CH 4 emissions, representing an unprecedented compilation of eddy covariance and chamber observations. We found a robust increasing trend of CH 4 emissions (+8.9%) with strong inter-annual variability. The majority of emission increases occurred in early summer (June and July) and were mainly driven by warming (52.3%) and ecosystem productivity (40.7%). Moreover, a 2 °C temperature anomaly in 2016 led to the highest recorded annual CH 4 emissions (22.3 Tg CH 4 yr −1 ) over this region, driven primarily by high emissions over Western Siberian lowlands. However, current-generation models from the Global Carbon Project failed to capture the emission magnitude and trend, and may bias the estimates in future wetland CH 4 emission driven by amplified Boreal–Arctic warming and greening.

54 ENVIRONMENTAL SCIENCES↗

Detection and attribution of long-term and fine-scale changes in spring phenology over urban areas: A case study in New York State

Spring phenology plays an essential role in climate change, terrestrial ecosystem, and public health. Field-based monitoring and understanding of changes in spring phenology for long periods and in large regions are challenging due to the limited in-site observations. Space-based remotely sensed observations offer great potentials for monitoring decadal spring phenology changes from regional to global scales. However, the coarse-scale remotely sensed observations are insufficient to capture fine-scale spring phenology dynamics, especially in urban areas, and this makes it challenging for understanding the combined effects of climate change and urbanization on spring phenology. We derived the start of phenology season (SOS) in New York State using 30 m Landsat observations from 1990 to 2015 to understand the impact of the environment and urbanization on SOS. The results show that SOS for different years reveals heterogeneous spatial distribution. Most regions of New York State have been experiencing significant spring phenology changes in form of earlier onset of vegetation greening, ranging from 0.2 to 0.6 day/year during 1990 to 2015, and this trend varies slightly with latitudes and urbanization levels. Further, spatial correlation analysis shows that the increase in temperature and urbanization could both promote the advancement of SOS. However, the effect of urbanization (partial correlation coefficient (R) ranges from −0.289 to −0.542) on SOS is greater than the effect of temperature (R ranges from 0.006 to −0.192). The study generates a high spatio-temporal resolution spring phenology dataset for ecological, environmental and public health studies, especially in urban areas, and reveals the importance of better accounting for the urbanization effects when quantifying the SOS dynamics in phenology models.

Landsat↗

End of Green Sahara amplified mid- to late Holocene megadroughts in mainland Southeast Asia

Between 5 and 4 thousand years ago, crippling megadroughts led to the disruption of ancient civilizations across parts of Africa and Asia, yet the extent of these climate extremes in mainland Southeast Asia (MSEA) has never been defined. This is despite archeological evidence showing a shift in human settlement patterns across the region during this period. We report evidence from stalagmite climate records indicating a major decrease of monsoon rainfall in MSEA during the mid- to late Holocene, coincident with African monsoon failure during the end of the Green Sahara. Through a set of modeling experiments, we show that reduced vegetation and increased dust loads during the Green Sahara termination shifted the Walker circulation eastward and cooled the Indian Ocean, causing a reduction in monsoon rainfall in MSEA. Our results indicate that vegetation-dust climate feedbacks from Sahara drying may have been the catalyst for societal shifts in MSEA via ocean-atmospheric teleconnections.

54 ENVIRONMENTAL SCIENCES↗

Model code and data: biomass allocation adjustments induced by elevated CO2 and warming in a C3 brackish marsh, 2017-2022, Maryland

This dataset and R script accompany the published paper Bruns et al. (2024) in Geophysical Research Letters. The data are from the first six years of a field manipulation of whole-ecosystem warming and elevated CO2 experiment (Salt Marsh Accretion Response to Temperature eXperiment, or SMARTX) in the Smithsonian's Global Change Research Wetland (GCReW), a brackish, microtidal wetland site on a subestuary of the Chesapeake Bay. These data were generated to understand how warming and elevated CO2 interact to structure ecosystem-level responses to global change, particularly in terms of carbon sequestration. The dataset covers 2017-2022 and includes peak annual above ground biomass, annual belowground fine root productivity, and porewater NH4 for each experimental plot. The overall experiment is replicated in two locations on the marsh, a lower elevation zone dominated the C3 sedge S. Americanus and a higher elevation plot dominated by the C4 species. This paper and its data release is only for the C3 plot. Variable descriptions for data file is available in variable_descriptions.pdf. The R script Bruns_et_al_2024_GRL_make_figures.Rmd contains model code and other scripts used to generate all paper figures.

54 ENVIRONMENTAL SCIENCES↗