USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31Linked to original sources

Efficient hydrogeological characterization of remote stream corridors using drones

This project demonstrates the successful use of small unoccupied aircraft system (sUASs) for hydrogeological characterization of a remote stream reach in a rugged mountain terrain. Thermal infrared, visual imagery, and derived digital surface models are used to inform conceptual models of groundwater/surface‐water exchange and efficiently geolocate zones of preferential groundwater discharge that can be quantified using various ground‐based methodology.

Hydrological Processes↗

Long-term streamflow trends in Hawai‘i and implications for native stream fauna

Climate change has fundamentally altered the water cycle in tropical islands, which is a critical driver of freshwater ecosystems. To examine how changes in streamflow regime have impacted habitat quality for native migratory aquatic species, we present a 50‐year (1967–2016) analysis of hydrologic records in 23 unregulated streams across the five largest Hawaiian Islands. For each stream, flow was separated into direct run‐off and baseflow and high‐ and low‐flow statistics (i.e., Q10 and Q90) with ecologically important hydrologic indices (e.g., frequency of flooding and low flow duration) derived. Using Mann–Kendall tests with a running trend analysis, we determined the persistence of streamflow trends through time. We analysed native stream fauna from ~400 sites, sampled from 1992 to 2007, to assess species richness among islands and streams. Declines in streamflow metrics indicated a general drying across the islands. In particular, significant declines in low flow conditions (baseflows), were experienced in 57% of streams, compared with a significant decline in storm flow conditions for 22% of streams. The running trend analysis indicated that many of the significant downward trends were not persistent through time but were only significant if recent decades (1987–2016) were included, with an average decline in baseflow and run‐off of 10.90% and 8.28% per decade, respectively. Streams that supported higher native species diversity were associated with moderate discharge and baseflow index, short duration of low flows, and negligible downward trends in flow. A significant decline in dry season flows (May–October) has led to an increase in the number of no‐flow days in drier areas, indicating that more streams may become intermittent, which has important implications for mauka to makai (mountain to ocean) hydrological connectivity and management of Hawai'i's native migratory freshwater fauna.

Hawaii↗

Sources of inherent infiltration variability in postwildfire soils

An automated disc infiltrometer was developed to improve the measurements of soil hydraulic properties (saturated hydraulic conductivity and sorptivity) of soils affected by wildfire. Guidelines are given for interpreting curves showing cumulative infiltration as a function of time measured by the autodisc. The autodisc was used to measure the variability of these soil hydraulic properties in three different sample sets: (a) a reference soil consisting of a nonrepellent, uniform, fine sand; (b) soils with the same soil textural classification derived from the same bedrock geology but having different initial burn severities; and (c) soils from different bedrock geology but having the same burn severity. The autodisc infiltrometer had greater sampling rates and volume resolution when compared with the visual minidisc infiltrometer from previous studies. There was no statistical difference in the mean values measured using the autodisc and visual minidisc, but the variability of the autodisc measurements was significantly less than the visual minidisc for a given set of samples. The greatest variability of soil hydraulic properties in reference samples with uniform particle size was attributed to different pore geometries (coefficient of variation [COV] = 0.28–0.34). Unburned field samples (same soil type) with heterogeneous particle sizes had greater variability (COV = 0.57–0.78) than the reference samples. However, this basic variability decreased or remained constant in these field samples as burn severity increased. Additional sources of variability (COV = 0.53–1.99) were attributed to multiple layers resulting from ash or sediment deposition. Results indicate that resolving differences in soil hydraulic properties from different sites requires more than the common 10 random samples because of the multiple sources of variability.

Hydrological Processes↗

Changes in event‐based streamflow magnitude and timing after suburban development with infiltration‐based stormwater management

Green stormwater infrastructure implementation in urban watersheds has outpaced our understanding of practice effectiveness on streamflow response to precipitation events. Long‐term monitoring of experimental urban watersheds in Clarksburg, Maryland, USA, provided an opportunity to examine changes in event‐based streamflow metrics in two treatment watersheds that transitioned from agriculture to suburban development with a high density of infiltration‐focused stormwater control measures (SCMs). Urban Treatment 1 has predominantly single family detached housing with 33% impervious cover and 126 SCMs. Urban Treatment 2 has a mix of single family detached and attached housing with 44% impervious cover and 219 SCMs. Differences in streamflow‐event magnitude and timing were assessed using a before‐after‐control‐reference‐impact design to compare urban treatment watersheds to a forested control and an urban control with detention‐focused SCMs. Streamflow and precipitation events were identified from 14 years of sub‐daily monitoring data with an automated approach to characterize peak streamflow, runoff yield, runoff ratio, streamflow duration, time to peak, rise rate, and precipitation depth for each event. Results indicated that streamflow magnitude and timing were altered by urbanization in the urban treatment watersheds, even with SCMs treating 100% of the impervious area. The largest hydrologic changes were observed in streamflow magnitude metrics, with greater hydrologic change in Urban Treatment 2 compared to Urban Treatment 1. While streamflow changes were observed in both urban treatment watersheds, SCMs were able to mitigate peak flows and runoff volumes compared to the urban control. The urban control had similar impervious cover to Urban Treatment 2, but Urban Treatment 2 had more than twice the precipitation depth needed to initiate a flow response and lower median peak flow and runoff yield for events less than 20 mm. Differences in impervious cover between the Urban Treatment watersheds appeared to be a large driver of differences in streamflow response, rather than SCM density. Overall, use of infiltration‐focused SCMs implemented at a watershed‐scale did provide enhanced attenuation of peak flow and runoff volumes compared to centralized‐detention SCMs.

Maryland↗

Quantitative guidance for efficient vertical flow measurements at the sediment-water interface using temperature-depth profiles

Upward discharge to surface water bodies can be quantified using analytical models based on temperature-depth (T-z) profiles. The use of sediment T-z profiles is attractive as discharge estimates can be obtained using point-in-time data that are collected inexpensively and rapidly. Previous studies have identified that T-z methods can only be applied at times of the year when there is significant difference between the streambed-water interface and deeper sediment temperatures (e.g., winter and summer). However, surface water temperatures also vary diurnally, and the influence of these variations on discharge estimates from T-z methods is poorly understood. For this study, synthetic T-z profiles were generated numerically using measured streambed interface temperature data to assess the influence of diurnal temperature variations on discharge estimation and provide insight into the suitable application of T-z methods. Results show that the time of day of data collection can have a substantial influence on vertical flux estimates using T-z methods. For low groundwater discharge fluxes (e.g. 0.1 m d-1), daily transience in streambed temperatures led to relatively large errors in estimated flow magnitude and direction. For higher discharge fluxes (1.5 m d-1), the influence of transient streambed temperatures on discharge estimates was strongly reduced. Discharge estimates from point-in-time T-z profiles were most accurate when the uppermost point in the T-z profile was near the bed interface daily mean (two time periods daily). Where temperature time series data are available, daily averaged T-z profiles can produce accurate discharge estimates across a wide range of discharge rates. Seasonality in shallow groundwater temperature generally had a negligible influence on vertical flow estimates. These findings can be used to plan field campaigns and provide guidance on the optimal application of T-z methods to quantify vertical groundwater discharge to surface water bodies.

Hydrological Processes↗

Chronic and episodic acidification of streams along the Appalachian Trail corridor, eastern United States

Acidic atmospheric deposition has adversely affected aquatic ecosystems globally. As emissions and deposition of sulfur (S) and nitrogen (N) have declined in recent decades across North America and Europe, ecosystem recovery is evident in many surface waters. However, persistent chronic and episodic acidification remain important concerns in vulnerable regions. We evaluated acidification in 269 headwater streams during 2010–2012 along the Appalachian Trail (AT) that transits several ecoregions and is located downwind of high levels of S and N emission sources. Discharge was estimated by matching sampled streams to those of a nearby gaged stream and assuming equivalent daily mean flow percentiles. Charge balance acid‐neutralizing capacity (ANC) values were adjusted to the 15th (Q15) and 85th flow percentiles (Q85) by applying the ANC/discharge slope among sample pairs collected at each stream. A site‐based approach was applied to streams sampled twice or more and a second regression‐based approach to streams sampled once to estimate episodic acidification magnitudes as the ANC difference from Q15 to Q85. Streams with ANC <0 μeq/L doubled from 16% to 32% as discharge increased from Q15 to Q85 according to the site‐based approach. The proportion of streams with ANC <0 μeq/L at low flow and high flow decreased from north to south. Base cation dilution explained the greatest amount of episodic acidification among streams and variation in sulfate (SO 4 2− ) concentrations was a secondary explanatory variable. Episodic SO 4 2− patterns varied geographically with dilution dominant in northern streams underlain by soils developed in glacial sediment and increased concentrations dominant in southern streams with older, highly weathered soils. Episodic acidification increased as low‐flow ANC increased, exceeding 90 μeq/L in 25% of streams. Episodic increases in ANC were the dominant pattern in streams with low‐flow ANC values <30 μeq/L. Chronic and episodic acidification remain an ecological concern among AT streams. The approach developed here could be applied to estimate the magnitude and extent of chronic and episodic acidification in other regions recovering from decreasing levels of atmospheric S and N deposition.

Connecticut, Georgia, Maine, Massachusetts, Maryla↗

Effect of an environmental flow on vegetation growth and health using ground and remote sensing metrics

Understanding the effectiveness of environmental flow deliveries along rivers requires monitoring vegetation. Monitoring data are often collected at multiple spatial scales. For riparian vegetation, optical remote sensing methods can estimate growth responses at the riparian corridor scale, and field‐based measures can quantify species composition; however, the extent to which these different measures are duplicative or complementary is important to understand when planning monitoring programmes with limited resources. In this study, we analysed riparian vegetation growth in the delta of the Colorado River in response to an experimental pulse flow. Our goal was to compare ground‐based measurements of vegetation structure and composition with satellite‐based Landsat radiometric variables, such as the normalized difference vegetation index (NDVI). We made this comparison in 21 transects following the delivery of 131.8 million cubic meters (mcm) of water in the stream channel during the spring of 2014 as a pulse flow and 38.4 mcm as base flows. Vegetation cover increased 14% and NDVI increased 0.02 (15%) by October 2015, and both variables returned to pre‐pulse flow values in October 2016. Observed changes in vegetation structure and composition did not persist after the second year. The highest increase in vegetation cover in October 2014 and October 2015 resulted from species that could respond rapidly to additional water such as reeds ( Arundo donax and Phragmites australis ), cattail ( Typha domingensis ), and herbaceous plants. Dominant shrubs, saltcedar ( Tamarix spp.) and arrowweed ( Pluchea sericea ), both indicative of nonrestored habitats showed variable increases in cover, and native trees ( Salicaceae family) presented low increases (1%). The strong NDVI–vegetation cover relationship indicates that NDVI is appropriate to detect changes at the riparian corridor scale but needs to be complemented with ground data to determine the contributions by different species to the observed trends.

Colorado River↗

Evaluation of hydrologic impact of an irrigation curtailment program in the Upper Klamath Lake Basin using Landsat satellite data

Upper Klamath Lake (UKL) is the source of the Klamath river that flows through southern Oregon and northern California. The UKL basin is home to two endangered species and provides water for 81,000+ ha (200,000+ acres) of irrigation on the United States Bureau of Reclamation (USBR) Klamath Project located downstream of the UKL basin. Irrigated agriculture also occurs along the tributaries to UKL. During 2013–2016, water right calls resulted in various levels of curtailment of irrigation diversions from the tributaries to UKL. However, information on the extent of curtailment, how much irrigation water was saved, and its impact on the UKL is unknown. In this study, we combined Landsat-based actual evapotranspiration (ETa) data obtained from the Operational Simplified Surface Energy Balance (SSEBop) model with gridded precipitation and USGS station discharge data to evaluate the hydrologic impact of the curtailment program. Analysis was performed for five base years (2004, 2006, 2008-2010) and four target years (2013-2016) over irrigated areas above UKL. Our results indicated that the impact of the curtailment program over the June to September time-period was highest during 2013 and declined in each of the following years. The total on-field water savings were approximately 60 hm3 in 2013 and 2014, 44 hm3 in 2015, and 32 hm3 in 2016. The instream water flow change or extra water available (EWA) were found at 92, 68, 45, and 26 hm3 respectively for 2013, 2014, 2015 and 2016. Most water savings came from pasture and wetlands. Alfalfa showed the most decline in water use among grain crops. The resulting EWA from the curtailment contributed to a maximum of 19% of the lake inflows and 50% of the lake volume. This study presents the use of Landsat-based ETa and other remote sensing datasets for evaluating water-related impacts of the irrigation curtailment program.

Upper Klamath Lake Basin↗

Hillslope groundwater discharges provide localized ecosystem buffers from regional PFAS contamination in a gaining coastal stream

Emerging groundwater contaminants such as per- and polyfluoroalkyl substances (PFAS) may impact surface-water quality and groundwater-dependent ecosystems of gaining streams. Although complex near-surface hydrogeology of stream corridors challenges sampling efforts, recent advances in heat tracing of discharge zones enable efficient and informed data collection. For this study we used a combination of streambed temperature push-probe and thermal infrared methods to guide a discharge-zone-oriented sample collection along approximately 6 km of a coastal trout stream on Cape Cod, MA where groundwater discharge constitutes approximately 95% of total streamflow. Eight surface-water locations and discharging groundwater from 24 streambed and bank seepages were analyzed for dissolved oxygen, specific conductance, stable water isotopes, and a range of PFAS compounds which are contaminants of emerging concern in aquatic environments. The results indicate a complex system of groundwater discharge source flowpaths, where the sum of concentrations of six PFAS compounds (Environmental Protection Agency third Unregulated Contaminant Monitoring Rule UCMR 3) showed a median concentration of 52 331 (SD) ng/L with two higher outliers and three discharges with non-detection of PFAS. Higher UCMR 3 PFAS concentration was related -0.66 (Spearman Rank, p<0.001) to discharging groundwater that showed an evaporative signature (deuterium excess), indicating flow through at least one upgradient kettle lake. Therefore, more regional groundwater flowpaths originating from outside the local river corridor tended to show higher PFAS concentrations as evaluated at their respective discharge zones. Conversely, UCMR 3 PFAS concentrations were typically low at discharges that did not indicate evaporation and were adjacent to steep hillslopes and, therefore, were classified as locally recharged groundwater. Previous research at this stream found that the native brook trout favor discharge points of groundwater recharged on local hillslopes for spawning, likely in response to generally higher levels of dissolved oxygen compared to discharge zones located further away from hillslopes. Our study shows that the trout may thereby be avoiding emerging contaminants such as PFAS in groundwater recharged farther from the stream.

Massachusetts↗

Vegetation‐groundwater dynamics at a former uranium mill site following invasion of a biocontrol agent: A time series analysis of Landsat normalized difference vegetation index data

Because groundwater recharge in dry regions is generally low, arid and semiarid environments have been considered well‐suited for long‐term isolation of hazardous materials (e.g., radioactive waste). In these dry regions, water lost (transpired) by plants and evaporated from the soil surface, collectively termed evapotranspiration (ET), is usually the primary discharge component in the water balance. Therefore, vegetation can potentially affect groundwater flow and contaminant transport at waste disposal sites. We studied vegetation health and ET dynamics at a Uranium Mill Tailings Radiation Control Act (UMTRCA) disposal site in Shiprock, New Mexico, where a floodplain alluvial aquifer was contaminated by mill effluent. Vegetation on the floodplain was predominantly deep‐rooted, non‐native tamarisk shrubs ( Tamarix sp.). After the introduction of the tamarisk beetle ( Diorhabda sp.) as a biocontrol agent, the health of the invasive tamarisk on the Shiprock floodplain declined. We used Landsat normalized difference vegetation index (NDVI) data to measure greenness and a remote sensing algorithm to estimate landscape‐scale ET along the floodplain of the UMTRCA site in Shiprock prior to (2000–2009) and after (2010–2018) beetle establishment. Using groundwater level data collected from 2011 to 2014, we also assessed the role of ET in explaining seasonal variations in depth to water of the floodplain. Growing season scaled NDVI decreased 30% ( p < .001), while ET decreased 26% from the pre‐ to post‐beetle period and seasonal ET estimates were significantly correlated with groundwater levels from 2011 to 2014 ( r 2 = .71; p = .009). Tamarisk greenness (a proxy for health) was significantly affected by Diorhabda but has partially recovered since 2012. Despite this, increased ET demand in the summer/fall period might reduce contaminant transport to the San Juan River during this period.

New Mexico↗

Stormwater control impacts on runoff volume and peak flow: A meta-analysis of watershed modelling studies

Decades of research has concluded that the percent of impervious surface cover in a watershed is strongly linked to negative impacts on urban stream health. Recently, there has been a push by municipalities to offset these effects by installing structural stormwater control measures (SCMs), which are landscape features designed to retain and reduce runoff to mitigate the effects of urbanisation on event hydrology. The goal of this study is to build generalisable relationships between the level of SCM implementation in urban watersheds and resulting changes to hydrology. A literature review of 185 peer‐reviewed studies of watershed‐scale SCM implementation across the globe was used to identify 52 modelling studies suitable for a meta‐analysis to build statistical relationships between SCM implementation and hydrologic change. Hydrologic change is quantified as the percent reduction in storm event runoff volume and peak flow between a watershed with SCMs relative to a (near) identical control watershed without SCMs. Results show that for each additional 1% of SCM‐mitigated impervious area in a watershed, there is an additional 0.43% reduction in runoff and a 0.60% reduction in peak flow. Values of SCM implementation required to produce a change in water quantity metrics were identified at varying levels of probability. For example, there is a 90% probability (high confidence) of at least a 1% reduction in peak flow with mitigation of 33% of impervious surfaces. However, as the reduction target increases or mitigated impervious surface decreases, the probability of reaching the reduction target also decreases. These relationships can be used by managers to plan SCM implementation at the watershed scale.

Hydrological Processes↗

Effect of spatial resolution of satellite images on estimating the greenness and evapotranspiration of urban green spaces

Urban green spaces (UGS), like most managed land covers, are getting progressively affected by water scarcity and drought. Preserving, restoring and expanding UGS require sustainable management of green and blue water resources to fulfil evapotranspiration (ET) demand for green plant cover. The heterogeneity of UGS with high variation in their microclimates and irrigation practices builds up the complexity of ET estimation. In oversized UGS, areas too large to be measured with in situ ET methods, remote sensing (RS) approaches of ET measurement have the potential to estimate the actual ET. Often in situ approaches are not feasible or too expensive. We studied the effects of spatial resolution using different satellite images, with high‐, medium‐ and coarse‐spatial resolutions, on the greenness and ET of UGS using Vegetation Indices (VIs) and VI‐based ET, over a 780‐ha urban park in Adelaide, Australia. We validated ET with the ground‐based ET method of Soil Water Balance. Three sets of imagery from WorldView2, Landsat and MODIS, and three VIs including the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Enhanced Vegetation Index 2 (EVI2), were used to assess long‐term changes of VIs and ET calculated from the different imagery acquired for this study (2011–2018). We found high correspondence between ET‐MODIS and ET‐Landsat (R2 > 0.99 for all VIs). Landsat‐VIs captured the seasonal changes of greenness better than MODIS‐VIs. We used artificial neural network (ANN) to relate the RS‐ET and ground data, and ET‐MODIS (EVI2) showed the highest correlation (R2 = 0.95 and MSE =0.01 for validation). We found a strong relationship between RS‐ET and in situ measurements, even though it was not explicable by simple regressions; black box models helped us to explore their correlation. The methodology used in this research makes a strong case for the value of remote sensing in estimating and managing ET of green spaces in water‐limited cities.

Adelaide↗

Parameter estimation for multiple post-wildfire hydrologic models

Predictions of post‐wildfire flooding and debris flows are needed, typically with short lead times. Measurements of soil‐hydraulic properties necessary for model parameterization are, however, seldom available. This study quantified soil‐hydraulic properties, soil‐water retention, and selected soil physical properties within the perimeter of the 2017 Thomas Fire in California. The Thomas Fire burn scar produced catastrophic debris flows in January 2018, highlighting the need for improved prediction capability. Soil‐hydraulic properties were also indirectly estimated using relations tied to soil‐water retention. These measurements and estimates are examined in the context of parameterizing post‐wildfire hydrologic models. Tension infiltrometer measurements showed significant decreases ( p < .05) in field‐saturated hydraulic conductivity ( K fs ) and sorptivity ( S ) in burned areas relative to unburned areas. Wildfire effects on soil water‐retention were dominated by significant decreases in saturated soil‐water content ( θ S ). The van Genuchten parameters α , N , and residual water content did not show significant wildfire effects. The impacts of the wildfire on hydraulic and physical soil properties were greatest in the top 1 cm, emphasizing that measurements of post‐fire soil properties should focus on the near‐surface. Reductions in K fs , θ s , and soil‐water retention in burned soils were attributed to fire‐induced decreases in soil structure evidenced by increases in dry bulk density. Sorptivity reductions in burned soils were attributed to increases in soil‐water repellency. Rapid post‐fire assessments of flash flood and debris flow hazards using physically‐based hydrologic models are facilitated by similarities between K fs , S , and the Green–Ampt wetting front potential ( ψ f ) with measurements at other southern CA burned sites. We suggest that ratios of burned to unburned K fs (0.37), S (0.36), and ψ f (0.66) could be used to scale unburned values for model parameterization. Alternatively, typical burned values ( K fs = 20 mm hr −1 ; S = 6 mm hr −0.5 ; ψ f = 1.6 mm) could be used for model parameterization.

California↗

Hillslopes in humid-tropical climates aren’t always wet: Implications for hydrologic response and landslide initiation in Puerto Rico, USA

The devastating impacts of the widespread flooding and landsliding in Puerto Rico following the September 2017 landfall of Hurricane Maria highlight the increasingly extreme atmospheric disturbances and enhanced hazard potential in mountainous humid‐tropical climate zones. Long‐standing conceptual models for hydrologically driven hazards in Puerto Rico posit that hillslope soils remain wet throughout the year, and therefore, that antecedent soil wetness imposes a negligible effect on hazard potential. Our post‐Maria in situ hillslope hydrologic observations, however, indicate that while some slopes remain wet throughout the year, others exhibit appreciable seasonal and intra‐storm subsurface drainage. Therefore, we evaluated the performance of hydro‐meteorological (soil wetness and rainfall) versus intensity‐duration (rainfall only) hillslope hydrologic response thresholds that identify the onset of positive pore‐water pressure, a predisposing factor for widespread slope instability in this region. Our analyses also consider the role of soil‐water storage and infiltration rates on runoff generation, which are relevant factors for flooding hazards. We found that the hydro‐meteorological thresholds outperformed intensity‐duration thresholds for a seasonally wet, coarse‐grained soil, although they did not outperform intensity‐duration thresholds for a perennially wet, fine‐grained soil. These end‐member soils types may also produce radically different stormflow responses, with subsurface flow being more common for the coarse‐grained soils underlain by intrusive rocks versus infiltration excess and/or saturation excess for the fine‐grained soils underlain by volcaniclastic rocks. We conclude that variability in soil‐hydraulic properties, as opposed to climate zone, is the dominant factor that controls runoff generation mechanisms and modulates the relative importance of antecedent soil wetness for our hillslope hydrologic response thresholds.

Puerto Rico↗

Flow‐ecology modelling to inform reservoir releases for riparian restoration and management

Linked hydrologic, hydraulic, and ecological models can facilitate planning and implementing water releases from reservoirs to achieve ecological objectives along rivers. We applied a flow‐ecology model, the Ecosystem Functions Model (HEC‐EFM), to the Bill Williams River in southwestern USA to estimate areas suitable for recruitment of riparian tree seedlings in the context of managing flow releases from a large dam for riparian restoration. Ecological variables in the model included timing of seed dispersal, tolerable rates of flow recession, and tolerable duration of inundation following germination and early seedling establishment for native Fremont cottonwood and Goodding's willow, and non‐native tamarisk. Hydrological variables included peak flow timing, rate of flow recession following the peak, and duration of inundation. A one‐dimensional hydraulic model was applied to estimate stage‐discharge relationships along ~58 river kilometres. We then used HEC‐EFM to apply relationships between seedling ecology and streamflow to link hydrological dynamics with ecological response. We developed and validated HEC‐EFM based on an examination of seedling recruitment following an experimental flow release from Alamo Dam in spring 2006. The model predicted the largest area of potential recruitment for cottonwood (280–481 ha), with smaller areas predicted for willow (174–188 ha) and tamarisk (59–60 ha). Correlations between observed and predicted patches with successful seedling recruitment for areas within 40 m of the main channel ranged from 0.66 to 0.94. Finally, we examined arrays of hydrographs to identify which are most conducive to seedling recruitment along the river, given different combinations of peak flow, recession rate, and water volume released. Similar application of this model could be useful for informing reservoir management in the context of riparian restoration along other rivers facing similar challenges.

Arizona↗

Climate sensitivity of water use by riparian woodlands at landscape scales

Semi‐arid riparian woodlands face threats from increasing extractive water demand and climate change in dryland landscapes worldwide. Improved landscape‐scale understanding of riparian woodland water use (evapotranspiration, ET) and its sensitivity to climate variables is needed to strategically manage water resources, as well as to create successful ecosystem conservation and restoration plans for potential climate futures. In this work, we assess the spatial and temporal variability of Cottonwood ( Populus fremontii )‐Willow ( Salix gooddingii ) riparian gallery woodland ET and its relationships to vegetation structure and climate variables for 80 km of the San Pedro River corridor in southeastern Arizona, USA, between 2014 and 2019. We use a novel combination of publicly available remote sensing, climate and hydrological datasets: cloud‐based Landsat thermal remote sensing data products for ET (Google Earth Engine EEFlux), Landsat multispectral imagery and field data‐based calibrations to vegetation structure (leaf‐area index, LAI), and open‐source climate and hydrological data. We show that at landscape scales, daily ET rates (6–10 mm day −1 ) and growing season ET totals (400–1,400 mm) matched rates of published field data, and modelled reach‐scale average LAI (0.80–1.70) matched lower ranges of published field data. Over 6 years, the spatial variability of total growing season ET (CV = 0.18) exceeded that of temporal variability (CV = 0.10), indicating the importance of reach‐scale vegetation and hydrological conditions for controlling ET dynamics. Responses of ET to climate differed between perennial and intermittent‐flow stream reaches. At perennial‐flow reaches, ET correlated significantly with temperature, whilst at intermittent‐flow sites ET correlated significantly with rainfall and stream discharge. Amongst reaches studied in detail, we found positive but differing logarithmic relationships between LAI and ET. By documenting patterns of high spatial variability of ET at basin scales, these results underscore the importance of accurately accounting for differences in woodland vegetation structure and hydrological conditions for assessing water‐use requirements. Results also suggest that the climate sensitivity of ET may be used as a remote indicator of subsurface water resources relative to vegetation demand, and an indicator for informing conservation management priorities.

Arizona↗