Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global climate 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,459 records · Page 81

A Variable-Resolution Stretched-Grid General Circulation Model and Data Assimilation System with Multiple Areas of Interest: Studying the Anomalous Regional Climate Events of 1998

The new stretched-grid design with multiple (four) areas of interest, one at each global quadrant, is implemented into both a stretched-grid GCM (general circulation model) and a stretched-grid data assimilation system (DAS). The four areas of interest include: the U.S./Northern Mexico, the El Nino area/Central South America, India/China, and the Eastern Indian Ocean/Australia. Both the stretched-grid GCM and DAS annual (November 1997 through December 1998) integrations are performed with 50 km regional resolution. The efficient regional down-scaling to mesoscales is obtained for each of the four areas of interest while the consistent interactions between regional and global scales and the high quality of global circulation, are preserved. This is the advantage of the stretched-grid approach. The global variable resolution DAS incorporating the stretched-grid GCM has been developed and tested as an efficient tool for producing regional analyses and diagnostics with enhanced mesoscale resolution. The anomalous regional climate events of 1998 that occurred over the U.S., Mexico, South America, China, India, African Sahel, and Australia are investigated in both simulation and data assimilation modes. Tree assimilated products are also used, along with gauge precipitation data, for validating the simulation results. The obtained results show that the stretched-grid GCM and DAS are capable of producing realistic high quality simulated and assimilated products at mesoscale resolution for regional climate studies and applications.

Fox-Rabinovitz, Michael S.↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

Measuring stomatal and guard cell metrics for plant physiology and growth using StoManager1

Abstract Automated guard cell detection and measurement are vital for understanding plant physiological performance and ecological functioning in global water and carbon cycles. Most current methods for measuring guard cells and stomata are laborious, time-consuming, prone to bias, and limited in scale. We developed StoManager1, a high-throughput tool utilizing geometrical, mathematical algorithms, and convolutional neural networks to automatically detect, count, and measure over 30 guard cell and stomatal metrics, including guard cell and stomatal area, length, width, stomatal aperture area/guard cell area, orientation, stomatal evenness, divergence, and aggregation index. Combined with leaf functional traits, some of these StoManager1-measured guard cell and stomatal metrics explained 90% and 82% of tree biomass and intrinsic water use efficiency (iWUE) variances in hardwoods, making them substantial factors in leaf physiology and tree growth. StoManager1 demonstrated exceptional precision and recall (mAP@0.5 over 0.96), effectively capturing diverse stomatal properties across over 100 species. StoManager1 facilitates the automation of measuring leaf stomatal and guard cells, enabling broader exploration of stomatal control in plant growth and adaptation to environmental stress and climate change. This has implications for global gross primary productivity (GPP) modeling and estimation, as integrating stomatal metrics can enhance predictions of plant growth and resource usage worldwide. Easily accessible open-source code and standalone Windows executable applications are available on a GitHub repository (https://github.com/JiaxinWang123/StoManager1) and Zenodo (https://doi.org/10.5281/zenodo.7686022).

Wang, Jiaxin (ORCID:0000000348085085)↗

Soil moisture and the persistence of North American drought

Numerical sensitivity experiments on the effects of soil moisture on North American summertime climate are performed using a 12-layer global atmospheric general circulation model. Consideration is given to the hypothesis that reduced soil moisture may induce and amplify warm, dry summers of midlatitude continental interiors. The simulations resemble the conditions of the summer of 1988, including an extensive drought over much of North America. It is found that a reduction in soil moisture leads to an increase in surface temperature, lower surface pressure, increased ridging aloft, and a northward shift of the jet stream. It is shown that low-level moisture advection from the Gulf of Mexico is important in the maintenance of persistent soil moisture deficits.

Oglesby, Robert J.↗

Ontolology Negotiation Between Scientific Archives

This paper describes an approach to ontology negotiation between information agents. Ontologies are declarative (data driven) expressions of an agent's "world": the objects, operations, facts, and rules that constitute the logical space within which an agent performs. Ontology negotiation enables agents to cooperate in performing a task, even if they are based on different ontologies. 'Me process allows agents to discover ontology conflicts and then, though incremental interpretation, clarification, and explanation, establish a common basis for communicating with each other. The need for ontology negotiation stems from the proliferation of information sources and of agents with widely varying specialty expertise. The unmanageability of massive amounts of web-based information is already becoming apparent. It is starting to have an impact on professions that rely on distributed archived information. If the expansion continues at its present rate without an ontology negotiation process being introduced, there will soon be no way to ensure the accuracy and completeness of information that scientists obtain from sources other than their own experiments. Ontology negotiation is becoming increasingly recognized as a crucial element of scalable agent technology. This is because agents, by their very nature, are supposed to operate with a fair amount of autonomy and independence from their end-users. Part of this independence is the ability to enlist other agents for help in performing a task (such as locating information on the web). The agents enlisted for help may be "owned" by a different end-user or organization (such as a document archive), and there is no guarantee that they will use the same terminology or understand the same concepts (objects, operators, theorems, rules) as the recruiting agent. For NASA, the need for ontology negotiation arises at the boundaries between scientific disciplines. For example: modeling the effects of global warming might involve knowledge about imaging, climate analysis, ecology, demographics, industrial economics, and biology. The need for ontology negotiation also arises at the boundaries between scientific programs. For example, a Principal Investigator may want to use information from a previous mission to complement downloads from the instruments currently deployed.

Bailin, Sidney C.↗

Geolocation Assessment Algorithm for CALIPSO Using Coastline Detection

Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) is a joint satellite mission between NASA and the French space agency CNES. The investigation will gather long-term, global cloud and aerosol optical and physical properties to improve climate models. The CALIPSO spacecraft is scheduled to launch in 2004 into a 98.2 inclination, 705 km circular orbit approximately 3 minutes behind the Aqua spacecraft. The payload consists of a two-wavelength polarization-sensitive lidar, and two passive imagers operating in the visible (0.645 mm) and infrared (8.7 - 12.0 mm) spectral regions. The imagers are nadir viewing and co-aligned with the lidar. Earth viewing measurements are geolocated to the Earth fixed coordinate system using satellite ephemeris, Earth rotation and geoid, and instrument pointing data. The coastline detection algorithm will assess the accuracy of the CALIPSO geolocation process by analyzing Wide Field Camera (WFC) visible ocean land boundaries. Processing space-time coincident MODIS and WFC scenes with the coastline algorithm will help verify the co-registration requirement with Moderate Resolution Imaging Spectrometer (MODIS) data. This paper quantifies the accuracy of the coastline geolocation assessment algorithm.

Currey, J. Chris↗

Far Infrared Measurements of Cirrus

Improved techniques for remote sensing of cirrus are needed to obtain global data for assessing the effect of cirrus in climate change models. Model calculations show that the far infrared/sub-millimeter spectral region is well suited for retrieving cirrus Ice Water Path and particle size parameters. Especially useful cirrus information is obtained at frequencies below 60 cm-1 where single particle scattering dominates over thermal emission for ice particles larger than about 50 m. Earth radiance spectra have been obtained for a range of cloud conditions using an aircraft-based Fourier transform spectrometer. The Far InfraRed Sensor for Cirrus (FIRSC) is a Martin-Puplett interferometer which incorporates a polarizer for the beamsplitter and can be operated in either intensity or linear polarization measurement mode. Two detector channels span 10 to 140 cm-1 with a spectral resolution of 0.1 cm-1; achieving a Noise Equivalent Temperature of approximately 1K at 30 cm-1 in a 4 sec scan. Examples are shown of measured and modeled Earth radiance for a range of cloud conditions from 1998 and 1999 flights.

Nolt, I. G.↗

International Satellite Cloud Climatology Project

The International Satellite Cloud Climatology Project (ISCCP) was established as the first project of the World Climate Research Programme (WCP-2) to collect and analyze satellite radiance measurements to infer the global distribution of cloud radiative properties and their diurnal and seasonal variations. The operational phase of ISCCP began in July 1983 and is currently planned to continue through June 2010. [Mission Objectives] To produce a global, reduced resolution, infrared and visible, calibrated and normalized radiance data set containing basic information on the radiative properties of the atmosphere from which cloud parameters can be derived. To stimulate and coordinate basic research on techniques for inferring the physical properties of clouds from the condensed radiance data set and to apply the resulting algorithms to derive and validate a global cloud climatology for improving the parameterization of clouds in climate models. To promote research using ISCCP data and contributing to improved understanding of the Earth's radiation budget (top of the atmosphere and surface) and hydrological cycle. [Temporal_Coverage: Start_Date=1983-07-01; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180].

Rossow, William B.↗

Performance Analysis of Scientific and Engineering Applications Using MPInside and TAU

In this paper, we present performance analysis of two NASA applications using performance tools like Tuning and Analysis Utilities (TAU) and SGI MPInside. MITgcmUV and OVERFLOW are two production-quality applications used extensively by scientists and engineers at NASA. MITgcmUV is a global ocean simulation model, developed by the Estimating the Circulation and Climate of the Ocean (ECCO) Consortium, for solving the fluid equations of motion using the hydrostatic approximation. OVERFLOW is a general-purpose Navier-Stokes solver for computational fluid dynamics (CFD) problems. Using these tools, we analyze the MPI functions (MPI_Sendrecv, MPI_Bcast, MPI_Reduce, MPI_Allreduce, MPI_Barrier, etc.) with respect to message size of each rank, time consumed by each function, and how ranks communicate. MPI communication is further analyzed by studying the performance of MPI functions used in these two applications as a function of message size and number of cores. Finally, we present the compute time, communication time, and I/O time as a function of the number of cores.

Saini, Subhash↗

High Spectral Resolution Lidar – Generation 2 (HSRL-2) Retrievals of Ocean Surface Wind Speed: Methodology and Evaluation

Ocean surface wind speed (i.e., wind speed 10 m above sea level) is a critical parameter used by atmospheric models to estimate the state of the marine atmospheric boundary layer (MABL). Accurate surface wind speed measurements in diverse locations are required to improve characterization of MABL dynamics and assess how models simulate large-scale phenomena related to climate change and global weather patterns. To provide these measurements, this study introduces and evaluates a new surface wind speed data product from the NASA Langley Research Center nadir-viewing High Spectral Resolution Lidar – generation 2 (HSRL-2) using data collected as part of the NASA Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) mission. The HSRL-2 can directly measure vertically resolved aerosol backscatter and extinction profiles without additional constraints or assumptions, enabling the instrument to accurately derive atmospheric attenuation and directly determine surface reflectance (i.e., surface backscatter). Also, the high horizontal spatial resolution of the HSRL-2 retrievals (0.5 s or ∼ 75 m along track) allows the instrument to probe the fine-scale spatial variability in surface wind speeds over time along the flight track and over breaks in broken cloud fields. A rigorous evaluation of these retrievals is performed by comparing coincident HSRL-2 and National Center for Atmospheric Research (NCAR) Airborne Vertical Atmosphere Profiling System (AVAPS) dropsonde data, owing to the joint deployment of these two instruments on the ACTIVATE King Air aircraft. These comparisons show correlations of 0.89, slopes of 1.04 and 1.17, and y intercepts of −0.13 and −1.05 m s−1 for linear and bisector regressions, respectively, and the overall accuracy is calculated to be 0.15 ± 1.80 m s−1. It is also shown that the dropsonde surface wind speed data most closely follow the HSRL-2 distribution of wave slope variance using the distribution proposed by Hu et al. (2008) rather than the ones proposed by Cox and Munk (1954) and Wu (1990) for surface wind speeds below 7 m s−1, with this category comprising most of the ACTIVATE data set. The retrievals are then evaluated separately for surface wind speeds below 7 m s−1 and between 7 and 13.3 m s−1 and show that the HSRL-2 retrieves surface wind speeds with a bias of ∼ 0.5 m s−1 and an error of ∼ 1.5 m s−1, a finding not apparent in the cumulative comparisons. Also, it is shown that the HSRL-2 retrievals are more accurate in the summer (−0.18 ± 1.52 m s−1) than in the winter (0.63 ± 2.07 m s−1), but the HSRL-2 is still able to make numerous (N=236) accurate retrievals in the winter. Overall, this study highlights the abilities and assesses the performance of the HSRL-2 surface wind speed retrievals, and it is hoped that further evaluation of these retrievals will be performed using other airborne and satellite data sets.

Sanja Dmitrovic↗

Model-Based Estimation of Sampling-Caused Uncertainty in Aerosol Remote Sensing for Climate Research Applications

To evaluate the effect of sampling frequency on the global monthly mean aerosol optical thickness (AOT), we use 6 years of geographical coordinates of Moderate Resolution Imaging Spectroradiometer (MODIS) L2 aerosol data, daily global aerosol fields generated by the Goddard Institute for Space Studies General Circulation Model and the chemical transport models Global Ozone Chemistry Aerosol Radiation and Transport, Spectral Radiationtransport Model for Aerosol Species and Transport Model 5, at a spatial resolution between 1.125 deg × 1.125 deg and 2 deg × 3◦: the analysis is restricted to 60 deg S-60 deg N geographical latitude. We found that, in general, the MODIS coverage causes an underestimate of the global mean AOT over the ocean. The long-term mean absolute monthly difference between all and dark target (DT) pixels was 0.01-0.02 over the ocean and 0.03-0.09 over the land, depending on the model dataset. Negative DT biases peak during boreal summers, reaching 0.07-0.12 (30-45% of the global long-term mean AOT). Addition of the Deep Blue pixels tempers the seasonal dependence of the DT biases and reduces the mean AOT difference over land by 0.01-0.02. These results provide a quantitative measure of the effect the pixel exclusion due to cloud contamination, ocean sun-glint and land type has on the MODIS estimates of the global monthly mean AOT. We also simulate global monthly mean AOT estimates from measurements provided by pixel-wide along-track instruments such as the Aerosol Polarimetry Sensor and the Cloud-Aerosol LiDAR with Orthogonal Polarization. We estimate the probable range of the global AOT standard error for an along-track sensor to be 0.0005-0.0015 (ocean) and 0.0029-0.01 (land) or 0.5-1.2% and 1.1-4% of the corresponding global means. These estimates represent errors due to sampling only and do not include potential retrieval errors. They are smaller than or comparable to the published estimate of 0.01 as being a climatologically significant change in the global mean AOT, suggesting that sampling density is unlikely to limit the use of such instruments for climate applications at least on a global, monthly scale.

long-term variability↗

The GGCMI Phase 2 emulators: global gridded crop model responses to changes in CO2, temperature, water, and nitrogen (version 1.0)

Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: that growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts.

James A Franke↗

On the Feasibility of Studying Shortwave Aerosol Radiative Forcing of Climate Using Dual-Wavelength Aerosol Backscatter Lidar

The current low confidence in the estimates of aerosol-induced perturbations of Earth's radiation balance is caused by the highly non-uniform compositional, spatial and temporal distributions of tropospheric aerosols on a global scale owing to their heterogeneous sources and short lifetimes. Nevertheless, recent studies have shown that the inclusion of aerosol effects in climate model calculations can improve agreement with observed spatial and temporal temperature distributions. In light of the short lifetimes of aerosols, determination of their global distribution with space-borne sensors seems to be a necessary approach. Until recently, satellite measurements of tropospheric aerosols have been approximate and did not provide the full set of information required to determine their radiative effects. With the advent of active aerosol remote sensing from space (e.g., PICASSO-CENA), the applicability fo lidar-derived aerosol 180 deg -backscatter data to radiative flux calculations and hence studies of aerosol effects on climate needs to be investigated.

Redemann, Jens↗

On a simplified solution of climate-carbon dynamics in idealized flat10MIP simulations

Abstract. Idealized experiments with coupled climate-carbon Earth system models (ESMs) provide a basis for understanding the response of the carbon cycle to external forcing and for quantifying climate-carbon feedbacks. Here, we analyze globally-averaged results from idealized esm-flat10 experiments and show that most models exhibit a quasi-linear relationship between cumulative carbon uptake on land and in the ocean during a period of constant fossil fuel emissions of 10 Pg C yr−1. We hypothesize that this relationship does not depend on emission pathways. Further, as a simplification, we quantify the relationship between cumulative ocean carbon uptake and changes in ocean heat content using a linear approximation. In this way, changes in oceanic heat content and atmospheric CO2 concentration become interdependent variables, reducing the coupled temperature-CO2 system to just one differential equation. The equation can be solved analytically or numerically for the atmospheric CO2 concentration as a function of fossil fuel emissions. This approach leads to a simplified description of global carbon and climate dynamics, which could be used for applications beyond existing analytical frameworks.

Brovkin, Victor↗

Assessing Effects of Climate and Technology Uncertainties in Large Natural Resource Allocation Problems

The productivity of the world's natural resources is critically dependent on a variety of highly uncertain factors, which obscure individual investors and governments that seek to make long-term, sometimes irreversible, investments in their exploration and utilization. These dynamic considerations are poorly represented in disaggregated resource models, as incorporating uncertainty into large-dimensional problems presents a challenging computational task. In this paper, we apply the SCEQ algorithm (Cai and Judd, 2023) to solve a large-scale dynamic stochastic global land resource use problem with stochastic crop yields due to adverse climate impacts and limits on further technological progress. For the same model parameters and bounded shocks, the range of land conversion is considerably smaller for the dynamic stochastic model than for deterministic scenario analysis.

numerical methods↗

Evaluating E3SM Global Storm‐Resolving Model Simulations of Deep Convection: Insights From DP‐SCREAM During TRACER

Global Storm-Resolving Models (GSRMs) are becoming increasingly vital for advancing climate modeling and improving the prediction of extreme weather events. Houston, a coastal region frequently affected by deep convective storms, offers an ideal setting to evaluate the ability of GSRMs to simulate deep convection. This study assesses the performance of the Doubly Periodic Simple Cloud-Resolving E3SM (Energy Exascale Earth System Model) Atmosphere Model (DP-SCREAM) using observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign. DP-SCREAM effectively reproduces the diurnal cycles of clouds and precipitation, demonstrating much greater skill than the E3SM single column model. The DP-SCREAM is demonstrated to be applicable to coastal regions, partially due to the forcing data sets already capturing the influence of breezes. DP-SCREAM also replicates biases persistent in the global version of SCREAM: the underrepresentation of boundary layer shallow clouds, a lack of mid-level congestus clouds, and the popcorn convection, characterized by small and disorganized convective cells generating the strongest precipitation. To investigate these issues, two sensitivity experiments were conducted: increasing the mixing length and scaling up the buoyancy flux within the Simplified Higher Order Closure scheme. Increasing the mixing length improved mid-level congestus representation and reduced unrealistic early morning fog occurrence. Enhancing buoyancy flux only marginally improved the bias of underproduced big convective cells. In conclusion, an additional resolution sensitivity test at 0.5 km grid spacing demonstrated that a refined horizontal resolution alone is insufficient to resolve these biases.

54 ENVIRONMENTAL SCIENCES↗