Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tropical variability”

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 793 records · Page 44

MetaFlux: Meta-learning global carbon fluxes from sparse spatiotemporal observations

We provide a global, long-term carbon flux dataset of gross primary production and ecosystem respiration generated using meta-learning, called MetaFlux. The idea behind meta-learning stems from the need to learn efficiently given sparse data by learning how to learn broad features across tasks to better infer other poorly sampled ones. Using meta-trained ensemble of deep models, we generate global carbon products on daily and monthly timescales at a 0.25-degree spatial resolution from 2001 to 2021, through a combination of reanalysis and remote-sensing products. Site-level validation finds that MetaFlux ensembles have lower validation error by 5–7% compared to their non-meta-trained counterparts. In addition, they are more robust to extreme observations, with 4–24% lower errors. We also checked for seasonality, interannual variability, and correlation to solar-induced fluorescence of the upscaled product and found that MetaFlux outperformed other machine-learning based carbon product, especially in the tropics and semi-arids by 10–40%. Overall, MetaFlux can be used to study a wide range of biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗

Global root traits (GRooT) database

Motivation: Trait data are fundamental to the quantitative description of plant form and function. Although root traits capture key dimensions related to plant responses to changing environmental conditions and effects on ecosystem processes, they have rarely been included in large-scale comparative studies and global models. For instance, root traits remain absent from nearly all studies that define the global spectrum of plant form and function. Thus, to overcome conceptual and methodological roadblocks preventing a widespread integration of root trait data into large-scale analyses we created the Global Root Trait (GRooT) Database. GRooT provides readyto- use data by combining the expertise of root ecologists with data mobilization and curation. Specifically, we (a) determined a set of core root traits relevant to the description of plant form and function based on an assessment by experts, (b) maximized species coverage through data standardization within and among traits, and (c) implemented data quality checks. Main types of variables contained: GRooT contains 114,222 trait records on 38 continuous root traits. Spatial location and grain: Global coverage with data from arid, continental, polar, temperate and tropical biomes. Data on root traits were derived from experimental studies and field studies. Time period and grain: Data were recorded between 1911 and 2019. Major taxa and level of measurement: GRooT includes root trait data for which taxonomic information is available. Trait records vary in their taxonomic resolution, with subspecies or varieties being the highest and genera the lowest taxonomic resolution available. It contains information for 184 subspecies or varieties, 6,214 species, 1,967 genera and 254 families. Owing to variation in data sources, trait records in the database include both individual observations and mean values. Software format: GRooT includes two csv files. A GitHub repository contains the csv files and a script in R to query the database.

59 BASIC BIOLOGICAL SCIENCES↗

Assessing Tropical Pacific-induced Predictability of Southern California Precipitation Using a Novel Multi-input Multi-output Autoencoder

We construct a novel Multi-Input Multi-Output Autoencoder-decoder (MIMO-AE) to capture the non-linear relationship of Southern California precipitation and tropical Pacific Ocean sea surface temperature. The MIMO-AE is trained on both monthly TP-SST and SC-PRECIP anomalies simultaneously. The co-variability of the two fields in the MIMO-AE shared nonlinear latent space can be condensed into an index, termed the MIMO-AE index. We use a transfer learning approach to train a MIMO-AE on the combined dataset of 100 years of output from a historical simulation with the Energy Exascale Earth Systems Model version 1 and a segment of observational data. We further use Long Short-Term Memory networks to assess sub-seasonal predictability of SC-PRECIP using the MIMO-AE index. We find that the MIMO-AE index provides enhanced predictability of SC-PRECIP for a lead-time of up-to four months as compared to Niño 3.4 index and the El Niño Southern Oscillation Longitudinal Index.

54 ENVIRONMENTAL SCIENCES↗

Changes in stratospheric ozone and temperature due to the eruptions of Mt. Pinatubo

The paper investigates the impact of the Mt. Pinatubo eruptions on the total column ozone measured from the Nimbus 8 TOMS and the NOAA-11 SBUV/2 spectrometers. Within a few months after the eruptions, the total column ozone decreased by 5-6 percent in the tropics, 3-4 percent at midlatitudes, and 6-9 percent at high latitudes in the Northern Hemisphere. After the effects of quasi-biennial oscillations and interannual variability are taken into account, the decrease in the column ozone attributed to volcanic eruptions at these latitudes may not be more than 2-4 percent. The most noticeable effect on the Pinatubo eruptions, as observed during the El Chichon period, is the breakdown of the phase relation between ozone and temperature. This is attributed to additional heating in the lower stratosphere caused by volcanic aerosols.

Chandra, S.↗

A New ML-Based Adaptive Thinning Methodology to Improve the Impact of AIRS and CrIS Assimilation on Global Tropical Cyclone Forecasts

This work builds on previous research performed by this team to improve the forecast of Tropical Cyclones (TCs) by assimilating AIRS and CrIS radiances into the NASA Global Earth Observing System (GEOS). Past published work demonstrated that the assimilation of radiances with variable density was beneficial to TC forecasting in the GEOS. In the previous setup, a fixed-size moving square named 'TC domain' was activated by the so-called TC-vitals, an international real-time message accessible to all NWP forecasting centers, that documents the existence of a TC, its estimated position, and its size. The information from TC-vitals activated a switch in the GEOS, which allowed to reduce the distance used for thinning AIRS and CrIS data inside a 15 degrees by 15 degrees moving TC domain centered on the storm, so that more data were assimilated in the vicinity of the TC during its lifetime. The methodology produced improved TC analyses and led to better forecasts, particularly related to intensity, without damaging the global forecast skill. In the new version, the adaptive thinning methodology is based on a machine-learning technique. The technique searches for TCs and creates TC masks by using cloud-top temperatures from all geostationary satellites without the need for additional information. It is being trained against the International Best Track Archive for Climate Stewardship (IBTrACS) data base. Once a TC mask is created, a switch identical to the one used in the previous adaptive thinning method is activated, allowing the GEOS to ingest more data in the TC-shaped size-changing domain that follows the storm. As of today, the team has been able to successfully assimilate data inside the ML-detected TC domains. Future work includes an improved capability of reducing false alarm rates (i.e., cloud systems that are erroneously labeled as TCs).

Oreste Reale↗

Leaf-cutting ant (Atta cephalotes) nests may be hotspots of methane and carbon dioxide emissions in tropical forests

Leaf-cutting ants of the genus Atta are widely distributed throughout the American tropics and subtropics and rival other herbivores in the consumption of surrounding foliage. Although numerous studies have been conducted on the role these insects play in herbivory and organic matter dynamics, only a handful of studies have examined their impacts on soil greenhouse gas emissions. Our study investigated fluxes of carbon dioxide (CO 2 ) and methane (CH 4 ) from three nests of Atta cephalotes using a portable greenhouse gas analyzer, and measured CO 2 and CH 4 emissions from soils containing nest holes that ranged 5.2–152.1 g CO 2 -C and –1.1 to 15,264.7 mg CH 4 -C m -2 day -1 , respectively. Fluxes of CO 2 and CH 4 were positively correlated above nest holes, but not in patches of soil away from leaf-cutting ant nests. Nearby non-nest soil emissions were significantly lower, ranging from 0.6 to 6.0 g CO 2 -C and –1.3 to 0.77 mg CH 4 -C m -2 day -1 . Fluxes of both gases among nests and among holes within a single nest were highly variable. This preliminary dataset is small in scale both temporarily and geographically, but the discovery of substantial greenhouse gas fluxes from Atta cephalotes nests may have important implications for carbon budgets of tropical and subtropical American forests. Further work will be necessary to determine the mechanisms behind enhanced greenhouse gas emissions from leaf-cutting ant nests, and how this may alter ecosystem-scale CO 2 emissions and CH 4 sink strength in tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗

Above-ground net primary productivity in regenerating seasonally dry tropical forest: Contributions of rainfall, forest age and soil

Identifying factors controlling forest productivity is critical to understanding forest-climate change feedbacks, modelling vegetation dynamics and carbon finance schemes. However, little research has focused on productivity in regenerating tropical forests which are expanding in their fraction of global area have an order of magnitude larger carbon uptake rates relative to older forest. Here We examined above-ground net primary productivity (ANPP) and its components (wood production and litterfall) over 10 years in forest plots that vary in successional age, soil characteristics and species composition using band dendrometers and litterfall traps in regenerating seasonally dry tropical forests in northwestern Costa Rica. We show that the components of ANPP are differentially driven by age and annual rainfall and that local soil variation is important. Total ANPP was explained by a combination of age, annual rainfall and soil variation. Wood production comprised 35% of ANPP on average across sites and years, and was explained by annual rainfall but not forest age. Conversely, litterfall increased with forest age and soil fertility yet was not affected by annual rainfall. In this region, edaphic variability is highly correlated with plant community composition. Thus, variation in ecosystem processes explained by soil may also be partially explained by species composition. These results suggest that future changes in annual rainfall can alter the secondary forest carbon sink, but this effect will be buffered by the litterfall flux which varies little among years. In determining the long-term strength of the secondary forest carbon sink, both rainfall and forest age will be critical variables to track. We also conclude that detailed understanding of local site variation in soils and plant community may be required to accurately predict the impact of changing rainfall on forest carbon uptake. Synthesis. We show that in seasonally dry tropical forest, annual rainfall has a positive relationship with the growth of above-ground woody tissues of trees and that droughts lead to significant reductions in above-ground productivity. These results provide evidence for climate change—carbon cycle feedbacks in the seasonal tropics and highlight the value of longitudinal data on forest regeneration.

54 ENVIRONMENTAL SCIENCES↗

The Roles of Climate Change and Climate Variability in the 2017 Atlantic Hurricane Season

The 2017 Atlantic hurricane season was extremely active with six major hurricanes, the third most on record. The sea-surface temperatures (SSTs) over the eastern Main Development Region (EMDR), where many tropical cyclones (TCs) developed during active months of August/September, were ~0.96 C above the 1901-2017 average (warmest on record): about ~0.42 C from a long-term upward trend and the rest (~80%) attributed to the Atlantic Meridional Mode (AMM). The contribution to the SST from the North Atlantic Oscillation (NAO) over the EMDR was a weak warming, while that from El Niño-Southern Oscillation (ENSO) was negligible. Nevertheless, ENSO, the NAO, and the AMM all contributed to favorable wind shear conditions, while the AMM also produced enhanced atmospheric instability. Compared with the strong hurricane years of 2005/2010, the ocean heat content (OHC) during 2017 was larger across the tropics, with higher SST anomalies over the EMDR and Caribbean Sea. On the other hand, the dynamical/thermodynamical atmospheric conditions, while favorable for enhanced TC activity, were less prominent than in 2005/2010 across the tropics. The results suggest that unusually warm SST in the EMDR together with the long fetch of the resulting storms in the presence of record-breaking OHC may be key factors in driving the strong TC activity in 2017.

SST↗

The Roles of Climate Change and Climate Variability in the 2017 Atlantic Hurricane Season

The 2017 Atlantic hurricane season was extremely active with six major hurricanes, the third most on record. The sea-surface temperatures (SSTs) over the eastern Main Development Region (EMDR), where many tropical cyclones (TCs) developed during active months of August/September, were ~0.96°C above the 1901-2017 average (warmest on record): about ~0.42°C from a long-term upward trend and the rest (~80%) attributed to the Atlantic Meridional Mode (AMM). The contribution to the SST from the North Atlantic Oscillation (NAO) over the EMDR was a weak warming, while that from El Nino Southern Oscillation (ENSO) was negligible. Nevertheless, ENSO, the NAO, and the AMM all contributed to favorable wind shear conditions, while the AMM also produced enhanced atmospheric instability. Compared with the strong hurricane years of 2005/2010, the ocean heat content (OHC) during 2017 was larger across the tropics, with higher SST anomalies over the EMDR and Caribbean Sea. On the other hand, the dynamical/thermodynamical atmospheric conditions, while favorable for enhanced TC activity, were less prominent than in 2005/2010 across the tropics. The results suggest that unusually warm SST in the EMDR together with the long fetch of the resulting storms in the presence of record-breaking OHC may be key factors in driving the strong TC activity in 2017.

Lim, Young-Kwon↗

The volcanic signal in Goddard Institute for Space Studies three-dimensional model simulations

Transient calculations of the Goddard Institute for Space Studies general circulation model for the climatic signal of volcanic eruptions are analyzed. By compositing the output for two different volcanoes for scenario A and five different volcanoes for scenario B, the natural variability is suppressed and the volcanic signals are extracted. Significant global mean surface air temperature cooling and precipitation reduction are found for several years following the eruptions, with larger changes in the Northern Hemisphere (NH) than in the Southern Hemisphere. The global-average temperature response lasts for more than four years, but the precipitation response disappears after three years. The largest cooling in the model occurs in the NH summer of the year after spring eruptions. Significant zonal-average temperature reductions begin in the tropics immediately after the eruptions and extend to 45 deg S - 45 deg N in the year after the eruptions. In the second year, cooling is still seen from 30 deg S to 30 deg N. Because of the low variability in this model as compared to the real world, these signals may appear more significant here than they would by attempting to isolate them using real data. The results suggest that volcanoes can enhance the drought in the Sahel. No evidence was found that stratospheric aerosols from the low-atitude volcanic eruptions can trigger El Nino-Southern Oscillation (ENSO) events in this model.

Robock, Alan↗

Fitting a Two-Component Scattering Model to Polarimetric SAR Data from Forests

Two simple scattering mechanisms are fitted to polarimetric synthetic aperture radar (SAR) observations of forests. The mechanisms are canopy scatter from a reciprocal medium with azimuthal symmetry and a ground scatter term that can represent double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants or Bragg scatter from a moderately rough surface, which is seen through a layer of vertically oriented scatterers. The model is shown to represent the behavior of polarimetric backscatter from a tropical forest and two temperate forest sites by applying it to data from the National Aeronautic and Space Agency/Jet Propulsion Laboratory's Airborne SAR (AIRSAR) system. Scattering contributions from the two basic scattering mechanisms are estimated for clusters of pixels in polarimetric SAR images. The solution involves the estimation of four parameters from four separate equations. This model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem. The model is used to develop an understanding of the ground-trunk double-bounce scattering that is present in the data, which is seen to vary considerably as a function of incidence angle. Two parameters in the model fit appear to exhibit sensitivity to vegetation canopy structure, which is worth further exploration. Results from the model fit for the ground scattering term are compared with estimates from a forward model and shown to be in good agreement. The behavior of the scattering from the ground-trunk interaction is consistent with the presence of a pseudo-Brewster angle effect for the air-trunk scattering interface. If the Brewster angle is known, it is possible to directly estimate the real part of the dielectric constant of the trunks, a key variable in forward modeling of backscatter from forests. It is also shown how, with a priori knowledge of the forest height, an estimate for the attenuation coefficient of the canopy can be obtained directly from the multi-incidence-angle polarimetric observations. This attenuation coefficient is another key variable in forward models and is generally related to the canopy density.

temperate forests↗

Regional Variations in the Diurnal Cycle of Tropical Precipitation as Represented by IMERG, ERA5, and Spaceborne Ku Radar

The diurnal cycle of precipitation is highly regional and is typically a product of multiple competing effects that can be highly localized. The diurnal cycle in high precipitation regions such as the Amazon and the Maritime Continent are of particular interest, especially due to the complex coastal effects which take place over the Maritime Continent. The high spatial and temporal resolution provided by the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (GPM) mission (IMERG) dataset, is used in this study to examine the fine-scale features of the diurnal cycle in these regions. Using an 18-year (2000 – 2018) record of IMERG precipitation observations, diurnal and semidiurnal phase and amplitude are calculated using a fast Fourier transform (FFT) method on precipitation averaged for each half-hour of the day at 0.1°x0.1° spatial resolution. We first introduce an objective method of identifying locations where the diurnal signals are robust and strong. Clear patterns of precipitation phase propagation with distance from shore are shown over both regions, with the diurnal phase and amplitude exhibiting a strong dependence on the distance from the coastline. Semidiurnal cycles are generally weaker than the diurnal cycle except in some isolated locations. Similar analysis is also conducted on the ERA5 reanalysis data in order to evaluate the model’s representation of the precipitation diurnal cycle. The model captures the broad scale patterns of diurnal variability but does not capture all the fine scale patterns nor the exact timing that is observed by IMERG. Comparisons are also made to a long record Ku radar dataset created by combining Tropical Rainfall Measuring Mission (TRMM) and GPM observations, thus providing an additional point of comparison for the timing of the ERA5 precipitation peak, since the timing precipitation can be different, even in between observational datasets.

L J Hayden↗

Carbon Dynamics in Vegetation and Soils

The overall goals of CD-08 team in Phase I were to quantify the contributions of different components of the carbon cycle to overall ecosystem carbon balance in Amazonian tropical forests and to undertake process studies at a number of sites along the eastern LBA transect to understand how and why these fluxes vary with site, season, and year. We divided this work into a number of specific tasks: (1) determining the average rate (and variability) of tree growth over the past 3 decades; (2) determining age demographics of tree populations, using radiocarbon to determine tree age; (3) assessing the rate of production and decomposition of dead wood debris; (4) determining turnover rates for organic matter in soils and the mean age of C respired from soil using radiocarbon measurements; and (5) comparing our results with models and constructing models to predict the potential of tropical forests to function as sources or sinks of C. This report summarizes the considerable progress made towards our original goals, which have led to increased understanding of the potential for central Amazon forests to act as sources or sinks of carbon with altered productivity. The overall picture of tropical forest C dynamics emerging from our Phase I studies suggests that the fraction of gross primary production allocated to growth in these forests is only 25-30%, as opposed to the 50% assumed by many ecosystem models. Consequent slow tree growth rates mean greater mean tree age for a given diameter, as reflected in our measurements and models of tree age. Radiocarbon measurements in leaf and root litter suggest that carbon stays in living tree biomass for several years up to a decade before being added to soils, where decomposition is rapid. The time lags predicted from 14C, when coupled with climate variation on similar time scales, can lead to significant interannual variation in net ecosystem C exchange.

Trumbore, Susan↗

Generating Vegetation Leaf Area Index Earth System Data Record from Multiple Sensors: Implementation, Analysis and Validation - Part 2

The evaluation of a new global monthly leaf area index (LAI) data set for the period July 1981 to December 2006 derived from AVHRR Normalized Difference Vegetation Index (NDVI) data is described. The physically based algorithm is detailed in the first of the two part series. Here, the implementation, production and evaluation of the data set are described. The data set is evaluated both by direct comparisons to ground data and indirectly through inter-comparisons with similar data sets. This indirect validation showed satisfactory agreement with existing LAI products, importantly MODIS, at a range of spatial scales, and significant correlations with key climate variables in areas where temperature and precipitation limit plant growth. The data set successfully reproduced well-documented spatio-temporal trends and inter-annual variations in vegetation activity in the northern latitudes and semi-arid tropics. Comparison with plot scale field measurements over homogeneous vegetation patches indicated a 7% underestimation when all major vegetation types are taken into account. The error in mean values obtained from distributions of AVHRR LAI and high-resolution field LAI maps for different biomes is within 0.5 LAI for six out of the ten selected sites. These validation exercises though limited by the amount of field data, and thus less than comprehensive, indicated satisfactory agreement between the LAI product and field measurements. Overall, the intercomparison with short-term LAI data sets, evaluation of long term trends with known variations in climate variables, and validation with field measurements together build confidence in the utility of this new 26 year LAI record for long term vegetation monitoring and modeling studies.

Ganguly, Sangram↗

Observations of the Global Characteristics and Regional Radiative Effects of Marine Cloud Liquid Water

The large-scale spatial distribution and temporal variability of cloud liquid water path (LWP) over the world's oceans and the relationship of cloud LWP to temperature and the radiation budget are investigated using recent satellite measurements from the Special Sensor Microwave/Imager (SSM/I), the Earth Radiation Budget Experiment (ERBE), and the International Satellite Cloud Climatology Project (ISCCP). Observations of cloud liquid water on a 2.5 deg x 2.5 deg and are used over a 53-month period beginning July 1987 and ending in December 1991. The highest values of cloud liquid water (greater than 0.13 kg/sq m) occur largely along principal routes of northern midlatitude storms and in areas dominated by tropical convection. The zonally averaged structure is distinctly trimodal, where maxima appear in the midlatitudes and near the equator. The average marine cloud LWP over the globe is estimated to be about 0.113 kg/sq m. Its highest seasonal variability is typically between 15% and 25% of the annual mean but in certain locations can exceed 30%. Comparisons of cloud LWP to temperature for low clouds during JJA and DJF of 1990 show significant positive correlations at colder temperatures and negative correlations at warmer temperatures. The correlations also exhibit strong seasonal and regional variation. Coincident and collocated observations of cloud LWP from the SSM/I and albedo measurements from the Earth Radiation Budget Satellite (ERBS) and the NOAA-10 satellite are compared for low clouds in the North Pacific and North Atlantic. The observed albedo-LWP relationships correspond reasonably well with theory, where the average cloud effective radius (r(sub e)) is 11.1 microns and the standard deviation is 5.2 microns. The large variability in the inferred values of r(sub e) suggests that other factors may be important in the albedo-LWP relationships. In terms of the effect of the LWP on the net cloud forcing, the authors find that a 0.05 kg/sq m increase in LWP (for LWP less than 0.2 kg/sq m) results in a -25 W/sq m change in the net cloud forcing at a solar zenith angle of 75 deg.

Greenwald, Thomas J.↗

Tropical Tropospheric Ozone Climatology: Approaches Based on SHADOZ Observations

The SHADOZ (Southern Hemisphere Additional Ozonesondes) ozone sounding network was initiated in 1998 to improve the coverage of tropical in-situ ozone measurements for satellite validation, algorithm development and related process studies. Over 2000 soundings have been archived at the central website, , for 12 stations that span the entire equatorial zone [Thompson et al., JGR, 108,8238, 2003]. The most striking features of tropospheric ozone profiles in SHADOZ are: (1) persistent longitudinal variability in tropospheric ozone profiles, with a 10-15 DU column-integrated difference between Atlantic and Pacific sites; (2) intense short-term variability triggered by changing meteorological conditions and advection of pollution. The implications of these results for profile climatologies and trends are described along with several approaches to classifying ozone profiles: 1) Seasonal means during MAM (March-April-May) and SON (September-October-November); 2) Maxima and minima, identified through correlation of TOMS-derived TTO (tropical tropospheric ozone) column depth with the sonde integrated tropospheric ozone column; and 3) Meteorological regimes, a technique that is effective in the subtropics where tropical and mid-latitude conditions alternate.

Thompson, Anne M.↗

Biodiversity and Stoichiometric Plasticity Increase Pico‐Phytoplankton Contributions to Marine Net Primary Productivity and the Biological Pump

Earth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low-latitude oligotrophic systems harboring pico-sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon-climate system. We have added three pico-phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico-eukaryote populations that contribute significantly to productivity and export within the subtropical to mid-latitude transition zone, with the western subtropical regions of each basin supporting the most P-poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C-rich/P-poor stoichiometry, preferentially inhabited by pico-eukaryotes and diatoms. Collectively, pico-phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico-phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.

Environmental Sciences & Ecology↗