The main objective of this study was to use an uncertainty version of a widely used monthly time step, semi-distributed model (the Pitman model) to explore the equifinalities in the way in which the main hydrological processes are simulated and any identifiable linkages with uncertainties in the available observational data. The study area is the Zambezi River basin and 18 gauged sub-basins have been included in the analyses. Unfortunately, it is not generally possible to quantify some of the observational uncertainties in such a data scarce area and mostly we are limited to identifying where these data are clearly deficient (i.e. erroneous or non-representative). The overall conclusion is that the equifinalities in the model are hugely dominant in terms of the uncertainties in the relative occurrence of different runoff generating processes, although water use uncertainties in the semi-arid parts of the basin can contribute to these uncertainties. The identification of landscape features that suggest the occurrence of saturation excess surface runoff provides some information to constrain the model. Improved independent estimates of groundwater recharge is also identified as a key source of observational data that would help a great deal in constraining the model parameter space and therefore reducing some of the model equifinality.

Snow is Earth’s most climatically sensitive land cover type. Air temperature and moisture availability are first-order controls on snowfall. Maximum snowfall occurs at temperatures very near 0°C, so even a slight increase in temperature will shift a snowy winter to one with midseason rainfall and melt events. Traditional snow metrics are not able to adequately capture the changing nature of snow cover. For example, April 1 snow water equivalent (SWE, the amount of water represented by the snowpack) is used as a streamflow predictor. Still, it cannot express the effects of midwinter melt events, now expected in warming snow climates. The multiple impacts of a changing snowpack require a suite of climate indicators derived from readily measured or modelled data that serve as proxies for relevant snow-related and climate-driven processes. Such indicators need to be simple enough to “tell the story” of snowpack changes over space and time, but not overly simplistic as to be conflated with other variables or, conversely, overly complicated in their interpretation. This paper describes a targeted set of spatially explicit, multi-temporal snow metrics for multiple sectors, stakeholders, and scientists. They include metrics based on satellite data from NASA’s Moderate Resolution Imaging Spectroradiometer, meteorological observations and snow data from ground-based stations, and climate model output. We describe and provide examples for Snow Cover Frequency (SCF), Snow Disappearance Date (SDD), snowstorm temperature (ST), At-Risk Snow (ARS), and Frequency of a Warm Winter (FWW).

The rapid expansion of impermeable surfaces in cities has a major impact on urban hydrology. Infiltration of rainwater is reduced and water runs off faster with higher runoff peaks. Urban trees as stormwater management tools are becoming more relevant to reduce flood risks in addition to other ecosystem services. An in-situ field experiment to measure throughfall on Norway maple ( Acer platanoides) and small-leaved lime ( Tilia cordata) was conducted to determine the interception of solitary urban trees with different degrees of surface sealing in the city of Freiburg, Germany. The relationships between rainfall characteristics, tree morphological traits, and the interception behavior were investigated with eight trees per species. 76 recorded rainfall events were evaluated from April to September 2021. Average interception values were higher for small-leaved lime (70.3 ± 6.6%) than for Norway maple (54.8 ± 10.3%) and hence much higher than in a typical forested environment. The average interception loss of all recorded events was 2.58 ± 0.60 mm for Norway maple and 3.73 ± 0.29 mm for small-leaved lime. For both tree species, significant linear correlations were found between the relative interception and other factors like rainfall depths, the leaf area index (LAI), and the plant area index (PAI) (adj.R 2 > 0.45). In contrast to Norway maple, small-leaved lime also showed significant relationships of several tree morphological parameters with the interception (adj.R 2 > 0.43). LAI, which also effects the interception, of both tree species significantly decreased with the degree of surface sealing. Our results provide a better understanding of the interception process of solitary trees for different urban sites and allows to parameterize interception based on measurable properties. However, further field experiments with various tree species need to be conducted to obtain a larger database for typical parameters in models and to support urban planners in managing stormwater runoff.

Groundwater-surface water interactions (GSI) connect rivers and streams with riparian areas and the adjacent aquifer. Although these interactions represent a substantial control of quantity and quality of both groundwater and surface water, knowledge on GSI along rivers at the regional scale, particularly inland waterways, is still limited. We investigated GSI along the river Moselle, an important federal inland waterway in Germany, by using radon and tritium to identify gaining and losing stream conditions, respectively. Gaining stream conditions were identified by continuously measuring radon along the river as part of boat surveys using a high spatial resolution (every two km) during intermediate (October 2020) and near low flow conditions (August/September 2021). Quarterly tritium inventories from 2017 to 2021 revealed losses of up to 27 % due to losing stream conditions at the upstream locations of damns (particularly near the hydroelectric power plant Lehmen) while gains (up to 51 %) likely triggered by a flood-induced mass transfer of water from the aquifer back into the river. Using radon mass balance modeling, good agreements of simulated versus measured radon data with respect to two groundwater end-member scenarios were obtained during intermediate flow (Spearman’s ρ: 0.97 and 0.99; MAE: 10.1 and 3.4 Bq l -1) and near low flow (Spearman’s ρ: 0.97 and 0.99; MAE: 11 and 6.5 Bq l -1). Important groundwater inflow was limited to the meander of Detzem, where cumulated groundwater inflow of about 19 m 3 s -1 (9.5 % of total discharge) and 4.2 m 3s -1 (3.8 % of total discharge) was simulated during intermediate and near low flow, respectively. However, these groundwater contributions were relatively low compared to alpine streams, for example. Finally, the study will help to better identify and quantify GSI at the regional scale and provide methodological guidance for future studies focusing on inland waterways.

A regional coupled approach to water cycle prediction is demonstrated for the 4-month period from November 2013 to February 2014 through analysis of precipitation, soil moisture, river flow and coastal ocean simulations produced by a km-scale atmosphere-land-ocean coupled system focussed on the United Kingdom (UK), running with horizontal grid spacing of around 1.5 km across all components. The Unified Model atmosphere component, in which convection is explicitly simulated, reproduces the observed UK rainfall accumulation (r2 of 0.62 for daily accumulation), but there is a notable bias in its distribution – too dry over western upland areas and too wet further east. The JULES land surface model soil moisture state is shown to be in broad agreement with a limited number of cosmic-ray neutron probe observations. A comparison of observed and simulated river flow shows the coupled system is useful for predicting broad scale features, such as distinguishing high and low flow regions and times during the period of interest but are shown to be less accurate than optimised hydrological models. The impact of simulated river discharge on NEMO model simulations of coastal ocean state is explored in the coupled system, with comparisons provided relative to experiments using climatological river input and no river input around the UK coasts. Results show that the freshwater flux around the UK contributes of order 0.2 psu to the mean surface salinity, and comparisons to profile observations give evidence of an improved vertical structure when applying simulated flows. This study represents a baseline assessment of the coupled system performance, with priorities for future model developments discussed.

Fully coupled atmospheric-hydrological models allow a more realistic representation of the land surface–boundary layer continuum, representing both high-resolution land-surface/subsurface water lateral redistribution and the related feedback towards the atmosphere. This study evaluates the potential contribution of the fully coupled approach in extended-range mesoscale hydrometeorological ensemble forecasts. Previous studies have shown, for deterministic simulations, that the effect of fully coupling for short-range forecasts is minor compared to other sources of uncertainty, however, it becomes not negligible when increasing the forecast period. Through a proof-of-concept consisting of an ensemble (50 members from the ECMWF Ensemble Prediction System) seven-days-in-advance forecast of a high impact event affecting the Calabrian peninsula (southern Italy, Mediterranean basin) on November 2019, the paper elucidates the extent to which the improved representation of the terrestrial water lateral transport in the Weather Research and Forecasting (WRF) – Hydro modeling system affects the ensemble water balance, focusing on the precipitation and the hydrological response, in terms of both soil moisture dynamics and streamflow in 14 catchments spanning over 42% of the region. The fully coupled approach caused an increase of surface soil moisture and latent heat flux from land in the days preceding the event, partially affecting the lower Planetary Boundary Layer. However, when shoreward moisture transport from surrounding sea rapidly increased becoming the dominant process, only a weak signature of soil moisture contribution could be detected, resulting in only slightly higher precipitation forecast and not clear variation trend of peak flow, even though the latter variable increased up to 10% in some catchments. Overall, this study highlighted a remarkable performance of the medium-range ensemble forecasts, suggesting a profitable use of the fully coupled approach for forecasting purposes in circumstances in which soil moisture dynamics is more relevant and needs to be better addressed.

Every application of soil erosion models brings the need of proper parametrization, i.e., finding physically or conceptually plausible parameter values that allow a model to reproduce measured values. No universal approach for model parametrization, calibration and validation exists, as it depends on the model, spatial and temporal resolution and the nature of the datasets used. We explored some existing options for parametrization, calibration and validation for erosion modelling exemplary with a specific dataset and modelling approach. A modified version of the Morgan-Morgan-Finney (MMF) model was selected, representing a balanced position between physically-based and empirical modelling approaches. The resulting calculator for soil erosion (CASE) model works in a spatially distributed way on the timescale of individual rainfall events. A dataset of 142 high-intensity rainfall experiments in Central Europe (AT, HU, IT, CZ), covering various slopes, soil types and experimental designs was used for calibration and validation with a modified Monte-Carlo approach. Subsequently, model parameter values were compared to parameter values obtained by alternative methods (measurements, pedotransfer functions, literature data). The model reproduced runoff and soil loss of the dataset in the validation setting with R 2 adj of 0.89 and 0.76, respectively. Satisfactory agreement for the water phase was found, with calibrated saturated hydraulic conductivity (k sat) values falling within the interquartile range of k sat predicted with 14 different PTFs, or being within one order of magnitude. The chosen approach also well reflected specific experimental setups contained in the dataset dealing with the effects of consecutive rainfall and different soil water conditions. For the sediment phase of the tested model agreement between calibrated cohesion, literature values and field measurements were only partially in line. For future applications of similar model applications or datasets, the obtained parameter combinations as well as the explored methods for deriving them may provide guidance.

Concentration-discharge (C-Q) relationships can provide insight into how catchments store and transport solutes, but analysis is often limited to long-term behavior assessed from infrequent grab samples. Increasing availability of high-frequency sensor data has shown that C-Q relationships can vary substantially across temporal scales, and in response to different hydrologic drivers. Here, we present four years of dissolved organic carbon (DOC) and nitrate-nitrogen (NO3-N) sensor data from a snowmelt- dominated catchment in the Rocky Mountains of Colorado. We assessed both the direction (enrichment vs. dilution) and hysteresis in C-Q relationships across a range of time scales, from interannual to sub-daily. Both solutes exhibited a seasonal flushing response, with concentrations initially increasing as solute stores are mobilized by the melt pulse, but then declining as these stores are depleted. The high-frequency data revealed that the seasonal melt pulse was composed of numerous individual daily melt pulses. The solute response to daily melt pulses was relatively chemostatic, suggesting mobilization and depletion to be progressive rather than episodic processes. In contrast, rainfall-induced pulses produced short-lived but substantial enrichment responses, suggesting they may activate alternative solute sources or transport pathways. The results clearly demonstrate that solute responses to individual events may differ significantly from the longer-term behavior these events combine to generate, something which only becomes apparent when high-frequency data are used.  

Irrigation activities are a major control on water movement and storage in irrigated river valleys in the Intermountain West, USA. Particularly in dry years, surface water diversions can deplete streams over the summer irrigation season, leading to more variable stream temperatures and increased risk for resident aquatic species. Cooler lateral inflows derived from irrigation activities can mitigate the impacts of depletion by buffering main channel stream temperatures. Given the increasing susceptibility of depleted streams to climate and land use changes, understanding stream temperature patterns and controls in these systems is critical. We used intensive field monitoring over three summers and thermal aerial imagery to characterize stream temperature patterns and irrigation influences in a 2.5 km reach of a small agricultural stream in northern Utah. Considering variable hydrology, weather, channel morphology, diversions, and lateral inflows we found stream temperatures to be relatively insensitive to flow depletion or lateral inflows in a wet year but very sensitive in drier years. Irrigation-related lateral inflows reduced longitudinal warming and diel variability during drier years and at times prevented temperatures from reaching stressful or lethal limits. Reaches with substantial lateral inflow contributions also had a greater areal proportion of low temperatures and spatial temperature diversity. These trends were enhanced by differences in channel morphology, with greater spatial and temporal variability in multi-thread than single-thread reaches. Study results highlight critical flow and weather conditions driving increased temperature variability that will likely become more extreme with additional climate change related reductions in baseflow. Regardless of the cause, this study highlights that decreased instream flows increase the importance of identifying, quantifying, and maintaining lateral inflows to maintain instream temperatures and preservation of these inflows should be considered in future water management decisions.

Informative/Abstract:Estimating streamflow is time and labour intensive due to the necessity of developing a rating curve. The development of a rating curve involves acquiring at least thirty in-field measurements of streamflow across a wide range of flow levels, which can be costly and impractical in remote regions with limited seasonal access. Here we showcase an automated system which accurately estimates streamflow multiple times each day, greatly facilitating the development of rating curves for remote or seasonally inaccessible sites. The system uses an emerging technique referred to as particle image velocimetry (PIV) to track the movement of objects and flow structure features on the mobile water surface to generate velocity vector grids. Velocity grids were used to calculate streamflow and facilitate the development of a rating curve. This represents the first use of an automated PIV system to estimate streamflow in small streams (< 5 m wide) and the first system to automatically distribute particles for facilitated PIV analysis.Keywords: Particle Image Velocimetry, Streamflow Monitoring, Automated Systems, Particle TracerFunding: This research was funded through the Ministry of Natural Resources and Forestry, the Canada-Ontario Agreement Fund, and the Queen Elizabeth II Graduate Scholarship in Science and Technology.

The streambed is the critical interface between the aquatic and terrestrial systems and hosts important biogeochemical hot spots within river corridors. Although the streambed characteristics are significantly different from those of its surrounding soil, the streambed itself has not been explicitly represented in watershed models. We developed an integrated hydrologic model that explicitly incorporated a streambed layer to examine the hydrological effects of streambed characteristics including hydraulic conductivity (K), layer thickness, and width on the exchange fluxes across the streambed as well as the streamflow at the watershed outlet. The numerical experiments were performed in the American River Watershed, a headwater, mountainous watershed within the Yakima River Basin in central Washington. Despite having a negligible effect on the watershed streamflow, an explicit representation of the streambed with distinctive properties dramatically changed the magnitude and variability of the exchange flux. In general, larger streambed K along with a thicker streambed layer induced larger exchange fluxes. The exchange flux was most sensitive to the streambed width or the mesh resolution of the streambed. A smaller streambed width (or a finer streambed resolution) increases exchange fluxes per unit area while reducing the overall exchange volumes across the entire streambed. The amount of baseflow decreased by 6% as the streambed width decreased from 250 m to 50 m. This finding is important because these hydrological changes may in turn affect the exchange of nutrients and contaminants between surface water and groundwater and the associated biogeochemical processes. Our work demonstrated the importance of representing streambed in fully distributed, process-based watershed models in better capturing the exchange flow dynamics in river corridors.

Snowmelt drives a large portion of streamflow in many mountain areas of the world. However, the water pathways since snow melts until water reaches the streams, and its associated transit time is still largely unknown. Such processes are important for drawing conclusions about the hydrological role of the upstream snowpack after melting. This work analyzes for first time the influence of snowmelt on spring streamflow in years of different snow accumulation and duration, in an alpine catchment of the central Spanish Pyrenees. A multi-approach research was performed, by combining the analysis of climatic, snow, streamflow, piezometric levels, water temperature, electrical conductivity and isotopic (δ 18O) data. Results show that snow played a preeminent role on the hydrological response of the catchment during spring. Liquid precipitation during the melting period also determined the shape of the spring hydrographs. When snow cover disappeared from the catchment, soil water storage and streamflow showed a sharp decline. Consequently, streamflow electrical conductivity, temperature and δ 18O showed a marked tipping point towards higher values. The fast hydrological response of the catchment to snow and meteorological fluctuations, as well as the marked diel fluctuations of streamflow δ 18O during the melting period, strongly suggests soil storage was small, leading to short meltwater transit times. As a consequence of this hydrological behavior, independently of the amount of snow accumulated and of melting date, summer streamflow remained always low, with small runoff peaks driven by rainfall events. The expected reduction of snow accumulation and duration in the area in a next future will bring an earlier snowmelt and rise of stream water temperature. However, given the low storage capacity of the catchment and the contribution of rainfall events to spring runoff, the annual water balance and the runoff seasonality of the catchment would not change drastically.

Evapotranspiration (ET) plays an important role in integrated water resource planning, development and management. This process is particularly relevant in semiarid regions. The aim of the present study is to compare the actual spatial and temporal evapotranspiration (ETa) patterns and temporal trends in two semiarid forests, one in Brazil (Aiuaba) and the other in Spain (Valladolid). We used the Surface Energy Balance Algorithm for Land (SEBAL) to assess the effect of climatic variation in both areas. In the Brazilian semiarid forest, Caatinga is the main vegetation, while it is stone pine in Spain. For this purpose, 69 Landsat-5 and 42 Landsat-8 images (1995 – 2019) were used. The Mann-Kendall test was applied to assess the occurrence of trends in precipitation, temperature and potential evapotranspiration data; and the Temporal Stability Index (TSI) to know which areas have greater seasonal ETa. The annual amplitude of the potential evapotranspiration (ET0) is the same in both areas, however, the Caatinga values are higher. In the Caatinga forest, when ET0 presents its highest values throughout the year, ETa presents the lowest, and vice versa. In the Pinares forest, ETa follows the ET0 dynamics during the year, and the difference between ET0 and ETa is maximum during the summer. The Caatinga forest showed a greater spatial variation of ETa than the Pinares forest as well as a greater extension with lower temporal stability of ETa than the Pinares forest. Both the Caatinga forest and the Pinares forest showed significant annual trends of increase for ET0 and ETa: 3.5 mm and 2.2 mm, and 7.0 mm and 3.9 mm, respectively.

Runoff and erosion can increase after wildfires, but little is known about the effects of wildfire plus post-fire salvage logging, or mitigating these effects. Past research has identified soil compaction and reduced surface cover as controls on runoff and erosion, but the relative contributions of these changes are not clear. Two years after high severity burning by the 2015 Valley Fire in California, replicated rainfall simulations were carried out in four soil conditions across compaction and cover factors: uncompacted/compacted by logging machinery and bare soil/60% wood slash-cover. Runoff after 71 mm of rainfall totaled 27 mm in the uncompacted bare plots and 39 mm in the compacted bare plots. Runoff in the slash-covered plots decreased by 50% and 33% as compared to the uncompacted and compacted bare plots, respectively, although none of the differences in runoff were significant. Rainsplash averaged 30 g for the bare plots, regardless of compaction, and decreased significantly by 70% on slash-covered plots. Sediment yield totaled 460 and 818 g m-2 for the uncompacted and compacted bare plots, respectively, and slash significantly reduced these amounts by 72% and 69%, respectively. Our results showed that post-fire soil erosion in high severity burned unlogged areas was still very high two years after the wildfire. The combination of wildfire and salvage logging doubled soil erosion by increases in both runoff amount and sediment concentration. Antecedent soil moisture (dry or wet) was the dominant factor for runoff, while surface cover was the dominant factor for erosion and sediment delivery. Covering the soil with slash reduced both runoff and erosion, suggesting this treatment would reduce long‐term sediment delivery from burned areas and skid trails. Saturated hydraulic conductivity (Ks) and interrill erodibility (Ki) calculated from these simulations confirmed previous research and will support modeling efforts related to wildfire and post-fire salvage logging.

Post-fire managers throughout the world use predictive models to estimate potential erosion risks to aid in evaluating downstream impacts of increased runoff and erosion, and to target critical areas within a fire for applying mitigation practices. Erosion prediction can be complicated by forest road networks. Using GIS technology and a soil erosion model, this study evaluated the effect of roads on erosion and sediment yield following a wildfire, and whether the predictive models were providing reasonable results. The GeoWEPP model was used to simulate onsite erosion and offsite sediment delivery before and after fire disturbance. A 2-m resolution DEM was used as the terrain layer. Erosion rates in excess of 4 Mg ha-1 yr-1 were predicted mainly from the moderate and high severity burn areas. Roads influenced both flow path and sub-catchment delineations, affecting the spatial distribution of sediment detachment and transport. Through that influence, roads tended to reduce estimated erosion on slopes below the roads, but road fillslopes and steep channels were areas of significant increases in erosion risks. Measured deposition amounts along roads and in sediment basins were similar to predicted amounts. The results confirm that road prisms, culverts and road ditches greatly influence sedimentation processes after wildfire, and they present opportunities to detain eroded sediments before they reach downstream water bodies.

It is widely recognized that multi-year drought can induce changes in catchment hydrological behaviors. However, at present, our understanding about multi-year, drought-induced changes in catchment hydrological behaviors and its driving factors at the process level is still very limited. This study proposed a new approach using a data assimilation technique with a process-based hydrological model to detect multi-year drought-induced changes in catchment hydrological behaviors and to identify driving factors for the changes in an unimpaired Australian catchment (Wee Jasper) which experienced prolonged drought from 1997 to 2009. Modelling experiments demonstrated that the multi-year drought caused a significant change in the catchment rainfall-runoff relationship, indicated by significant step changes in the estimated time-variant hydrological parameters SC (indicating catchment active water storage capacity) and C (reflecting catchment evapotranspiration dynamics), whose average values increased 23.4% and 10.2%, respectively, due to drought. The change in the rainfall-runoff relationship identified by the data assimilation method is consistent with that arrived at by a statistical examination. The proposed method provides insights about the drivers of the changes in the rainfall-runoff relationship at the processes level. Declining groundwater and deep soil moisture depleted by persistent evapotranspiration of deep-rooted woody vegetation during drought are the main driving factors for the catchment behaviors change in the Wee Jasper catchment. The new method proposed in this study was found to be an effective technique for detecting both the change of hydrological behaviors induced by prolonged drought and its driving factors at the process level.

Increased urbanization, coupled with the projected impacts of climatic change, mandates further evaluation of the impact of urban development on water flow paths to guide sustainable land use planning. Though the general urbanization impacts of increased storm runoff peaks and reduced baseflows are well known; how the complex, non-stationary interaction of the dominant water fluxes within dynamic urban water stores sustain streamflow regimes over longer periods of time are less well quantified. In particular, there is a challenge in how hydrological modelling should integrate the juxtaposition of rapid and slower flow pathways of the urban ‘karst’ landscape and different approaches need evaluation. In this context, we utilized hydrological and water stable isotope datasets within a modelling framework that combined the commonly used HEC urban runoff model along with a simple hydrological tracer module and transit time modelling to evaluate the spatial and temporal variation of water flow paths and ages within a heavily urbanized 217km 2 catchment in Berlin, Germany. Deeper groundwater was the primary flow component within less urbanized regions of the catchments, with increased direct runoff and shallow subsurface contributions in more urbanized areas near the catchment outlet. The addition of wastewater effluent in the mid-reaches of the catchment was the dominant water supply to the lower stream, and sustained baseflows during the summer months. Water ages from each modelling approach imitated flow contributions and opportunity for mixing with subsurface storage; with older water and lower young water contributions in less urbanized sub-catchments and younger water and higher young water contributions in more urbanized regions. The results form a first step towards more integrated modelling tools for similar peri-urban catchments, given the potential limitations of more simple model frameworks. The results have broader implications for assessing the uncertainty in evaluating urban impacts on hydrological function under environmental change.

Wildfires are a cause of soil water repellency (hydrophobicity), which reduces infiltration while increasing erosion and flooding from post-fire rainfall. Post-fire soil water repellency degrades over time, often in response to repeated wetting and drying of the soil. However, in mountainous fire-prone forests such as those in the Western USA, the fire season often terminates in a cold and wet winter, during which soils not only wet and dry, but also freeze and thaw. Little is know about the effect of repeated freezing and thawing of soil on the breakdown of post-fire hydrophobicity. This study characterized the changes in hydrophobicity of Sierra Nevada mountain soils exposed to different combinations of wet-dry and freeze-thaw cycling. Following each cycle, hydrophobicity was measured using the Molarity of Ethanol test. Hydrophobicity declined similarly across all experiments that included a wetting cycle. Repeated freezing and thawing of dry soil did not degrade soil water repellency. Total soil organic matter content was not different between soils of contrasting hydrophobicity. Macroscopic changes such as fissures and cracks were observed to form as soil hydrophobicity decayed. Microscopic changes revealed by scanning electron microscope imagery suggest different levels of soil aggregation occurred in samples with distinct hydrophobicities, although the size of aggregates was not clearly correlated to the change in water repellency due to wet-dry and freeze-thaw cycling. A nine year climate and soil moisture record from Providence Critical Zone Observatory was combined with the laboratory results to estimate that hydrophobicity would persist an average of 144 days post-fire at this well-characterized, typical mid-elevation Sierra Nevada site. Most of the breakdown in soil water repellency (79%) under these climate conditions would be attributable to freeze-thaw cycling, underscoring the importance of this process in soil recovery from fire in the Sierra Nevada.