Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snowfall rate”

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 55 records · Page 3

Arctic Soil Patterns Analogous to Fluid Instabilities: Supporting Data

This dataset characterizes solifluction lobe morphology and spatial patterns using pre-existing LiDAR-derived digital elevation models of 25 sites across Norway with accompanying long term climate data for each site. Data were collected as part of an effort to better understand controls on the formation of solifluction patterns and to test the idea that they are analogous to fluid instabilities. We also provide soil velocity profiles and estimates of effective viscosity from across the world, drawn from literature. They were collected to improve our understanding of the rheology of soliflucting soil. See this article for more information on the theoretical motivation behind this dataset see "Arctic soil patterns analogous to fluid instabilities" (Glade et al., 2021). Data files arranged in a hierarchy and include image files *.tif and *.png, GIS shapefiles and geopackages (*.gpkg), *.csv (with same file as *.xlsx), and *.py (Python scripts readable with a text editor). Files also bundled into *.zip files. Note (2021-10-20): unit corrections made on two files: RR.csv and snowfall.csv. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Precipitation observed by the airborne advanced microwave moisture sounder (AMMS)

Airborne radiometric measurements at frequencies near 92 GHz and 183 GHz were conducted over two precipitation events near Wallops Island, Virginia during February, 1986. The measured brightness temperatures are compared with those from calculations to estimate the snowfall and rainfall rates for both events. The estimated rates over water surface are within a factor of two of those derived from the concurrent measurements by the SPANDAR radar at the NASA Wallops Flight Facility. These estimated rates, however, suggest that both snowfall events are light and close to the threshold of radiometric detection especially over land surface. Observations of additional snowfall events with lower frequency channels are needed to demonstrate the approach and to account for cloud effects.

Chang, A. T. C.↗

Role of Snowfall Versus Air Temperatures for Greenland Ice Sheet Melt-Albedo Feedbacks

The Greenland Ice Sheet is a leading contributor to global sea-level rise because climate warming has enhanced surface meltwater runoff. Melt rates are particularly sensitive to air temperatures due to feedbacks with albedo. The primary melt-albedo feedback, fluctuation of seasonal snowlines, however, is determined not only by melt but also by antecedent snowfall which could delay the onset of dark glacier ice exposure. Here we investigate the role of snowfall versus air temperatures on ice sheet melt-albedo feedbacks using satellite remote sensing and atmospheric reanalysis data. We find several lines of evidence that snowline fluctuations are driven primarily by air temperatures and that snowfall is a secondary control. First, standardized linear regressions indicate that the timing of glacier ice exposure is nearly twice as sensitive to air temperatures than antecedent snowfall. Second, in 74% of the ablation zone by area, winter snowfall rates are not significantly correlated with winter air temperatures. This relationship implies that ice sheet melt due to climate warming is unlikely to be compensated by higher snowfall rates in the ablation zone. Third, we find no significant change in snowfall rates in the ablation zone during our 1981–2021 study period. Our findings demonstrate that snowfall is unlikely to reduce future ice sheet melt and that ice sheet meltwater runoff should be accurately predicted by air temperatures. Although given the importance of melt-albedo feedbacks, ice sheet models that parameterize albedo or are coupled with regional climate models are likely to provide the most accurate projections of mass loss.

J. C. Ryan↗

Cold-Season Precipitation Sensitivity to Microphysical Parameterizations: Hydrologic Evaluations Leveraging Snow Lidar Datasets

Abstract Cloud microphysical processes are an important facet of atmospheric modeling, as they can control the initiation and rates of snowfall. Thus, parameterizations of these processes have important implications for modeling seasonal snow accumulation. We conduct experiments with the Weather Research and Forecasting (WRF V4.3.3) Model using three different microphysics parameterizations, including a sophisticated new scheme (ISHMAEL). Simulations are conducted for two cold seasons (2018 and 2019) centered on the Colorado Rockies’ ∼750-km 2 East River watershed. Precipitation efficiencies are quantified using a drying-ratio mass budget approach and point evaluations are performed against three NRCS SNOTEL stations. Precipitation and meteorological outputs from each are used to force a land surface model (Noah-MP) so that peak snow accumulation can be compared against airborne snow lidar products. We find that microphysical parameterization choice alone has a modest impact on total precipitation on the order of ±3% watershed-wide, and as high as 15% for certain regions, similar to other studies comparing the same parameterizations. Precipitation biases evaluated against SNOTEL are 15% ± 13%. WRF Noah-MP configurations produced snow water equivalents with good correlations with airborne lidar products at a 1-km spatial resolution: Pearson’s r values of 0.9, RMSEs between 8 and 17 cm, and percent biases of 3%–15%. Noah-MP with precipitation from the PRISM geostatistical precipitation product leads to a peak SWE underestimation of 32% in both years examined, and a weaker spatial correlation than the WRF configurations. We fall short of identifying a clearly superior microphysical parameterization but conclude that snow lidar is a valuable nontraditional indicator of model performance.

54 ENVIRONMENTAL SCIENCES↗

Microphysical Properties of Snow and Their Link to Ze-S Relations During BAECC 2014

This study uses snow events from the Biogenic Aerosols-Effects on Clouds and Climate (BAECC) 2014 campaign to investigate the connection between properties of snow and radar observations. The general hydrodynamic theory is applied to video-disdrometer measurements to retrieve masses of falling ice particles. Errors associated with the observation geometry and the measured particle size distribution (PSD) are addressed by devising a simple correction procedure. The value of the correction factor is determined by comparison of the retrieved precipitation accumulation with weighing-gauge measurements. Derived mass-dimensional relations are represented in the power-law form m = amDbm. It is shown that the retrieved prefactor am and exponent bm react to changes in prevailing microphysical processes. From the derived microphysical properties, event-specific relations between the equivalent reflectivity factor Ze and snowfall precipitation rate S (Ze = 5 azsSbzs ) are determined. For the studied events, the prefactor of the Ze-S relation varied between 53 and 782 and the exponent was in the range of 1.19-1.61. The dependence of the factors azs and bzs on the m(D) relation and PSD are investigated. The exponent of the Ze-S relation mainly depends on the exponent of the m(D) relation, whereas the prefactor azs depends on both the intercept parameterN0 of the PSDand the prefactors of them(D) andy(D) relations. Changes in azs for a given N0 are shown to be linked to changes in liquid water path, which can be considered to be a proxy for degree of riming.

von Lerber, Annakaisa↗

Estimation of Seasonal Snow Water Equivalent Using Landsat Observations

This work presents a methodology for estimating seasonal snow water equivalent (SWE) from the use of remotely sensed Visible and Near Infrared observations from the Landsat mission. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data over the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlations with in situ SWE observations of greater than 0.95. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (32 years), and accuracy provide a unique dataset for investigating snow processes. In particular, this presentation illustrates how the reanalysis dataset was used to provide climatology of the seasonal snowfall accumulation rates, distributions, and variability over the last three decades.

SWE↗

A Comprehensive Northern Hemisphere Particle Microphysics Data Set From the Precipitation Imaging Package

Microphysical observations of precipitating particles are critical data sources for numerical weather prediction models and remote sensing retrieval algorithms. However, obtaining coherent data sets of particle microphysics is challenging as they are often unindexed, distributed across disparate institutions, and have not undergone a uniform quality control process. This work introduces a unified, comprehensive Northern Hemisphere particle microphysical data set from the National Aeronautics and Space Administration precipitation imaging package (PIP), accessible in a standardized data format and stored in a centralized, public repository. Data is collected from 10 measurement sites spanning 34° latitude (37°N–71°N) over 10 years (2014–2023), which comprise a set of 1,070,000 precipitating minutes. The provided data set includes measurements of a suite of microphysical attributes for both rain and snow, including distributions of particle size, vertical velocity, and effective density, along with higher-order products including an approximation of volume-weighted equivalent particle densities, liquid equivalent snowfall, and rainfall rate estimates. The data underwent a rigorous standardization and quality assurance process to filter out erroneous observations to produce a self-describing, scalable, and achievable data set. Case study analyses demonstrate the capabilities of the data set in identifying physical processes like precipitation phase-changes at high temporal resolution. Bulk precipitation characteristics from a multi-site intercomparison also highlight distinct microphysical properties unique to each location. This curated PIP data set is a robust database of high-quality particle microphysical observations for constraining future precipitation retrieval algorithms, and offers new insights toward better understanding regional and seasonal differences in bulk precipitation characteristics.

54 ENVIRONMENTAL SCIENCES↗

Passive Microwave Signatures and Retrieval of High-Latitude Snowfall Over Open Oceans and Sea Ice: Insights From Coincidences of GPM and CloudSat Satellites

This article studies changes in microwave signals of oceanic snowfall in response to the formation of snow-covered sea ice using active and passive coincident data from the radar and radiometer onboard the CloudSat and the global precipitation measurement satellites. Using reanalysis data of liquid and ice water path as well as satellite retrievals of sea ice snow-cover depth, spectral regions are determined over which the snowfall signatures are likely to be obscured or falsely detected. Relying on an a priori database populated with the active–passive coincidences, a Bayesian snowfall retrieval algorithm is presented that links a k -nearest neighbor matching with the inverse Gaussian estimator used in the Goddard profiling algorithm. Without relying on any ancillary data of air temperature, the results demonstrate that over open oceans (sea ice), we can passively retrieve the CloudSat active snowfalls with a true positive rate of 92 (85%) and the root mean squared error of 0.24 (0.15) mm per h.

Sajad Vahedizade↗

Evaluating the Performance of Single and Double Moment Microphysics Schemes During a Synoptic-Scale Snowfall Event

Increases in computing resources have allowed for the utilization of high-resolution weather forecast models capable of resolving cloud microphysical and precipitation processes among varying numbers of hydrometeor categories. Several microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, ranging from single-moment predictions of precipitation content to double-moment predictions that include a prediction of particle number concentrations. Each scheme incorporates several assumptions related to the size distribution, shape, and fall speed relationships of ice crystals in order to simulate cold-cloud processes and resulting precipitation. Field campaign data offer a means of evaluating the assumptions present within each scheme. The Canadian CloudSat/CALIPSO Validation Project (C3VP) represented collaboration among the CloudSat, CALIPSO, and NASA Global Precipitation Measurement mission communities, to observe cold season precipitation processes relevant to forecast model evaluation and the eventual development of satellite retrievals of cloud properties and precipitation rates. During the C3VP campaign, widespread snowfall occurred on 22 January 2007, sampled by aircraft and surface instrumentation that provided particle size distributions, ice water content, and fall speed estimations along with traditional surface measurements of temperature and precipitation. In this study, four single-moment and two double-moment microphysics schemes were utilized to generate hypothetical WRF forecasts of the event, with C3VP data used in evaluation of their varying assumptions. Schemes that incorporate flexibility in size distribution parameters and density assumptions are shown to be preferable to fixed constants, and that a double-moment representation of the snow category may be beneficial when representing the effects of aggregation. These results may guide forecast centers in optimal configurations of their forecast models for winter weather and identify best practices present within these various schemes.

Molthan, Andrew L.↗

Modelling Greenland Outlet Glaciers

The objective of this project was to develop simple yet realistic models of Greenland outlet glaciers to better understand ongoing changes and to identify possible causes for these changes. Several approaches can be taken to evaluate the interaction between climate forcing and ice dynamics, and the consequent ice-sheet response, which may involve changes in flow style. To evaluate the icesheet response to mass-balance forcing, Van der Veen (Journal of Geophysical Research, in press) makes the assumption that this response can be considered a perturbation on the reference state and may be evaluated separately from how this reference state evolves over time. Mass-balance forcing has an immediate effect on the ice sheet. Initially, the rate of thickness change as compared to the reference state equals the perturbation in snowfall or ablation. If the forcing persists, the ice sheet responds dynamically, adjusting the rate at which ice is evacuated from the interior to the margins, to achieve a new equilibrium. For large ice sheets, this dynamic adjustment may last for thousands of years, with the magnitude of change decreasing steadily over time as a new equilibrium is approached. This response can be described using kinematic wave theory. This theory, modified to pertain to Greenland drainage basins, was used to evaluate possible ice-sheet responses to perturbations in surface mass balance. The reference state is defined based on measurements along the central flowline of Petermann Glacier in north-west Greenland, and perturbations on this state considered. The advantage of this approach is that the particulars of the dynamical flow regime need not be explicitly known but are incorporated through the parameterization of the reference ice flux or longitudinal velocity profile. The results of the kinematic wave model indicate that significant rates of thickness change can occur immediately after the prescribed change in surface mass balance but adjustments in flow rapidly diminish these rates to a few cm/yr at most. The time scale for adjustment is of the order of a thousand years or so.

vanderVeen, Cornelis↗

Polarization Lidar Liquid Cloud Detection Algorithm for Winter Mountain Storms

We have collected an extensive polarization lidar dataset from elevated sites in the Tushar Mountains of Utah in support of winter storm cloud seeding research and experiments. Our truck-mounted ruby lidar collected zenith, dual-polarization lidar data through a roof window equipped with a wiper system to prevent snowfall accumulation. Lidar returns were collected at a rate of one shot every 1 to 5 min during declared storm periods over the 1985 and 1987 mid-Jan. to mid-Mar. Field seasons. The mid-barrier remote sensor field site was located at 2.57 km MSL. Of chief interest to weather modification efforts are the heights of supercooled liquid water (SLW) clouds, which must be known to assess their 'seedability' (i.e., temperature and height suitability for artificially increasing snowfall). We are currently re-examining out entire dataset to determine the climatological properties of SLW clouds in winter storms using an autonomous computer algorithm.

Sassen, Kenneth↗

Detection Thresholds of Falling Snow From Satellite-Borne Active and Passive Sensors

There is an increased interest in detecting and estimating the amount of falling snow reaching the Earths surface in order to fully capture the global atmospheric water cycle. An initial step toward global spaceborne falling snow algorithms for current and future missions includes determining the thresholds of detection for various active and passive sensor channel configurations and falling snow events over land surfaces and lakes. In this paper, cloud resolving model simulations of lake effect and synoptic snow events were used to determine the minimum amount of snow (threshold) that could be detected by the following instruments: the W-band radar of CloudSat, Global Precipitation Measurement (GPM) Dual-Frequency Precipitation Radar (DPR)Ku- and Ka-bands, and the GPM Microwave Imager. Eleven different nonspherical snowflake shapes were used in the analysis. Notable results include the following: 1) The W-band radar has detection thresholds more than an order of magnitude lower than the future GPM radars; 2) the cloud structure macrophysics influences the thresholds of detection for passive channels (e.g., snow events with larger ice water paths and thicker clouds are easier to detect); 3) the snowflake microphysics (mainly shape and density)plays a large role in the detection threshold for active and passive instruments; 4) with reasonable assumptions, the passive 166-GHz channel has detection threshold values comparable to those of the GPM DPR Ku- and Ka-band radars with approximately 0.05 g *m(exp -3) detected at the surface, or an approximately 0.5-1.0-mm * h(exp -1) melted snow rate. This paper provides information on the light snowfall events missed by the sensors and not captured in global estimates.

remote sensing↗

Distributed Disdrometer and Rain Gauge Measurement Infrastructure Developed for GPM Ground Validation

Global Precipitation Mission (GPM)retrieval algorithm validation requires datasets characterizing the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (FOV;<10 km). Collection of this data provides a means to assess retrieval errors related to beam filling and algorithm PSD assumptions. Hence, GPM Ground Validation is developing a deployable network of precipitation gauges and disdrometers to provide fine-scale measurements of PSD and precipitation accumulation variability. These observations will be combined with dual-frequency, polarimetric, and profiling radar data in a bootstrapping fashion to extend validated PSD measurements to a large coverage domain. Accordingly, a total of 24 Parsivel disdrometers(PD), 5 3rd-generation 2D Video Disdrometers (2DVD), 70 tipping bucket rain gauges (TBRG),9 weighing gauges, 7 Hot-Plate precipitation sensors (HP), and 3 Micro Rain Radars (MRR) have been procured. In liquid precipitation the suite of TBRG, PD and 2DVD instruments will quantify a broad spectrum of rain rate and PSD variability at sub-kilometer scales. In the envisioned network configuration 5 2DVDs will act as reference points for 16 collocated PD and TBRG measurements. We find that PD measurements provide similar measures of the rain PSD as observed with collocated 2DVDs (e.g., D0, Nw) for rain rates less than 15 mm/hr. For heavier rain rates we will rely on 2DVDs for PSD information. For snowfall we will combine point-redundant observations of SWER distributed over three or more locations within a FOV. Each location will contain at least one fenced weighing gauge, one HP, two PDs, and a 2DVD. MRRs will also be located at each site to extend the measurement to the column. By collecting SWER measurements using different instrument types that employ different measurement techniques our objective is to separate measurement uncertainty from natural variability in SWER and PSD. As demonstrated using C3VP polarimetric radar, gauge, and 2DVD/PD datasets these measurements can be combined to bootstrap an area wide SWER estimate via constrained modification of density-diameter and radar reflectivity-snowfall relationships. These data will be combined with snowpack, airborne microphysics, radar, radiometer, and tropospheric sounding data to refine GPM snowfall retrievals. The gauge and disdrometer instruments are being developed to operate autonomously when necessary using solar power and wireless communications. These systems will be deployed in numerous field campaigns through 2016. Planned deployment of these systems include field campaigns in Finland (2010), Oklahoma (2011), Canada (2012) and North Carolina (2013). GPM will also deploy 20 pairs of TBRGs within a 25 km2 region along the Virginia coast under NASA NPOL radar coverage in order to quantify errors in point-area rainfall measurements.

Petersen, Walter A.↗

GPCP Version 3.2 Products and Results

The Global Precipitation Climatology Project (GPCP) products address the need for long-term precipitation products that emphasize homogeneity, following Climate Data Record (CDR) principles. The new-generation Version 3.2 provides key improvements over the operational Version 2.3 such as: finer spatial resolution of 0.5°x0.5°; wider geosynchronous infrared estimation (58°N-S) upgraded with the PERSIANN-CDR algorithm; upgraded retrievals from selected passive microwave sensors (GPROF algorithm) that calibrate the IR input; revised intercalibrations of TOVS and AIRS data (used at high latitudes); climatologies based on CloudSat, TRMM, and GPM to provide overall calibration by modern satellite estimates; the latest Global Precipitation Climatology Centre (GPCC) precipitation gauge analyses over land areas; regional modifications to the gauge undercatch correction; and IMERG half-hourly data input to the Daily V3.2 product. We will show sample analyses that demonstrate aspects of the Version 3.2 precipitation record, such as the global climatology, the time series for global land and ocean total precipitation and snowfall, and the time series of tropical land and ocean daily precipitation rate histograms. For selected analyses we will show improvements in both the Monthly and Daily products in Version 3.2 compared to the operational Version 2.3. In particular, the climatological zonal profile of precipitation in the Southern Ocean, extending south of 40°S, improves a suspected artifact in V2.3. Similarly, the Daily histograms over ocean in Version 3.2 lack the jump in the predecessor Version 1.3 Daily over ocean at the start of 2009, although a smaller jump is introduced in June 2014. The presentation will conclude with a prospectus for the future satellites/sensors and community datasets necessary to continue computation of a consistent CDR product on the one hand, while also potentially contributing to improvements in the historical record.

Global Precipitation Measurement↗

Mapping of ice layer extent and snow accumulation in the percolation zone of the Greenland ice sheet

The Greenland ice sheet underwent record extensive melt in 2002 and prolonged melt in 2003. The severe melting created a significant and extensive ice layer over the Greenland ice sheet. An innovative approach is developed to detect the ice layer formation using data acquired by the SeaWinds scatterometer on the QuikSCAT satellite. QuikSCAT backscatter together with in situ data from automatic weather stations of the Greenland Climate Network are used to map the extent of ice layer formation. The results reveal areas of extensive ice layer formed by the 2002 melt, which is consistent with the maximum melt extent in 2002. Moreover, during freezing seasons, QuikSCAT data show a linear decrease in backscatter (in decibels or dB) that is related to the amount of snow accumulation in the ice layer formation region. This snow accumulation signature is caused by the attenuation of radar waves in the snow layer, accumulating since the last major melt event, whose thickness appears as an exponential function in relation to the backscatter signature. We use the Greenland Climate Network data to calibrate the QuikSCAT accumulation rate in order to estimate and map snow accumulation. QuikSCAT results capture the extreme snowfall in mid-April 2003, which deposited more than 0.5 m of snow in a day as measured by the automated weather station at the NASA South East site. Large-scale QuikSCAT results show an anomalous increase of snow accumulation over the southeast region of Greenland during the 2002-2003 freezing season.

snow↗

Cold Climate Degradation: An Analysis of Double-Axis Tracked, E-W Vertical, and Fixed-Tilt Photovoltaic Deployments in Alaska

As countries around the world transition towards renewable energy, there is increasing interest in using photovoltaic (PV) technologies to help decarbonize remote northern communities due to their scalability and affordability. However, a major barrier towards large-scale adoption of PV in cold climates is performance uncertainty under extreme environmental conditions including snowfall, freeze-thaw cycles, and high wind loads. Existing literature on PV degradation rates in the North is relatively limited, with published degradation rates varying between -0.2%/year (Sweden) to -1.3%/year (Scotland). At this workshop, we will present preliminary results on the long-term performance of two diverse photovoltaic sites located in Fairbanks, Alaska at 64.8 degrees N: a monofacial Al-BSF double-axis tracking site maintained by the Cold Climate Housing Research Center (CCHRC), and a bifacial PERC/SHJ E-W vertical and south-facing fixed-tilt site maintained by the Alaska Center Energy and Power (ACEP). CCHRC data has been collected over a period of 15 years, while ACEP site data has been collected over 4 years. Using the degradation analysis tool, RdTools, we will present annual system degradation rates, seasonal performance ratio, and identify potential cold climate failure mechanisms for commercially available PV technologies. This analysis will add to existing literature by directly comparing the performance of multiple PV configurations in Alaska.

bifacial↗

Tracking snowmelt during hydrological surface processes using a distributed hydrological model in a mesoscale basin on the Tibetan Plateau

We report that mountain snowpack is an important water resource for the high altitude and latitude regions where the terrain is complex. However, the snowmelt pathway and its actual contribution to streamflow and soil moisture are rarely reported and remain unclear in such regions. To fill in this knowledge gap, we incorporate a snowmelt pathway tracking algorithm to a high-resolution physics-based distributed-hydrology-soil-vegetation model (DHSVM), to track snowmelt movement and quantify snowmelt contributions in the surface hydrologic processes. A simple reservoir operation scheme is also incorporated in the model. The modified model is applied to a dammed meso-scale watershed in the northeast Tibetan Plateau, China to study the snow and reservoir effects. The results show that annual snow contribution to soil moisture (SC-SM) and snow contribution to streamflow (SC-S) significantly decrease over 1965-2019. At a monthly scale, SC-SM has the largest amplitude at the top soil layer and its peak in the deeper layer lags behind the upper layer, and mean monthly SC-S at all stations show bimodal distributions corresponding to snowfall season. Reservoir regulation has subtle impacts (≤2.0%) on SC-S. If the current climate change rate continues, monthly and annual streamflow at the outlet will decrease primarily due to snowpack reduction. To mitigate climate change impacts, better water resource management is needed in this watershed.

54 ENVIRONMENTAL SCIENCES↗