Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Groundwater Flux”

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 109 records · Page 6

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗

Remote semi-continuous flow rate logging seepage meter

The movement of groundwater and its associated solutes from upland regions was implicated in the degradation of receiving surface water bodies. Current efforts to directly measure this influx of water incorporate manually operated seepage meters which are hindered by severe limitations. A prototype seepage meter was developed by NASA Langley Research Center and Virginia Polytechnic Institute and State University that will allow for the semi-continuous collection and data logging of seepage flux across the sediment water interface. The meter is designed to operate at depths to 40 meters, and alleviate or minimize all disadvantages associated with traditional methods while remaining cost effective. The unit was designed to operate independently for time periods on the order of weeks with adjustable sample sequences depending upon hydrologic conditions. When used in conjunction with commercially available pressure transducers, this seepage meter allows for correlations to be made between groundwater discharge and tidal/sea state conditions in coastal areas. Field data from the Chesapeake Bay and Florida Bay systems are presented.

Reay, William G.↗

Riverbed Temperature and 4D ERT Monitoring Reveals Heterogenous Horizontal and Vertical Groundwater-Surface Water Exchange Flows Under Dynamic Stage Conditions

Groundwater surface water exchange plays a critical role in physical, biological, and geochemical function of coastal and riverine systems. Observing exchange flow behavior in heterogeneous systems is a primary challenge, particularly when flows are governed by dynamic river stage or tidal variations. In this paper we demonstrate a novel application of time-lapse 3D electrical resistivity tomography and temperature monitoring where an array of thermistors installed beneath a riverbed double as resistivity electrodes. We use the array to monitor stage driven exchange flows over a 6-day period in a dynamic, stage-driven high order stream. We present a method for addressing the otherwise confounding effects of the moving river-surface boundary on the raw resistivity data, thereby enabling successful tomographic imaging. Temperature time-series at each thermistor location and time-lapse 3D images of changes in bulk electrical conductivity together provide a detailed description of exchange dynamics over a 10-meter by 45-meter section of the riverbed, to a depth of approximately 5 m. Results reveal highly variable flux behavior throughout the monitoring domain including both horizontal and vertical exchange flows.

54 ENVIRONMENTAL SCIENCES↗

Vadose zone flushing of fertilizer tracked by isotopes of water and nitrate

Abstract A substantial fraction of nitrogen (N) fertilizer applied in agricultural systems is not incorporated into crops and moves below the rooting zone as nitrate (NO 3 − ). Understanding mechanisms for soil N retention below the rooting zone and leaching to groundwater is essential for our ability to track the fate of added N. We used dual stable isotopes of nitrate (δ 15 N–NO 3 − and δ 18 O–NO 3 − ) and water (δ 18 O–H 2 O and δ 2 H–H 2 O) to understand the mechanisms driving nitrate leaching at three depths (0.8, 1.5, and 3.0 m) of an irrigated corn field sampled every 2 weeks from 2016 to 2020 in the southern Willamette Valley, Oregon, USA. Distinct periods of high nitrate concentrations with lower δ 15 N–NO 3 − values indicated that a portion of that nitrate was from recent fertilizer applications. We used a mixing model to quantify nitrate fluxes associated with recently added fertilizer N versus older, legacy soil N during these “fertilizer signal periods.” Nitrate leached below 3.0 m in these periods made up a larger proportion of the total N leached at that depth (∼52%) versus the two shallower depths (∼13%–16%), indicating preferential movement of recently applied fertilizer N through the deep soil into groundwater. Further, N associated with recent fertilizer additions leached more easily when compared to remobilized legacy N. A high volume of fall and winter precipitation may push residual fertilizer N to depth, potentially posing a larger threat to groundwater than legacy N. Optimizing fertilizer N additions could minimize fertilizer losses and reduce nitrate leaching to groundwater.

Agriculture↗

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗

Mitigating drought mortality by incorporating topography into variable forest thinning strategies

Abstract Drought-induced productivity reductions and tree mortality have been increasing in recent decades in forests around the globe. Developing adaptation strategies hinges on an adequate understanding of the mechanisms governing the drought vulnerability of forest stands. Prescribed reduction in stand density has been used as a management tool to reduce water stress and wildfire risk, but the processes that modulate fine-scale variations in plant water supply and water demand are largely missing in ecosystem models. We used an ecohydrological model that couples plant hydraulics with groundwater hydrology to examine how within-stand variations in tree spatial arrangements and topography might mitigate forest vulnerability to drought at individual-tree and stand scales. Our results demonstrated thinning generally ameliorated plant hydraulic stress and improved carbon and water fluxes of the remaining trees, although the effectiveness varied by climate and topography. Variable thinning that adjusted thinning intensity based on topography-mediated water availability achieved higher stand productivity and lower mortality risk, compared to evenly-spaced thinning at comparable intensities. The results from numerical experiments provided mechanistic evidence that topography mediates the effectiveness of thinning and highlighted the need for an explicit consideration of within-stand heterogeneity in trees and abiotic environments when designing forest thinning to mitigate drought impacts.

Environmental Sciences & Ecology↗

Carbon Dioxide Production in Bedrock Beneath Soils Substantially Contributes to Forest Carbon Cycling

Abstract Soils are widely considered the primary terrestrial organic matter pool mediating carbon transactions with the atmosphere and groundwater. Because soils are both a host and a product of rhizosphere activity, they are thought to mark the location where photosynthetic fixation of carbon dioxide (CO 2 ) is balanced by the oxidation of organic matter. However, in many terrestrial environments, the rhizosphere extends below soils and into fractured bedrock, and it is unknown if the resulting biological and hydrologic dynamics in bedrock have a significant impact on carbon cycling. Here we show substantial production of CO 2 in weathered bedrock at 4–8 m below the thin soils (<0.5 m thick) of a Northern California forest using innovative monitoring technology for sampling gases and water in fractured rock. The deep CO 2 production supports a persistent upward flux of CO 2(g) year‐round from bedrock to soil, constituting between 2% and 29% of the average daily CO 2 efflux from the land surface. When water is rapidly transported across the fractured bedrock vadose zone, nearly 50% of the CO 2 produced in the bedrock dissolves into water, promoting water‐rock interaction and export of dissolved inorganic carbon (DIC) from the unsaturated zone to groundwater, constituting as much as 80% of the DIC exiting the hillslope. Such CO 2 production in weathered bedrock is subject to unique moisture, temperature, biological, and mineralogical conditions which are currently missing from terrestrial carbon models.

54 ENVIRONMENTAL SCIENCES↗

DPFEHM: a differentiable subsurface physics simulator

The Earth’s subsurface is a key resource that provides energy via fossil fuels and geothermal sources, stores drinking water, and is used in the fight against climate change via carbon sequestration. Simulating the physical processes that occur in the Earth’s subsurface with computers enables better use of this resource. DPFEHM is a Julia package that includes computer models with a focus on the Earth’s subsurface, especially fluid flow, which is critical for the aforementioned applications. DPFEHM is able to solve the groundwater flow equations (single phase flow), Richards equation (air/water), the advection-dispersion equation, and the 2d wave equation. One of the key features of DPFEHM is that it supports automatic differentiation, so it can be integrated into machine learning workflows using frameworks such as Flux or PyTorch. The automatic differentiation capabilities give it the same performance as adjoint methods.

54 ENVIRONMENTAL SCIENCES↗

Stable isotopes of carbon dioxide (CO2), dissolved inorganic carbon (DIC), water and concentrations of greenhouse gases in a mountainous watershed (Crested Butte, Colorado)

Temperature and precipitation variability in response to climate change affects water cycling of watersheds and can potentially impact water quality, water availability, elemental and molecular fluxes, and biogeochemical processes. Here we report the measurements of light stable isotopes and greenhouse gases analysis to a large multidisciplinary project addressing watershed function. Water hydrogen (H) and oxygen (O) stable isotopes of snow and water saturated soil were measured to constrain source contributions to groundwater and streamflow. In an effort to better constrain existing nitrogen and carbon cycling models in high-altitude ecosystems, carbon dioxide (CO2), nitrous oxide (N2O), methane (CH4) greenhouse gases concentrations and stable isotope signatures of CO2, dissolved inorganic carbon (DIC), and water were measured. This dataset includes one CSV file and one Excel file containing data including stable isotopes of carbon dioxide (CO2), dissolved inorganic carbon (DIC), water and concentrations of greenhouse gases. The CSV and Excel files contain the same information/data but have included both versions for accessibility.

54 ENVIRONMENTAL SCIENCES↗

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Underground sources of drinking water chemistry changes in response to potential CO 2 leakage

The purpose of this study was to quantify changes to underground sources of drinking water (USDW) quality in response to potential CO 2 leakage from geologic CO 2 sequestration (GCS) reservoirs. Here, the developed a framework of combined laboratory experiments and reactive transport simulations and used this framework to evaluate the Ogallala aquifer overlying the Farnsworth Unit (FWU), an active GCS site, as a case study. Using chemical reaction parameters obtained from laboratory experiments and numerical simulations, site-specific mechanisms of CO 2 -water-sediment interactions at the USDW aquifer were interpreted. Long-term risks of potential CO 2 leakage were then evaluated with field-scale numerical models using the regional hydrogeological characteristics and reaction parameters obtained from our experiments and simulations. Results suggest that carbonate mineral impurity and cation exchange are key mechanisms for interactions between CO 2 and the aquifer sediment. Additionally, for a large leakage rate of 0.1 % injection from one leaky well, the leakage plume might impact an area of 300 m in diameter and significantly affect the local water quality by changing pH and cation concentrations (e.g., Zn, Ba and Sr). After leakage ceases, the zone of impacted fluids would not migrate significantly in subsequent decades due to a low regional groundwater flowrate (for this case study). The relatively small area of impact might not be detected in a monitoring well given the broader spacing in a typical field scenario. Effective early leakage detection may require additional tools, e.g., borehole CO 2 movement, four-dimensional seismicity, CO 2 soil flux, samples from deeper aquifers, etc., to ensure effective leakage detection and long-term safety of GCS projects.

58 GEOSCIENCES↗

AmeriFlux US-CMW Charleston Mesquite Woodland

This is the AmeriFlux version of the carbon flux data for the site US-CMW Charleston Mesquite Woodland. Site Description - Site is a riparian mesquite woodland located along an old alluvial terrace along the banks of the San Pedro River in southeastern Arizona. Site climate is semiarid with monsoonal (Jul-Sep) rainfall. The mesquite trees access groundwater from an alluvial aquifer recharged elsewhere. Site water balance is not closed and ET greatly exceeds precipitation. The understory grasses and forbs are rainfall dependent. This groundwater-dependent ecosystem is far more ebulliant and productive than the upland vegetation around it (shrublands and grasslands).

Scott, Russell↗

AmeriFlux US-LS1 San Pedro River Lewis Springs Sacaton Grassland

This is the AmeriFlux version of the carbon flux data for the site US-LS1 San Pedro River Lewis Springs Sacaton Grassland. Site Description - This is a riparian grassland located along an old alluvial terrace along the banks of the entrenched San Pedro River in southeastern Arizona. Site climate is semiarid with monsoonal (Jul-Sep) rainfall. The grasses access groundwater (depth of 2-4 m)from an alluvial aquifer recharged elsewhere. Site water balance is not closed and ET greatly exceeds precipitation. The perennial bunchgrasses use both rainfall, groundwater and occasional flood water. This groundwater-using ecosystem is far more ebulliant and productive than the upland vegetation around it (shrublands and grasslands) outside of the San Pedro River corridor.

Scott, Russell↗

AmeriFlux US-LS2 San Pedro River Lewis Springs Savanna

This is the AmeriFlux version of the carbon flux data for the site US-LS2 San Pedro River Lewis Springs Savanna. Site Description - A riparian savanna located on an old alluvial terrace about 0.5 km from the banks of the entrenched San Pedro River in southeastern Arizona USA. Site climate is semiarid with monsoonal (Jul-Sep) rainfall. The trees access groundwater (depth of ~7 m)from an alluvial aquifer recharged elsewhere. Therefore, site water balance is not closed and ET greatly exceeds precipitation. The perennial bunchgrasses use primarily rainfall. This groundwater-using ecosystem is far more ebulliant and productive than the upland vegetation around it (shrublands and grasslands) outside of the San Pedro River corridor.

Scott, Russell↗

From legacy contamination to watershed systems science: a review of scientific insights and technologies developed through DOE-supported research in water and energy security

Abstract Water resources, including groundwater and prominent rivers worldwide, are under duress because of excessive contaminant and nutrient loads. To help mitigate this problem, the United States Department of Energy (DOE) has supported research since the late 1980s to improve our fundamental knowledge of processes that could be used to help clean up challenging subsurface problems. Problems of interest have included subsurface radioactive waste, heavy metals, and metalloids (e.g. uranium, mercury, arsenic). Research efforts have provided insights into detailed groundwater biogeochemical process coupling and the resulting geochemical exports of metals and nutrients to surrounding environments. Recently, an increased focus has been placed on constraining the exchanges and fates of carbon and nitrogen within and across bedrock to canopy compartments of a watershed and in river–floodplain settings, because of their important role in driving biogeochemical interactions with contaminants and the potential of increased fluxes under changing precipitation regimes, including extreme events. While reviewing the extensive research that has been conducted at DOE’s representative sites and testbeds (such as the Oyster Site in Virginia, Savannah River Site in South Carolina, Oak Ridge Reservation in Tennessee, Hanford in Washington, Nevada National Security Site in Nevada, Riverton in Wyoming, and Rifle and East River in Colorado), this review paper explores the nature and distribution of contaminants in the surface and shallow subsurface (i.e. the critical zone) and their interactions with carbon and nitrogen dynamics. We also describe state-of-the-art, scale-aware characterization approaches and models developed to predict contaminant fate and transport. The models take advantage of DOE leadership-class high-performance computers and are beginning to incorporate artificial intelligence approaches to tackle the extreme diversity of hydro-biogeochemical processes and measurements. Recognizing that the insights and capability developments are potentially transferable to many other sites, we also explore the scientific implications of these advances and recommend future research directions.

54 ENVIRONMENTAL SCIENCES↗

AmeriFlux FLUXNET-1F US-CMW Charleston Mesquite Woodland

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-CMW Charleston Mesquite Woodland. This is the FLUXNET version of the carbon flux data for the site US-CMW Charleston Mesquite Woodland produced by applying the standard ONEFlux (1F) software. Site Description - Site is a riparian mesquite woodland located along an old alluvial terrace along the banks of the San Pedro River in southeastern Arizona. Site climate is semiarid with monsoonal (Jul-Sep) rainfall. The mesquite trees access groundwater from an alluvial aquifer recharged elsewhere. Site water balance is not closed and ET greatly exceeds precipitation. The understory grasses and forbs are rainfall dependent. This groundwater-dependent ecosystem is far more ebulliant and productive than the upland vegetation around it (shrublands and grasslands).

Scott, Russell [USDA-ARS]↗

Fate of Water Pumped from Underground and Contributions to Sea Level Rise

The contributions from terrestrial water sources to sea-level rise, other than ice caps and glaciers, are highly uncertain and heavily debated1-5. Recent assessments indicate that groundwater depletion (GWD) may become the most important positive terrestrial contribution6-10 over the next 50 years, probably equal in magnitude to the current contributions from glaciers and ice caps6. However, the existing estimates assume that nearly 100% of groundwater extracted eventually ends up in the oceans. Owing to limited knowledge of the pathways and mechanisms governing the ultimate fate of pumped groundwater, the relative fraction of global GWD that contributes to sea-level rise remains unknown. Here, using a coupled climate-hydrological model11,12 simulation, we show that only 80% of GWDends up in the ocean. An increase in runo to the ocean accounts for roughly two-thirds, whereas the remainder results from the enhanced net flux of precipitation minus evaporation over the ocean, due to increased atmospheric vapour transport from the land to the ocean. The contribution of GWD to global sea-level rise amounted to 0.02 (+/- 0.004)mm yr(sup-1) in 1900 and increased to 0.27 (+/- 0.04)mm yr(sup-1) in 2000. This indicates that existing studies have substantially overestimated the contribution of GWD to global sea-level rise by a cumulative amount of at least 10 mm during the twentieth century and early twenty-first century. With other terrestrial water contributions included, we estimate the net terrestrial water contribution during the period 1993-2010 to be +0.12 +/-0.04)mm yr(sup-1), suggesting that the net terrestrialwater contribution reported in the IPCC Fifth Assessment Report report is probably overestimated by a factor of three.

water↗

Meteorological, reference evapotranspiration and estimated transpiration data, July 2020-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of meteorological data (air temperature, relative humidity, barometric pressure, precipitation, wind speed, soil heat flux, incoming shortwave radiation, incoming longwave radiation, outgoing shortwave radiation and outgoing longwave radiation) from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package also includes reference evapotranspiration calculated via the Penman-Monteith equation and an estimate of willow (Salix spp.) evapotranspiration using the method of Irmak et al., J. Hydrology, 2013. Due to two malfunctioning sensors (incoming shortwave radiation from Oct 2020-May 2021 and incoming longwave radiation from May 2021-March 2022), observed net radiation and calculated evapotranspiration are only included from July-October 2020. However, users may be able to extend calculated evapotranspiration using estimates of these two parameters from other sources. Both sensors have been replaced and data beyond 2021 be uploaded to this repository when it is available.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 3 "*data.csv" files, two for meteorological data for 2020 and 2021 and one for reference evapotranspiration and transpiration data from 2020 ("SLT_OBJ2_PenmanMonteith_2020_data.csv"). Each "*data.csv" file of meteorological data has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding to each measured parameter. We also include a Jupyter notebook ("PenmanMonteithET_Calculations.ipynb") that details the assumptions and equations used to calculate daily and hourly reference evapotranspiration and daily willow transpiration from raw meteorological data.

54 ENVIRONMENTAL SCIENCES↗