Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial 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 775 records · Page 43

Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie

Ecosystems at coastal terrestrial–aquatic interfaces play a significant role in global biogeochemical cycles. In this study, we aimed to characterize coastal wetlands with particular focus on the co-variability between plant dynamics, topography, soil, and other environmental factors. We proposed a functional zonation approach based on machine learning clustering to identify the spatial regions, i.e., zones that capture these co-varied properties. This approach was applied to publicly available datasets along Lake Erie, in the Great Lakes Region. We investigated the heterogeneity of coastal ecosystem structures as a function of along-shore distance and transverse distance, based on the spatial data layers, including topography, wetland vegetation cover, and the time series of Landsat’s enhanced vegetation index (EVI) between 1990 and 2020. Results showed that the topographic metrics (elevation and slope), soil texture, and plant productivity influence the spatial distribution of wetland land-covers (emergent and phragmites). These results highlight a natural organization along the transverse axis, where the elevation and the EVI increase further away from the coastline. In addition, the clustering analysis allowed us to identify regions with distinct environmental characteristics, as well as the ones that are more sensitive to interannual lake-level variations.

47 OTHER INSTRUMENTATION↗

Understanding the Spatiotemporal Variability of Tropical Orographic Rainfall Using Convective Plume Buoyancy

Abstract Mechanical forcing by orography affects precipitating convection across many tropical regions, but controls on the intensity and horizontal extent of the orographic precipitation peak and rain shadow remain poorly understood. A recent theory explains this control of precipitation as arising from modulation of lower-tropospheric temperature and moisture by orographic mechanical forcing, setting the distribution of convective rainfall by controlling parcel buoyancy. Using satellite and reanalysis data, we evaluate this theory by investigating spatiotemporal precipitation variations in six mountainous tropical regions spanning South and Southeast Asia, and the Maritime Continent. We show that a strong relationship holds in these regions between daily precipitation and a measure of convective plume buoyancy. This measure depends on boundary layer thermodynamic properties and lower-free-tropospheric moisture and temperature. Consistent with the theory, temporal variations in lower-free-tropospheric temperature are primarily modulated by orographic mechanical lifting through changes in cross-slope wind speed. However, winds directed along background horizontal moisture gradients also influence lower-tropospheric moisture variations in some regions. The buoyancy measure is also shown to explain many aspects of the spatial patterns of precipitation. Finally, we present a linear model with two horizontal dimensions that combines mountain wave dynamics with a linearized closure exploiting the relationship between precipitation and plume buoyancy. In some regions, this model skillfully captures the spatial structure and intensity of rainfall; it underestimates rainfall in regions where time-mean ascent in large-scale convergence zones shapes lower-tropospheric humidity. Overall, these results provide new understanding of fundamental processes controlling subseasonal and spatial variations in tropical orographic precipitation.

54 ENVIRONMENTAL SCIENCES↗

Visualization and Quantification of Wind Induced Variability in Hydrogen Clouds Following Releases of Liquid Hydrogen: Preprint

Well characterized experimental data for consequence model validation is important in progressing the use of liquid hydrogen as an energy carrier. In 2019, the Health and Safety Executive (HSE) undertook a series of liquid hydrogen dispersion and combustion experiments as a part of the Pre-normative Research into the Safe Use of Liquid Hydrogen (PRESLHY) project. In partnership between the National Renewable Energy Laboratory (NREL) and HSE, time and spatially varying hydrogen concentration measurements were made in 25 dispersion experiments and 23 congested ignition experiments associated with PRESLHY WP3 and WP5, respectively. These measurements were undertaken using the hydrogen wide area monitoring system developed by NREL. During the 23 congested ignition experiments, high variability was observed in the measured explosion severity during experiments with similar initial conditions. This led to the conclusion that wind, including localized gusts, had a large influence on the dispersion of the hydrogen, and therefore the quantity of hydrogen that was present in the congested region of the explosions. Using the hydrogen concentration measurements taken immediately prior to ignition, the hydrogen clouds were visualized in an attempt to rationalize the variability in overpressure between the tests. Gaussian process regression was applied to quantify the variability of the measured hydrogen concentrations. This analysis could also be used to guide modifications in experimental designs for future research on hydrogen combustion behavior.

HSR&D↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗

Deep Learning and Photogrammetric Reconstruction for Automated Crack Detection and Dimensional Measurement in Mining Operations

Surface crack detection and dimensional measurement at active mining sites present significant safety and operational challenges. Manual inspection methods are labor-intensive, spatially incomplete, and expose personnel to hazardous environments, while existing automated approaches have been developed primarily for concrete civil infrastructure and have not been validated on the complex, variable surfaces characteristic of mining environments. This dissertation presents an automated pipeline that integrates deep learning semantic segmentation with Structure-from-Motion photogrammetry to detect surface cracks and measure their aperture, length, and vertical displacement from standard RGB imagery acquired during routine Uncrewed Aerial Vehicle (UAV) survey operations, without requiring additional sensor hardware or manual measurement. The pipeline combines a U-Net architecture with an EfficientNet-B0 encoder, pretrained on the SDNET2018 concrete crack dataset and fine-tuned on a mining-specific dataset spanning laboratory concrete specimens, coal refuse impoundment embankments, and post-blast limestone quarry benches. Photogrammetric reconstruction is performed using COLMAP Structure-from-Motion and Multi-View Stereo, with crack segmentation masks projected into the reconstructed point cloud to enable three-dimensional vertical displacement measurement through local plane fitting and bimodal surface detection. The pipeline was validated across 36 controlled laboratory specimens at three imaging distances and four vertical displacement levels, achieving aperture measurement RMSE of 0.047 cm and R² of 0.954, and vertical displacement RMSE of 0.140 cm and R² of 0.966, against independent caliper measurements. Field application at a coal refuse impoundment in southwestern Pennsylvania detected 71 crack components across the embankment crest, with a dominant longitudinal crack exhibiting aperture values reaching 28 cm and a 95th percentile vertical displacement of 35.53 cm, consistent in magnitude and spatial distribution with simultaneously acquired LiDAR-derived estimates. Application across four post-blast limestone quarry bench datasets in California successfully characterized blast-induced fracture networks at ground sampling distances ranging from 0.59 to 1.23 cm/pixel, with detected crack geometries physically consistent with observable surface conditions at each site. The results demonstrate that deep learning-based crack detection and photogrammetric measurement can be integrated into routine UAV inspection workflows at mining sites, providing repeatable, scalable, and quantitative crack characterization across surface types, crack scales, and displacement magnitudes not previously addressed in the literature. The pipeline requires no dedicated surveying equipment beyond the UAV platforms already deployed at mine sites for survey and monitoring purposes, supporting practical adoption within existing operational workflows.

Crack detection, Dimensional Measurement↗

High-Resolution Meteorology with Climate Change Impacts from Global Climate Model Data Using Generative Machine Learning

As renewable energy generation increases, the impacts of weather and climate on energy generation and demand become critical to the reliability of the energy system. However, these impacts are often overlooked. Global climate models (GCMs) can be used to understand possible changes to our climate, but their coarse resolution makes them difficult to use in energy system modelling. Here we present open-source generative machine learning methods that produce meteorological data at a nominal spatial resolution of 4 km at an hourly frequency based on inputs from 100 km daily-average GCM data. These methods run 40 times faster than traditional downscaling methods and produce data that have high-resolution spatial and temporal attributes similar to historical datasets. We demonstrate that these methods can be used to downscale projected changes in wind, solar and temperature variables across multiple GCMs including projections for more frequent low-wind and high-temperature events in the Eastern United States.

climate change↗

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗

Understanding Air-Sea Feedbacks to the MJO through Process Evaluation of Observations and Models (Final Report)

The goal of our project was to enhance our understanding of the role of coupling in the simulation of the Madden Julian Oscillation (MJO) in coupled climate models. We developed diagnostics of precipitation and Sea Surface Temperature variability in climate models and investigated how they are related to metrics of MJO simulation performance. We found a relationship between models with a high spatial coherence of precipitation and a good propagation of the MJO through the Maritime Continent region, but no relationship between metrics of precipitation organisation and intensity. Analysis of coupled and uncoupled model simulations from the Sixth Coupled Model Intercomparison Project (CMIP6) shows that coupling systematically improved the MJO simulation in both metrics that we investigated. We analysed the role of background moisture gradients in the MJO simulation performance following Lim et al. (2018) and found that although the model to model variation in moisture gradients to the west of Indonesia was related to model MJO performance, in general no single moisture gradient was as good at explaining a large fraction of MJO performance across the models as compared to a linear combination of those metrics, which was able to explain a significant fraction of the variance in the MJO performance. Those relationships were much stronger for coupled simulations than uncoupled simulations and stronger for those metrics that measured propagation of the MJO across the Maritime Continent than for overall MJO activity.

54 ENVIRONMENTAL SCIENCES↗

A high accuracy/resolution spectral element/Fourier–Galerkin method for the simulation of shoaling non-linear internal waves and turbulence in long domains with variable bathymetry

A high-order hybrid continuous-Galerkin numerical method, designed for the simulation of non-linear, non -hydrostatic internal waves and turbulence in long computational domains with complex bathymetry, is presented. The spatial discretization in the non-periodic wave-propagating directions, utilizes the nodal spectral element method. Such a high-order element-based discretization allows the highly accurate representation of complex domain geometry along with the flexibility of concentrating resolution in areas of interest. Under the assumption of the normal-to-isobath propagation of non-linear internal waves, a third periodic direction is incorporated via a Fourier-Galerkin discretization. The distinct non-hydrostatic nature of non-linear internal waves and, any instabilities and turbulence therein, necessitates the numerically challenging solution of the pressure Poisson problem. A defining feature of this work is the application of a domain decomposition approach, combined with block-Jacobi/deflation-based preconditioning to the pressure Poisson problem. Such a combined approach is particularly suitable for the long high aspect-ratio complex domains of interest and enables the efficient high-accuracy reproduction of the non-hydrostatic dynamics of non-linear internal waves. Implementation details are also described in the context of the stability of the solver and its parallelization strategy. A series of benchmarks of increasing complexity demonstrate the robustness of the flow solver. The benchmarks culminate with the three-dimensional simulation of a convectively breaking mode-one non-linear internal wave over a realistic South-China-Sea bathymetric transect and background current/stratification profiles.

Deflation↗

Impacts of climate change on estuarine stratification and implications for hypoxia within a shallow subtropical system

Vertical density stratification often plays an important role in the formation and expansion of coastal hypoxic zones through its effect on near-bed circulation and vertical oxygen flux. However, the impact of future climate change on estuarine circulation is widely unknown. Here, we developed and calibrated a three-dimensional hydrodynamic model for Pensacola Bay, a shallow subtropical estuary in the northeastern Gulf of Mexico. Model simulations based on years 2013–2017 were applied to examine changes in salinity, temperature, and density under future climate scenarios, including increased radiative forcing (IR) and temperature (T), increased freshwater discharge (D), sea level rise (SLR), and wind intensification (W). Simulations showed that the impacts of climate change on modeled state variables varied over time with external forcing conditions. The model demonstrated the potential for sea level rise and increased freshwater discharge to episodically increase vertical density gradients in the Bay. However, increased wind forcing destabilized vertical gradients, at times reducing the spatial extent and duration of stable stratification. For time periods with low freshwater discharge, moderate increases in wind speed (10%) can destabilize density gradients strengthened by increased discharge (10%) and sea level rise (0.48 m). In contrast, destruction of strong density gradients that form near the mid-Bay channel following flood events requires stronger wind forcing. Furthermore, these results highlight the importance of considering natural variability in freshwater and wind forcing, as well as local phenomena that are generally unresolved by global climate models.

54 ENVIRONMENTAL SCIENCES↗

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Likeness: a toolkit for connecting the social fabric of place to human dynamics

The ability to produce richly-attributed synthetic populations is key for understanding human dynamics, responding to emergencies, and preparing for future events, all while protecting individual privacy. The Likeness toolkit accomplishes these goals with a suite of Python packages: pymedm/pymedm\_legacy, livelike, and actlike. This production process is initialized in pymedm (or pymedm\_legacy) that utilizes census microdata records as the foundation on which disaggregated spatial allocation matrices are built. The next step, performed by livelike, is the generation of a fully autonomous agent population attributed with hundreds of demographic census variables. The agent population synthesized in livelike is then attributed with residential coordinates in actlike based on block assignment and, finally, allocated to an optimal daytime activity location via the street network. We present a case study in Knox County, Tennessee, synthesizing 30 populations of public K–12 school students \& teachers and allocating them to schools. Validation of our results shows they are highly promising by replicating reported school enrollment and teacher capacity with a high degree of fidelity.

Tuccillo, Joe↗

Barriers and variable spacing enhance convective cooling and increase power output in solar PV plants

When the temperature of solar photovoltaic (PV) modules rises, efficiency drops and module degradation accelerates. Thus, it is beneficial to reduce module operating temperatures. Previous studies of solar power plants have illustrated that incoming flow characteristics, turbulent mixing, and array geometry can strongly impact convective cooling, as measured by the convective heat transfer coefficient h. In the fields of heat transfer and plant canopy flow, previous work has shown that system-scale arrangement modifications—e.g., variable spacing, barriers, or windbreaks—can passively alter the flow, enhance turbulent mixing, and influence convection. However, researchers have not yet explored how variable spacing or barriers might enhance convective cooling in solar power plants. Here, high-resolution large-eddy simulations model the air flow and heat transfer through solar power plant arrangements modified with missing modules and barrier walls. We then perform a control volume analysis to evaluate the net heat flux and compute h, which quantifies the influence of these spatial modifications on convective cooling and, thus, module temperature and power output. Installing barrier walls yields the greatest improvements, increasing h by 3.4%, reducing module temperature by an estimated 2.5 °C, and boosting power output by an estimated 1.4% on average. These findings indicate that incorporating variable spacing or barrier-type elements into PV plant designs can reduce module temperature and, thus, improve PV performance and service life.

14 SOLAR ENERGY↗

A spatially regularized detector for emergent/re-emergent disease outbreaks

Early detection of outbreaks caused by emergent pathogens, using epidemiological surveillance data i.e., daily case counts, is difficult. This is because the data tend to be noisy during the early epoch of the outbreak. In contrast, the spread-rate of the disease tends to be well-behaved, as it depends only on the mixing patterns of the population and the characteristics of the pathogen, neither of which behave erratically in space-time. In this report, we explore whether the spread-rate can be used for epidemiological surveillance, conditional on case count data. Estimating the spread-rate from case count data allows us to exploit exogenous information, e.g., incubation period distributions etc., which can considerably smooth out any erratic temporal behavior. Further, epidemiological dynamics are spatially correlated, and if case counts are available for multiple areal units e.g., counties, these correlations could potentially be used to suppress noise in the early epoch data. These exogenous information and structure are not exploited by conventional syndromic surveillance detectors to extract a well-behaved latent variable for monitoring purposes. The technical challenge lies in the estimation of the spread-rate field jointly over a collection of areal units; further, the spread-rate varies over time. We develop a method based on mean-field variational inference to approximately estimate the spread-rate field, using a Gaussian Random Field Model for spatial regularization. The method is tested on the estimation of spread-rate in the thirty-three counties of New Mexico and detect the arrival of the Fall 2020 COVID-19 wave in September 2020. We find that the method is scalable, but underestimates the uncertainty in the estimated spread-rate field. We detect the arrival of the Fall 2020 wave a week ahead of conventional syndromic surveillance algorithms, but our simplistic detection algorithm, based on simple anomaly detection, suffers from a high false positive rate, similar to conventional detectors.

59 BASIC BIOLOGICAL SCIENCES↗

Small-Scale Irrigation: Improving Food Security under Changing Climate and Water Resource Conditions in Ethiopia

We develop a new systems modeling tool that integrates knowledge from hydrology, agriculture, and economics to understand the effect of small-scale irrigation on food security and groundwater sustainability in Ethiopia. Irrigation is an effective tool to mitigate climate impacts and improve agricultural yields. Small-scale irrigation, such as decentralized groundwater irrigation, is well suited for developing countries where smallholder farming communities are widely dispersed and can only afford small infrastructure investment. We study the underlying interdependencies between food and water systems in Ethiopia, where small-holder agriculture is the foundation of the nation’s economy and climate variability has led to great challenges to its food security. Our coupled market and crop model with groundwater module captures the interdependencies of climate, water availability (including irrigation), crop yield, farmland allocation, crop production, transport and consumption based on a system approach across multiple spatial scales. We study the implication of small-scale irrigation to Ethiopia’s food security and water resource conditions as a “what-if” question by comparing an irrigation scenario to the calibrated baseline in 2015, a year of significant drought and crop failure over a large portion of Ethiopia. Our model offers fresh insights into geographic disparities in outcomes that are driven by baseline climate variability, soil fertility, and market conditions. In general, we find that small-scale irrigation can potentially improve food security through increases in food consumption, but it requires policy support to direct the increases of production to domestic consumption while maintaining a sustainable groundwater condition. By using Ethiopia as an example, we show the strength of our model to study how water infrastructure resources support critical functions and service in water and food systems.

Zhang, Ying↗

Quantifying Variability and Controls of Riverine Dissolved Organic Carbon Exported to Arctic Coastal Margins of North America (Final Report)

This project involved implementation and application of a coupled permafrost hydrology and dissolved organic carbon process models to investigate how spatial and seasonal variations in terrestrial hydrology and soil freeze/thaw dynamics influence the mobilization, loading, and export of organic carbon to the stream network for selected arctic basins across northern Alaska and northwest Canada. The model simulations were constrained by detailed observations of in- stream chemistry, soil active layer profile moisture and temperature dynamics, streamflow, soil carbon inventories and satellite microwave remote sensing based assessments of surface soil freeze-thaw dynamics. We developed and applied the numerical modeling and data analysis, incorporating observed data for calibration and validation, in order to investigate the terrestrial hydrology, permafrost dynamics, and associated DOC production and loading to rivers across a region encompassing watersheds draining to the coast. The project produced six publications and three datasets archived in public repositories.

54 ENVIRONMENTAL SCIENCES↗