Redes neuronales de convolución profundas para la regionalización estadística de proyecciones de cambio climático

RESUMEN Las proyecciones climáticas a escala local y/o regional son muy demandadas por diversos sectores socioeconómicos para elaborar sus planes de adaptación y mitigación al cambio climático. Sin embargo, los modelos climáticos globales actuales presentan una resolución espacial muy baja, lo que dificulta enormemente la elaboración de este tipo de estudios. Una manera de aumentar esta resolución es establecer relaciones estadísticas entre la variable local de interés (por ejemplo la temperatura y/o precipitación en una localidad dada) y un conjunto de variables de larga escala (por ejemplo, geopotencial y/o vientos en distintos niveles verticales) dadas por los modelos climáticos. En particular, en esta Tesis se explora la idoneidad de las redes neuronales de convolución (CNN) como método de downscaling estadístico para generar proyecciones de cambio climático a alta resolución sobre Europa. Para ello se evalúa primero la capacidad de estos modelos para reproducir la variabilidad local de precipitación y de temperatura en un período histórico reciente, comparándolas contra otros métodos estadísticos de referencia. A continuación, se analiza la idoneidad de estos modelos para regionalizar las proyecciones climáticas en el futuro (hasta el año 2100). Además, se desarrollan diversos estudios de interpretabilidad sobre redes neuronales para ganar confianza y conocimiento sobre el uso de este tipo de técnicas para aplicaciones climáticas, puesto que a menudo son rechazadas por ser consideradas “cajas negras”.ABSTRACT Regional climate projections are very demanded by different socioeconomics sectors to elaborate their adaptation and mitigation plans to climate change. Nevertheless, the state-of-the-art Global Glimate Models (GCMs) present very coarse spatial resolutions what limits their use in most of practical applications and impact studies. One way to increase this limited spatial resolution is to establish empirical/statistical functions which link the local variable of interest (e.g. temperature and/or precipitation at a given site) with a set of large-scale atmospheric variables (e.g. geopotential and/or winds at different vertical levels), which are typically well-reproduced by GCMs. In this context, this Thesis explores the suitability of deep learning, and in particular modern Convolutional Neural Networks (CNNs), as statistical downscaling techniques to produce regional climate change projections over Europe. To achieve this ambitious goal, the capacity of CNNs to reproduce the local variability of precipitation and temperature fields in present climate conditions is first assessed by comparing their performance with that from a set of traditional, benchmark statistical methods. Subsequently, their suitability to produce plausible future (up to 2100) high-resolution scenarios is put to the test by comparing their projected signals of change with those given by a set of state-of-the-art GCMs from CMIP5 and Regional Climate Models (RCMs) from the flagship EURO-CORDEX initiative. Also, a variety of interpretability techniques are also carried out to gain confidence and knowledge on the use of CNNs for climate applications, which have typically discarded until now for being considered as "black-boxes"

Deep Ensembles to Improve Uncertainty Quantification of Statistical Downscaling Models under Climate Change Conditions

Recently, deep learning has emerged as a promising tool for statistical downscaling, the set of methods for generating high-resolution climate fields from coarse low-resolution variables. Nevertheless, their ability to generalize to climate change conditions remains questionable, mainly due to the stationarity assumption. We propose deep ensembles as a simple method to improve the uncertainty quantification of statistical downscaling models. By better capturing uncertainty, statistical downscaling models allow for superior planning against extreme weather events, a source of various negative social and economic impacts. Since no observational future data exists, we rely on a pseudo reality experiment to assess the suitability of deep ensembles for quantifying the uncertainty of climate change projections. Deep ensembles allow for a better risk assessment, highly demanded by sectoral applications to tackle climate change.Comment: Accepted at the ICLR 2023 Tackling Climate Change with Machine Learning Worksho

On the use of Deep Generative Models for Perfect Prognosis Climate Downscaling

Deep Learning has recently emerged as a perfect prognosis downscaling technique to compute high-resolution fields from large-scale coarse atmospheric data. Despite their promising results to reproduce the observed local variability, they are based on the estimation of independent distributions at each location, which leads to deficient spatial structures, especially when downscaling precipitation. This study proposes the use of generative models to improve the spatial consistency of the high-resolution fields, very demanded by some sectoral applications (e.g., hydrology) to tackle climate change.Comment: Accepted at the NeurIPS 2021 Tackling Climate Change with Machine Learning Worksho

Configuration and intercomparison of deep learning neural models for statistical downscaling

Deep learning techniques (in particular convolutional neural networks, CNNs) have recently emerged as a promising approach for statistical downscaling due to their ability to learn spatial features from huge spatiotemporal datasets. However, existing studies are based on complex models, applied to particular case studies and using simple validation frameworks, which makes a proper assessment of the (possible) added value offered by these techniques difficult. As a result, these models are usually seen as black boxes, generating distrust among the climate community, particularly in climate change applications. In this paper we undertake a comprehensive assessment of deep learning techniques for continental-scale statistical downscaling, building on the VALUE validation framework. In particular, different CNN models of increasing complexity are applied to downscale temperature and precipitation over Europe, comparing them with a few standard benchmark methods from VALUE (linear and generalized linear models) which have been traditionally used for this purpose. Besides analyzing the adequacy of different components and topologies, we also focus on their extrapolation capability, a critical point for their potential application in climate change studies. To do this, we use a warm test period as a surrogate for possible future climate conditions. Our results show that, while the added value of CNNs is mostly limited to the reproduction of extremes for temperature, these techniques do outperform the classic ones in the case of precipitation for most aspects considered. This overall good performance, together with the fact that they can be suitably applied to large regions (e.g., continents) without worrying about the spatial features being considered as predictors, can foster the use of statistical approaches in international initiatives such as Coordinated Regional Climate Downscaling Experiment (CORDEX).The authors acknowledge the funding provided by the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER). They also acknowledge the E-OBS dataset from the EU-FP6 project UERRA (http://www.uerra.eu, last access: 23 April 2020) and the Copernicus Climate Change Service, and the data providers in the ECA&D project (https://www.ecad.eu, last access: 23 April 2020)

On the suitability of deep convolutional neural networks for continental-wide downscaling of climate change projections

In a recent paper, Baño-Medina et al. (Confguration and Intercomparison of deep learning neural models for statistical downscaling. preprint, 2019) assessed the suitability of deep convolutional neural networks (CNNs) for downscaling of temperature and precipitation over Europe using large-scale 'perfect' reanalysis predictors. They compared the results provided by CNNs with those obtained from a set of standard methods which have been traditionally used for downscaling purposes (linear and generalized linear models), concluding that CNNs are well suited for continental-wide applications. That analysis is extended here by assessing the suitability of CNNs for downscaling future climate change projections using Global Climate Model (GCM) outputs as predictors. This is particularly relevant for this type of 'black-box' models, whose results cannot be easily explained based on physical reasons and could potentially lead to implausible downscaled projections due to uncontrolled extrapolation artifacts. Based on this premise, we analyze in this work the two key assumptions that are made in perfect prognosis downscaling: (1) the predictors chosen to build the statistical model should be well reproduced by GCMs and (2) the statistical model should be able to reliably extrapolate out of sample (climate change) conditions. As a first step to test the suitability of these models, the latter assumption is assessed here by analyzing how the CNNs afect the raw GCM climate change signal (defned as the diference, or delta, between future and historical climate). Our results show that, as compared to well-established generalized linear models (GLMs), CNNs yield smaller departures from the raw GCM outputs for the end of century, resulting in more plausible downscaling results for climate change applications. Moreover, as a consequence of the automatic treatment of spatial features, CNNs are also found to provide more spatially homogeneous downscaled patterns than GLMs.The authors acknowledge partial support from the ATLAS project, funded by the Spanish Research Program (PID2019-111481RB-I00). Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature

Downscaling multi-model climate projection ensembles with deep learning (DeepESD): contribution to CORDEX EUR-44

Deep learning (DL) has recently emerged as an innovative tool to downscale climate variables from large-scale atmospheric fields under the perfect-prognosis (PP) approach. Different convolutional neural networks (CNNs) have been applied under present-day conditions with promising results, but little is known about their suitability for extrapolating future climate change conditions. Here, we analyze this problem from a multi-model perspective, developing and evaluating an ensemble of CNN-based downscaled projections (hereafter DeepESD) for temperature and precipitation over the European EUR-44i (0.5º) domain, based on eight global circulation models (GCMs) from the Coupled Model Intercomparison Project Phase 5 (CMIP5). To our knowledge, this is the first time that CNNs have been used to produce downscaled multi-model ensembles based on the perfect-prognosis approach, allowing us to quantify inter-model uncertainty in climate change signals. The results are compared with those corresponding to an EUR-44 ensemble of regional climate models (RCMs) showing that DeepESD reduces distributional biases in the historical period. Moreover, the resulting climate change signals are broadly comparable to those obtained with the RCMs, with similar spatial structures. As for the uncertainty of the climate change signal (measured on the basis of inter-model spread), DeepESD preserves the uncertainty for temperature and results in a reduced uncertainty for precipitation. To facilitate further studies of this downscaling approach, we follow FAIR principles and make publicly available the code (a Jupyter notebook) and the DeepESD dataset. In particular, DeepESD is published at the Earth System Grid Federation (ESGF), as the first continental-wide PP dataset contributing to CORDEX (EUR-44).This research has been supported by the Spanish Government (MCIN/AEI /10.13039/501100011033) through project CORDyS (grant no. PID2020-116595RB-I00)

Statistical downscaling with the downscaleR package (v3.1.0): contribution to the VALUE intercomparison experiment

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.We thank the European Union Cooperation in Science and Technology (EU COST) Action ES1102 VALUE (http://www.value-cost.eu) for making publicly available the data used in this article and the tools implementing the comprehensive set of validation measures and indices. We also thank the THREDDS Data Server (TDS) software developed by UCAR/Unidata (https://doi.org/10.5065/D6N014KG, Unidata, 2006) and all R developers and their supporting community for providing free software facilitating open science. We acknowledge the World Climate Research Program’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the EC-EARTH Consortium for producing and making available their model output used in this paper. For CMIP the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led the development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We are very grateful to the two anonymous referees participating in the interactive discussion for their insightful comments, helping us to considerably improve the original paper. Financial support. The authors acknowledge partial funding from the MULTI-SDM project (MINECO/FEDER, CGL2015-66583-R) and from the project INDECIS, part of the European Research Area for Climate Services Consortium (ERA4CS) with co-funding by the uropean Union (grant no. 690462)

A comparison of statistical downscaling techniques for daily precipitation: results from the CORDEX flagship pilot study in South America

Southeast South America (SESA) is one of the regions of the planet where extreme precipitation events occur and have high impact on human activities. These extreme events result from the complex interactions of a broad range of scales, therefore their study, modelling and projections in a changing climate continue to be a challenging task. The CORDEX Flagship Pilot Study in South America (FPSSESA) addresses this topic in order to advance in the understanding and modelling of extreme precipitation events based on coordinated experiments using different downscaling approaches. In this work we present the results from the collaborative action to intercompare different statistical downscaling techniques in simulating daily precipitation in SESA with special focus on extremes. To this end, seven statistical downscaling models based on the regression and analog families were evaluated over SESA. The sensitivity to the different predictor and predictand datasets were tested using two reanalyses (ECMWF ERA-Interim and Japanese 55-year Reanalysis JRA-55) and two daily precipitation (station data and MSWEP) datasets. The models were calibrated and cross-validated during the 1979-2009 period and also evaluated in the independent warm season of 2009-2010. This season, with record of extreme precipitation events, is the target season chosen in the FPS-SESA to perform the dynamical downscaling simulations as well, and therefore it allows for comparisons between both approaches. The results show that the methods are more skillful when combined predictors including circulation variables at middle levels and local humidity at low levels of the atmosphere are considered. The performance of the models is also sensitive to reanalysis choice. The methods show overall good performance in simulating daily precipitation characteristics over the region, but no single model performs best over all validation metrics and aspects evaluated.Fil: Bettolli, Maria Laura. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Ciencias de la Atmósfera y los Océanos; ArgentinaFil: Gutiérrez Llorente, José Manuel. Universidad de Cantabria; EspañaFil: Iturbide, Maialen. Universidad de Cantabria; EspañaFil: Baño Medina, Jorge. Universidad de Cantabria; EspañaFil: Huth, Radan. Karlova Univerzita (cuni); República ChecaFil: Solman, Silvina Alicia. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Fernández, Jesús. Universidad de Cantabria; EspañaFil: da Rocha, Rosmeri Porfirio. Universidade de Sao Paulo; BrasilFil: Llopart, Marta. Universidad Estadual de Sao Paulo; BrasilFil: Lavín Gullón, Álvaro. Universidad de Cantabria; EspañaFil: Coppola, Erika. The Abdus Salam; ItaliaFil: Chou, Sin Chan. Centro de Previsao de Tempo e Estudos Climaticos; BrasilFil: Doyle, Moira Evelina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Olmo, Matías Ezequiel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Ciencias de la Atmósfera y los Océanos; ArgentinaFil: Feijoó, Martín. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaInternational Conference on Regional ClimateBeijingChinaWorld Climate Research ProgrammeInstituto Sueco de Meteorología e Hidrologí

Implementation of FAIR principles in the IPCC: the WGI AR6 Atlas repository

The Sixth Assessment Report (AR6) of the Intergovernmental Panel on Climate Change (IPCC) has adopted the FAIR Guiding Principles. We present the Atlas chapter of Working Group I (WGI) as a test case. We describe the application of the FAIR principles in the Atlas, the challenges faced during its implementation, and those that remain for the future. We introduce the open source repository resulting from this process, including coding (e.g., annotated Jupyter notebooks), data provenance, and some aggregated datasets used in some figures in the Atlas chapter and its interactive companion (the Interactive Atlas), open to scrutiny by the scientific community and the general public. We describe the informal pilot review conducted on this repository to gather recommendations that led to significant improvements. Finally, a working example illustrates the re-use of the repository resources to produce customized regional information, extending the Interactive Atlas products and running the code interactively in a web browser using Jupyter notebooks.Peer reviewe