skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, July 12 until 2:00 AM ET on Saturday, July 13 due to maintenance. We apologize for the inconvenience.

Title: Optimizing Automated Kriging to Improve Spatial Interpolation of Monthly Rainfall over Complex Terrain

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
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Hydrometeorology
Page Range / eLocation ID:
p. 561-572
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Patterns ofδ18O andδ2H in Earth's precipitation provide essential scientific data for use in hydrological, climatological, ecological and forensic research. Insufficient global spatial data coverage promulgated the use of gridded datasets employing geostatistical techniques (isoscapes) for spatiotemporally coherent isotope predictions. Cluster‐based isoscape regionalization combines the advantages of local or regional prediction calibrations into a global framework. Here we present a revision of a Regionalized Cluster‐Based Water Isotope Prediction model (RCWIP2) incorporating new isotope data having extensive spatial coverage and a wider array of predictor variables combined with high‐resolution gridded climatic data. We introduced coupling ofδ18O andδ2H (e.g.,d‐excess constrained) in the model predictions to prevent runaway isoscapes when each isotope is modelled separately and cross‐checked observed versus modelledd‐excess values. We improved model error quantification by adopting full uncertainty propagation in all calculations. RCWIP2 improved the RMSE over previous isoscape models by ca. 0.3 ‰ forδ18O and 2.5 ‰ forδ2H with an uncertainty <1.0 ‰ forδ18O and < 8 ‰ forδ2H for most regions of the world. The determination of the relative importance of each predictor variable in each ecoclimatic zone is a new approach to identify previously unrecognized climatic drivers on mean annual precipitationδ18O andδ2H. The improved RCWIP2 isoscape grids and maps (season, monthly, annual, regional) are available for download at

    more » « less
  2. Abstract

    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
  3. Abstract

    High temporal and spatial resolution precipitation datasets are essential for hydrological and flood modeling to assist water resource management and emergency responses, particularly for small watersheds, such as those in Hawai‘i in the United States. Unfortunately, fine temporal (subdaily) and spatial (<1 km) resolutions of rainfall datasets are not always readily available for applications. Radar provides indirect measurements of the rain rate over a large spatial extent with a reasonable temporal resolution, while rain gauges provide “ground truth.” There are potential advantages to combining the two, which have not been fully explored in tropical islands. In this study, we applied kriging with external drift (KED) to integrate hourly gauge and radar rainfall into a 250 m × 250 m gridded dataset for the tropical island of O‘ahu. The results were validated with leave-one-out cross validation for 18 severe storm events, including five different storm types (e.g., tropical cyclone, cold front, upper-level trough, kona low, and a mix of upper-level trough and kona low), and different rainfall structures (e.g., stratiform and convective). KED-merged rainfall estimates outperformed both the radar-only and gauge-only datasets by 1) reducing the error from radar rainfall and 2) improving the underestimation issues from gauge rainfall, especially during convective rainfall. We confirmed the KED method can be used to merge radar with gauge data to generate reliable rainfall estimates, particularly for storm events, on mountainous tropical islands. In addition, KED rainfall estimates were consistently more accurate in depicting spatial distribution and maximum rainfall value within various storm types and rainfall structures.

    Significance Statement

    The results of this study show the effectiveness of utilizing kriging with external drift (KED) in merging gauge and radar rainfall data to produce highly accurate, reliable rainfall estimates in mountainous tropical regions, such as O‘ahu. The validated KED dataset, with its high temporal and spatial resolutions, offers a valuable resource for various types of rainfall-related research, particularly for extreme weather response and rainfall intensity analyses in Hawai’i. Our findings improve the accuracy of rainfall estimates and contribute to a deeper understanding of the performance of various rainfall estimation methods under different storm types and rainfall structures in a mountainous tropical setting.

    more » « less
  4. Abstract

    The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB), funded by the US National Science Foundation, National Institutes of Health, and Department of Energy, has served structural biologists and Protein Data Bank (PDB) data consumers worldwide since 1999. RCSB PDB, a founding member of the Worldwide Protein Data Bank (wwPDB) partnership, is the US data center for the global PDB archive housing biomolecular structure data. RCSB PDB is also responsible for the security of PDB data, as the wwPDB‐designated Archive Keeper. Annually, RCSB PDB serves tens of thousands of three‐dimensional (3D) macromolecular structure data depositors (using macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro‐electron diffraction) from all inhabited continents. RCSB PDB makes PDB data available from its research‐focusedRCSB.orgweb portal at no charge and without usage restrictions to millions of PDB data consumers working in every nation and territory worldwide. In addition, RCSB PDB operates an outreach and educationPDB101.RCSB.orgweb portal that was used by more than 800,000 educators, students, and members of the public during calendar year 2020. This invited Tools Issue contribution describes (i) how the archive is growing and evolving as new experimental methods generate ever larger and more complex biomolecular structures; (ii) the importance of data standards and data remediation in effective management of the archive and facile integration with more than 50 external data resources; and (iii) new tools and features for 3D structure analysis and visualization made available during the past yearviatheRCSB.orgweb portal.

    more » « less
  5. This paper discusses the design and implementation of the Hawai‘i Rainfall Analysis and Mapping Application (HI-RAMA) decision support tool. HI-RAMA provides researchers and community stakeholders interactive access to and visualization of hosted historical and near-real-time monthly rainfall maps and aggregated rainfall station observational data for the State of Hawai‘i. The University of Hawai‘i Information Technology Services Cyberinfrastructure team in partnership with members of the Hawai‘i Established Program to Stimulate Competitive Research (EPSCoR) ‘Ike Wai project team developed this application as part of the ‘Ike Wai Gateway to support water sustainability research for the state of Hawai‘i. This tool is designed to provide user-friendly access to information that can reveal the impacts of climate changes related to precipitation so users can make data-driven decisions. 
    more » « less