scholarly journals Non-linear statistical downscaling of present and LGM precipitation and temperatures over Europe

2007 ◽  
Vol 3 (4) ◽  
pp. 899-933 ◽  
Author(s):  
M. Vrac ◽  
D. Paillard ◽  
P. Naveau

Abstract. The needs of small-scale climate information have become prevalent to study the impacts of future climate change as well as for paleoclimate researches where the reconstructions from proxies are obviously local. In this study we develop a non-linear statistical downscaling method to generate local temperatures and precipitation values from large-scale variables (e.g. Global Circulation Model – GCM – outputs), through Generalized Additive Models (GAMs) calibrated on the present Western Europe climate. First, various monthly GAMs (i.e. one model for each month) are tested for preliminary analysis. Then, annual GAMs (i.e. one model for the 12 months altogether) are developed and tailored for two sets of predictors (geographical and physical) to downscale local temperatures and precipitation. As an evaluation of our approach under large-scale conditions different from present Western Europe, projections are realized (1) for present North America and Northern Europe and compared to local observations (spatial test); and (2) for the Last Glacial Maximum (LGM) period, and compared to local reconstructions and GCMs outputs (temporal test). In general, both spatial and temporal evaluations indicate that the GAMs are flexible and efficient tools to capture and downscale non-linearities between large- and local-scale variables. More precisely, the results emphasize that, while physical predictors alone are not capable of downscaling realistic values when applied to climate strongly different from the one used for calibration, the inclusion of geographical-type variables – such as altitude, advective continentality and W-slope – into GAM predictors brings robustness and improvement to the method and its local projections.

2007 ◽  
Vol 3 (4) ◽  
pp. 669-682 ◽  
Author(s):  
M. Vrac ◽  
P. Marbaix ◽  
D. Paillard ◽  
P. Naveau

Abstract. Local-scale climate information is increasingly needed for the study of past, present and future climate changes. In this study we develop a non-linear statistical downscaling method to generate local temperatures and precipitation values from large-scale variables of a Earth System Model of Intermediate Complexity (here CLIMBER). Our statistical downscaling scheme is based on the concept of Generalized Additive Models (GAMs), capturing non-linearities via non-parametric techniques. Our GAMs are calibrated on the present Western Europe climate. For this region, annual GAMs (i.e. models based on 12 monthly values per location) are fitted by combining two types of large-scale explanatory variables: geographical (e.g. topographical information) and physical (i.e. entirely simulated by the CLIMBER model). To evaluate the adequacy of the non-linear transfer functions fitted on the present Western European climate, they are applied to different spatial and temporal large-scale conditions. Local projections for present North America and Northern Europe climates are obtained and compared to local observations. This partially addresses the issue of spatial robustness of our transfer functions by answering the question "does our statistical model remain valid when applied to large-scale climate conditions from a region different from the one used for calibration?". To asses their temporal performances, local projections for the Last Glacial Maximum period are derived and compared to local reconstructions and General Circulation Model outputs. Our downscaling methodology performs adequately for the Western Europe climate. Concerning the spatial and temporal evaluations, it does not behave as well for Northern America and Northern Europe climates because the calibration domain may be too different from the targeted regions. The physical explanatory variables alone are not capable of downscaling realistic values. However, the inclusion of geographical-type variables – such as altitude, advective continentality and moutains effect on wind (W–slope) – as GAM explanatory variables clearly improves our local projections.


2010 ◽  
Vol 4 (4) ◽  
pp. 2233-2275 ◽  
Author(s):  
G. Levavasseur ◽  
M. Vrac ◽  
D. M. Roche ◽  
D. Paillard ◽  
A. Martin ◽  
...  

Abstract. We quantify the agreement between permafrost distributions from PMIP2 (Paleoclimate Modeling Intercomparison Project) climate models and permafrost data. We evaluate the ability of several climate models to represent permafrost and assess the inter-variability between them. Studying an heterogeneous variable such as permafrost implies to conduct analysis at a smaller spatial scale compared with climate models resolution. Our approach consists in applying statistical downscaling methods (SDMs) on large- or regional-scale atmospheric variables provided by climate models, leading to local permafrost modelling. Among the SDMs, we first choose a transfer function approach based on Generalized Additive Models (GAMs) to produce high-resolution climatology of surface air temperature (SAT). Then, we define permafrost distribution over Eurasia by SAT conditions. In a first validation step on present climate (CTRL period), GAM shows some limitations with non-systemic improvements in comparison with the large-scale fields. So, we develop an alternative method of statistical downscaling based on a stochastic generator approach through a Multinomial Logistic Regression (MLR), which directly models the probabilities of local permafrost indices. The obtained permafrost distributions appear in a better agreement with data. In both cases, the provided local information reduces the inter-variability between climate models. Nevertheless, this also proves that a simple relationship between permafrost and the SAT only is not always sufficient to represent local permafrost. Finally, we apply each method on a very different climate, the Last Glacial Maximum (LGM) time period, in order to quantify the ability of climate models to represent LGM permafrost. Our SDMs do not significantly improve permafrost distribution and do not reduce the inter-variability between climate models, at this period. We show that LGM permafrost distribution from climate models strongly depends on large-scale SAT. The differences with LGM data, larger than in the CTRL period, reduce the contribution of downscaling and depend on several factors deserving further studies.


2004 ◽  
Vol 35 (3) ◽  
pp. 261-278 ◽  
Author(s):  
Maj-Lena Linderson ◽  
Christine Achberger ◽  
Deliang Chen

Statistical downscaling models for precipitation in Scania, southern Sweden, have been developed and applied to calculate the changes in the future Scanian precipitation climate due to projected changes in the atmospheric composition. The models are based on multiple linear regression, linking large-scale predictors at monthly time resolution to regional statistics of daily precipitation on a monthly basis. To account for spatial precipitation variability within the area, the precipitation statistics were derived for different regions in Scania. The final downscaling models, developed for different regions and seasons, use atmospheric circulation, large-scale humidity and precipitation as predictors. Among the precipitation statistics examined, only the models for estimating the mean precipitation and the frequency of wet days were skilful. Based on the Canadian Global Circulation Model 1 (CGCM1), a future scenario of these two statistics was created. The downscaled scenario shows a significant increase of the annual mean precipitation by about 10% and a slight decrease in the frequency of wet days, indicating an increase in the precipitation amounts as well as in the precipitation intensity. The main increase of precipitation amounts and intensity occur during winter, while the summer precipitation amounts decrease slightly. The seasonal changes found in precipitation are likely attributed to changes in the westerly flow of the atmospheric circulation.


2011 ◽  
Vol 7 (4) ◽  
pp. 1225-1246 ◽  
Author(s):  
G. Levavasseur ◽  
M. Vrac ◽  
D. M. Roche ◽  
D. Paillard ◽  
A. Martin ◽  
...  

Abstract. We quantify the agreement between permafrost distributions from PMIP2 (Paleoclimate Modeling Intercomparison Project) climate models and permafrost data. We evaluate the ability of several climate models to represent permafrost and assess the variability between their results. Studying a heterogeneous variable such as permafrost implies conducting analysis at a smaller spatial scale compared with climate models resolution. Our approach consists of applying statistical downscaling methods (SDMs) on large- or regional-scale atmospheric variables provided by climate models, leading to local-scale permafrost modelling. Among the SDMs, we first choose a transfer function approach based on Generalized Additive Models (GAMs) to produce high-resolution climatology of air temperature at the surface. Then we define permafrost distribution over Eurasia by air temperature conditions. In a first validation step on present climate (CTRL period), this method shows some limitations with non-systematic improvements in comparison with the large-scale fields. So, we develop an alternative method of statistical downscaling based on a Multinomial Logistic GAM (ML-GAM), which directly predicts the occurrence probabilities of local-scale permafrost. The obtained permafrost distributions appear in a better agreement with CTRL data. In average for the nine PMIP2 models, we measure a global agreement with CTRL permafrost data that is better when using ML-GAM than when applying the GAM method with air temperature conditions. In both cases, the provided local information reduces the variability between climate models results. This also confirms that a simple relationship between permafrost and the air temperature only is not always sufficient to represent local-scale permafrost. Finally, we apply each method on a very different climate, the Last Glacial Maximum (LGM) time period, in order to quantify the ability of climate models to represent LGM permafrost. The prediction of the SDMs (GAM and ML-GAM) is not significantly in better agreement with LGM permafrost data than large-scale fields. At the LGM, both methods do not reduce the variability between climate models results. We show that LGM permafrost distribution from climate models strongly depends on large-scale air temperature at the surface. LGM simulations from climate models lead to larger differences with LGM data than in the CTRL period. These differences reduce the contribution of downscaling.


Mathematics ◽  
2021 ◽  
Vol 9 (4) ◽  
pp. 299
Author(s):  
Jaime Pinilla ◽  
Miguel Negrín

The interrupted time series analysis is a quasi-experimental design used to evaluate the effectiveness of an intervention. Segmented linear regression models have been the most used models to carry out this analysis. However, they assume a linear trend that may not be appropriate in many situations. In this paper, we show how generalized additive models (GAMs), a non-parametric regression-based method, can be useful to accommodate nonlinear trends. An analysis with simulated data is carried out to assess the performance of both models. Data were simulated from linear and non-linear (quadratic and cubic) functions. The results of this analysis show how GAMs improve on segmented linear regression models when the trend is non-linear, but they also show a good performance when the trend is linear. A real-life application where the impact of the 2012 Spanish cost-sharing reforms on pharmaceutical prescription is also analyzed. Seasonality and an indicator variable for the stockpiling effect are included as explanatory variables. The segmented linear regression model shows good fit of the data. However, the GAM concludes that the hypothesis of linear trend is rejected. The estimated level shift is similar for both models but the cumulative absolute effect on the number of prescriptions is lower in GAM.


2019 ◽  
Vol 4 (3) ◽  
pp. 194-207
Author(s):  
PIET GELEYNS

The Hoge Kempen rural industrial transition landscape: a layered landscape of Outstanding Universal Value? Up until the beginning of the 20th century, the eastern part of the Belgian province of Limburg was a sparsely populated and not very productive part of the country. The dominating heathland was maintained with sheep, which were an essential part of a small-scale extensive farming system. This all changed when coal was discovered in 1901. Seven large coalmines were established in a few decades, each one employing thousands of coal-miners. This also meant that entire new garden cities were built, to house the coal-miners and their families. The confrontation between the small-scale traditional land-use and the new large-scale industrial developments defines the landscape up to today. The scale and the force of the turnover are considered unprecedented for Western Europe, which is why it is being presented by Belgium for inclusion in the World Heritage List.


2007 ◽  
Vol 4 (5) ◽  
pp. 3413-3440 ◽  
Author(s):  
E. P. Maurer ◽  
H. G. Hidalgo

Abstract. Downscaling of climate model data is essential to most impact analysis. We compare two methods of statistical downscaling to produce continuous, gridded time series of precipitation and surface air temperature at a 1/8-degree (approximately 140 km² per grid cell) resolution over the western U.S. We use NCEP/NCAR Reanalysis data from 1950–1999 as a surrogate General Circulation Model (GCM). The two methods included are constructed analogues (CA) and a bias correction and spatial downscaling (BCSD), both of which have been shown to be skillful in different settings, and BCSD has been used extensively in hydrologic impact analysis. Both methods use the coarse scale Reanalysis fields of precipitation and temperature as predictors of the corresponding fine scale fields. CA downscales daily large-scale data directly and BCSD downscales monthly data, with a random resampling technique to generate daily values. The methods produce comparable skill in producing downscaled, gridded fields of precipitation and temperatures at a monthly and seasonal level. For daily precipitation, both methods exhibit some skill in reproducing both observed wet and dry extremes and the difference between the methods is not significant, reflecting the general low skill in daily precipitation variability in the reanalysis data. For low temperature extremes, the CA method produces greater downscaling skill than BCSD for fall and winter seasons. For high temperature extremes, CA demonstrates higher skill than BCSD in summer. We find that the choice of most appropriate downscaling technique depends on the variables, seasons, and regions of interest, on the availability of daily data, and whether the day to day correspondence of weather from the GCM needs to be reproduced for some applications. The ability to produce skillful downscaled daily data depends primarily on the ability of the climate model to show daily skill.


2018 ◽  
Vol 15 (144) ◽  
pp. 20180327 ◽  
Author(s):  
Mitchell B. Lyons ◽  
Charlotte H. Mills ◽  
Christopher E. Gordon ◽  
Mike Letnic

Vegetation cover is fundamental in the formation and maintenance of geomorphological features in dune systems. In arid Australia, increased woody shrub cover has been linked to removal of the apex predator (Dingoes, Canis dingo ) via subsequent trophic cascades. We ask whether this increase in shrubs can be linked to altered physical characteristics of the dunes. We used drone-based remote sensing to measure shrub density and construct three-dimensional models of dune morphology. Dunes had significantly different physical characteristics either side of the ‘dingo-proof fence’, inside which dingoes are systematically eradicated and shrub density is higher over vast spatial extents. Generalized additive models revealed that dunes with increased shrub density were higher, differently shaped and more variable in height profile. We propose that low shrub density induces aeolian and sedimentary processes that result in greater surface erosion and sediment transport, whereas high shrub density promotes dune stability. We speculate that increased vegetation cover acts to push dunes towards an alternate stable state, where climatic variation no longer has a significant effect on their morphodynamic state within the bi-stable state model. Our study provides evidence that anthropogenically induced trophic cascades can indirectly lead to large-scale changes in landscape geomorphology.


2004 ◽  
Vol 132 (6) ◽  
pp. 1065-1071 ◽  
Author(s):  
S. BROOKER ◽  
N. B. KABATEREINE ◽  
E. M. TUKAHEBWA ◽  
F. KAZIBWE

The spatial epidemiology of intestinal nematodes in Uganda was investigated using generalized additive models and geostatistical methods. The prevalence of Ascaris lumbricoides and Trichuris trichiura was unevenly distributed in the country with prevalence greatest in southwest Uganda whereas hookworm was more homogeneously distributed. A. lumbricoides and T. Trichiura prevalence were nonlinearly related to satellite sensor-based estimates of land surface temperature; hookworm was nonlinearly associated with rainfall. Semivariogram analysis indicated that T. trichiura prevalence exhibited no spatial structure and that A. lumbricoides exhibited some spatial dependency at small spatial distances, once large-scale, mainly environmental, trends had been removed. In contrast, there was much more spatial structure in hookworm prevalence although the underlying factors are at present unclear. The implications of the results are discussed in relation to parasite spatial epidemiology and the prediction of infection distributions.


Sign in / Sign up

Export Citation Format

Share Document