skip to main content


Title: Trend Analyses of Baseflow and BFI for Undisturbed Watersheds in Michigan—Constraints from Multi-Objective Optimization
Documenting how ground- and surface water systems respond to climate change is crucial to understanding water resources, particularly in the U.S. Great Lakes region, where drastic temperature and precipitation changes are observed. This study presents baseflow and baseflow index (BFI) trend analyses for 10 undisturbed watersheds in Michigan using (1) multi-objective optimization (MOO) and (2) modified Mann–Kendall (MK) tests corrected for short-term autocorrelation (STA). Results indicate a variability in mean baseflow (0.09–8.70 m3/s) and BFI (67.9–89.7%) that complicates regional-scale extrapolations of groundwater recharge. Long-term (>60 years) MK trend tests indicate a significant control of total precipitation (P) and snow- to rainfall transitions on baseflow and BFI. In the Lower Peninsula Rifle River watershed, increasing P and a transition from snow- to rainfall has increased baseflow at a lower rate than streamflow; an overall pattern that may contribute to documented flood frequency increases. In the Upper Peninsula Ford River watershed, decreasing P and a transition from rain- to snowfall had no significant effects on baseflow and BFI. Our results highlight the value of an objectively constrained BFI parameter for shorter-term (<50 years) hydrologic trend analysis because of a lower STA susceptibility.  more » « less
Award ID(s):
1936671
NSF-PAR ID:
10218532
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Water
Volume:
13
Issue:
4
ISSN:
2073-4441
Page Range / eLocation ID:
564
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Snow surveys in two Arctic watersheds located in Alaska, USA, provide 32 years of spatially distributed snow water equivalent (SWE) and snow depth observations. Annual snow surveys from the Imnavait Creek (20,036 measurements from 1985 to 2017) and Upper Kuparuk River (5,804 measurements from 1997 to 2017) watersheds were conducted to capture end‐of‐winter snow accumulation. The average end‐of‐winter SWE in the Upper Kuparuk River watershed (102 ± 29 mm) is consistently less than the Imnavait Creek watershed (130 ± 34 mm) during the common period of record (1997–2017). The average end‐of‐winter SWE in both watersheds indicates a positive trend. Comparison of SWE records with cumulative solid precipitation measured at the Imnaviat [sic] SNOTEL site highlights the undercatch of gauge precipitation and difference in long‐term trends. In this paper, we present a historic overview of data collection, discuss data accuracy, and point out advantages and limitations associated with ground‐based snow measurements in remote Arctic locations. As new methods and techniques of measuring SWE and solid precipitation become available, the presented data set will provide a historic perspective for new observations and will quantitatively relate current or future snow conditions to those that have occurred since the late twentieth century.

     
    more » « less
  2. Abstract Some of the largest climatic changes in the Arctic have been observed in Alaska and the surrounding marginal seas. Near-surface air temperature (T2m), precipitation ( P ), snowfall, and sea ice changes have been previously documented, often in disparate studies. Here, we provide an updated, long-term trend analysis (1957–2021; n = 65 years) of such parameters in ERA5, NOAA U.S. Climate Gridded Dataset (NClimGrid), NOAA National Centers for Environmental Information (NCEI) Alaska climate division, and composite sea ice products preceding the upcoming Fifth National Climate Assessment (NCA5) and other near-future climate reports. In the past half century, annual T2m has broadly increased across Alaska, and during winter, spring, and autumn on the North Slope and North Panhandle (T2m > 0.50°C decade −1 ). Precipitation has also increased across climate divisions and appears strongly interrelated with temperature–sea ice feedbacks on the North Slope, specifically with increased (decreased) open water (sea ice extent). Snowfall equivalent (SFE) has decreased in autumn and spring, perhaps aligned with a regime transition of snow to rain, while winter SFE has broadly increased across the state. Sea ice decline and melt-season lengthening also have a pronounced signal around Alaska, with the largest trends in these parameters found in the Beaufort Sea. Alaska’s climatic changes are also placed in context against regional and contiguous U.S. air temperature trends and show ∼50% greater warming in Alaska relative to the lower-48 states. Alaska T2m increases also exceed those of any contiguous U.S. subregion, positioning Alaska at the forefront of U.S. climate warming. Significance Statement This study produces an updated, long-term trend analysis (1957–2021) of key Alaska climate parameters, including air temperature, precipitation (including snowfall equivalent), and sea ice, to inform upcoming climate assessment reports, including the Fifth National Climate Assessment (NCA5) scheduled for publication in 2023. Key findings include widespread annual and seasonal warming with increased precipitation across much of the state. Winter snowfall has broadly increased, but spring and autumn snowfalls have decreased as rainfall increased. Autumn warming and precipitation increases over the North Slope, in particular, appear related to decreased sea ice coverage in the Beaufort Sea and Chukchi Seas. These trends may result from interrelated processes that accelerate Alaska climate changes relative to those of the contiguous United States. 
    more » « less
  3. Abstract

    For evaluating the climatic and landscape controls on long‐term baseflow, baseflow index (BFI, defined as the ratio of baseflow to streamflow) and baseflow coefficient (BFC, defined as the ratio of baseflow to precipitation) are formulated as functions of climate aridity index, storage capacity index (defined as the ratio of average soil water storage capacity to precipitation), and a shape parameter for the spatial variability of storage capacity. The derivation is based on the two‐stage partitioning framework and a cumulative distribution function for storage capacity. Storage capacity has a larger impact on BFI than on BFC. When storage capacity index is smaller than 1, BFI is less sensitive to storage capacity index in arid regions compared to that in humid regions; whereas, when storage capacity index is larger than 1, BFI is less sensitive to storage capacity index in humid regions. The impact of storage capacity index on BFC is only significant in humid regions. The shape parameter plays an important role on fast flow generation at the first‐stage partitioning in humid regions and baseflow generation at the second‐stage partitioning in arid regions. The derived formulae were applied to more than 400 catchments where storage capacity index was found to follow a logarithmic function with climate aridity index. The role of climate forcings at finer timescales on baseflow were quantified, indicating that seasonality in climate forcings has a significant control especially on BFI.

     
    more » « less
  4. null (Ed.)
    Precipitation occurs in two basic forms defined as liquid state and solid state. Different from rain-fed watershed, modeling snow processes is of vital importance in snow-dominated watersheds. The seasonal snowpack is a natural water reservoir, which stores snow water in winter and releases it in spring and summer. The warmer climate in recent decades has led to earlier snowmelt, a decline in snowpack, and change in the seasonality of river flows. The Soil and Water Assessment Tool (SWAT) could be applied in the snow-influenced watershed because of its ability to simultaneously predict the streamflow generated from rainfall and from the melting of snow. The choice of parameters, reference data, and calibration strategy could significantly affect the SWAT model calibration outcome and further affect the prediction accuracy. In this study, SWAT models are implemented in four upland watersheds in the Tulare Lake Basin (TLB) located across the Southern Sierra Nevada Mountains. Three calibration scenarios considering different calibration parameters and reference datasets are applied to investigate the impact of the Parallel Energy Balance Model (ParBal) snow reconstruction data and snow parameters on the streamflow and snow water-equivalent (SWE) prediction accuracy. In addition, the watershed parameters and lapse rate parameters-led equifinality is also evaluated. The results indicate that calibration of the SWAT model with respect to both streamflow and SWE reference data could improve the model SWE prediction reliability in general. Comparatively, the streamflow predictions are not significantly affected by differently lumped calibration schemes. The default snow parameter values capture the extreme high flows better than the other two calibration scenarios, whereas there is no remarkable difference among the three calibration schemes for capturing the extreme low flows. The watershed and lapse rate parameters-induced equifinality affects the flow prediction more (Nash-Sutcliffe Efficiency (NSE) varies between 0.2–0.3) than the SWE prediction (NSE varies less than 0.1). This study points out the remote-sensing-based SWE reconstruction product as a promising alternative choice for model calibration in ungauged snow-influenced watersheds. The streamflow-reconstructed SWE bi-objective calibrated model could improve the prediction reliability of surface water supply change for the downstream agricultural region under the changing climate. 
    more » « less
  5. Abstract

    Stable isotope variations are extremely useful for flow partitioning within the hydrologic cycle but remain poorly understood throughout the tropics, particularly in watersheds with rapidly infiltrating soils, such as Andisols in Central America. This study examines the fluctuations of stable isotope ratios (δ18O and δ2H) in the hydrologic components of a tropical coffee agroforestry watershed (~1 km2) with Andisol soils in Costa Rica. Samples were collected in precipitation, groundwater, springs, and stream water over 2 years. The local meteoric water line for the study site was δ2H = 8.5 δ18O + 18.02 (r2 = 0.97,n = 198). The isotope ratios in precipitation exhibited an enriched trend during the dry season and a notable depletion at the beginning of the wet season. The δ18O compositions in groundwater (average = −6.4‰, σ = 0.7) and stream water (average = −6.7‰, σ = 0.6) were relatively stable over time, and both components exhibited more enriched values in 2013, which was the drier year. No strong correlation was observed between the isotope ratios and the precipitation amount at the event or daily time‐step, but a correlation was observed on a monthly scale. Stream water and base flow hydrograph separations based on isotope end‐member estimations showed that pre‐event water originating from base flow was prevalent. However, isotope data indicate that event water originating from springs appears to have been the primary driver of initial rises in stream flow and peak flows. These results indicate that isotope sampling improves the understanding of water balance components, even in a tropical humid location, where significant variations in rainfall challenge current modelling efforts. Further research using fine‐scale hydrometric and isotopic data would enhance understanding the processes driving spring flow generation in watersheds.

     
    more » « less