Intermittent snow depth observations can be leveraged with data assimilation (DA) to improve model estimates of snow water equivalent (SWE) at the point scale. A key consideration for scaling a DA system to the basin scale is its performance at locations with forest cover – where canopy-snow interactions affect snow accumulation and melt, yet are difficult to model and parameterize. We implement a particle filter (PF) assimilation technique to assimilate intermittent depth observations into the Flexible Snow Model (FSM2), and validate the output against snow density and SWE measurements across paired forest and open sites, at two locations with different climates and forest structures. Assimilation reduces depth error by 70-90%, density error by 5-30%, and SWE error by 50-70% at forest locations (relative to control model runs) and brings errors in-line with adjacent open sites. The PF correctly simulates the seasonal evolution of the snowpack under forest canopy, including cases where interception lowers SWE in the forest during accumulation, and shading reduces melt during the ablation season (relative to open sites). The snow model outputs are sensitive to canopy-related parameters, but DA reduces the range in depth and SWE estimates resulting from spatial variations or uncertainties in these parameters by more than 50%. The results demonstrate that the challenge of accurately measuring, estimating, or calibrating canopy-related parameters is reduced when snow depth observations are assimilated to improve SWE estimates.</p>
more »
« less
A Meteorology and Snow Data Set From Adjacent Forested and Meadow Sites at Crested Butte, CO, USA
Abstract We present meteorology and snow observation data collected at sites in the southwestern Colorado Rocky Mountains (USA) over three consecutive water years with different amounts of snow water equivalent (SWE) accumulation: A year with above average SWE (2019), a year with average SWE (2020), and a year with below average SWE (2021). This data set is distinguished by its emphasis on paired open‐forest sites in a continental snow climate. Approximately once a month during February–May, we collected data from 15 to 20 snow pits and took 8 to 19 snow depth transects. Our sampling sites were in open and adjacent forested areas at 3,100 m and in a lower elevation aspen (3,035 m) and higher elevation conifer stand (3,395 m). In total, we recorded 270 individual snow pit density and temperature profiles and over 4,000 snow depth measurements. These data are complimented by continuous meteorological measurements from two weather stations: One in the open and one in the adjacent forest. Meteorology data—including incoming shortwave and longwave radiation, outgoing shortwave radiation, relative humidity, wind speed, snow depth, and air and infrared surface temperature—were quality controlled and the forcing data were gap‐filled. These data are available to download from Bonner, Smyth, et al. (2022) athttps://doi.org/10.5281/zenodo.6618553, at three levels of processing, including a level with downscaled, adjusted precipitation based on data assimilation using observed snow depth and a process‐based snow model. We demonstrate the utility of these data with a modeling experiment that explores open‐forest differences and identifies opportunities for improvements in model representation.
more »
« less
- Award ID(s):
- 1761441
- PAR ID:
- 10375840
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 9
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This dataset contains meteorology and snow observation data collected at sites in the southwestern Colorado Rocky Mountains during water years 2019-2021. Data collection had an emphasis on paired open-forest sites and included three forested elevations. In total, we present 270 snow pit observations, 4,019 snow depth measurements, and three years of meteorological forcing from two weather stations (one in a meadow, the other in an adjacent forest). The dataset is described in a forthcoming publication of the same name: A meteorology and snow dataset from adjacent forested and meadow sites at Crested Butte, CO, USA</em> (Bonner et al., 2022).</p> All snow observation and meteorological forcing data are available as both .nc and .mat files. Additionally, original digitized copies of snow pit observations are provided as .gsheet/.xlxs files.</p> This dataset will continue to be updated, via this repository, as additional years of data are collected.</p>more » « less
-
The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) produced a wealth of observational data along the drift of the R/V Polarstern in the Arctic Ocean from October 2019 to September 2020. These data can further process-level understanding and improvements in models. However, the observational records contain temporal gaps and are provided in different formats. One goal of the MOSAiC Single Column Model Working Group (MSCMWG: https://mosaic-expedition.org/science/cross-cutting_groups/) is to provide consistently-formatted, gap-filled, merged datasets representing the conditions at the MOSAiC Central Observatory (the intensively studied region within a few km of R/V Polarstern) that are suitable for driving models on this spatial domain (e.g., single column models, large eddy simulations, etc). The MSCMWG is an open group, please contact the dataset creators if you would like to contribute to future versions of these merged datasets (including new variables). This dataset contains version 1 of these merged datasets, and comprises the variables necessary to force a single column ice model (e.g., Icepack: https://zenodo.org/doi/10.5281/zenodo.1213462). The atmospheric variables are primarily derived from Met City (~66 percent (%) of record, https://doi.org/10.18739/A2PV6B83F), with temporal gaps filled by bias and advection corrected data from Atmospheric Surface Flux Stations ( https://doi.org/10.18739/A2XD0R00S, https://doi.org/10.18739/A25X25F0P, https://doi.org/10.18739/A2FF3M18K). Some residual gaps in shortwave radiation were filled with ARM ship-board radiometer data. Three different options for snowfall precipitation rate (prsn) are provided, based on in-situ observations that precipitation greatly exceeded accumulation on level ice, and accumulation rates varied on different ice types. MOSAiC_kazr_snow_MDF_20191005_20201001.nc uses 'snowfall_rate1' derived from the vertically-pointing, ka-band radar on the vessel (https://doi.org/10.5439/1853942). MOSAiC_Raphael_snow_fyi_MDF_20191005_20201001.nc and MOSAiC_Raphael_snow_syi_MDF_20191005_20201001.nc use snow accumulation measurements from manual mass balance sites (https://doi.org/10.18739/A2NK36626) to derived a pseudo-precipitation. MOSAiC_Raphael_snow_fyi_MDF_20191005_20201001.nc is based on the First Year Ice (fyi) sites. MOSAiC_Raphael_snow_syi_MDF_20191005_20201001.nc is based on the Second Year Ice (syi) sites. The other atmospheric variables for these files are identical. Oceanic variables are in MOSAiC_ocn_MDF_20191006_20200919.nc and are derived from https://doi.org/10.18739/A21J9790B. The data are netCDF files formatted according to the Merged Data File format (https://doi.org/10.5194/egusphere-2023-2413, https://gitlab.com/mdf-makers/mdf-toolkit). The code 'recipes' that were used to produce these data are available at: https://doi.org/10.5281/zenodo.10819497. If you use these datasets, please also cite the appropriate publications: Meteorological variables (excluding precipitation): Cox et al., 2023 (https://doi.org/10.1038/s41597-023-02415-5) Oceanographic variables: Schulz et al., 2023 (https://doi.org/10.31223/X5TT2W) KAZR-derived precipitation: Matrosov et al., 2022 (https://doi.org/10.1525/elementa.2021.00101) Accumulation-derived pseudo-precipitation: Raphael et al., in review. The following are known issues that will be addressed in future dataset releases: 1. Residual gaps occupy approximately 20% of the data record (see addendum) 2. Some transitions to shiprad downwelling shortwave are unreasonable abrupt 3. MDF format does not currently include a field for point-by-point data source Addendum: For atmospheric variables, below indicates the percentage sourced from each dataset (and the amount missing a.k.a NaN) Air Temperature metcity 0.661943 NaN 0.193333 asfs30 0.134910 asfs40 0.008607 asfs50 0.001207 Specific Humidity metcity 0.658890 NaN 0.196298 asfs40 0.008695 Wind Velocity metcity 0.666334 NaN 0.255003 asfs30 0.068828 asfs40 0.008630 asfs50 0.001205 Downwelling Longwave metcity 0.549417 asfs30 0.241502 NaN 0.209081 Downwelling Shortwave metcity 0.674166 NaN 0.158814 asfs30 0.140794 shipradS1 0.026226 Note that the 21 day gap from the end of Central Observatory 2 to the start of Central Observatory 3 occupies 5.8% of the record.more » « less
-
Abstract Understanding how the presence of a forest canopy influences the underlying snowpack is critical to making accurate model predictions of bulk snow density and snow water equivalent (SWE). To investigate the relative importance of forest processes on snow density and SWE, we applied the SUMMA model at three sites representing diverse snow climates in Colorado (USA), Oregon (USA), and Alberta (Canada) for 5 years. First, control simulations were run for open and forest sites. Comparisons to observations showed the uncalibrated model with NLDAS‐2 forcing performed reasonably. Then, experiments were completed to isolate how forest processes affected modelled snowpack density and SWE, including: (1) mass reduction due to interception loss, (2) changes in the phase and amount of water delivered from the canopy to the underlying snow, (3) varying new snow density from reduced wind speed, and (4) modification of incoming longwave and shortwave radiation. Delivery effects (2) increased forest snowpack density relative to open areas, often more than 30%. Mass effects (1) and wind effects (3) decreased forest snowpack density, but generally by less than 6%. The radiation experiment (4) yielded negligible to positive effects (i.e., 0%–10%) on snowpack density. Delivery effects on density were greatest at the warmest times in the season and at the warmest site (Oregon): higher temperatures increased interception and melted intercepted snow, which then dripped to the underlying snowpack. In contrast, mass effects and radiation effects were shown to have the greatest impact on forest‐to‐open SWE differences, yielding differences greater than 30%. The study highlights the importance of delivery effects in models and the need for new types of observations to characterize how canopies influence the flux of water to the snow surface.more » « less
-
Abstract Type Ia supernova explosions (SN Ia) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (withZbetween 11 and 20) and iron group elements—for example, about 70% of the solar iron is expected to be made by SN Ia. In this work, we calculate complete abundance yields for 39 models of SN Ia explosions, based on three progenitors—a 1.4M⊙deflagration detonation model, a 1.0M⊙double detonation model, and a 0.8M⊙double detonation model—and 13 metallicities, with22Ne mass fractions of 0, 1 × 10−7, 1 × 10−6, 1 × 10−5, 1 × 10−4, 1 × 10−3, 2 × 10−3, 5 × 10−3, 1 × 10−2, 1.4 × 10−2, 5 × 10−2, and 0.1, respectively. Nucleosynthesis calculations are done using the NuGrid suite of codes, using a consistent nuclear reaction network between the models. Complete tables with yields and production factors are provided online at Zenodo:Yields (https://doi.org/10.5281/zenodo.8060323). We discuss the main properties of our yields in light of the present understanding of SN Ia nucleosynthesis, depending on different progenitor mass and composition. Finally, we compare our results with a number of relevant models from the literature.more » « less
An official website of the United States government
