Search

Hot, cold and dry meteorological extremes are often linked with severe impacts on the public health, agricultural, energy and environmental sectors. Skillful predictions of such extremes could therefore enable stakeholders to better plan and adapt to future impacts of these events. The intensity, duration and frequency of such extremes are affected by anthropogenic climate change and modulated by different modes of climate variability. Here we use a large multi-model ensemble from the Coupled Model Intercomparison Project Phase 6 and constrain these simulations by sub-selecting those members whose global SST anomaly patterns are most similar to observations at a given point in time, thereby phasing in the decadal climate variability with observations. Hot and cold extremes are skillfully predicted over most of the globe, with also a widespread added value from using the constrained ensemble compared to the unconstrained full CMIP6 ensemble. On the other hand, dry extremes show skill only in some regions with results sensitive to the index used. Still, we find skillful predictions and added skill for dry extremes in some regions such as western north America, southern central and eastern Europe, southeastern Australia, southern Africa and the Arabian peninsula. We also find that the added skill in the constrained ensemble is due to a combination of improved multi-decadal variations in phase with observed climate extremes and improved representation of long-term changes. Our results demonstrate that constraining decadal variability in climate projections can provide improved estimates of temperature extremes and drought in the next twenty years, which can inform targeted adaptation strategies to near-term climate change., This research has been partly supported by the Horizon2020 LANDMARC project (grant agreement No. 869367) and the Horizon Europe ASPECT project (grant number 101081460). PDL has received funding from the European Union’s Horizon Europe Research and Innovation Programme under grant agreement No 101059659. CDT acknowledges financial support from the Spanish Ministry for Science and Innovation (FPI PRE2019-509 08864 financed by MCIN/AEI/10.13039/501100011033 and by FSE invierte en tu futuro). MGD is grateful for support by the AXA Research Fund. The authors are further grateful for the support by the Department of Research and Universities of the Government of Catalonia to the Climate Variability and Change Research Group (Code: 2021 SGR 00786)., Peer Reviewed, Postprint (author's final draft)

Future precipitation changes are typically estimated from climate model simulations, while the credibility of such projections needs to be assessed by their ability to capture observed precipitation changes. Here we evaluate how skillfully historical climate simulations contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6) capture observed changes in mean and extreme precipitation. We find that CMIP6 historical simulations skillfully represent observed precipitation changes over large parts of Europe, Asia, northeastern North America, parts of South America and western Australia, whereas a lack of skill is apparent in western North America and parts of Africa. In particular in regions with moderate skill the availability of very large ensembles can be beneficial to improve the simulation accuracy. CMIP6 simulations are regionally skillful where they capture observed (positive or negative) trends, whereas a lack of skill is found in regions characterized by negative observed precipitation trends where CMIP6 simulates increases., We are grateful for support by the Departament de Recerca i Universitats de la Generalitat de Catalunya for the Climate Variability and Change (CVC) Research Group (Reference: 2021 SGR 00786), and research funding by the Horizon 2020 LANDMARC project (grant agreement no. 869367), the Horizon Europe ASPECT project (Grant 101081460), and the AXA Research Fund. CDT acknowledges financial support from the Spanish Ministry for Science and Innovation (FPI PRE2019–509 08864 financed by MCIN/AEI/http://doi.org/10.13039/501100011033). PDL received funding from the Horizon Europe Research and Innovation Programme, Grant 101059659. We thank the climate modeling groups contributing to CMIP6 for producing and making available their model output. We are grateful to Margarida Samsó and Pierre–Antoine Bretonnière for downloading, formatting and managing the large data sets of climate simulations and observations used in this study., Peer Reviewed, Postprint (published version)

While carbon dioxide is the main cause for global warming, modeling short-lived climate forcers (SLCFs) such as methane, ozone, and particles in the Arctic allows us to simulate near-term climate and health impacts for a sensitive, pristine region that is warming at 3 times the global rate. Atmospheric modeling is critical for understanding the long-range transport of pollutants to the Arctic, as well as the abundance and distribution of SLCFs throughout the Arctic atmosphere. Modeling is also used as a tool to determine SLCF impacts on climate and health in the present and in future emissions scenarios. In this study, we evaluate 18 state-of-the-art atmospheric and Earth system models by assessing their representation of Arctic and Northern Hemisphere atmospheric SLCF distributions, considering a wide range of different chemical species (methane, tropospheric ozone and its precursors, black carbon, sulfate, organic aerosol, and particulate matter) and multiple observational datasets. Model simulations over 4 years (2008–2009 and 2014–2015) conducted for the 2022 Arctic Monitoring and Assessment Programme (AMAP) SLCF assessment report are thoroughly evaluated against satellite, ground, ship, and aircraft-based observations. The annual means, seasonal cycles, and 3-D distributions of SLCFs were evaluated using several metrics, such as absolute and percent model biases and correlation coefficients. The results show a large range in model performance, with no one particular model or model type performing well for all regions and all SLCF species. The multi-model mean (mmm) was able to represent the general features of SLCFs in the Arctic and had the best overall performance. For the SLCFs with the greatest radiative impact (CH4, O3, BC, and SO), the mmm was within ±25 % of the measurements across the Northern Hemisphere. Therefore, we recommend a multi-model ensemble be used for simulating climate and health impacts of SLCFs. Of the SLCFs in our study, model biases wer, Assessments from the Russian ship-based campaign were performed with the support of RFBR project no. 20-55-12001 and according to the development program of the Interdisciplinary Scientific and Educational School of M.V. Lomonosov Moscow State University “Future Planet and Global Environmental Change”. Development of the methodology for aethalometric data treatment was supported by RSF project no. 19-77-30004. The BC observations on R/V Mirai were supported by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan (Arctic Challenge for Sustainability (ArCS) project). Contributions by SMHI were funded by the Swedish Environmental Protection Agency under contract NV-03174-20 and the Swedish Climate and Clean Air Research program (SCAC) as well as partly by the Swedish National Space Board (NORD-SLCP, grant agreement ID: 94/16) and the EU Horizon 2020 project Integrated Arctic Observing System (INTAROS, grant agreement ID: 727890). Work on ACE-FTS analysis was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC). Julia Schmale received funding from the Swiss National Science Foundation (project no. 200021_188478). Duncan Watson-Parris received funding from NERC projects NE/P013406/1 (A-CURE) and NE/S005390/1 (ACRUISE) as well as funding from the European Union's Horizon 2020 research and innovation program iMIRACLI under Marie Skłodowska-Curie grant agreement no. 860100. LATMOS has been supported by the EU iCUPE (Integrating and Comprehensive Understanding on Polar Environments) project (grant agreement no. 689443) under the European Network for Observing our Changing Planet (ERA-Planet), as well as access to IDRIS HPC resources (GENCI allocation A009017141) and the IPSL mesoscale computing center (CICLAD: Calcul Intensif pour le CLimat, l’Atmosphère et la Dynamique) for model simulations. Naga Oshima was supported by the Japan Society for the Promotion of Science KAKENHI (grant nos. JP18H03363, JP18H05292, and, Peer Reviewed, "Article signat per més de 50 autors/es: Cynthia H. Whaley, Rashed Mahmood, Knut von Salzen, Barbara Winter, Sabine Eckhardt, Stephen Arnold, Stephen Beagley, Silvia Becagli, Rong-You Chien, Jesper Christensen, Sujay Manish Damani, Xinyi Dong, Konstantinos Eleftheriadis, Nikolaos Evangeliou, Gregory Faluvegi, Mark Flanner, Joshua S. Fu, Michael Gauss, Fabio Giardi, Wanmin Gong, Jens Liengaard Hjorth, Lin Huang, Ulas Im, Yugo Kanaya, Srinath Krishnan, Zbigniew Klimont, Thomas Kühn, Joakim Langner, Kathy S. Law, Louis Marelle, Andreas Massling, Dirk Olivié, Tatsuo Onishi, Naga Oshima, Yiran Peng, David A. Plummer, Olga Popovicheva, Luca Pozzoli, Jean-Christophe Raut, Maria Sand, Laura N. Saunders, Julia Schmale, Sangeeta Sharma, Ragnhild Bieltvedt Skeie, Henrik Skov, Fumikazu Taketani, Manu A. Thomas, Rita Traversi, Kostas Tsigaridis, Svetlana Tsyro, Steven Turnock, Vito Vitale, Kaley A. Walker, Minqi Wang, Duncan Watson-Parris, and Tahya Weiss-Gibbons ", Postprint (published version)

Model initialization is a matter of transferring the observed information available at the start of a forecast to the model. An optimal initialization is generally recognized to be able to improve climate predictions up to a few years ahead. However, systematic errors in models make the initialization process challenging. When the observed information is transferred to the model at the initialization time, the discrepancy between the observed and model mean climate causes the drift of the prediction toward the model-biased attractor. Although such drifts can be generally accounted for with a posteriori bias correction techniques, the bias evolving along the prediction might affect the variability that we aim at predicting, and disentangling the small magnitude of the climate signal from the initial drift to be removed represents a challenge. In this study, we present an innovative initialization technique that aims at reducing the initial drift by performing a quantile matching between the observed state at the initialization time and the model state distribution. The adjusted initial state belongs to the model attractor and the observed variability amplitude is scaled toward the model one. Multi-annual climate predictions integrated for 5 years and run with the EC-Earth3 Global Coupled Model have been initialized with this novel methodology, and their prediction skill has been compared with the non-initialized historical simulations from CMIP6 and with the same decadal prediction system but based on full-field initialization. We perform a skill assessment of the surface temperature, the heat content in the ocean upper layers, the sea level pressure, and the barotropic ocean circulation. The added value of the quantile matching initialization is shown in the North Atlantic subpolar region and over the North Pacific surface temperature as well as for the ocean heat content up to 5 years. Improvements are also found in the predictive skill of the Atlantic Meridional O, This study was supported by the project LISTEN funded by the European Commission Horizon 2020 Marie Skłodowska-Curie Actions - IF (GA 799930). The authors thankfully acknowledge the computer resources from the ECMWF special project INCIPIT (spitvolp) and the technical assistance provided by ECMWF and BSC. The climate simulations have been performed using Autosubmit workflow manager (Manubens-Gil et al., 2016). The authors thank Paolo Davini for providing the restart files of the historical simulation used to implement the quantile matching. We acknowledge Saskia Loosveldt and the earthdiags and ESMValTool suite developers, as well as the startR and s2dverification (Manubens et al., 2018) software packages developers, as these tools were used to postprocess, analyze, and visualize the results presented in this work. PO acknowledges support by the Spanish Ministry of Economy, Industry and Competitiveness through the Ramon y Cajal grant (RYC-2017-22772)., Peer Reviewed, Postprint (published version)

The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990-2014) and future (2015-2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (> 60 degrees N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Inter-comparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO42-), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO(4)(2-)burdens decrease significantly in all simulations by 10 %-60% following the reductions of 7 %-78% in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030-2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol-radiation interactions (RFARI) of -0.39 +/- 0.01Wm(-2), which is -0.08Wm(-2) larger than the 1990-2010 mean forcing (-0.32Wm(-2)), of which -0.24 +/- 0.01Wm(-2) was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of --0.35 to -0.40Wm(-2) for the same period, which is -0.01 to -0.06Wm(-2) larger than the 199

The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60∘ N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO2−4), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO2−4 burdens decrease significantly in all simulations by 10 %–60 % following the reductions of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of −0.39±0.01 W m−2, which is −0.08 W m−2 larger than the 1990–2010 mean forcing (−0.32 W m−2), of which −0.24±0.01 W m−2 was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W m−2 for the same period, which is −0.01 to −0.06 W m−2 larger than the 1990–2010 mean forcing of, This research has been supported by the Aarhus University Interdisciplinary Centre for Climate Change (iClimate) OH fund (no. 2020-0162731), the FREYA project funded by the Nordic Council of Ministers (grant agreement nos. MST-227-00036 and MFVM-2019-13476), and the EVAM-SLCF funded by the Danish Environmental Agency (grant agreement no. MST-112-00298). Kostas Tsigaridis and Gregory Faluvegi thank the NASA Modeling, Analysis and Prediction program (MAP) for support. Zbigniew Klimont was financially supported by the EU-funded Action on Black Carbon in the Arctic (EUA-BCA) under the EU Partnership Instrument. Joshua P. French was partially supported by NSF award 1915277., Peer Reviewed, Postprint (published version)

Observations facilitate model evaluation and provide constraints that are relevant to future predictions and projections. Constraints for uninitialized projections are generally based on model performance in simulating climatology and climate change. For initialized predictions, skill scores over the hindcast period provide insight into the relative performance of models, and the value of initialization as compared to projections. Predictions and projections combined can, in principle, provide seamless decadal to multi-decadal climate information. For that, though, the role of observations in skill estimates and constraints needs to be understood in order to use both consistently across the prediction and projection time horizons. This paper discusses the challenges in doing so, illustrated by examples of state-of-the-art methods for predicting and projecting changes in European climate. It discusses constraints across prediction and projection methods, their interpretation, and the metrics that drive them such as process accuracy, accurate trends or high signal-to-noise ratio. We also discuss the potential to combine constraints to arrive at more reliable climate prediction systems from years to decades. To illustrate constraints on projections, we discuss their use in the UK's climate prediction system UKCP18, the case of model performance weights obtained from the Climate model Weighting by Independence and Performance (ClimWIP) method, and the estimated magnitude of the forced signal in observations from detection and attribution. For initialized predictions, skill scores are used to evaluate which models perform well, what might contribute to this performance, and how skill may vary over time. Skill estimates also vary with different phases of climate variability and climatic conditions, and are influenced by the presence of external forcing. This complicates the systematic use of observational constraints. Furthermore, we illustrate that sub-selecting simulati, All authors were supported by the EUCP project funded by the European Commission's Horizon 2020 programme, Grant Agreement number 776613. JM was also supported by the french ANR MOPGA project ARCHANGE and by the EU-H2020 Blue Action (GA 727852) and 4C projects (GA 821003). MGD also received funding by the Spanish Ministry for the Economy, Industry and Competitiveness grant reference RYC-2017-22964., Peer Reviewed, "Article signat per 14 autors/es: Gabriele C. Hegerl, Andrew P. Ballinger, Ben B. B. Booth, Leonard F. Borchert, Lukas Brunner, Markus G. Donat, Francisco J. Doblas-Reyes, Glen R. Harris, Jason Lowe, Rashed Mahmood, Juliette Mignot, James M. Murphy, Didier Swingedouw and Antje Weisheimer", Postprint (published version)

Low clouds persist in the summer Arctic with important consequences for the radiation budget. In this study, we simulate the linear relationship between liquid water content (LWC) and cloud droplet number concentration (CDNC) observed during an aircraft campaign based out of Resolute Bay, Canada, conducted as part of the Network on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments study in July 2014. Using a single-column model, we find that autoconversion can explain the observed linear relationship between LWC and CDNC. Of the three autoconversion schemes we examined, the scheme using continuous drizzle (Khairoutdinov and Kogan, 2000) appears to best reproduce the observed linearity in the tenuous cloud regime (Mauritsen et al., 2011), while a scheme with a threshold for rain (Liu and Daum, 2004) best reproduces the linearity at higher CDNC. An offline version of the radiative transfer model used in the Canadian Atmospheric Model version 4.3 is used to compare the radiative effects of the modelled and observed clouds. We find that there is no significant difference in the upward longwave cloud radiative effect at the top of the atmosphere from the three autoconversion schemes (p=0.05) but that all three schemes differ at p=0.05 from the calculations based on observations. In contrast, the downward longwave and shortwave cloud radiative effect at the surface for the Wood (2005b) and Khairoutdinov and Kogan (2000) schemes do not differ significantly (p=0.05) from the observation-based radiative calculations, while the Liu and Daum (2004) scheme differs significantly from the observation-based calculation for the downward shortwave but not the downward longwave fluxes., This research has been supported by the Natural Sciences and Engineering Research Council of Canada (Discovery Grants RGPIN-2014-05173 and RGPIN 155649) and the Marine Environmental Observation, Prediction and Response Network (MEOPAR), which is a federally funded Networks of Centres of Excellence (NCE) (EC1-RC-DAL)., Peer Reviewed, Postprint (published version)

Motivated by the need to predict how the Arctic atmosphere will change in a warming world, this article summarizes recent advances made by the research consortium NETCARE (Network on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments) that contribute to our fundamental understanding of Arctic aerosol particles as they relate to climate forcing. The overall goal of NETCARE research has been to use an interdisciplinary approach encompassing extensive field observations and a range of chemical transport, earth system, and biogeochemical models. Several major findings and advances have emerged from NETCARE since its formation in 2013. (1) Unexpectedly high summertime dimethyl sulfide (DMS) levels were identified in ocean water (up to 75 nM) and the overlying atmosphere (up to 1 ppbv) in the Canadian Arctic Archipelago (CAA). Furthermore, melt ponds, which are widely prevalent, were identified as an important DMS source (with DMS concentrations of up to 6 nM and a potential contribution to atmospheric DMS of 20 % in the study area). (2) Evidence of widespread particle nucleation and growth in the marine boundary layer was found in the CAA in the summertime, with these events observed on 41 % of days in a 2016 cruise. As well, at Alert, Nunavut, particles that are newly formed and grown under conditions of minimal anthropogenic influence during the months of July and August are estimated to contribute 20 % to 80 % of the 30–50 nm particle number density. DMS-oxidation-driven nucleation is facilitated by the presence of atmospheric ammonia arising from seabird-colony emissions, and potentially also from coastal regions, tundra, and biomass burning. Via accumulation of secondary organic aerosol (SOA), a significant fraction of the new particles grow to sizes that are active in cloud droplet formation. Although the gaseous precursors to Arctic marine SOA remain poorly defined, the measured levels of common continental SOA precursors (isoprene