The Reynolds Creek Long-Term Dataset: A long-term meteorological dataset derived from station observations in the Reynolds Creek Experimental Watershed

<h2 dir="ltr"><b>Abstract</b></h2><p dir="ltr">A forty-year gridded meteorological forcing dataset spanning the water years 1984 to 2023 (October 1<sup>st</sup> to September 30<sup>th</sup>) has been compiled for the Reynolds Creek Experimental Watershed (RCEW) in southwest Idaho, USA. This Reynold Creek Long-Term (RCLT) dataset consists of hourly, 10-meter resolution grids of air temperature, vapor pressure, precipitation mass and phase, incoming shortwave and longwave radiation, and wind speed and direction. These variables are foundational for many ecological and hydrological Land Surface Models (LSMs) used in research and operational applications and have been interpolated and calculated from hourly measurements from the dense meteorological station network within the RCEW. The elevation relief of the 240 km<sup>2</sup> RCEW spans the historical winter rain-to-snow transition, so a research application example is presented in which we show how the snow-dominated area of the basin has evolved over the forty-year data record. This dataset, stored in cloud-optimized Zarr format, enables future model development, benchmarking, and uncertainty analyses of existing models, independent validation of gridded atmospheric reanalysis datasets, and novel investigations of hydroclimatic variability across snow-dominated semi-arid environments.</p><h2 dir="ltr"><b>1 Site Description</b></h2><p dir="ltr">In 1960, the United States Congress allocated funding for an experimental research watershed to advance hydrologic research in western U.S. rangelands. Since then, the USDA Agricultural Research Service (ARS) Northwest Watershed Research Center (NWRC) has managed the scientific infrastructure in the Reynolds Creek Experimental Watershed (RCEW) in Southwest Idaho (43.205°, -116.75°). At approximately 240 km<sup>2</sup>, the RCEW is characteristic of semiarid snow-dominated environments found throughout the Great Basin and spans an elevation gradient of 1100 to 2244 meters above sea level. Characteristic of mountain climatology, a significant elevational and directional precipitation gradient exists due to the prevailing northeast-trending storms during the winter and spring, when most of the annual precipitation occurs (Hanson, 2001). Since 1984, average annual precipitation ranged from 228 mm at site RC.057 site in the low northeast elevations to 1,086 mm at site RC.163 in the highest southwest elevations of the watershed. The NWRC has a long history of publishing station-based hydrometeorological datasets, including Slaughter et al. (2001), Reba et al. (2011) and Godsey et al. (2018). This dataset updates the gridded temperature, humidity, and precipitation dataset reported by Kormos et al. (2018), includes additional wind and radiation data, and appends nine more recent years of conditions (2015-2023) that have experienced a wider range of weather variability (Monteiro & Morin, 2023).</p><h2 dir="ltr"><b>2 Instrumentation and Variable Distribution</b></h2><p dir="ltr">The Spatial Modeling for Resources Framework (SMRF; Havens et al. (2017)) was employed to distribute each of the ten land-surface meteorological variables to a 10-meter regular grid. Each forcing variable has either been empirically derived or directly interpolated from hourly station measurements across the catchment domain. Owing to the considerable length of time encompassed by the dataset, many different sensors have been deployed in the watershed over the forty-year data record with differing levels of accuracy, which are not reported here. Invalid data were preliminarily removed for all measured variables besides precipitation (which necessitated a unique approach detailed in subsection 2.6), and temporal interpolation was performed for data gaps of two hours or less. Gaps lasting longer than two hours were left empty with the foreknowledge that spatial interpolation from nearby sites in the high-density network would act as data surrogates. It is worth noting that the NWRC has always employed a full-time staff of technicians tasked with the calibration and servicing of each sensor deployed in the RCEW. The gridded interpolation methods for all ten modeled variables from the six measured variables are described in Hedrick et al. (2018) and further elaborated upon in the following subsections.</p><h3 dir="ltr"><b>3.1 Air Temperature</b></h3><p dir="ltr">Hourly measurements of air temperature (T<sub>a</sub>) were made at 35 individual sites over the 40-year period (Figure 2), with the number of sites significantly increasing after water year 2000. Measurements are currently made using various incarnations of the widely used Vaisala HMP series of temperature and humidity sensors with ventilated radiation shields.</p><p dir="ltr">A modified inverse distance weighting (IDW) approach was used to distribute T<sub>a</sub> across the 10-meter grid. In this process, the elevational trend is calculated at each time step, constrained to be negative due to the general relationship between elevation and T<sub>a</sub>, then subtracted from the station measurements to produce a temperature residual. These residuals are distributed using standard IDW and added to each grid cell’s position on the elevation gradient slope line. This approach for distributing T<sub>a</sub> improves station representativeness for areas with complex local topography.</p><h3 dir="ltr"><b>2.2 Vapor Pressure</b></h3><p dir="ltr">Gridded actual vapor pressure (e<sub>a</sub>) values were interpolated from measurements of relative humidity (RH) at the same 35 sites as T<sub>a</sub> over the 40-year data record, using the same Vaisala HMP instruments referred to in subsection 2.1. The empirical Tetens equation was used for deriving e<sub>a</sub> from RH and T<sub>a</sub>:</p><p dir="ltr">e<sub>a </sub>= RH × 0.6108 × e ^ (17.3 × T<sub>a </sub>/ 237.3 + T<sub>a</sub>). (1)</p><p dir="ltr">Station-derived e<sub>a</sub> was then distributed to the 10-meter grid using the same modified IDW detrending approach described in subsection 2.1. Dew point temperature was also calculated for the wet bulb temperature calculation (subsection 2.6.2 below) but was not stored in the RCLT dataset because it can be calculated from the vapor pressure and air temperature using the existing empirical relationships (e.g., the Clausius-Clapeyron equation).</p><h3 dir="ltr"><b>2.3 Wind</b></h3><p dir="ltr">Wind speed (u<sub>s</sub>) and direction (u<sub>dir</sub>) were measured at a total of 29 sites over the data record, though only three sites were available prior to 1994 and four sites prior to 2002. The sparseness of wind measurements through the early years of this dataset is likely a source of uncertainty in the distributed wind grids for that period, though we should note that the pre-2002 measurements captured the full elevation gradient in RCEW at low (RC.076), mid (RC.127), and high elevation sites (RC.176).</p><p dir="ltr">Station measurements of wind were distributed to the 10-meter grid using the maximum upwind slope (maxus) terrain parameter described in Winstral, et al. (2002) and Winstral et al. (2009). In short, the underlying digital elevation model (DEM) is used to calculate a maxus value (in degrees) over a user-defined upwind distance (here 300 meters) for all possible upwind directions (0˚ to 360˚) in 5˚ increments. The resulting 72 layers of maxus grids are stored in a lookup library. Then, for each station the measured wind speed is adjusted to simulate what the wind speed would have been on a flat surface (‘flatwind’) using the maxus value for the measured wind direction at the site. Once the adjusted ‘flatwind’ speeds and wind direction components have been distributed across the entire grid using standard IDW, the distributed wind directions are used to find the maxus value for each grid cell and the distributed ‘flatwind’ speeds are converted back to actual wind speeds.</p><p dir="ltr">For the gridded dataset, wind speed and direction were converted into U- and V-components to match the conventions of NWP models such as the WRF and HRRR models. The U-component represents the East-West wind speed, with positive values indicating wind out of the west, while the V-component represents the North-South wind speed, positive value indicating wind out of the south.</p><h3 dir="ltr"><b>2.4 Shortwave Radiation</b></h3><p dir="ltr">Gridded incoming shortwave radiation (S<sub>in</sub>) has been measured at 23 sites across the RCEW but cannot be directly spatially interpolated from measurements due to the complex terrain and the variable vegetation canopy present across the catchment. Instead, a three-step process produced the hourly gridded S<sub>in</sub> product.</p><ol><li>Hourly clear sky atmospheric S<sub>in</sub> was modeled, then corrected for surrounding terrain in each 10-meter grid cell following Dozier (1980) and Dubayah (1994), respectively.</li><li>Station measurements of S<sub>in</sub> were divided by clear sky radiation to derive a cloud factor (C<sub>fac</sub>) at each station pixel (C<sub>fac</sub> =1 represents cloud-free conditions), which was then distributed using standard IDW across the domain.</li><li>Canopy-corrected S<sub>in</sub> values were estimated using empirical relationships presented by Link and Marks (1999), where direct beam shortwave radiation under canopy (R<sub>b</sub>) can be represented by:</li></ol><p dir="ltr">R<sub>b </sub>= S<sub>b,in</sub> × e ^ (-μh / cos(⁡θ)).</p><p dir="ltr">In this equation, S<sub>b,in</sub> is the above canopy cloud corrected direct beam radiation, μ is a canopy extinction coefficient, h is the height of the canopy, and θ is the solar zenith angle. Diffuse shortwave radiation under canopy (R<sub>d</sub>) is computed by adjusting the above canopy cloud corrected diffuse radiation (S<sub>d,in</sub>) by the canopy optical transmissivity (τ):</p><p dir="ltr">R<sub>d</sub>= τ × S<sub>d,in</sub>.</p><p dir="ltr">The terrain-, cloud-, and can...

Tags: earth-sciences, environmental-sciences, hydrology, land-surface-models, meteorology, snow-hydrology, spatial-data-and-applications
Organization: Department of Agriculture