Improving Madden–Julian Oscillation Simulation in Atmospheric General Circulation Models by Coupling with Snow–Ice–Thermocline One-dimensional Ocean Model

A one-column turbulent kinetic energy–type ocean mixed-layer model Snow–Ice–Thermocline (SIT) when coupled with three atmospheric general circulation models (AGCMs) to yielded superior Madden–Julian Oscillation (MJO) simulation. SIT is designed to have fine layers similar to those observed near the ocean surface and therefore can realistically simulate the diurnal warm layer and cool skin. This refined discretization of the near ocean surface in SIT provides accurate sea surface temperature (SST) simulation, thus facilitating realistic air–sea interaction. Coupling SIT with European Centre Hamburg Model, Version 5 (ECHAM5); Community Atmosphere Model, Version 5 (CAM5); and High Resolution Atmospheric Model (HiRAM) significantly improved MJO simulation in three coupled AGCMs compared with the AGCM driven with prescribed SST. This study suggests two major improvements to the coupling process. First, during the preconditioning phase of MJO over Maritime Continent (MC), the over underestimated surface latent heat bias in AGCMs can be corrected. Second, during the phase of strongest convection over MC, the change of the intraseasonal circulation in the meridional circulation is the dominant factor in the coupled simulations relative to the uncoupled experiments. The study results indicate that a fine vertical resolution near the surface, which better captures temperature variations in the upper few meters of the ocean, considerably improves different models with different configurations and physical parameterization schemes; this could be an essential factor for accurate MJO simulation.

Specified dynamics scheme impacts on wave-mean flow dynamics, convection, and tracer transport in CESM2 (WACCM6)

Specified dynamics schemes are ubiquitous modeling tools for isolating the roles of dynamics and transport on chemical weather and climate. They typically constrain the circulation of a chemistry–climate model to the circulation in a reanalysis product through linear relaxation. However, recent studies suggest that these schemes create a divergence in chemical climate and the meridional circulation between models and do not accurately reproduce trends in the circulation. In this study we perform a systematic assessment of the specified dynamics scheme in the Community Earth System Model version 2, Whole Atmosphere Community Climate Model version 6 (CESM2 (WACCM6)), which proactively nudges the circulation toward the reference meteorology. Specified dynamics experiments are performed over a wide range of nudging timescales and reference meteorology frequencies, with the model's circulation nudged to its own free-running output – a clean test of the specified dynamics scheme. Errors in the circulation scale robustly and inversely with meteorology frequency and have little dependence on the nudging timescale. However, the circulation strength and errors in tracers, tracer transport, and convective mass flux scale robustly and inversely with the nudging timescale. A 12 to 24 h nudging timescale at the highest possible reference meteorology frequency minimizes errors in tracers, clouds, and the circulation, even up to the practical limit of one reference meteorology update every time step. The residual circulation and eddy mixing integrate tracer errors and accumulate them at the end of their characteristic transport pathways, leading to elevated error in the upper troposphere and lower stratosphere and in the polar stratosphere. Even in the most ideal case, there are non-negligible errors in tracers introduced by the nudging scheme. Future development of more sophisticated nudging schemes may be necessary for further progress.

The dependence of mean climate state on shortwave absorption by water vapor

State-of-the-art climate models exhibit significant spread in the climatological value of atmospheric shortwave absorption (SWA). This study investigates both the possible causes and climatic impacts of this SWA inter-model spread. The inter-model spread of global-mean SWA largely originates from the inter-model difference in water vapor shortwave absorptivity. Hence, we alter the water vapor shortwave absorptivity in the Community Earth System Model, version 1, with Atmosphere Model, version 4 (CESM1-CAM4). Increasing the water vapor shortwave absorptivity leads to a reduction in global-mean precipitation and a La Niña-like cooling over the tropical Pacific. The global-mean atmospheric energy budget suggests that the precipitation is suppressed as a way to compensate for the increased SWA. The precipitation reduction is driven by the weakened surface winds, stabilized planetary boundary layer, and surface cooling. The La Niña-like cooling over the tropical Pacific is attributed to the zonal asymmetry of climatological evaporative damping efficiency and the low cloud enhancement over the eastern basin. Complementary fixed SSTs simulations suggest that the latter is more fundamental and that it primarily arises from atmospheric processes. Consistent with our experiments, the CMIP5/6 models with a higher global-mean SWA tend to exhibit the tropical Pacific toward a more La Niña-like mean state, highlighting the possible role of water vapor shortwave absorptivity for shaping the mean-state climate patterns.

Description and evaluation of the community aerosol dynamics model MAFOR v2.0

Numerical models are needed for evaluating aerosol processes in the atmosphere in state-of-the-art chemical transport models, urban-scale dispersion models and climatic models. This article describes a publicly available aerosol dynamics model MAFOR (Multicomponent Aerosol FORmation model; version 2.0); we address the main structure of the model, including the types of operation and the treatments of the aerosol processes. The main advantage of MAFOR v2.0 is the consistent treatment of both the mass- and number-based concentrations of particulate matter. An evaluation of the model is also presented, against a high-resolution observational dataset in a street canyon located in the centre of Helsinki (Finland) during an afternoon traffic rush hour on 13 December 2010. The experimental data included measurements at different locations in the street canyon of ultrafine particles, black carbon, and fine particulate mass PM1. This evaluation has also included an intercomparison with the corresponding predictions of two other prominent aerosol dynamics models, AEROFOR and SALSA. All three models fairly well simulated the decrease of the measured total particle number concentrations with increasing distance from the vehicular emission source. The MAFOR model reproduced the evolution of the observed particle number size distributions more accurately than the other two models. The MAFOR model also predicted the variation of the concentration of PM1 better than the SALSA model. We also analysed the relative importance of various aerosol processes based on the predictions of the three models. As expected, atmospheric dilution dominated over other processes; dry deposition was the second most significant process. Numerical sensitivity tests with the MAFOR model revealed that the uncertainties associated with the properties of the condensing organic vapours affected only the size range of particles smaller than 10 nm in diameter. These uncertainties do not therefore affect significantly the predictions of the whole of the number size distribution and the total number concentration. The MAFOR model version 2 is well documented and versatile to use, providing a range of alternative parametrizations for various aerosol processes. The model includes an efficient numerical integration of particle number and mass concentrations, an operator-splitting of processes, and the use of a fixed sectional method. The model could be used as a module in various atmospheric and climatic models.

Explicit silicate cycling in the Kiel Marine Biogeochemistry Model version 3 (KMBM3) embedded in the UVic ESCM version 2.9

We describe and test a new model of biological marine silicate cycling, implemented in the Kiel Marine Biogeochemical Model version 3 (KMBM3), embedded in the University of Victoria Earth System Climate Model (UVic ESCM) version 2.9. This new model adds diatoms, which are a key component of the biological carbon pump, to an existing ecosystem model. This new model combines previously published parameterizations of a diatom functional type, opal production and export with a novel, temperature-dependent dissolution scheme. Modelled steady-state biogeochemical rates, carbon and nutrient distributions are similar to those found in previous model versions. The new model performs well against independent ocean biogeochemical indicators and captures the large-scale features of the marine silica cycle to a degree comparable to similar Earth system models. Furthermore, it is computationally efficient, allowing both fully coupled, long-timescale transient simulations and “offline” transport matrix spinups. We assess the fully coupled model against modern ocean observations, the historical record starting from 1960 and a business-as-usual atmospheric CO2 forcing to the year 2300. The model simulates a global decline in net primary production (NPP) of 1.4 % having occurred since the 1960s, with the strongest declines in the tropics, northern midlatitudes and Southern Ocean. The simulated global decline in NPP reverses after the year 2100 (forced by the extended RCP8.5 CO2 concentration scenario), and NPP returns to 98 % of the pre-industrial rate by 2300. This recovery is dominated by increasing primary production in the Southern Ocean, mostly by calcifying phytoplankton. Large increases in calcifying phytoplankton in the Southern Ocean offset a decline in the low latitudes, producing a global net calcite export in 2300 that varies only slightly from pre-industrial rates. Diatom distribution moves southward in our simulations, following the receding Antarctic ice front, but diatoms are outcompeted by calcifiers across most of their pre-industrial Southern Ocean habitat. Global opal export production thus drops to 75 % of its pre-industrial value by 2300. Model nutrients such as phosphate, silicate and nitrate build up along the Southern Ocean particle export pathway, but dissolved iron (for which ocean sources are held constant) increases in the upper ocean. This different behaviour of iron is attributed to a reduction of low-latitude NPP (and consequently, a reduction in both uptake and export and particle, including calcite scavenging), an increase in seawater temperatures (raising the solubility of particulate iron) and stratification that “traps” the iron near the surface. These results are meant to serve as a baseline for sensitivity assessments to be undertaken with this model in the future.

The Summer Asia-North America Teleconnection and its Modulation by ENSO in Community Atmosphere Model, Version 5 (CAM5)

Seasonal forecasts of summer continental United States (CONUS) rainfall have relatively low skill, partly due to a lack of consensus about its sources of predictability. The East Asian monsoon (EAM) can excite a cross-Pacific Rossby wave train, also known as the Asia-North America (ANA) teleconnection. In this study, we analyze the ANA teleconnection in observations and model simulations from the Community Atmospheric Model, version 5 (CAM5), comparing experiments with prescribed climatological SSTs and prescribed observed SSTs. Observations indicate a statistically significant relationship between a strong EAM and increased probability of positive precipitation anomalies over the U.S. west coast and the Plains-Midwest. The ANA teleconnection and CONUS rainfall patterns are improved in the CAM5 experiment with prescribed observed SSTs, suggesting that SST variability is necessary to simulate this teleconnection over CONUS. We find distinct ANA patterns between ENSO phases, with the La Niña-related patterns in CAM5_obsSST disagreeing with observations. Using linear steady-state quasi-geostrophic theory, we conclude that incorrect EAM forcing location greatly contributed to CAM5 biases, and jet stream disparities explained the ENSO-related biases. Finally, we compared EAM forcing experiments with different mean states using a simple dry nonlinear atmospheric general circulation model. Overall, the ANA pattern over CONUS and its modulation by ENSO forcing are well described by dry dynamics on seasonal-to-interannual timescales, including the constructive (destructive) interference between El Niño (La Niña) modulation and the ANA patterns over CONUS.

Evaluation de deux méthodes de correction de biais des sorties de modèles climatiques régionaux Cordex-Africa pour la prévision des pluies : cas du bassin côtier oranais

Ce travail vise à évaluer les pluies simulées issues des sorties de modèles climatiques régionaux Cordex-Africa dans le bassin côtier oranais en Algérie. Pour cela les simulations du modèle RCA4 (Rossby Centre Atmosphere model, version 4) forcé par deux modèles de circulation globale (MPI-ESM-LR et CNRM-CM5) sous deux scenarios de forçages radiatifs «Representative Concentration Pathways» (RCPs) RCP 4.5 et RCP 8.5 sont comparées aux pluies observées au niveau de cinq stations pluviométriques, au cours de la période de contrôle 1981–2005 à l'échelle mensuelle. Les données futures simulées sont ensuite corrigées à l'aide de deux méthodes de correction de biais, à savoir, la méthode quantile-quantile et la méthode Delta, afin de mieux analyser leur évolution au cours de la période de projection 2075–2099. Les coefficients d'échange estimés au cours de la période 2075–2099 montrent que les simulations corrigées par la méthode Delta sont moins biaisées que les simulations corrigées par la méthode quantile-quantile. L'analyse de l'évolution future des pluies met en évidence une réduction de −12 % à −38 % d'ici la fin du 21ème siècle selon le RCP 4.5. Cette réduction qui est encore plus importante selon le scénario pessimiste RCP 8.5, risque d'affecter la disponibilité des ressources en eau dans la région qui a connu par le passé une période de sècheresse sévère et persistante. Enfin, cette étude peut être utilisée comme outil d'aide à la décision destiné aux parties prenantes de la gestion intégrée des ressources en eau et de l'agriculture. Néanmoins, pour une meilleure appréciation des impacts socio-économiques, une étude plus approfondie en considérant plusieurs modèles climatiques et d'autres paramètres climatiques, est recommandée.

Decadal climate predictions with the Canadian Earth System Model version 5 (CanESM5)

The Canadian Earth System Model version 5 (CanESM5) developed at Environment and Climate Change Canada's Canadian Centre for Climate Modelling and Analysis (CCCma) is participating in phase 6 of the Coupled Model Intercomparison Project (CMIP6). A 40-member ensemble of CanESM5 retrospective decadal forecasts (or hindcasts) is integrated for 10 years from realistic initial states once a year during 1961 to the present using prescribed external forcing. The results are part of CCCma's contribution to the Decadal Climate Prediction Project (DCPP) component of CMIP6. This paper evaluates CanESM5 large ensemble decadal hindcasts against observational benchmarks and against historical climate simulations initialized from pre-industrial control run states. The focus is on the evaluation of the potential predictability and actual skill of annual and multi-year averages of key oceanic and atmospheric fields at regional and global scales. The impact of initialization on prediction skill is quantified from the hindcasts decomposition into uninitialized and initialized components. The dependence of potential and actual skill on ensemble size is examined. CanESM5 decadal hindcasts skillfully predict upper-ocean states and surface climate with a significant impact from initialization that depend on climate variable, forecast range, and geographic location. Deficiencies in the skill of North Atlantic surface climate are identified and potential causes discussed. The inclusion of biogeochemical modules in CanESM5 enables the prediction of carbon cycle variables which are shown to be potentially skillful on decadal timescales, with a strong long-lasting impact from initialization on skill in the ocean and a moderate short-lived impact on land.

The Use of Remote Sensing-Based ET Estimates to Improve Global Hydrological Simulations in the Community Land Model Version 5.0

Terrestrial evapotranspiration (ET) is a critical component of water and energy cycles, and improving global land evapotranspiration is one of the challenging works in the development of land surface models (LSMs). In this study, we apply a bias correction approach into the Community Land Model version 5.0 (CLM5) globally by utilizing the remote sensing-based ET dataset. Results reveal that the correction approach can alleviate both overestimation and underestimation of ET by CLM5 over the globe. The adjustment to overestimation is generally effective, whereas the effectiveness for underestimation is determined by the ET regime, namely water-limited or energy-limited. In the areas with abundant precipitation, the underestimation is effectively corrected by increasing ET without the water supply limit. In areas with rare precipitation, however, increasing ET is limited by water supply, which leads to an undesirable correction effect. Compared with the ET simulated by CLM5, the bias correction approach can reduce the global-averaged relative bias (RB) and the root mean square error (RMSE) by 51.8% and 65.9% against Global Land Evaporation Amsterdam Model (GLEAM) ET data, respectively. Meanwhile, the correlation coefficient (CC) can also be improved from 0.93 to 0.98. Continentally, the most substantial ET improvement occurs in Asia, with the RB and RMSE decreased by 69.7% (from 7.04% to 2.14%) and 70.2% (from 0.312 mm day−1 to 0.093 mm day−1, equivalent to from 114 mm year−1 to 34 mm year−1), and the CC increased from 0.92 to 0.99, respectively. Consequently, benefiting from the improvement of ET, the simulations of runoff and soil moisture are also improved over the globe and each of the six continents, and the improvement varies with region. This study demonstrates that the use of satellite-based ET products is beneficial to hydrological simulations in land surface models over the globe.