Use of variability modes to evaluate AR4 climate models over the Euro-Atlantic region

This paper analyzes the ability of the multi-model simulations from the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) to simulate the main leading modes of variability over the Euro-Atlantic region in winter: the North-Atlantic Oscillation (NAO), the Scandinavian mode (SCAND), the East/Atlantic Oscillation (EA) and the East Atlantic/Western Russia mode (EA/WR). These modes of variability have been evaluated both spatially, by analyzing the intensity and location of their anomaly centres, as well as temporally, by focusing on the probability density functions and e-folding time scales. The choice of variability modes as a tool for climate model assessment can be justified by the fact that modes of variability determine local climatic conditions and their likely change may have important implications for future climate changes. It is found that all the models considered are able to simulate reasonably well these four variability modes, the SCAND being the mode which is best spatially simulated. From a temporal point of view the NAO and SCAND modes are the best simulated. UKMO-HadGEM1 and CGCM3.1(T63) are the models best at reproducing spatial characteristics, whereas CCSM3 and CGCM3.1(T63) are the best ones with regard to the temporal features. GISS-AOM is the model showing the worst performance, in terms of both spatial and temporal features. These results may bring new insight into the selection and use of specific models to simulate Euro-Atlantic climate, with some models being clearly more successful in simulating patterns of temporal and spatial variability than others.     

Climate change in terms of modes of atmospheric variability and circulation regimes over southern South America

The simulated low-frequency variability patterns of the atmospheric circulation, ranging from interannual to interdecadal timescales, are studied in an area encompassing southern South America. The experiment is a transient simulation performed with the IPSL CCM2 coupled global model, in which the greenhouse forcing is continuously increasing. The main modes of low-frequency variability are found to remain stationary throughout the simulation, suggesting they depend more on the internal dynamics of the atmospheric flow than on its external forcing. Inspection of the circulation regimes that represent the more recurrent patterns at interannual and interdecadal timescales showed that climate change manifests itself as a change in regime population, suggesting that the negative phase of the Antarctic Oscillation-like pattern becomes more frequented in a climate change scenario. Changes of regime occurrence are superimposed to a positive trend whose spatial pattern is reminiscent of the structure of the Antarctic Oscillation-mode of variability. Moreover, it resembles the spatial patterns of those regimes that show a significant change in population. The change in regime frequencies of the circulation patterns of low-frequency variability are in opposite phase with respect to the trend, thus, the behaviour of these patterns of variability, superimposed to a changing mean state, modulates the climate change signal. The analysis of the high frequencies, in terms of recurrent patterns representing intraseasonal and synoptic-scale of variability, shows no significant changes in regime characteristics, concerning both spatial and temporal behaviour.     

A method for regional climate change detection using smooth temporal patterns

This paper introduces an original method for climate change detection, called temporal optimal detection method. The method consists in searching for a smooth temporal pattern in the observations. This pattern can be either the response of the climate system to a specific forcing or to a combination of forcings. Many characteristics of this new method are different from those of the classical “optimal fingerprint” method. It allows to infer the spatial distribution of the detected signal, without providing any spatial guess pattern. The spatial properties of the internal climate variability doesn’t need to be estimated either. The estimation of such quantities being very challenging at regional scale, the proposed method is particularly well-suited for such scale. The efficiency of the method is illustrated by applying it on real homogenized datasets of temperatures and precipitation over France. A multimodel detection is performed in both cases, using an ensemble of atmosphere-ocean general circulation models for estimating the temporal patterns. Regarding temperatures, new results are highlighted, especially by showing that a change is detected even after removing the uniform part of the warming. The sensitivity of the method is discussed in this case, relatively to the computation of the temporal patterns and to the choice of the model. The method also allows to detect a climate change signal in precipitation. This change impacts the spatial distribution of the precipitation more than the mean over the domain. The ability of the method to provide an estimate of the spatial distribution of the change following the prescribed temporal patterns is also illustrated.     

Assessment of a climate model to reproduce rainfall variability and extremes over Southern Africa

It is increasingly accepted that any possible climate change will not only have an influence on mean climate but may also significantly alter climatic variability. A change in the distribution and magnitude of extreme rainfall events (associated with changing variability), such as droughts or flooding, may have a far greater impact on human and natural systems than a changing mean. This issue is of particular importance for environmentally vulnerable regions such as southern Africa. The sub-continent is considered especially vulnerable to and ill-equipped (in terms of adaptation) for extreme events, due to a number of factors including extensive poverty, famine, disease and political instability. Rainfall variability and the identification of rainfall extremes is a function of scale, so high spatial and temporal resolution data are preferred to identify extreme events and accurately predict future variability. The majority of previous climate model verification studies have compared model output with observational data at monthly timescales. In this research, the assessment of ability of a state of the art climate model to simulate climate at daily timescales is carried out using satellite-derived rainfall data from the Microwave Infrared Rainfall Algorithm (MIRA). This dataset covers the period from 1993 to 2002 and the whole of southern Africa at a spatial resolution of 0.1° longitude/latitude. This paper concentrates primarily on the ability of the model to simulate the spatial and temporal patterns of present-day rainfall variability over southern Africa and is not intended to discuss possible future changes in climate as these have been documented elsewhere. Simulations of current climate from the UK Meteorological Office Hadley Centre’s climate model, in both regional and global mode, are firstly compared to the MIRA dataset at daily timescales. Secondly, the ability of the model to reproduce daily rainfall extremes is assessed, again by a comparison with extremes from the MIRA dataset. The results suggest that the model reproduces the number and spatial distribution of rainfall extremes with some accuracy, but that mean rainfall and rainfall variability is under-estimated (over-estimated) over wet (dry) regions of southern Africa.     

Hydroclimate risk to economic growth in sub-Saharan Africa

In order to plan strategies for adaptation to climate change, the current effects of climate on economic growth need to be understood. This study reviews evidence of climate effects on economic growth and presents original analysis of the effect in Sub-Saharan Africa (SSA). Case studies from the literature demonstrate that historically, climate has had significant and negative effects on household income, agricultural productivity and economic growth in SSA. This study focuses on the effects hydroclimatic variability on economic growth in the countries of SSA. We utilize a new national level precipitation statistic that incorporates spatial and temporal variability within each country. Country level economic growth statistics are analyzed in panel regressions. Persistent negative precipitation anomalies (drought) are found to be the most significant climate influence on GDP per capita growth. Temperature and precipitation variability show significant effects in some cases. Results imply the consideration of hydroclimatic risks, namely drought, may be the priority concern for adaptation to a changing climate for Sub-Saharan Africa. This conclusion is contrary to the premise of many climate change impact assessments that focus on temperature increases as the primary concern.     

A multi-proxy reconstruction of spatial and temporal variations in Asian summer temperatures over the last millennium

To investigate climate variability in Asia during the last millennium, the spatial and temporal evolution of summer (June–July–August; JJA) temperature in eastern and south-central Asia is reconstructed using multi-proxy records and the regularized expectation maximization (RegEM) algorithm with truncated total least squares (TTLS), under a point-by-point regression (PPR) framework. The temperature index reconstructions show that the late 20th century was the warmest period in Asia over the past millennium. The temperature field reconstructions illustrate that temperatures in central, eastern, and southern China during the 11th and 13th centuries, and in western Asia during the 12th century, were significantly higher than those in other regions, and comparable to levels in the 20th century. Except for the most recent warming, all identified warm events showed distinct regional expressions and none were uniform over the entire reconstruction area. The main finding of the study is that spatial temperature patterns have, on centennial time-scales, varied greatly over the last millennium. Moreover, seven climate model simulations, from the Coupled Model Intercomparison Project Phase 5 (CMIP5), over the same region of Asia, are all consistent with the temperature index reconstruction at the 99 % confidence level. Only spatial temperature patterns extracted as the first empirical orthogonal function (EOF) from the GISS-E2-R and MPI-ESM-P model simulations are significant and consistent with the temperature field reconstruction over the past millennium in Asia at the 90 % confidence level. This indicates that both the reconstruction and the simulations depict the temporal climate variability well over the past millennium. However, the spatial simulation or reconstruction capability of climate variability over the past millennium could be still limited. For reconstruction, some grid points do not pass validation tests and reveal the need for more proxies with high temporal resolution, accurate dating, and sensitive temperature signals, especially in central Asia and before AD 1400.     

USE OF A STOCHASTIC WEATHER GENERATOR IN THE DEVELOPMENT OF CLIMATE CHANGE SCENARIOS

Climate change scenarios with a high spatial and temporal resolution are required in the evaluation of the effects of climate change on agricultural potential and agricultural risk. Such scenarios should reproduce changes in mean weather characteristics as well as incorporate the changes in climate variability indicated by the global climate model (GCM) used. Recent work on the sensitivity of crop models and climatic extremes has clearly demonstrated that changes in variability can have more profound effects on crop yield and on the probability of extreme weather events than simple changes in the mean values. The construction of climate change scenarios based on spatial regression downscaling and on the use of a local stochastic weather generator is described. Regression downscaling translated the coarse resolution GCM grid-box predictions of climate change to site-specific values. These values were then used to perturb the parameters of the stochastic weather generator in order to simulate site-specific daily weather data. This approach permits the incorporation of changes in the mean and variability of climate in a consistent and computationally inexpensive way. The stochastic weather generator used in this study, LARS-WG, has been validated across Europe and has been shown to perform well in the simulation of different weather statistics, including those climatic extremes relevant to agriculture. The importance of downscaling and the incorporation of climate variability are demonstrated at two European sites where climate change scenarios were constructed using the UK Met. Office high resolution GCM equilibrium and transient experiments.     

THE IMPACT OF CLIMATE CHANGE ON THE SNOW COVER PATTERN IN ESTONIA

The paper deals with problems of temporal and spatial variability of snow cover duration, of correlation between snow cover and winter mean air temperature patterns and of the impact of climate change on the snow cover pattern in Estonia. Snow cover fields are presented in form of IDRISI raster images. Snow cover duration measured at ca 100 stations and observation points have been interpolated into raster cells. On the base of time series of raster images, a map of mean territorial distribution of snow cover duration is calculated. Estonia is characterized by a great spatial variability of snow cover mostly caused by the influence of the Baltic Sea. General regularities of snow cover pattern are determined. A 104-year time series of spatial mean values of snow cover duration is composed and analyzed. A decreasing trend and periodical fluctuations have detected. Standardized principal component analysis is used for the time series of IDRISI raster images. It enables to study the influence of different factors on the formation of snow cover fields and territorial extent of coherent fluctuations. Correlation between snow cover duration and winter mean air temperature fields is analyzed. A spatial regression model is created for estimation of the influence of climate change on snow cover pattern in Estonia. Using incremental climate change scenarios (2 °C, 4 °C and 6 °C of warming in winter) mean decrease of snow cover duration in different regions in Estonia is calculated. According to results of model calculation, the highest decrease of snow cover duration will be take place on islands and in the coastal region of West Estonia. A permanent snow cover may not form at all. In the areas with maximum snow cover duration in North-East and South-East Estonia, that decrease should be much lower.     

A European pattern climatology 1766–2000

Using monthly independently reconstructed gridded European fields for the 500 hPa geopotential height, temperature, and precipitation covering the last 235 years we investigate the temporal and spatial evolution of these key climate variables and assess the leading combined patterns of climate variability. Seasonal European temperatures show a positive trend mainly over the last 40 years with absolute highest values since 1766. Precipitation indicates no clear trend. Spatial correlation technique reveals that winter, spring, and autumn covariability between European temperature and precipitation is mainly influenced by advective processes, whereas during summer convection plays the dominant role. Empirical Orthogonal Function analysis is applied to the combined fields of pressure, temperature, and precipitation. The dominant patterns of climate variability for winter, spring, and autumn resemble the North Atlantic Oscillation and show a distinct positive trend during the past 40 years for winter and spring. A positive trend is also detected for summer pattern 2, which reflects an increased influence of the Azores High towards central Europe and the Mediterranean coinciding with warm and dry conditions. The question to which extent these recent trends in European climate patterns can be explained by internal variability or are a result of radiative forcing is answered using cross wavelets on an annual basis. Natural radiative forcing (solar and volcanic) has no imprint on annual European climate patterns. Connections to CO₂ forcing are only detected at the margins of the wavelets where edge effects are apparent and hence one has to be cautious in a further interpretation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Influences of Climate Change and Its Interannual Variability on Surface Energy Fluxes from 1948 to 2000

Understanding changes in land surface processes over the past several decades requires knowledge of trends and interannual variability in surface energy fluxes in response to climate change. In our study, the Community Land Model version 3.5 (CLM3.5), driven by the latest updated hybrid reanalysis-observational surface climate data from Princeton University, is used to obtain global distributions of surface energy fluxes during 1948 to 2000. Based on the climate data and simulation results, long-term trends and interannual variability (IAV) of both climatic variables and surface energy fluxes for this span of 50+ years are derived and analyzed. Regions with strong long-term trends and large IAV for both climatic variables and surface energy fluxes are identified. These analyses reveal seasonal variations in the spatial patterns of climate and surface fluxes; however, spatial patterns in trends and IAV for surface energy fluxes over the past ～50 years do not fully correspond to those for climatic variables, indicating complex responses of land surfaces to changes in the climatic forcings.     

Climate variability and trends in downscaled high-resolution simulations and projections over Metropolitan France

In order to fulfill the society demand for climate information at the spatial scale allowing impact studies, long-term high-resolution climate simulations are produced, over an area covering metropolitan France. One of the major goals of this article is to investigate whether such simulations appropriately simulate the spatial and temporal variability of the current climate, using two simulation chains. These start from the global IPSL-CM4 climate model, using two regional models (LMDz and MM5) at moderate resolution (15–20 km), followed with a statistical downscaling method in order to reach a target resolution of 8 km. The statistical downscaling technique includes a non-parametric method that corrects the distribution by using high-resolution analyses over France. First the uncorrected simulations are evaluated against a set of high-resolution analyses, with a focus on temperature and precipitation. Uncorrected downscaled temperatures suffer from a cold bias that is present in the global model as well. Precipitations biases have a season- and model-dependent behavior. Dynamical models overestimate rainfall but with different patterns and amplitude, but both have underestimations in the South-Eastern area (Cevennes mountains) in winter. A variance decomposition shows that uncorrected simulations fairly well capture observed variances from inter-annual to high-frequency intra-seasonal time scales. After correction, distributions match with analyses by construction, but it is shown that spatial coherence, persistence properties of warm, cold and dry episodes also match to a certain extent. Another aim of the article is to describe the changes for future climate obtained using these simulations under Scenario A1B. Results are presented on the changes between current and mid-term future (2021–2050) averages and variability over France. Interestingly, even though the same global climate model is used at the boundaries, regional climate change responses from the two models significantly differ.     

在气候影响研究中引入随机天气发生器的方法和不确定性

通过采用不同的随机天气发生器生成一定气候背景下各种气候变率情景，许多学者在最近的研究中已经认识到气候变率对农作物生长发育影响的重要性。传统的气候影响评估方法直接以大气环流模式的模拟试验结果作为未来气候情景，这样不可能理解如上的重要性。本文着重评述将随机天气发生器应用于气候变化影响研究的一般方法框架，以及作者的具体个例研究方法。文中最后分析了目前该领域研究中还存在的一些不确定性。 在当前的气候变化影响研究中，有不同的方法用来研制一种称为WGEN的典型随机天气发生器的参数化方案及其随机试验方法。不同的研究者也有不同的参数调控方法。通常的思路是通过气候控制试验和2×CO2试验之间的气候变量平均值和方差的变化来扰动随机天气发生器的参数，以生成未来逐日气候变化情景。本文作者根据短期气候预测模式的输出产品建立了一套WGEN的参数化方案及其随机试验方法，并且在时间和空间两个尺度上检验和评估了此参数化方案下WGEN的模拟能力。另外，作者由未来降水的变化，调试随机天气发生器参数，生成了气候变率变化情景。这些参数调节可以产生各种不同类型和定性大小的气候变率变化，用于气候影响评估的敏感性分析。通过如上方法，作为一个个例，文中评估了未来气候变率变化     

Assessment of different metrics for physical climate feedbacks

We quantify the feedbacks from the physical climate system on the radiative forcing for idealized climate simulations using four different methods. The results differ between the methods and differences are largest for the cloud feedback. The spatial and temporal variability of each feedback is used to estimate the averaging scale necessary to satisfy the feedback concept of one constant global mean value. We find that the year-to-year variability, combined with the methodological differences, in estimates of the feedback strength from a single model is comparable to the model-to-model spread in feedback strength of the CMIP3 ensemble. The strongest spatial and temporal variability is in the short-wave component of the cloud feedback. In our simulations, where many sources of natural variability are neglected, long-term averages are necessary to get reliable feedback estimates. Considering the large natural variability and relatively small forcing present in the real world, as compared to the forcing imposed by doubling CO₂ concentrations in the simulations, implies that using observations to constrain feedbacks is a challenging task and requires reliable long-term measurements.     

Multivariate and multiscale dependence in the global climate system revealed through complex networks

A systematic characterization of multivariate dependence at multiple spatio-temporal scales is critical to understanding climate system dynamics and improving predictive ability from models and data. However, dependence structures in climate are complex due to nonlinear dynamical generating processes, long-range spatial and long-memory temporal relationships, as well as low-frequency variability. Here we utilize complex networks to explore dependence in climate data. Specifically, networks constructed from reanalysis-based atmospheric variables over oceans and partitioned with community detection methods demonstrate the potential to capture regional and global dependence structures within and among climate variables. Proximity-based dependence as well as long-range spatial relationships are examined along with their evolution over time, yielding new insights on ocean meteorology. The tools are implicitly validated by confirming conceptual understanding about aggregate correlations and teleconnections. Our results also suggest a close similarity of observed dependence patterns in relative humidity and horizontal wind speed over oceans. In addition, updraft velocity, which relates to convective activity over the oceans, exhibits short spatiotemporal decorrelation scales but long-range dependence over time. The multivariate and multi-scale dependence patterns broadly persist over multiple time windows. Our findings motivate further investigations of dependence structures among observations, reanalysis and model-simulated data to enhance process understanding, assess model reliability and improve regional climate predictions.     

Temperature and Precipitation Changes in China During the Holocene

We review here proxy records of temperature and precipitation in China during the Holocene, especially the last two millennia. The quality of proxy data, methodology of reconstruction, and uncertainties in reconstruction were emphasized in comparing different temperature and precipitation reconstruction and clarifying temporal and spatial patterns of temperature and precipitation during the Holocene. The Holocene climate was generally warm and wet. The warmest period occurred in 9.6-6.2 cal ka BP, whereas a period of maximum monsoon precipitation started at about 11.0 cal ka BP and lasted until about 8.0-5.0 cal ka BP. There were a series of millennial-scale cold or dry events superimposed on the general trend of climate changes. During past two millennia, a warming trend in the 20th century was clearly detected, but the warming magnitude was smaller than the maximum level of the Medieval Warm Period and the Middle Holocene. Cold conditions occurred over the whole of China during the Little Ice Age （AD 1400-AD 1900）, but the warming of the Medieval Warm Period （AD 900-AD 1300） was not distinct in China, especially west China. The spatial pattern of precipitation showed significant regional differences in China, especially east China. The modern warm period has lasted 20 years from 1987 to 2006. Bi-decadal oscillation in precipitation variability was apparent over China during the 20th century. Solar activity and volcanic eruptions both were major forcings governing the climate variability during the last millennium.     

Natural variability of summer rainfall over China in HadCM3

Summer rainfall over China has shown decadal variability in the past half century, which has resulted in major north–south shifts in rainfall with important implications for flooding and water resource management. This study has demonstrated how multi-century climate model simulations can be used to explore interdecadal natural variability in the climate system in order to address important questions around recent changes in Chinese summer rainfall, and whether or not anthropogenic climate change is playing a role. Using a 1,000-year simulation of HadCM3 with constant pre-industrial external forcing, the dominant modes of total and interdecadal natural variability in Chinese summer rainfall have been analysed. It has been shown that these modes are comparable in magnitude and in temporal and spatial characteristics to those observed in the latter part of the twentieth century. However, despite 1,000 years of model simulation it has not been possible to demonstrate that these modes are related to similar variations in the global circulation and surface temperature forcing occurring during the latter half of the twentieth century. This may be in part due to model biases. Consequently, recent changes in the spatial distribution of Chinese summer rainfall cannot be attributed solely to natural variability, nor has it been possible to eliminate the likelihood that anthropogenic climate change has been the driving factor. It is more likely that both play a role.     

集合样本数影响月平均500 hPa高度场预测技巧分析

气候系统模式输出结果是当前开展气候预测业务的重要参考依据之一,如何提高气候系统模式输出结果的可信度是改进气候业务预测能力的关键之一。利用1999—2010年NCEP CFSv2模式每日四次预测未来45天的回算数据,分析了集合样本数对模式预测能力的影响。分析结果表明,模式对月平均500 hPa位势高度的预测技巧在热带地区较高,而中高纬度地区较低;模式对500 hPa位势高度时间异常的预测能力优于空间异常。无论是空间异常还是时间异常,随着模式超前时间的增加,预测技巧均逐渐降低,但是在不同区域和不同月份,模式预测技巧随超前时间的变化存在差异。此外,模式预测技巧存在非常大的年际变率。增加集合样本数,对不同月份和不同起报时间预测技巧的稳定度和预测技巧值均有明显正效果,特别是对亚洲中纬度地区改善度较大。增加集合样本数也可以在一定程度上降低模式预测技巧年际变率。集合样本数增加对于500 hPa位势高度空间异常的改进优于时间异常。      

Comparing the spatial structure of variability in two datasets against each other on the basis of EOF-modes

In analysis of climate variability or change it is often of interest how the spatial structure in modes of variability in two datasets differ from each other, e.g. between past and future climate or between models and observations. Often such analysis is based on Empirical Orthogonal Function (EOF) analysis or other simple indices of large-scale spatial structures. The present analysis lays out a concept on how two datasets of multivariate climate variability can be compared against each other on basis of EOF analysis and how the differences in the multivariate spatial structure between the two datasets can be quantified in terms of explained variance in the leading spatial patterns. It is also illustrated how the patterns of largest differences between the two datasets can be defined and interpreted. We illustrate this method on the basis of several well-defined artificial examples and by comparing our approach with examples of climate change studies from the literature. These literature examples include analysis of changes in the modes of variability under climate change for the sea level pressure (SLP) of the North Atlantic and Europe, the SLP of the Southern Hemisphere, the surface temperature of the Northern Hemisphere, the sea surface temperature of the North Pacific and for precipitation in the tropical Indo-Pacific.     

Climatic drivers of tree growth at tree line in Southwest Yukon change over time and vary between landscapes

Growth of trees at their altitudinal and latitudinal range limits is expected to increase as climate warms, but trees often exhibit unexplained spatial and temporal variation in climate-growth responses, particularly in alpine regions. Until this variability is explained, predictions of future tree growth are unlikely to be accurate. We sampled Picea glauca (white spruce) growing at forest and tree line on north and south aspects in two mountain ranges of southwest Yukon to determine how and why ring-width patterns vary between topographic settings, and over time. We used multivariate statistical analysis to characterize variation in ring-width patterns between topographic factors and time periods, and calculated correlations between ring-width indices and climate variables to explain this variation. Ring-width patterns varied more between mountain ranges than elevations or aspects, particularly in recent decades when ring-widths increased in one mountain range but not the other. Growth responses to summer temperature were notably weaker during warmer time periods, but growth was not positively correlated to summer precipitation, suggesting trees may not be suffering from temperature-induced drought stress. Rather, ring-width indices began responding positively to spring snow depth after 1976. We conclude that tree growth is unlikely to increase in synchrony with rising air temperatures across subarctic tree lines in southwest Yukon. Instead, they may decline in areas that are prone to thin snowpacks or rapid spring runoff due to the negative influence warming springs will have on snow depth and, consequently, early growing season soil moisture.     

New precipitation and temperature grids for northern Patagonia: Advances in relation to global climate grids

Climate data of mean monthly temperature and total monthly precipitation compiled from different sources in northern Patagonia were interpolated to 20-km resolution grids over the period 1997–2010. This northern Patagonian climate grid (NPCG) improves upon previous gridded products in terms of its spatial resolution and number of contributing stations, since it incorporates 218 and 114 precipitation and temperature records, respectively. A geostatistical method using surface elevation from a Digital Elevation Model (DEM) as the ancillary variable was used to interpolate station data into even spaced points. The maps provided by NPCG are consistent with the broad spatial and temporal patterns of the northern Patagonian climate, showing a comprehensive representation of the latitudinal and altitudinal gradients in temperature and precipitation, as well as their related patterns of seasonality and continentality. We compared the performance of NPCG and various other datasets available to the climate community for northern Patagonia. The grids used for the comparison included those of the Global Precipitation Climatology Project, ERAInterim, Climate Research Unit (University of East Anglia), and University of Delaware. Based on three statistics that quantitatively assess the spatial coherence of gridded data against available observations (bias, MAE, and RMSE), NPCG outperforms other global grids. NPCG represents a useful tool for understanding climate variability in northern Patagonia and a valuable input for regional models of hydrological and ecological processes. Its resolution is optimal for validating data from the general circulation models and working with raster data derived from remote sensing, such as vegetation indices.