Search for: All records

Methods that combine in-situ measurements, statistical methods, and model simulations with remotely sensed data provide a pathway for improving the robustness of surface flux products. For this research, we acquired eddy-covariance fluxes along a five-tower transect in a snowy montane forest over three consecutive winters to characterize near-field variability of the subcanopy environment. The novel experiment design enabled discriminating near-field evaposublimation sources. Boosted regression trees reveal that the predictive capacity of state variables change with season and storm cycle frequency. High rates of post-storm evaposublimation of canopy-intercepted snow at this site were constrained by short residence time of snow in the canopy due to throughfall and melt. The snow melt-out date for open vs. closed canopy conditions depended on total snowfall accumulation. Compared with low accumulation years, the snow melt-out date under the dense canopy during the high accumulation winter was later than for the open area, as shading became more important later in the season. The field experiments informed an environmental response function that was used to integrate ERA5-Land latent heat flux data at 20-km nominal resolution with USFS Tree Canopy Fraction data at 30-m resolution and showed near-field flux variability that was not resolved in model simulations. Previous evaposublimation results from experiments in alpine and subalpine environments do not directly translate to a montane forest due to differences in process rates. 

To examine spatial and temporal scales of katabatic flow, a distributed temperature sensing (DTS) optical fiber was deployed 2 km down a mild slope irregularly interrupted by small-scale drainage features as part of the Mountain Terrain Atmospheric Modeling and Observation (MATERHORN) experiment conducted at the U.S. Army Dugway Proving Ground, Utah. The fiber was suspended at two heights near the surface, enabling measurement of variations in lapse rate near the surface at meter-scale spatial resolution with 1-min temporal resolution. Experimental results derived from the DTS and tower-mounted instrumentation indicate that airflow through small-scale drainage features regulated the local cooling rate whereas topographic slope and distance along the drainage strongly influenced the larger-scale cooling rate. Empirical results indicate that local cooling rate decays exponentially after local sunset and basin-wide cooling rate decreases linearly with time. The difference in the functional form for cooling rate between local and basin-wide scales suggests that small-scale features have faster timescales that manifests most strongly shortly after local sunset. More generally, partitioning drainage flow by scale provides insight and a methodology for improved understanding of drainage flow in complex terrain. 

Abstract Integrated quadrant analysis is a novel technique to identify and to characterize the trajectory and strength of turbulent coherent structures in the atmospheric surface layer. By integrating the three-dimensional velocity field characterized by traditional quadrant analysis with respect to time, the trajectory history of individual coherent structures can be preserved with Eulerian turbulence measurements. We develop a method to identify the ejection phase of coherent structures based on turbulence kinetic energy (TKE). Identifying coherent structures within a time series using TKE performs better than identifying them with the streamwise and vertical velocity components because some coherent structures are dominated by the cross-stream velocity component as they pass the sensor. By combining this identification method with the integrated quadrant analysis, one can animate or plot the trajectory of individual coherent structures from high-frequency velocity measurements. This procedure links a coherent ejection with the subsequent sweep and quiescent period in time to visualize and quantify the strength and the duration of a coherent structure. We develop and verify the method of integrated quadrant analysis with data from two field studies: the Eclipse Boundary Layer Experiment (EBLE) in Corvallis, Oregon in August 2017 (grass field) and the Vertical Cherry Array Experiment (VACE) in Linden, California in November 2019 (cherry orchard). The combined TKE identification method and integrated quadrant analysis are promising additions to conditional sampling techniques and coherent structure characterization because the identify coherent structures and couple the sweep and ejection components in space. In an orchard (VACE), integrated quadrant analysis verifies each coherent structure is dominated by a sweep. Conversely, above the roughness sublayer (EBLE), each coherent structure is dominated by an ejection. 

Abstract Predicting winter flooding is critical to protecting people and securing water resources in California’s Sierra Nevada. Rain-on-snow (ROS) events are a common cause of widespread flooding and are expected to increase in both frequency and magnitude with anthropogenic climate change in this region. ROS flood severity depends on terrestrial water input (TWI), the sum of rain and snowmelt that reaches the land surface. However, an incomplete understanding of the processes that control the flow and refreezing of liquid water in the snowpack limits flood prediction by operational and research models. We examine how antecedent snowpack conditions alter TWI during 71 ROS events between water years 1981 and 2019. Observations across a 500-m elevation gradient from the Independence Creek catchment were input into SNOWPACK, a one-dimensional, physically based snow model, initiated with the Richards equation and calibrated with collocated snow pillow observations. We compare observed “historical” and “scenario” ROS events, where we hold meteorologic conditions constant but vary snowpack conditions. Snowpack variables include cold content, snow density, liquid water content, and snow water equivalent. Results indicate that historical events with TWI > rain are associated with the largest observed streamflows. A multiple linear regression analysis of scenario events suggests that TWI is sensitive to interactions between snow density and cold content, with denser (>0.30 g cm⁻³) and colder (<−0.3 MJ of cold content) snowpacks retaining >50 mm of TWI. These results highlight the importance of hydraulic limitations in dense snowpacks and energy limitations in warm snowpacks for retaining liquid water that would otherwise be available as TWI for flooding. Significance StatementThe purpose of this study is to understand how the snowpack modulates quantities of water that reach the land surface during rain-on-snow (ROS) events. While the amount of near-term storm rainfall is reasonably predicted by meteorologists, major floods associated with ROS are more difficult to predict and are expected to increase in frequency. Our key findings are that liquid water inputs to the land surface vary with snowpack characteristics, and although many hydrologic models incorporate snowpack cold content and density to some degree, the complexity of ROS events justifies the need for additional observations to improve operational forecasting model results. Our findings suggest additional comparisons between existing forecasting models and those that physically represent the snowpack, as well as field-based observations of cold content and density and liquid water content, would be useful follow-up investigations.