Origin of the warm eastern tropical Atlantic SST bias in a climate model

The substantial warm sea surface temperature bias in the eastern Tropical Atlantic reported in most CMIP5 climate simulations with various models, in particular along the coast of Namibia and Angola, remains an issue in more recent and CMIP6-ready versions of climate models such as EC-Earth. A complete and original set of experiments with EC-Earth3.1 is performed to investigate the causes and mechanisms responsible for the emergence and persistence of this bias. The fully-developed bias is studied in a historical experiment that has reached quasi-equilibrium, while retrospective prediction experiments are used to highlight the development/growth from an observed initial state. Prediction experiments are performed at both low and high resolution to assess the possible dependence of the bias on horizontal resolution. Standalone experiments with the ocean and the atmosphere components of EC-Earth are also analyzed to separate the respective contributions of the ocean and atmosphere to the development of the bias. EC-Earth3.1 exhibits a bias similar to that reported in most climate models that took part in CMIP5. The magnitude of this bias, however, is weaker than most CMIP5 models by few degrees. Increased horizontal resolution only leads to a minor reduction of the bias in EC-Earth. The warm SST bias is found to be the result of an excessive solar absorption in the ocean mixed layer, which can be linked to the excessive solar insolation due to unrealistically low cloud cover, and the absence of spatial and temporal variability of the biological productivity in the ocean component. The warm SST bias is further linked to deficient turbulent vertical mixing of cold water to the mixed layer. Our study points at a need for better representation of clouds in the vicinity of eastern boundaries in atmosphere models, and better representation of solar penetration and turbulent mixing in the ocean models in order to eliminate the Tropical Atlantic biases.We would like to acknowledge the anonymous reviewer who provided constructive comments that led to a considerable improvement of the manuscript. We would also like to thank Aurore Voldoire and Anna-Lena Deppenmeier for the useful discussion, and Yann Planton for providing the code for implementing the tendency diagnostics in NEMO. This research has received funding from the EU Seventh Framework Programme FP7 (2007–2013) under grant agreements 308378 (SPECS), 603521 (PREFACE) and the Horizon 2020 EU program under grand agreements 641727 (PRIMAVERA). We acknowledge RES and ECMWF for awarding access to supercomput- ing facilities in the Barcelona Supercomputing Center in Spain and the ECMWF Supercomputing Center in the UK, through the HiResClim and SPESICCF projects, recpectively. We acknowledge the work of the developers of the s2dverification R-based package (http://cran.r-project. org/web/packages/s2dverification/index.html). The visualization of some of the figures was done with the NCAR Command Language (NCL, Version 6.3.0, 2016, Boulder, Colorado: UCAR/NCAR/CISL/TDD, http://dx.doi.org/10.5065/D6WD3XH5).Peer ReviewedPostprint (author's final draft

Mozambique Channel Eddies in GCMs: A question of resolution and slippage

Hydrographic observations in the 21st century have shown that the flow within the MozambiqueChannel is best described by a series of large poleward-propagating anticyclonic eddies, rather than, aspreviously thought, a continuous intense western boundary current. The portrayal of this region in various runs of the NEMO 75-level model is found to vary between those two descriptions depending upon the resolution used and the implementation of the model's lateral boundary conditions. In a comparison of 1/4 ? resolution runs, the change of these conditions from free-slip to no-slip leads to the mean southward flow moving further offshore, with greater variability in the zonal and meridional velocities as the flow organises itself into eddies, and a reduction in total transport. If a realization of a model is unable to get these aspects of the physical flow correct, then this will significantly reduce its ability to show a realistic biological signal or long-term response to climate change. Further south, beyond Durban, the application of no-slip conditions similarly causes the mean Agulhas Current to lie further offshore, making it much more able to simulate Natal Pulses.<br/

Review and assessment of latent and sensible heat flux accuracy over the global oceans

For over a decade, several research groups have been developing air-sea heat flux information over the global ocean, including latent (LHF) and sensible (SHF) heat fluxes over the global ocean. This paper aims to provide new insight into the quality and error characteristics of turbulent heat flux estimates at various spatial and temporal scales (from daily upwards). The study is performed within the European Space Agency (ESA) Ocean Heat Flux (OHF) project. One of the main objectives of the OHF project is to meet the recommendations and requirements expressed by various international programs such as the World Research Climate Program (WCRP) and Climate and Ocean Variability, Predictability, and Change (CLIVAR), recognizing the need for better characterization of existing flux errors with respect to the input bulk variables (e.g. surface wind, air and sea surface temperatures, air and surface specific humidities), and to the atmospheric and oceanic conditions (e.g. wind conditions and sea state). The analysis is based on the use of daily averaged LHF and SHF and the asso- ciated bulk variables derived from major satellite-based and atmospheric reanalysis products. Inter-comparisons of heat flux products indicate that all of them exhibit similar space and time patterns. However, they also reveal significant differences in magnitude in some specific regions such as the western ocean boundaries during the Northern Hemisphere winter season, and the high southern latitudes. The differences tend to be closely related to large differences in surface wind speed and/or specific air humidity (for LHF) and to air and sea temperature differences (for SHF). Further quality investigations are performed through comprehensive comparisons with daily-averaged LHF and SHF estimated from moorings. The resulting statistics are used to assess the error of each OHF product. Consideration of error correlation between products and observations (e.g., by their assimilation) is also given. This reveals generally high noise variance in all products and a weak signal in common with in situ observations, with some products only slightly better than others. The OHF LHF and SHF products, and their associated error characteristics, are used to compute daily OHF multiproduct-ensemble (OHF/MPE) estimates of LHF and SHF over the ice-free global ocean on a 0.25° × 0.25° grid. The accuracy of this heat multiproduct, determined from comparisons with mooring data, is greater than for any individual product. It is used as a reference for the anomaly characterization of each individual OHF product

Impact of a new anisotropic rheology on simulations of Arctic sea ice

new rheology that explicitly accounts for the subcontinuum anisotropy of the sea ice cover is implemented into the Los Alamos sea ice model. This is in contrast to all models of sea ice included in global circulation models that use an isotropic rheology. The model contains one new prognostic variable, the local structure tensor, which quantifies the degree of anisotropy of the sea ice, and two parameters that set the time scale of the evolution of this tensor. The anisotropic rheology provides a subcontinuum description of the mechanical behavior of sea ice and accounts for a continuum scale stress with large shear to compression ratio and tensile stress component. Results over the Arctic of a stand-alone version of the model are presented and anisotropic model sensitivity runs are compared with a reference elasto-visco-plastic simulation. Under realistic forcing sea ice quickly becomes highly anisotropic over large length scales, as is observed from satellite imagery. The influence of the new rheology on the state and dynamics of the sea ice cover is discussed. Our reference anisotropic run reveals that the new rheology leads to a substantial change of the spatial distribution of ice thickness and ice drift relative to the reference standard visco-plastic isotropic run, with ice thickness regionally increased by more than 1 m, and ice speed reduced by up to 50%

Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean–atmosphere models

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Light, C., Arbic, B., Martin, P., Brodeau, L., Farrar, J., Griffies, S., Kirtman, B., Laurindo, L., Menemenlis, D., Molod, A., Nelson, A., Nyadjro, E., O’Rourke, A., Shriver, J., Siqueira, L., Small, R., & Strobach, E. Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean–atmosphere models. Climate Dynamics, (2022): 1–27, https://doi.org/10.1007/s00382-022-06257-6.High-frequency precipitation variance is calculated in 12 different free-running (non-data-assimilative) coupled high resolution atmosphere–ocean model simulations, an assimilative coupled atmosphere–ocean weather forecast model, and an assimilative reanalysis. The results are compared with results from satellite estimates of precipitation and rain gauge observations. An analysis of irregular sub-daily fluctuations, which was applied by Covey et al. (Geophys Res Lett 45:12514–12522, 2018. https://doi.org/10.1029/2018GL078926) to satellite products and low-resolution climate models, is applied here to rain gauges and higher-resolution models. In contrast to lower-resolution climate simulations, which Covey et al. (2018) found to be lacking with respect to variance in irregular sub-daily fluctuations, the highest-resolution simulations examined here display an irregular sub-daily fluctuation variance that lies closer to that found in satellite products. Most of the simulations used here cannot be analyzed via the Covey et al. (2018) technique, because they do not output precipitation at sub-daily intervals. Thus the remainder of the paper focuses on frequency power spectral density of precipitation and on cumulative distribution functions over time scales (2–100 days) that are still relatively “high-frequency” in the context of climate modeling. Refined atmospheric or oceanic model grid spacing is generally found to increase high-frequency precipitation variance in simulations, approaching the values derived from observations. Mesoscale-eddy-rich ocean simulations significantly increase precipitation variance only when the atmosphere grid spacing is sufficiently fine (< 0.5°). Despite the improvements noted above, all of the simulations examined here suffer from the “drizzle effect”, in which precipitation is not temporally intermittent to the extent found in observations.Support for CXL’s effort on this project was provided by a Research Experiences for Undergraduates (REU) supplement for National Science Foundation (NSF) grant OCE-1851164 to BKA, which also provided partial support for PEM. In addition, BKA acknowledges NSF grant OCE-1351837, which provided partial support for AKO, Office of Naval Research grant N00014-19-1-2712 and NASA grants NNX17AH55G, which also provided partial support for ADN, and 80NSSC20K1135. JTF’s participation, and the SPURS-II buoy data, were funded by NASA grants 80NSSC18K1494 and NNX15AG20G

Eddy-Permitting Ocean Circulation Hindcasts of Past Decades

International audienc

Impact of melt ponds on Arctic sea ice simulations from 1990 to 2007

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2007 and 2011 and climate models did not predict this decline. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds form on Arctic sea ice during the melting season and their presence affects the heat and mass balances of the ice cover, mainly by decreasing the value of the surface albedo by up to 20%. We have developed a melt pond model suitable for forecasting the presence of melt ponds based on sea ice conditions. This model has been incorporated into the Los Alamos CICE sea ice model, the sea ice component of several IPCC climate models. Simulations for the period 1990 to 2007 are in good agreement with observed ice concentration. In comparison to simulations without ponds, the September ice volume is nearly 40% lower. Sensitivity studies within the range of uncertainty reveal that, of the parameters pertinent to the present melt pond parameterization and for our prescribed atmospheric and oceanic forcing, variations of optical properties and the amount of snowfall have the strongest impact on sea ice extent and volume. We conclude that melt ponds will play an increasingly important role in the melting of the Arctic ice cover and their incorporation in the sea ice component of Global Circulation Models is essential for accurate future sea ice forecasts

Oceanic hindcast simulations at high resolution suggest that the Atlantic MOC is bistable

All climate models predict a freshening of the North Atlantic at high latitude that may induce an abrupt change of the Atlantic Meridional Overturning Circulation (hereafter AMOC) if it resides in the bistable regime, where both a strong and a weak state coexist. The latter remains uncertain as there is no consensus among observations and ocean reanalyses, where the AMOC is bistable, versus most climate models that reproduce a mono-stable strong AMOC. A series of four hindcast simulations of the global ocean at 1/12° resolution, which is presently unique, are used to diagnose freshwater transport by the AMOC in the South Atlantic, an indicator of AMOC bistability. In all simulations, the AMOC resides in the bistable regime: it exports freshwater southward in the South Atlantic, implying a positive salt advection feedback that would act to amplify a decreasing trend in subarctic deep water formation as projected in climate scenarios