Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ensemble models”

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 1,315 records · Page 73

Fouling modeling and prediction approach for heat exchangers using deep learning

In this article, we develop a generalized and scalable statistical model for accurate prediction of fouling resistance using commonly measured parameters of industrial heat exchangers. This prediction model is based on deep learning where a scalable algorithmic architecture learns non-linear functional relationships between a set of target and predictor variables from large number of training samples. Here, the efficacy of this modeling approach is demonstrated for predicting fouling in an analytically modeled cross-flow heat exchanger, designed for waste heat recovery from flue-gas using room temperature water. The performance results of the trained models demonstrate that the mean absolute prediction errors are under 10 –4 KW –1 for flue-gas side, water side and overall fouling resistances. The coefficients of determination (R 2 ), which characterize the goodness of fit between the predictions and observed data, are over 99%. Even under varying levels of measurement noise in the inputs, we demonstrate that predictions over an ensemble of multiple neural networks achieves better accuracy and robustness to noise. We find that the proposed deep-learning fouling prediction framework learns to follow heat exchanger flow and heat transfer physics, which we confirm using locally interpretable model agnostic explanations around randomly selected operating points. Overall, we provide a robust algorithmic framework for fouling prediction that can be generalized and scaled to various types of industrial heat exchangers.

42 ENGINEERING↗

The design and performance of a real-time self excited vocoder

This paper is concerned with a generic class of predictive speech coders that includes the newly proposed Self Excited Vocoder (SEV) and the well known Code-Excited Linear Predictive Coder (CELPC). All members of this class form an excitation sequence for a linear predictive model filter using the same general model for the excitation signal. The general excitation model is based on a block coding technique where each sequence is drawn from an ensemble of sequences. This paper reports on two developments related to this general model. The first development is a new type of excitation ensemble that can in general be populated by many different types of sequences. The second development is a means of populating this new type of ensemble based on a vector quantizer design procedure using a new distortion measure.

Rose, Richard C.↗

Process-oriented stochastic perturbations applied to the parametrization of turbulence and shallow convection for ensemble prediction

Weather forecasting nowadays often requires some estimation of uncertainties associated with the output of meteorological models, in order to better inform decision making, especially in the context of intense weather events. Ensemble prediction systems provide such information through sets of possible scenarios which are designed to represent various uncertainty sources, including model uncertainties. A wide variety of methods have been proposed to estimate model uncertainties, among which perturbation methods targeting uncertain processes are a promising research field. Here, we focus on the representation of small-scale variability by process-oriented perturbation schemes applied to two key physical processes, namely turbulence and shallow convection. The perturbations are applied to a single-column version of the convection-permitting AROME model, in three idealized boundary-layer cases. Large-eddy simulations (LESs) of the same cases serve as a reference for the subgrid variability that has to be represented, and the results are also compared to those given by the Stochastically Perturbed Parametrization Tendencies (SPPT) method, which is a method commonly used by weather forecast centres to represent model uncertainty. The spread produced by our process-oriented perturbations of turbulence and shallow convection does not represent all the small-scale variability implied by the LESs for temperature and humidity. However, it is of a similar order of magnitude for the wind, thanks to perturbations generated by the stochastic turbulence scheme. The dispersion is structurally different from what is obtained with SPPT. It is non-negligible in the lower levels, where SPPT perturbations are usually suppressed because of numerical instabilities, indicating a possible complementarity between the schemes.

54 ENVIRONMENTAL SCIENCES↗

Benchmarking the performance of uncertainty quantification methods for neural network-based interatomic potentials

Machine-learned interatomic potentials (ML-IAPs) continue to gain popularity as accurate, computationally efficient replacements for traditional, physics-based interatomic potentials and expensive ab initio methods. Uncertainty quantification (UQ) of ML-IAPs is a growing area of research as UQ is critical in many applications of IAPs, such as developing curated datasets, active learning-based data augmentation, self-improving models, and estimating the uncertainty of molecular dynamics simulations. In this paper, we construct and benchmark a series of different neural network potentials (NNPs) with varying network architectures to determine the performance of these models with respect to both the mean and uncertainty calibration error. Each NNP method is specifically designed to predict either epistemic or aleatoric uncertainty with particular focus on the differences in behavior between the epistemic and aleatoric uncertainty estimates. We benchmark these methods using multiple datasets common in the ML-IAP literature. The results show that the aleatoric uncertainty from single-shot model architectures is a competitive alternative to ensemble-based epistemic uncertainty predictions in regions of sufficient data-density. However, in regions where the representative data is sparse, aleatoric uncertainty models tend to overpredict and epistemic methods tend to underpredict the actual model error. We conclude that the type of UQ is crucial when discussing performance of probabilistic model results as different methods have different performance characteristics depending on the regime in which they are evaluated. Therefore, the type of UQ method should be carefully evaluated against both the data characteristics and requirements for the intended application.

97 MATHEMATICS AND COMPUTING↗

Automatic Determination of the Conic Coronal Mass Ejection Model Parameters

Characterization of the three-dimensional structure of solar transients using incomplete plane of sky data is a difficult problem whose solutions have potential for societal benefit in terms of space weather applications. In this paper transients are characterized in three dimensions by means of conic coronal mass ejection (CME) approximation. A novel method for the automatic determination of cone model parameters from observed halo CMEs is introduced. The method uses both standard image processing techniques to extract the CME mass from white-light coronagraph images and a novel inversion routine providing the final cone parameters. A bootstrap technique is used to provide model parameter distributions. When combined with heliospheric modeling, the cone model parameter distributions will provide direct means for ensemble predictions of transient propagation in the heliosphere. An initial validation of the automatic method is carried by comparison to manually determined cone model parameters. It is shown using 14 halo CME events that there is reasonable agreement, especially between the heliocentric locations of the cones derived with the two methods. It is argued that both the heliocentric locations and the opening half-angles of the automatically determined cones may be more realistic than those obtained from the manual analysis

Pulkkinen, A.↗

Northern Winter Climate Change: Assessment of Uncertainty in CMIP5 Projections Related to Stratosphere-Troposphere Coupling

Future changes in the stratospheric circulation could have an important impact on northern winter tropospheric climate change, given that sea level pressure (SLP) responds not only to tropospheric circulation variations but also to vertically coherent variations in troposphere-stratosphere circulation. Here we assess northern winter stratospheric change and its potential to influence surface climate change in the Coupled Model Intercomparison Project-Phase 5 (CMIP5) multimodel ensemble. In the stratosphere at high latitudes, an easterly change in zonally averaged zonal wind is found for the majority of the CMIP5 models, under the Representative Concentration Pathway 8.5 scenario. Comparable results are also found in the 1% CO2 increase per year projections, indicating that the stratospheric easterly change is common feature in future climate projections. This stratospheric wind change, however, shows a significant spread among the models. By using linear regression, we quantify the impact of tropical upper troposphere warming, polar amplification, and the stratospheric wind change on SLP. We find that the intermodel spread in stratospheric wind change contributes substantially to the intermodel spread in Arctic SLP change. The role of the stratosphere in determining part of the spread in SLP change is supported by the fact that the SLP change lags the stratospheric zonally averaged wind change. Taken together, these findings provide further support for the importance of simulating the coupling between the stratosphere and the troposphere, to narrow the uncertainty in the future projection of tropospheric circulation changes.

climate change↗

Attribution of the role of climate change in the forest fires in Sweden 2018

In this study, we analyse the role of climate change in the forestfires that raged through large parts of Sweden in the summer of 2018 from ameteorological perspective. This is done by studying the Canadian FireWeather Index (FWI) based on sub-daily data, both in reanalysis data sets(ERA-Interim, ERA5, the Japanese 55 year Reanalysis, JRA-55, and Modern-Era Retrospective analysis for Research and Applications version 2, MERRA-2) and three large-ensemble climatemodels (EC-Earth, weather@home, W@H, and Community Earth System Model, CESM) simulations. The FWI, based on reanalysis, correlates well with the observed burnt area in summer (r=0.6 to 0.8). We find that the maximum FWI in July 2018 had return times of ~24 years (90 % CI, confidence interval, > 10 years) for southern and northern Sweden. Furthermore, we find a negative trend of the FWI for southern Sweden over the 1979 to 2017 time period in the reanalyses, yielding a non-significant reduced probability of such an event. However, the short observational record, large uncertainty between the reanalysis products and large natural variability of the FWI give a large confidence interval around this number that easily includes no change, so we cannot draw robust conclusions from reanalysis data. The three large-ensembles with climate models point to a roughly 1.1 (0.9 to 1.4) times increased probability (non-significant) for such events in the current climate relative to preindustrial climate. For a future climate (2 °C warming), we find a roughly 2 (1.5 to 3) times increased probability for such events relative to the preindustrial climate. The increased fire weather risk is mainly attributed to the increase intemperature. The other main factor, i.e. precipitation during summer months, isprojected to increase for northern Sweden and decrease for southern Sweden. We, however, do not find a clear change in prolonged dry periods in summer months that could explain the increased fire weather risk in the climate models. In summary, we find a (non-significant) reduced probability of such eventsbased on reanalyses, a small (non-significant) increased probability due to global warming up to now and a more robust (significant) increase in the risk for such events in the future based on the climate models.

54 ENVIRONMENTAL SCIENCES↗

Projected changes in extreme streamflow and inland flooding in the mid-21st century over Northeastern United States using ensemble WRF-Hydro simulations

Study region: Northeastern United States (NEUS). Study focus: We investigate the potential impacts of climate change on precipitation, streamflow, and inland flooding in the NEUS during the mid-21st century. Dynamically downscaled climate projections from three global climate models for historical (1995-2004) and future (2045-2054) periods under business-as-usual scenarios were used to force the hydrologic model WRF-Hydro at 200-meter resolution and create ensemble hydrologic simulations. Additionally, an extreme value model was developed to project the risks associated with low-frequency hydrologic events. New hydrological insights for the region: Results from four major watersheds indicate a significantly wetter regime in winter months and potential drier conditions during late spring to early summer. Discharges in fall are projected to decrease in the northern watersheds and increase toward the south. Extreme flow and water depths resulting from extreme inland flooding are projected to increase by 5-20% and > 100%, respectively. The extent of the total flooded area is likely to be 20% greater by the mid-century. These increased risks can be attributed to (i) an approximate 25% increase in decadal mean and > 40% increase in decadal extreme precipitation intensity, (ii) up to 30% lower snow availability and 5-25% higher evapotranspiration throughout the year, and (iii) a projected 5% increase in soil moisture in all seasons except summer. Furthermore, rapid snow melting in winter will likely cause an earlier peak flow in the rivers.

54 ENVIRONMENTAL SCIENCES↗

Ice Sheet and Ice Shelf Freshwater Forcing of the NASA GISS ModelE GCM

Recent mass loss from ice sheets and ice shelves is now persistent and prolonged enough that it impacts downstream oceanographic conditions. To demonstrate this, we use an ensemble of coupled GISS-E2.1-G (a CMIP6 model) simulations forced with historical estimates of anomalous freshwater, in addition to other climate forcings, from 1990 through 2019. In this ensemble there are detectable differences in zonal-mean sea surface temperatures (SST) and sea ice in the Southern Ocean, and in regional sea level around Antarctica and in the western North Atlantic. These impacts mostly improve the model's representation of historical changes, including reversing the forced trends in Southern Ocean surface temperature and Antarctic sea ice. The changes in SST may have implications for estimates of the SST pattern effect on climate sensitivity and for cloud feedbacks. We conclude that the changes are sufficiently large that model groups should strive to include more accurate estimates of these drivers in all-forcing historical simulations in future coupled model intercomparisons.

Ice sheet↗

Performance of Two Cloud-Radiation Parameterization Schemes in the Finite Volume General Circulation Model for Anomalously Wet May and June 2003 Over the Continental United States and Amazonia

An objective assessment of the impact of a new cloud scheme, called Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRAS) (together with its radiation modules), on the finite volume general circulation model (fvGCM) was made with a set of ensemble forecasts that invoke performance evaluation over both weather and climate timescales. The performance of McRAS (and its radiation modules) was compared with that of the National Center for Atmospheric Research Community Climate Model (NCAR CCM3) cloud scheme (with its NCAR physics radiation). We specifically chose the boreal summer months of May and June 2003, which were characterized by an anomalously wet eastern half of the continental United States as well as northern regions of Amazonia. The evaluation employed an ensemble of 70 daily 10-day forecasts covering the 61 days of the study period. Each forecast was started from the analyzed initial state of the atmosphere and spun-up soil moisture from the first-day forecasts with the model. Monthly statistics of these forecasts with up to 10-day lead time provided a robust estimate of the behavior of the simulated monthly rainfall anomalies. Patterns of simulated versus observed rainfall, 500-hPa heights, and top-of-the-atmosphere net radiation were recast into regional anomaly correlations. The correlations were compared among the simulations with each of the schemes. The results show that fvGCM with McRAS and its radiation package performed discernibly better than the original fvGCM with CCM3 cloud physics plus its radiation package. The McRAS cloud scheme also showed a reasonably positive response to the observed sea surface temperature on mean monthly rainfall fields at different time leads. This analysis represents a method for helpful systematic evaluation prior to selection of a new scheme in a global model.

Sud, Y. C.↗

Feature‐shared adaptive‐boost deep learning for invasiveness classification of pulmonary subsolid nodules in CT images

Purpose In clinical practice, invasiveness is an important reference indicator for differentiating the malignant degree of subsolid pulmonary nodules. These nodules can be classified as atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA), or invasive adenocarcinoma (IAC). The automatic determination of a nodule's invasiveness based on chest CT scans can guide treatment planning. However, it is challenging, owing to the insufficiency of training data and their interclass similarity and intraclass variation. To address these challenges, we propose a two‐stage deep learning strategy for this task: prior‐feature learning followed by adaptive‐boost deep learning. Methods The adaptive‐boost deep learning is proposed to train a strong classifier for invasiveness classification of subsolid nodules in chest CT images, using multiple 3D convolutional neural network (CNN)‐based weak classifiers. Because ensembles of multiple deep 3D CNN models have a huge number of parameters and require large computing resources along with more training and testing time, the prior‐feature learning is proposed to reduce the computations by sharing the CNN layers between all weak classifiers. Using this strategy, all weak classifiers can be integrated into a single network. Results Tenfold cross validation of binary classification was conducted on a total of 1357 nodules, including 765 noninvasive (AAH and AIS) and 592 invasive nodules (MIA and IAC). Ablation experimental results indicated that the proposed binary classifier achieved an accuracy of with an AUC of 81.3 . These results are superior compared to those achieved by three experienced chest imaging specialists who achieved an accuracy of , , and , respectively. About 200 additional nodules were also collected. These nodules covered 50 cases for each category (AAH, AIS, MIA, and IAC, respectively). Both binary and multiple classifications were performed on these data and the results demonstrated that the proposed method definitely achieves better performance than the performance achieved by nonensemble deep learning methods. Conclusions It can be concluded that the proposed adaptive‐boost deep learning can significantly improve the performance of invasiveness classification of pulmonary subsolid nodules in CT images, while the prior‐feature learning can significantly reduce the total size of deep models. The promising results on clinical data show that the trained models can be used as an effective lung cancer screening tool in hospitals. Moreover, the proposed strategy can be easily extended to other similar classification tasks in 3D medical images.

Wang, Jun↗

Three different directions in which the European Union could replace Russian natural gas

Russia's invasion of Ukraine fuelled an energy crisis, which considerably impacted Europe given its heavy reliance on Russian natural gas imports. This study uses an ensemble of four global integrated assessment models, which are further soft-linked to two sectoral models, and explores the synergies and trade-offs among three approaches to living without Russian gas in Europe: (a) replacing with other gas imports, (b) boosting domestic energy production, and (c) reducing demand and accelerating energy efficiency. We find that substituting Russian gas from other trade partners would miss an opportunity to accelerate decarbonisation in end-use sectors while risking further fossil-fuel lock-ins, despite featuring the lowest gas price spikes and potentially reducing heating costs for end-users in the near term. Boosting domestic, primarily renewable, energy production on the other hand would instead require considerable investments, potentially burdening consumers. Energy demand reductions, however, could offer considerable space for further emissions cuts at the lowest power-sector investment costs; nonetheless, an energy efficiency-driven strategy would also risk relocation of energy-intensive industries, an aspect of increasing relevance to EU policymakers.

03 NATURAL GAS↗

Hydrological projections over the Upper Indus Basin at 1.5 °C and 2.0 °C temperature increase

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.

54 ENVIRONMENTAL SCIENCES↗

Smaller Is Better: The Case for Lower-Order Iodoplumbate Species Dominating MAPbI 3 /Dimethylformamide Solutions

Here, using complementary experimental measurements and computational predictions of spectroscopic measurements (EXAFS, XANES, and UV–vis), we have determined the identity of the most stable iodoplumbate species in dilute lead halide perovskite precursor solutions. We have determined which species are most likely to be thermodynamically stable compared to others that are unstable or metastable. Condensed phase ab initio models were constructed, and the resulting ensembles were used to directly compare the computed signals to the experimental results of the EXAFS, XANES, and UV–vis spectra of PbI 2 :MAI in DMF. The results of this study suggest that only Pb 2+ , PbI + , and PbI 2 are dominant in the dilute lead perovskite precursor solutions as thermodynamically stable entities. Our interpretation of the relative stability of iodoplumbate species in solution, based on an analysis of EXAFS and XANES spectra, provides critically important new insight into the species most likely to be responsible for crystal nucleation and growth in these materials. This insight will have a significant consequence on the broad scientific community and will necessitate the reinterpretation of peaks in the UV–vis spectra of lead halide perovskite precursor solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superlattices are Greener on the Other Side: How Light Transforms Self-Assembled Mixed Halide Perovskite Nanocrystals

Perovskite nanocrystal superlattices (NC SLs) are the nearest real-world approximations to monodisperse NC ensembles. NC SLs thus represent ideal model systems for evaluating the optical and structural stability of CsPb(I 1– x Br x ) 3 NCs at a macroscopic level. Here, photoinduced changes to CsPb(I 1– x Br x ) 3 NC SLs (0 < x < 1.0) are probed via in situ photoluminescence, X-ray diffraction, and electron microscopy. We find that prolonged (~10–20 h) ultraviolet–visible irradiation causes irreversible PL blueshifts, photobrightening, and crystal structure contractions. These changes stem from gradual photoinduced I 2 sublimation, which transforms CsPb(I 1– x Br x ) 3 into CsPbBr 3 . Despite eliminating half of the initial halides from individual CsPb(I 0.53 Br 0.47 ) 3 particles, NCs within SLs remarkably preserve their initial crystallinity, cuboidal shapes, edge lengths, and size distributions. This work illustrates compositional control toward generating precisely engineered perovskite NC SLs. It also highlights iodide photo-oxidation as a hurdle that must be overcome if mixed halide perovskite nanomaterials are to be applied beyond fundamental studies.

36 MATERIALS SCIENCE↗

Uncertainty in Reconstructing Paleo‐Elevation of the Antarctic Ice Sheet From Temperature‐Sensitive Ice Core Records

Paleotemperature reconstructions from ice cores are mixed signals of changes in climate and ice‐surface elevation. A common, temperature‐based paleoaltimetry method suggests these signals can be disentangled by comparing two proxy locations with similar climates. The difference between the records is assumed to be due to elevation, which is estimated by scaling the temperature difference by a lapse rate. We investigate the uncertainty associated with this approach using a case study of the Antarctic Ice Sheet during the Last Glacial Maximum. From an ensemble of climate simulations, we extract modeled temperatures at locations of real ice cores. We find uncertainty on the order of hundreds of meters that results from spatial heterogeneity in non‐adiabatic temperature change, which itself stems in part from elevation‐induced atmospheric circulation change. Our findings suggest that caution is needed when interpreting temperature‐based paleoaltimetry results for ice sheets.

58 GEOSCIENCES↗

The Potential Benefits of Handling Mixture Statistics via a Bi-Gaussian EnKF: Tests With All-Sky Satellite Infrared Radiances

The meteorological characteristics of cloudy atmospheric columns can be very different from their clear counterparts. Thus, when a forecast ensemble is uncertain about the presence/absence of clouds at a specific atmospheric column (i.e., some members are clear while others are cloudy), that column's ensemble statistics will contain a mixture of clear and cloudy statistics. Such mixtures are inconsistent with the ensemble data assimilation algorithms currently used in numerical weather prediction. Hence, ensemble data assimilation algorithms that can handle such mixtures can potentially outperform currently used algorithms. In this study, we demonstrate the potential benefits of addressing such mixtures through a bi-Gaussian extension of the ensemble Kalman filter (BGEnKF). The BGEnKF is compared against the commonly used ensemble Kalman filter (EnKF) using perfect model observing system simulated experiments (OSSEs) with a realistic weather model (the Weather Research and Forecast model). Synthetic all-sky infrared radiance observations are assimilated in this study. In these OSSEs, the BGEnKF outperforms the EnKF in terms of the horizontal wind components, temperature, specific humidity, and simulated upper tropospheric water vapor channel infrared brightness temperatures. This study is one of the first to demonstrate the potential of a Gaussian mixture model EnKF with a realistic weather model. Our results thus motivate future research toward improving numerical Earth system predictions though explicitly handling mixture statistics.

54 ENVIRONMENTAL SCIENCES↗

Skillful Multi‐Month Predictions of Ecosystem Stressors in the Surface and Subsurface Ocean

Abstract Anthropogenic carbon emissions and associated climate change are driving rapid warming, acidification, and deoxygenation in the ocean, which increasingly stress marine ecosystems. On top of long‐term trends, short term variability of marine stressors can have major implications for marine ecosystems and their management. As such, there is a growing need for predictions of marine ecosystem stressors on monthly, seasonal, and multi‐month timescales. Previous studies have demonstrated the ability to make reliable predictions of the surface ocean physical and biogeochemical state months to years in advance, but few studies have investigated forecast skill of multiple stressors simultaneously or assessed the forecast skill below the surface. Here, we use the Community Earth System Model (CESM) Seasonal to Multiyear Large Ensemble (SMYLE) along with novel observation‐based biogeochemical and physical products to quantify the predictive skill of dissolved inorganic carbon (DIC), dissolved oxygen, and temperature in the surface and subsurface ocean. CESM SMYLE demonstrates high physical and biogeochemical predictive skill multiple months in advance in key oceanic regions and frequently outperforms persistence forecasts. We find up to 10 months of skillful forecasts, with particularly high skill in the Northeast Pacific (Gulf of Alaska and California Current Large Marine Ecosystems) for temperature, surface DIC, and subsurface oxygen. Our findings suggest that dynamical marine ecosystem prediction could support actionable advice for decision making.

Environmental Sciences & Ecology↗