Influence of variations in penetrating solar radiation on the diurnal and intraseasonal structure of the oceanic boundary layer, The

2010 Summer.Includes bibliographical references.The upper portion of the ocean is fairly well mixed and turbulent. The turbulence within the ocean boundary layer (OBL) is regulated by many mechanisms. One process that is receiving a renewed interest is the effect of penetrating component of surface shortwave radiation on ocean dynamics. The influence of solar radiation has been parameterized in two ways. A limited set of models force all the incoming solar radiation to be absorbed in the top model layer. The second parameterization assumes that the irradiance (light) at a given level follows a multiple term exponential. Most commonly it is assumed that shortwave radiation is absorbed in two bands: visible and near infrared. The strength of the infrared absorption is assumed to be fixed. For the visible band, absorption depends on water clarity. Until recently, water clarity could take six different values (Jerlov water types). On climate scales, spatial and temporal variations in water clarity, based on surface chlorophyll, have a strong impact on the simulated ocean temperature, salinity, and momentum. For example, the sea surface temperature (SST) in the cold tongue is reduced. In addition, the strength of the Walker circulation is increased. However, this response is not consistent among different models and parameterizations. When chlorophyll is predicted, the influence of vertically variable water clarity on the thermodynamic and dynamic fields of the ocean can be examined. Studies that have incorporated an ecosystem model find minimal changes relative to using observed surface chlorophyll. Previous research has focused on longer climate time scales and most models do not consider vertical variations in water clarity. In this study the response of the ocean to diurnal and intraseasonal variations of water clarity is examined. The sensitivity to vertical variations in water clarity is also considered. To study the impact of variable solar radiation a model that accurately represents upper ocean physics is required. A new ocean mixing model is proposed that addresses some of the known deficiencies in previous models. The new model predicts entrainment based on turbulence at the OBL base, unlike other ocean models. An over prediction of the vertical heat flux in previous mixed layer models is avoided. The model framework discussed can be easily extended to any coordinate system. Further, this model can be coupled to an ocean biological model, which would determine the water clarity with depth, in a natural way. An evaluation of the new model against observations and a newly developed vector vorticity large eddy simulation (LES) model has shown that the new model preforms as well or better than previous OBL models in certain circumstances. This is especially with low vertical resolution. Since this version of the new model is local, it does not perform as well in pure convective simulations as OBL models with non-local forcing In this new model and K-Profile Parameterization (KPP), the temperature and velocity is very sensitive to variations in water clarity. Trapping more heat near the surface increases the temperature near the surface and confines daytime momentum input to a shallow layer. In addition, the depth of the thermocline is reduced as water clarity decreases. The simulated temperature and velocity fields are insensitive to subsurface variations in water clarity. The responses of the new model and KPP are similar when the turbidity of the column is taken as the near surface average. Two-dimensional simulations examining the influence of spatially variable turbidity lead to a slightly deeper thermocline and weaker near surface velocity relative to simulations with a zonally constant water clarity. It is found that models must allow solar radiation to penetrate beyond the top model level. Further, water clarity should be diagnosed from observed or predicted surface chlorophyll instead of the six Jerlov water types

Comparing ocean surface boundary vertical mixing schemes including langmuir turbulence

Six recent Langmuir turbulence parameterization schemes and five traditional schemes are implemented in a common single‐column modeling framework and consistently compared. These schemes are tested in scenarios versus matched large eddy simulations, across the globe with realistic forcing (JRA55‐do, WAVEWATCH‐III simulated waves) and initial conditions (Argo), and under realistic conditions as observed at ocean moorings. Traditional non‐Langmuir schemes systematically underpredict large eddy simulation vertical mixing under weak convective forcing, while Langmuir schemes vary in accuracy. Under global, realistic forcing Langmuir schemes produce 6% (−1% to 14% for 90% confidence) or 5.2 m (−0.2 m to 17.4 m for 90% confidence) deeper monthly mean mixed layer depths than their non‐Langmuir counterparts, with the greatest differences in extratropical regions, especially the Southern Ocean in austral summer. Discrepancies among Langmuir schemes are large (15% in mixed layer depth standard deviation over the mean): largest under wave‐driven turbulence with stabilizing buoyancy forcing, next largest under strongly wave‐driven conditions with weak buoyancy forcing, and agreeing during strong convective forcing. Non‐Langmuir schemes disagree with each other to a lesser extent, with a similar ordering. Langmuir discrepancies obscure a cross‐scheme estimate of the Langmuir effect magnitude under realistic forcing, highlighting limited understanding and numerical deficiencies. Maps of the regions and seasons where the greatest discrepancies occur are provided to guide further studies and observations

Resolving and parameterising the ocean mesoscale in earth system models

Purpose of Review. Assessment of the impact of ocean resolution in Earth System models on the mean state, variability, and future projections and discussion of prospects for improved parameterisations to represent the ocean mesoscale. Recent Findings. The majority of centres participating in CMIP6 employ ocean components with resolutions of about 1 degree in their full Earth Systemmodels (eddy-parameterising models). In contrast, there are alsomodels submitted toCMIP6 (both DECK and HighResMIP) that employ ocean components of approximately 1/4 degree and 1/10 degree (eddy-present and eddy-rich models). Evidence to date suggests that whether the ocean mesoscale is explicitly represented or parameterised affects not only the mean state of the ocean but also the climate variability and the future climate response, particularly in terms of the Atlantic meridional overturning circulation (AMOC) and the Southern Ocean. Recent developments in scale-aware parameterisations of the mesoscale are being developed and will be included in future Earth System models. Summary. Although the choice of ocean resolution in Earth System models will always be limited by computational considerations, for the foreseeable future, this choice is likely to affect projections of climate variability and change as well as other aspects of the Earth System. Future Earth System models will be able to choose increased ocean resolution and/or improved parameterisation of processes to capture physical processes with greater fidelity

The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model’s strong aerosol-related effective radiative forcing (ERFari+aci = -1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).Plain Language SummaryThe U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1’s capabilities are demonstrated by performing a set of standardized simulation experiments described by the Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima protocol at standard horizontal spatial resolution of approximately 1° latitude and longitude. The model reproduces global and regional climate features well compared to observations. Simulated warming between 1850 and 2015 matches observations, but the model is too cold by about 0.5 °C between 1960 and 1990 and later warms at a rate greater than observed. A thermodynamic analysis of the model’s response to greenhouse gas and aerosol radiative affects may explain the reasons for the discrepancy.Key PointsThis work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System ModelThe performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 yearsE3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m2)Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/1/jame20860_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/2/jame20860.pd

Characterizing Tropical Cyclones in the Energy Exascale Earth System Model Version 1

Abstract In this study, we analyze the realism with which tropical cyclones (TCs) are simulated in the fully coupled low‐ and high‐resolution Energy Exascale Earth System Model (E3SM) version 1, with a focus on the latter. Compared to the low‐resolution (grid spacing of ∼1°), the representation of TCs improves considerably in the high‐resolution configuration (grid spacing of ∼0.25°). Significant improvements are found in the global TC frequency, TC lifetime maximum intensities, and the relative distribution of TCs among the different basins. However, at both resolutions, spurious TC activity is found in some basins, notably in the subtropical regions. Contrasting the simulated large‐scale TC environment with observations reveals that the model environment is unrealistically conducive for TC development in those regions. Further analysis indicates that these biases are likely related to those in thermodynamic potential intensity, caused by systematic SST biases, and vertical wind shear in the coupled model. TC‐ocean interaction is also examined in the high‐resolution configuration of the model. The salient features of the ocean's response to TC‐induced mixing and the ocean's impact on TC intensification are well‐reproduced. Finally, an evaluation of the influence of El Niño Southern Oscillation (ENSO) on TCs in the high‐resolution configuration of the model reveals that the ENSO‐TC relationship in the model has the right sign and is significant for the North Atlantic and Northwest Pacific, albeit weaker than in observations. In summary, the high‐resolution configuration of the E3SM model simulates TC activity reasonably and hence could be a useful tool for TC‐related research

The Community ocean Vertical Mixing (CVMix) project

Community Ocean Vertical Mixing (CVMix)} is a software package that aims to provide transparent, robust, flexible, well documented, shared Fortran source code for use in parameterizing vertical mixing processes in numerical ocean models. The project is focused on developing software for a consensus of closures that return a vertical diffusivity, viscosity, and possibly a non-local transport, with each quantity dependent on the tracer or velocity being mixed. CVMix modules are written as kernals designed for use in a variety of Fortran ocean model codes such as MPAS-ocean, MOM, and POP. CVMix modules use MKS units and expect the same for input and output. Code development occurs within a community of scientists and engineers who make use of CVMix modules for a variety of ocean codes. CVMix modules are freely distributed to the open source community under GPLv2 using an open source methodology

A global perspective on Langmuir turbulence in the ocean surface boundary layer

The turbulent mixing in thin ocean surface boundary layers (OSBL), which occupy the upper 100m or so of the ocean, control the exchange of heat and trace gases between the atmosphere and ocean. Here we show that current parameterizations of this turbulent mixing lead to systematic and substantial errors in the depth of the OSBL in global climate models, which then leads to biases in sea surface temperature. One reason, we argue, is that current parameterizations are missing key surface-wave processes that force Langmuir turbulence that deepens the OSBL more rapidly than steady wind forcing. Scaling arguments are presented to identify two dimensionless parameters that measure the importance of wave forcing against wind forcing, and against buoyancy forcing. A global perspective on the occurrence of wave-forced turbulence is developed using re-analysis data to compute these parameters globally. The diagnostic study developed here suggests that turbulent energy available for mixing the OSBL is under-estimated without forcing by surface waves. Wave-forcing and hence Langmuir turbulence could be important over wide areas of the ocean and in all seasons in the Southern Ocean. We conclude that surface-wave-forced Langmuir turbulence is an important process in the OSBL that requires parameterization

The DOE E3SM Coupled Model Version 1: Description and Results at High Resolution

This study provides an overview of the coupled high-resolution Version 1 of the Energy Exascale Earth System Model (E3SMv1) and documents the characteristics of a 50-year-long high-resolution control simulation with time-invariant 1950 forcings following the HighResMIP protocol. In terms of global root-mean-squared error metrics, this high-resolution simulation is generally superior to results from the low-resolution configuration of E3SMv1 (due to resolution, tuning changes, and possibly initialization procedure) and compares favorably to models in the CMIP5 ensemble. Ocean and sea ice simulation is particularly improved, due to better resolution of bathymetry, the ability to capture more variability and extremes in winds and currents, and the ability to resolve mesoscale ocean eddies. The largest improvement in this regard is an ice-free Labrador Sea, which is a major problem at low resolution. Interestingly, several features found to improve with resolution in previous studies are insensitive to resolution or even degrade in E3SMv1. Most notable in this regard are warm bias and associated stratocumulus deficiency in eastern subtropical oceans and lack of improvement in El Niño. Another major finding of this study is that resolution increase had negligible impact on climate sensitivity (measured by net feedback determined through uniform +4K prescribed sea surface temperature increase) and aerosol sensitivity. Cloud response to resolution increase consisted of very minor decrease at all levels. Large-scale patterns of precipitation bias were also relatively unaffected by grid spacing

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations - a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100¯km atmosphere and land, a 0.125 river-routing model, and 14→60¯km ocean and sea ice. By design, the computational cost of NARRM is 1/43× of the uniform low-resolution (LR) model at 100¯km but only 1/4¯10¯%-20¯% of a globally uniform high-resolution model at 25¯km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air-sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land-atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations