ALBERT

Description: A scalable second-order analytical orbit propagator programme based on modern and classical perturbation methods is being developed. As a first step in the validation and verification of part of our orbit propagator programme, we only consider the perturbation produced by zonal harmonic coefficients in the Earth’s gravity potential, so that it is possible to analyze the behaviour of the mathematical expressions involved in Delaunay normalization and the Krylov-Bogoliubov-Mitropolsky method in depth and determine their limits.

Description: We present a new approach in astrodynamics and celestial mechanics fields, called hybrid perturbation theory. A hybrid perturbation theory combines an integrating technique, general perturbation theory or special perturbation theory or semianalytical method, with a forecasting technique, statistical time series model or computational intelligence method. This combination permits an increase in the accuracy of the integrating technique, through the modeling of higher-order terms and other external forces not considered in the integrating technique. In this paper, neural networks have been used as time series forecasters in order to help two economic general perturbation theories describe the motion of an orbiter only perturbed by the Earth’s oblateness.

Description: Deprit’s method has been revisited in order to take advantage of certain arbitrariness arising when the inverse of the Lie operator is applied to obtain the generating function of the Lie transform. This arbitrariness is intrinsic to all perturbation techniques and can be used to demonstrate the equivalence among different perturbation methods, to remove terms from the generating function of the Lie transform, or to eliminate several angles simultaneously in the case of having a degenerate Hamiltonian.

Description: This is the second in a pair of papers in which the performance of Statistical Downscaling Methods (SDMs) is critically re-assessed with respect to their robust applicability in climate change studies. Whereas Part I focused on temperatures (Gutiérrez et al. 2013), the present manuscript deals with precipitation and considers an ensemble of twelve SDMs from the analog, weather typing and regression families. First, the performance of the methods is cross-validated considering reanalysis predictors, screening different geographical domains and predictor sets. Standard accuracy and distributional similarity scores, and a test for extrapolation capability are considered. The results are highly dependent on the predictor sets, with optimum configurations including information from middle tropospheric humidity. Second, a reduced ensemble of good performing SDMs is applied to four GCMs to properly assess the uncertainty of downscaled future climate projections. The results are compared with an ensemble of Regional Climate Models (RCMs) produced in the ENSEMBLES project. Generally, the mean signal is similar with both methodologies (with the exception of Summer, which is drier for the RCMs) but the uncertainty (spread) is larger for the SDM ensemble. Finally, the spread contribution of the GCM and SDM-derived components is assessed using a simple analysis of variance previously applied to the RCMs, obtaining larger interaction terms. Results show that the main contributor to the spread is the choice of the GCM, although the SDM dominates the uncertainty in some cases during Autumn and Summer due to the diverging projections from different families.

Description: The increasing demand for high-resolution climate information has attracted growing attention to statistical downscaling (SDS) methods, due in part to their relative advantages and merits as compared to dynamical approaches (based on regional climate model simulations), such as their much lower computational cost and their fitness for purpose for many local-scale applications. As a result, a plethora of SDS methods is nowadays available to climate scientists, which has motivated recent efforts for their comprehensive evaluation, like the VALUE initiative (http://www.value-cost.eu, last access: 29 March 2020). The systematic intercomparison of a large number of SDS techniques undertaken in VALUE, many of them independently developed by different authors and modeling centers in a variety of languages/environments, has shown a compelling need for new tools allowing for their application within an integrated framework. In this regard, downscaleR is an R package for statistical downscaling of climate information which covers the most popular approaches (model output statistics – including the so-called “bias correction” methods – and perfect prognosis) and state-of-the-art techniques. It has been conceived to work primarily with daily data and can be used in the framework of both seasonal forecasting and climate change studies. Its full integration within the climate4R framework (Iturbide et al., 2019) makes possible the development of end-to-end downscaling applications, from data retrieval to model building, validation, and prediction, bringing to climate scientists and practitioners a unique comprehensive framework for SDS model development. In this article the main features of downscaleR are showcased through the replication of some of the results obtained in VALUE, placing an emphasis on the most technically complex stages of perfect-prognosis model calibration (predictor screening, cross-validation, and model selection) that are accomplished through simple commands allowing for extremely flexible model tuning, tailored to the needs of users requiring an easy interface for different levels of experimental complexity. As part of the open-source climate4R framework, downscaleR is freely available and the necessary data and R scripts to fully replicate the experiments included in this paper are also provided as a companion notebook.

Description: Clouds play an important role in the climate system: (1) cooling Earth by reflecting incoming sunlight to space and (2) warming Earth by reducing thermal energy loss to space. Cloud radiative effects are especially important in polar regions and have the potential to significantly alter the impact of sea ice decline on the surface radiation budget. Using CERES (Clouds and the Earth's Radiant Energy System) data and 32 CMIP5 (Coupled Model Intercomparison Project) climate models, we quantify the influence of polar clouds on the radiative impact of polar sea ice variability. Our results show that the cloud short-wave cooling effect strongly influences the impact of sea ice variability on the surface radiation budget and does so in a counter-intuitive manner over the polar seas: years with less sea ice and a larger net surface radiative flux show a more negative cloud radiative effect. Our results indicate that 66±2% of this change in the net cloud radiative effect is due to the reduction in surface albedo and that the remaining 34±1 % is due to an increase in cloud cover and optical thickness. The overall cloud radiative damping effect is 56±2 % over the Antarctic and 47±3 % over the Arctic. Thus, present-day cloud properties significantly reduce the net radiative impact of sea ice loss on the Arctic and Antarctic surface radiation budgets. As a result, climate models must accurately represent present-day polar cloud properties in order to capture the surface radiation budget impact of polar sea ice loss and thus the surface albedo feedback.

Description: The performance of statistical downscaling (SD) techniques is critically reassessed with respect to their robust applicability in climate change studies. To this end, in addition to standard accuracy measures and distributional similarity scores, the authors estimate the robustness of the methods under warming climate conditions working with anomalous warm historical periods. This validation framework is applied to intercompare the performances of 12 different SD methods (from the analog, weather typing, and regression families) for downscaling minimum and maximum temperatures in Spain. First, a calibration of these methods is performed in terms of both geographical domains and predictor sets; the results are highly dependent on the latter, with optimum predictor sets including near-surface temperature data (in particular 2-m temperature), which appropriately discriminate cold episodes related to temperature inversion in the lower troposphere. Although regression methods perform best in terms of correlation, analog and weather generator approaches are more appropriate for reproducing the observed distributions, especially in case of wintertime minimum temperature. However, the latter two families significantly underestimate the temperature anomalies of the warm periods considered in this work. This underestimation is found to be critical when considering the warming signal in the late twenty-first century as given by a global climate model [the ECHAM5–Max Planck Institute (MPI) model]. In this case, the different downscaling methods provide warming values with differences in the range of 1°C, in agreement with the robustness significance values. Therefore, the proposed test is a promising technique for detecting lack of robustness in statistical downscaling methods applied in climate change studies.

Description: This work shows that local-scale climate projections obtained by means of statistical downscaling are sensitive to the choice of reanalysis used for calibration. To this aim, a generalized linear model (GLM) approach is applied to downscale daily precipitation in the Philippines. First, the GLMs are trained and tested separately with two distinct reanalyses (ERA-Interim and JRA-25) using a cross-validation scheme over the period 1981–2000. When the observed and downscaled time series are compared, the attained performance is found to be sensitive to the reanalysis considered if climate change signal–bearing variables (temperature and/or specific humidity) are included in the predictor field. Moreover, performance differences are shown to be in correspondence with the disagreement found between the raw predictors from the two reanalyses. Second, the regression coefficients calibrated either with ERA-Interim or JRA-25 are subsequently applied to the output of a global climate model (MPI-ECHAM5) in order to assess the sensitivity of local-scale climate change projections (up to 2100) to reanalysis choice. In this case, the differences detected in present climate conditions are considerably amplified, leading to “delta-change” estimates differing by up to 35% (on average for the entire country) depending on the reanalysis used for calibration. Therefore, reanalysis choice is an important contributor to the uncertainty of local-scale climate change projections and, consequently, should be treated with as much care as other better-known sources of uncertainty (e.g., the choice of the GCM and/or downscaling method). Implications of the results for the entire tropics, as well as for the model output statistics downscaling approach are also briefly discussed.

Description: Limited area models (LAMs) are widely used tools to downscale the wind speed forecasts issued by general circulation models. However, only a few studies have systematically analyzed the value added by the LAMs to the coarser-resolution-model wind. The goal of the present work is to investigate how added value depends on the resolution of the driving global model. With this aim, the Weather Research and Forecasting (WRF) Model was used to downscale three different global datasets (GFS, ERA-Interim, and NCEP–NCAR) to a 9-km-resolution grid for a 1-yr period. Model results were compared with a large set of surface observations, including land station and offshore buoy data. Substantial biases were found at this resolution over mountainous terrain, and a slight modification to the subgrid orographic drag parameterization was introduced to alleviate the problem. It was found that, at this resolution, WRF is able to produce significant added value with respect to the NCEP–NCAR reanalysis and ERA-Interim but only a small amount of added value with respect to GFS forecasts. Results suggest that, as model resolution increases, traditional skill scores tend to saturate. Thus, adding value to high-resolution global models becomes significantly more difficult.

Description: The increasing demand for high-resolution climate information has attracted a growing attention for statistical downscaling methods (SD), due in part to their relative advantages and merits as compared to dynamical approaches (based on regional climate model simulations), such as their much lower computational cost and their fitness-for-purpose for many local-scale applications. As a result, a plethora of SD methods is nowadays available for climate scientists, which has motivated recent efforts for their comprehensive evaluation, like the VALUE Project (http://www.value-cost.eu). The systematic intercomparison of a large number of SD techniques undertaken in VALUE, many of them independently developed by different authors and modeling centers in a variety of languages/environments, has shown a compelling need for new tools allowing for their application within an integrated framework. With this regard, downscaleR is an R package for statistical downscaling of climate information which covers the most popular approaches (Model Output Statistics – including the so called 'bias correction' methods – and Perfect Prognosis) and state-of-the-art techniques. It has been conceived to work primarily with daily data and can be used in the framework of both seasonal forecasting and climate change studies. Its full integration within the climate4R framework (Iturbide et al. 2019) makes possible the development of end-to-end downscaling applications, from data retrieval to model building, validation and prediction, bringing to climate scientists and practitioners a unique comprehensive framework for SD model development. In this article the main features of downscaleR are showcased through the replication of some of the results obtained in the VALUE Project, making an emphasis in the most technically complex stages of perfect-prog model calibration (predictor screening, cross-validation and model selection) that are accomplished through simple commands allowing for extremely flexible model tuning, tailored to the needs of users requiring an easy interface for different levels of experimental complexity. As part of the open-source climate4R framework, downscaleR is freely available and the necessary data and R scripts to fully replicate the experiments included in this paper are also provided as a companion notebook.