Spatial interpolation of large climate data sets using bivariate thin plate smoothing splines

2006 ◽  
Vol 21 (12) ◽  
pp. 1684-1694 ◽  
Author(s):  
P.A. Hancock ◽  
M.F. Hutchinson
2010 ◽  
Vol 27 ◽  
pp. 91-98 ◽  
Author(s):  
S. van der Heijden ◽  
U. Haberlandt

Abstract. For ecohydrological modeling climate variables are needed on subbasin basis. Since they usually originate from point measurements spatial interpolation is required during preprocessing. Different interpolation methods yield data of varying quality, which can strongly influence modeling results. Four interpolation methods to be compared were selected: nearest neighbour, inverse distance, ordinary kriging, and kriging with external drift (Goovaerts, 1997). This study presents three strategies to evaluate the influence of the interpolation method on the modeling results of discharge and nitrate load in the river in a mesoscale river catchment (~1000 km2) using the Soil and Water Assessment Tool (SWAT, Neitsch et al., 2005) model: I. Automated calibration of the model with a mixed climate data set and consecutive application of the four interpolated data sets. II. Consecutive automated calibration of the model with each of the four climate data sets. III. Random generation of 1000 model parameter sets and consecutive application of the four interpolated climate data sets on each of the 1000 realisations, evaluating the number of realisations above a certain quality criterion threshold. Results show that strategies I and II are not suitable for evaluation of the quality of the interpolated data. Strategy III however proves a significant influence of the interpolation method on nitrate modeling. A rank order from the simplest to the most sophisticated method is visible, with kriging with external drift (KED) outperforming all others. Responsible for this behaviour is the variable temperature, which benefits most from more sophisticated methods and at the same time is the main driving force for the nitrate cycle. The missing influence of the interpolation methods on discharge modeling is explained by a much higher measuring network density for precipitation than for all other climate variables.


2021 ◽  
Author(s):  
Jouke de Baar ◽  
Gerard van der Schrier ◽  
Irene Garcia-Marti ◽  
Else van den Besselaar

<p><strong>Objective</strong></p><p>The purpose of the European Copernicus Climate Change Service (C3S) is to support society by providing information about the past, present and future climate. For the service related to <em>in-situ</em> observations, one of the objectives is to provide high-resolution (0.1x0.1 and 0.25x0.25 degrees) gridded wind speed fields. The gridded wind fields are based on ECA&D daily average station observations for the period 1970-2020.</p><p><strong>Research question</strong> </p><p>We address the following research questions: [1] How efficiently can we provide the gridded wind fields as a statistically reliable ensemble, in order to represent the uncertainty of the gridding? [2] How efficiently can we exploit high-resolution geographical auxiliary variables (e.g. digital elevation model, terrain roughness) to augment the station data from a sparse network, in order to provide gridded wind fields with high-resolution local features?</p><p><strong>Approach</strong></p><p>In our analysis, we apply greedy forward selection linear regression (FSLR) to include the high-resolution effects of the auxiliary variables on monthly-mean data. These data provide a ‘background’ for the daily estimates. We apply cross-validation to avoid FSLR over-fitting and use full-cycle bootstrapping to create FSLR ensemble members. Then, we apply Gaussian process regression (GPR) to regress the daily anomalies. We consider the effect of the spatial distribution of station locations on the GPR gridding uncertainty.</p><p>The goal of this work is to produce several decades of daily gridded wind fields, hence, computational efficiency is of utmost importance. We alleviate the computational cost of the FSLR and GPR analyses by incorporating greedy algorithms and sparse matrix algebra in the analyses.</p><p><strong>Novelty</strong>   </p><p>The gridded wind fields are calculated as a statistical ensemble of realizations. In the present analysis, the ensemble spread is based on uncertainties arising from the auxiliary variables as well as from the spatial distribution of stations.</p><p>Cross-validation is used to tune the GPR hyper parameters. Where conventional GPR hyperparameter tuning aims at an optimal prediction of the gridded mean, instead, we tune the GPR hyperparameters for optimal prediction of the gridded ensemble spread.</p><p>Building on our experience with providing similar gridded climate data sets, this set of gridded wind fields is a novel addition to the E-OBS climate data sets.</p>


Climate ◽  
2019 ◽  
Vol 7 (5) ◽  
pp. 68 ◽  
Author(s):  
Flora Gofa ◽  
Anna Mamara ◽  
Manolis Anadranistakis ◽  
Helena Flocas

The creation of realistic gridded precipitation fields improves our understanding of the observed climate and is necessary for validating climate model output for a wide range of applications. The challenge in trying to represent the highly variable nature of precipitation is to overcome the lack of density of observations in both time and space. Data sets of mean monthly and annual precipitations were developed for Greece in gridded format with an analysis of 30 arcsec (∼800 m) based on data from 1971 to 2000. One hundred and fifty-seven surface stations from two different observation networks were used to cover a satisfactory range of elevations. Station data were homogenized and subjected to quality control to represent changes in meteorological conditions rather than changes in the conditions under which the observations were made. The Meteorological Interpolation based on Surface Homogenized Data Basis (MISH) interpolation method was used to develop data sets that reproduce, as closely as possible, the spatial climate patterns over the region of interest. The main geophysical factors considered for the interpolation of mean monthly precipitation fields were elevation, latitude, incoming solar irradiance, Euclidian distance from the coastline, and land-to-sea percentage. Low precipitation interpolation uncertainties estimated with the cross-validation method provided confidence in the interpolation method. The resulting high-resolution maps give an overall realistic representation of precipitation, especially in fall and winter, with a clear longitudinal dependence on precipitation decreasing from western to eastern continental Greece.


2017 ◽  
Vol 17 (23) ◽  
pp. 14593-14629 ◽  
Author(s):  
Craig S. Long ◽  
Masatomo Fujiwara ◽  
Sean Davis ◽  
Daniel M. Mitchell ◽  
Corwin J. Wright

Abstract. Two of the most basic parameters generated from a reanalysis are temperature and winds. Temperatures in the reanalyses are derived from conventional (surface and balloon), aircraft, and satellite observations. Winds are observed by conventional systems, cloud tracked, and derived from height fields, which are in turn derived from the vertical temperature structure. In this paper we evaluate as part of the SPARC Reanalysis Intercomparison Project (S-RIP) the temperature and wind structure of all the recent and past reanalyses. This evaluation is mainly among the reanalyses themselves, but comparisons against independent observations, such as HIRDLS and COSMIC temperatures, are also presented. This evaluation uses monthly mean and 2.5° zonal mean data sets and spans the satellite era from 1979–2014. There is very good agreement in temperature seasonally and latitudinally among the more recent reanalyses (CFSR, MERRA, ERA-Interim, JRA-55, and MERRA-2) between the surface and 10 hPa. At lower pressures there is increased variance among these reanalyses that changes with season and latitude. This variance also changes during the time span of these reanalyses with greater variance during the TOVS period (1979–1998) and less variance afterward in the ATOVS period (1999–2014). There is a distinct change in the temperature structure in the middle and upper stratosphere during this transition from TOVS to ATOVS systems. Zonal winds are in greater agreement than temperatures and this agreement extends to lower pressures than the temperatures. Older reanalyses (NCEP/NCAR, NCEP/DOE, ERA-40, JRA-25) have larger temperature and zonal wind disagreement from the more recent reanalyses. All reanalyses to date have issues analysing the quasi-biennial oscillation (QBO) winds. Comparisons with Singapore QBO winds show disagreement in the amplitude of the westerly and easterly anomalies. The disagreement with Singapore winds improves with the transition from TOVS to ATOVS observations. Temperature bias characteristics determined via comparisons with a reanalysis ensemble mean (MERRA, ERA-Interim, JRA-55) are similarly observed when compared with Aura HIRDLS and Aura MLS observations. There is good agreement among the NOAA TLS, SSU1, and SSU2 Climate Data Records and layer mean temperatures from the more recent reanalyses. Caution is advised for using reanalysis temperatures for trend detection and anomalies from a long climatology period as the quality and character of reanalyses may have changed over time.


Author(s):  
Jeffrey Sukharev ◽  
Chaoli Wang ◽  
Kwan-Liu Ma ◽  
Andrew T. Wittenberg

Sign in / Sign up

Export Citation Format

Share Document