skip to main content

Title: Constructing retrospective gridded daily weather data for agro‐hydrological applications in Oklahoma

Regional, automated meteorological networks, such as the Oklahoma Mesonet can potentially provide high quality forcing data for generating gridded surfaces, but proven methods of interpolating weather variables between the station locations are needed. We compared two interpolation methods, ordinary kriging (OK) and empirical Bayesian kriging (EBK), with and without using long‐term climate imprints (CI), for creating spatially continuous, daily weather datasets. Daily meteorological variables (maximum and minimum temperature, solar radiation, and precipitation) from the Oklahoma Mesonet for the period 1997–2014 were interpolated using geoprocessing tools in ArcGIS. Cross‐validation was used for evaluation of interpolation methods, with 90% of sites chosen randomly for the training set and the remaining 10% left for validation. For all interpolation approaches, cross‐validation showed coefficient of determination (R2) values of .99 and .98 for daily maximum and minimum air temperatures, with mean absolute error (MAE) ranging from ±0.45–0.50 °C for maximum temperature and ±0.77–0.80 °C for minimum temperature. Likewise, for daily solar radiation,R2values of .94 and .93 showed overall good prediction accuracy with MAE values 1.00 and 1.01 MJ m–2 d–1for EBK and OK, respectively. However, for rainfall, all methods yieldedR2values ≤.67, suggesting a need for more effective interpolation method. Based on its lower computational time and lower input data requirement, OK appears preferable to the other approaches tested here to provide the daily weather data for gridded models in Oklahoma and other regions with similar monitoring networks.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Agrosystems, Geosciences & Environment
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Gridded monthly rainfall estimates can be used for a number of research applications, including hydrologic modeling and weather forecasting. Automated interpolation algorithms, such as the “autoKrige” function in R, can produce gridded rainfall estimates that validate well but produce unrealistic spatial patterns. In this work, an optimized geostatistical kriging approach is used to interpolate relative rainfall anomalies, which are then combined with long-term means to develop the gridded estimates. The optimization consists of the following: 1) determining the most appropriate offset (constant) to use when log-transforming data; 2) eliminating poor quality data prior to interpolation; 3) detecting erroneous maps using a machine learning algorithm; and 4) selecting the most appropriate parameterization scheme for fitting the model used in the interpolation. Results of this effort include a 30-yr (1990–2019), high-resolution (250-m) gridded monthly rainfall time series for the state of Hawai‘i. Leave-one-out cross validation (LOOCV) is performed using an extensive network of 622 observation stations. LOOCV results are in good agreement with observations (R2= 0.78; MAE = 55 mm month−1; 1.4%); however, predictions can underestimate high rainfall observations (bias = 34 mm month−1; −1%) due to a well-known smoothing effect that occurs with kriging. This research highlights the fact that validation statistics should not be the sole source of error assessment and that default parameterizations for automated interpolation may need to be modified to produce realistic gridded rainfall surfaces. Data products can be accessed through the Hawai‘i Data Climate Portal (HCDP;

    Significance Statement

    A new method is developed to map rainfall in Hawai‘i using an optimized geostatistical kriging approach. A machine learning technique is used to detect erroneous rainfall maps and several conditions are implemented to select the optimal parameterization scheme for fitting the model used in the kriging interpolation. A key finding is that optimization of the interpolation approach is necessary because maps may validate well but have unrealistic spatial patterns. This approach demonstrates how, with a moderate amount of data, a low-level machine learning algorithm can be trained to evaluate and classify an unrealistic map output.

    more » « less
  2. Abstract

    Improved understanding of the effects of meteorological conditions on the transmission of SARS-CoV-2, the causative agent for COVID-19 disease, is needed. Here, we estimate the relationship between air temperature, specific humidity, and ultraviolet radiation and SARS-CoV-2 transmission in 2669 U.S. counties with abundant reported cases from March 15 to December 31, 2020. Specifically, we quantify the associations of daily mean temperature, specific humidity, and ultraviolet radiation with daily estimates of the SARS-CoV-2 reproduction number (Rt) and calculate the fraction ofRtattributable to these meteorological conditions. Lower air temperature (within the 20–40 °C range), lower specific humidity, and lower ultraviolet radiation were significantly associated with increasedRt. The fraction ofRtattributable to temperature, specific humidity, and ultraviolet radiation were 3.73% (95% empirical confidence interval [eCI]: 3.66–3.76%), 9.35% (95% eCI: 9.27–9.39%), and 4.44% (95% eCI: 4.38–4.47%), respectively. In total, 17.5% ofRtwas attributable to meteorological factors. The fractions attributable to meteorological factors generally were higher in northern counties than in southern counties. Our findings indicate that cold and dry weather and low levels of ultraviolet radiation are moderately associated with increased SARS-CoV-2 transmissibility, with humidity playing the largest role.

    more » « less
  3. Abstract

    Gridded air temperature data are required in various fields such as ecological modeling, weather forecasting, and surface energy balance assessment. In this work, a piecewise multiple linear regression model is used to produce high‐resolution (250 m) daily maximum (Tmax), minimum (Tmin), and mean (Tmean) near‐surface air temperature maps for the State of Hawaiʻi for a 32‐year period (1990–2021). Multiple meteorological and geographical variables such as the elevation, daily rainfall, coastal distance index, leaf area index, albedo, topographic position index, and wind speed are independently tested to determine the most well‐suited predictor variables for optimal model performance. During the mapping process, input data scarcity is addressed first by gap‐filling critical stations at high elevation using a predetermined linear relationship with other strongly‐correlated stations, and second, by supplementing the training dataset with station data from neighboring islands. Despite the numerous covariates physically linked to temperature, the most parsimonious model selection uses elevation as its sole predictor, and the inclusion of the additional variables results in increased cross‐validation errors. The mean absolute error of resultant estimatedTmaxandTminmaps over the Hawaiian Islands from 1990 to 2021 is 1.7°C and 1.3°C, respectively. Corresponding bias values are 0.01°C and −0.13°C, respectively for the same variables. Overall, the results show the proposed methodology can robustly generate daily air temperature maps from point‐scale measurements over complex topography.

    more » « less
  4. Abstract

    Soil moisture spatial patterns with length scales of 1‐100 km influence hydrological, ecological, and agricultural processes, but the footprint or support volume of existing monitoring systems, for example, satellite‐based radiometers and sparse in situ monitoring networks, is often either too large or too small to effectively observe these mesoscale patterns. This measurement scale gap hinders our understanding of soil water processes and complicates calibration and validation of hydrologic models and soil moisture satellites. One possible solution is to utilize geostatistical techniques that have proven effective for mapping static patterns in soil properties. The objective of this study was to determine how effectively dynamic, mesoscale soil moisture patterns can be mapped by applying regression kriging to the data from a sparse, large‐scale in situ network. The fully automated system developed here uses several data sets: daily soil moisture measurements from the Oklahoma Mesonet, sand content estimates from the Natural Resource Conservation Service Soil Survey Geographic Database, and an antecedent precipitation index computed from National Weather Service multisensor precipitation estimates. A multiple linear regression model is fitted daily to the observed data, and the residuals of that model are used in a semivariogram estimation and kriging routine to produce daily statewide maps of soil moisture at 5‐, 25‐, and 60‐cm depths at 800‐m resolution. During over 3 years of operation, this mapping system has revealed complex, dynamic, and depth‐specific mesoscale patterns, reflecting the shifting influences of both soil texture and precipitation, with a mean absolute error of ≤0.0576 cm3/cm3across all three depths.

    more » « less
  5. Daily summary data at NPP P-TOBO met station. Average/maximum/minimum air temperature; average/maximum relative humidity and wind speed and average wind direction; solar radiation; albedo. These are measured and calculated based on 1-second scan rate of all sensors located at an automated meteorological station installed at Jornada LTER NPP P-TOBO site. Wind speed is measured at 75 cm, 150 cm, and 300 cm, wind direction at approximately 3m, and air temperature and relative humidity at approximate 2.5m. Solar radiation is measured at 3m. This climate station is operated by the Jornada LTER Program. This is an ONGOING dataset. 
    more » « less