Peak refreezing in the Greenland firn layer under future warming scenarios

Firn (compressed snow) covers approximately 90% of the Greenland ice sheet (GrIS) and currently retains about half of rain and meltwater through refreezing, reducing runoff and subsequent mass loss. The loss of firn could mark a tipping point for sustained GrIS mass loss, since decades to centuries of cold summers would be required to rebuild the firn buffer. Here we estimate the warming required for GrIS firn to reach peak refreezing, using 51 climate simulations statistically downscaled to 1 km resolution, that project the long-term firn layer evolution under multiple emission scenarios (1850–2300). We predict that refreezing stabilises under low warming scenarios, whereas under extreme warming, refreezing could peak and permanently decline starting in southwest Greenland by 2100, and further expanding GrIS-wide in the early 22nd century. After passing this peak, the GrIS contribution to global sea level rise would increase over twenty-fold compared to the last three decades

Understanding advances in the simulation of intraseasonal variability in the ECMWF model. Part II: the application of process-based diagnostics

In Part I of this study it was shown that moving from a moisture-convergent- to a relative-humidity-dependent organized entrainment rate in the formulation for deep convection was responsible for significant advances in the simulation of the Madden – Julian Oscillation (MJO) in the ECMWF model. However, the application of traditional MJO diagnostics were not adequate to understand why changing the control on convection had such a pronounced impact on the representation of the MJO. In this study a set of process-based diagnostics are applied to the hindcast experiments described in Part I to identify the physical mechanisms responsible for the advances in MJO simulation. Increasing the sensitivity of the deep convection scheme to environmental moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid troposphere. Due to the modified precipitation – moisture relationship more moisture is able to build up, which effectively preconditions the tropical atmosphere for the t ransition t o d eep convection. R esults from this study suggest that a tropospheric moisture control on convection is key to simulating the interaction between the convective heating and the large-scale wave forcing associated with the MJO

Finite departure from convective quasi-equilibrium: periodic cycle and discharge-recharge mechanism

A simple self–contained theory is proposed for describing life cycles of convective systems as a discharge–recharge process. A closed description is derived for the dynamics of an ensemble of convective plumes based on an energy cycle. The system consists of prognostic equations for the cloud work function and the convective kinetic energy. The system can be closed by intro ducing a functional relationship between the convective kinetic energy and the cloud–base mass flux. The behaviour of this system is considered under a bulk simplification. Previous cloud–resolving mo delling as well as bulk statistical theories for ensemble convective systems suggest that a plausible relationship would be to assume that the convective kinetic energy is linearly proportional to the cloud–base mass flux. As a result, the system reduces to a nonlinear dynamical system with two dependent variables, the cloud–base mass flux and the cloud work function. The fully nonlinear solution of this system always represents a periodic cycle regardless of the initial condition under constant large–scale forcing. Importantly, the inclusion of energy dissipation in this model does not in itself lead the system to an equilibrium

Description and Demonstration of the Coupled Community Earth System Model v2-Community Ice Sheet Model v2 (CESM2-CISM2)

Earth system/ice-sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea-level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy-balance-based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice-sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions. </div

Vertical structure and physical processes of the Madden-Julian oscillation: synthesis and summary

The "Vertical structure and physical processes of the Madden-Julian oscillation (MJO)" project comprises three experiments, designed to evaluate comprehensively the heating, moistening and momentum associated with tropical convection in general circulation models (GCMs). We consider here only those GCMs that performed all experiments. Some models display relatively higher or lower MJO fidelity in both initialized hindcasts and climate simulations, while others show considerable variations in fidelity between experiments. Fidelity in hindcasts and climate simulations are not meaningfully correlated. The analysis of each experiment led to the development of process-oriented diagnostics, some of which distinguished between GCMs with higher or lower fidelity in that experiment. We select the most discriminating diagnostics and apply them to data from all experiments, where possible, to determine if correlations with MJO fidelity hold across scales and GCM states. While normalized gross moist stability had a small but statistically significant correlation with MJO fidelity in climate simulations, we find no link with fidelity in medium-range hindcasts. Similarly, there is no association between timestep-to-timestep rainfall variability, identified from short hindcasts, and fidelity in medium-range hindcasts or climate simulations. Two metrics that relate precipitation to free-tropospheric moisture--the relative humidity for extreme daily precipitation, and variations in the height and amplitude of moistening with rain rate--successfully distinguish between higher- and lower-fidelity GCMs in hindcasts and climate simulations. To improve the MJO, developers should focus on relationships between convection and both total moisture and its rate of change. We conclude by offering recommendations for further experiments

Application of MJO Simulation Diagnostics to Climate Models

The ability of eight climate models to simulate the Madden-Julian oscillation (MJO) is examined using diagnostics developed by the U.S. Climate Variability and Predictability (CLIVAR) MJO Working Group. Although the MJO signal has been extracted throughout the annual cycle, this study focuses on the boreal winter (November-April) behavior. Initially, maps of the mean state and variance and equatorial space-time spectra of 850-hPa zonal wind and precipitation are compared with observations. Models best represent the intraseasonal space-time spectral peak in the zonal wind compared to that of precipitation. Using the phase-space representation of the multivariate principal components (PCs), the life cycle properties of the simulated MJOs are extracted, including the ability to represent how the MJO evolves from a given subphase and the associated decay time scales. On average, the MJO decay (e-folding) time scale for all models is shorter (~20- 29 days) than observations (~31 days). All models are able to produce a leading pair of multivariate principal components that represents eastward propagation of intraseasonal wind and precipitation anomalies, although the fraction of the variance is smaller than observed for all models. In some cases, the dominant time scale of these PCs is outside of the 30-80-day band. Several key variables associated with the model's MJO are investigated, including the surface latent heat flux, boundary layer (925 hPa) moisture convergence, and the vertical structure of moisture. Low-level moisture convergence ahead (east) of convection is associated with eastward propagation in most of the models. A few models are also able to simulate the gradual moistening of the lower troposphere that precedes observed MJO convection, as well as the observed geographical difference in the vertical structure of moisture associated with the MJO. The dependence of rainfall on lower tropospheric relative humidity and the fraction of rainfall that is stratiform are also discussed, including implications these diagnostics have for MJO simulation. Based on having the most realistic intraseasonal multivariate empirical orthogonal functions, principal component power spectra, equatorial eastward propagating outgoing longwave radiation (OLR), latent heat flux, low-level moisture convergence signals, and vertical structure of moisture over the Eastern Hemisphere, the superparameterized Community Atmosphere Model (SPCAM) and the ECHAM4/Ocean Isopycnal Model (OPYC) show the best skill at representing the MJO.open1149

Does convective aggregation need to be represented in cumulus parameterizations?

Tropical deep convection exhibits a variety of levels of aggregation over a wide range of scales. Based on a multisatellite analysis, the present study shows at mesoscale that different levels of aggregation are statistically associated with differing large-scale atmospheric states, despite similar convective intensity and large-scale forcings. The more aggregated the convection, the dryer and less cloudy the atmosphere, the stronger the outgoing longwave radiation, and the lower the planetary albedo. This suggests that mesoscale convective aggregation has the potential to affect couplings between moisture and convection and between convection, radiation, and large-scale ascent. In so doing, aggregation may play a role in phenomena such as “hot spots” or the Madden-Julian Oscillation. These findings support the need for the representation of mesoscale organization in cumulus parameterizations; most parameterizations used in current climate models lack any such representation. The ability of a cloud system-resolving model to reproduce observed relationships suggests that such models may be useful to guide attempts at parameterizations of convective aggregation

Using a case-study approach to improve the Madden-Julian oscillation in the Hadley Centre model

In its default configuration, the Hadley Centre climate model (GA2.0) simulates roughly one-half the observed level of Madden–Julian oscillation activity, with MJO events often lasting fewer than seven days. We use initialised, climate-resolution hindcasts to examine the sensitivity of the GA2.0 MJO to a range of changes in sub-grid parameterisations and model configurations. All 22 changes are tested for two cases during the Years of Tropical Convection. Improved skill comes only from (a) disabling vertical momentum transport by convection and (b) increasing mixing entrainment and detrainment for deep and mid-level convection. These changes are subsequently tested in a further 14 hindcast cases; only (b) consistently improves MJO skill, from 12 to 22 days. In a 20-year integration, (b) produces near-observed levels of MJO activity, but propagation through the Maritime Continent remains weak. With default settings, GA2.0 produces precipitation too readily, even in anomalously dry columns. Implementing (b) decreases the efficiency of convection, permitting instability to build during the suppressed MJO phase and producing a more favourable environment for the active phase. The distribution of daily rain rates is more consistent with satellite data; default entrainment produces 6–12 mm/day too frequently. These results are consistent with recent studies showing that greater sensitivity of convection to moisture improves the representation of the MJO

Vertical structure and physical processes of the Madden-Julian oscillation: Biases and uncertainties at short range

An analysis of diabatic heating and moistening processes from 12-36 hour lead time forecasts from 12 Global Circulation Models are presented as part of the "Vertical structure and physical processes of the Madden-Julian Oscillation (MJO)" project. A lead time of 12-36 hours is chosen to constrain the large scale dynamics and thermodynamics to be close to observations while avoiding being too close to the initial spin-up for the models as they adjust to being driven from the YOTC analysis. A comparison of the vertical velocity and rainfall with the observations and YOTC analysis suggests that the phases of convection associated with the MJO are constrained in most models at this lead time although the rainfall in the suppressed phase is typically overestimated. Although the large scale dynamics is reasonably constrained, moistening and heating profiles have large inter-model spread. In particular, there are large spreads in convective heating and moistening at mid-levels during the transition to active convection. Radiative heating and cloud parameters have the largest relative spread across models at upper levels during the active phase. A detailed analysis of time step behaviour shows that some models show strong intermittency in rainfall and differences in the precipitation and dynamics relationship between models. The wealth of model outputs archived during this project is a very valuable resource for model developers beyond the study of the MJO. In addition, the findings of this study can inform the design of process model experiments, and inform the priorities for field experiments and future observing systems

Temporal Accumulation and Decision Processes in the Duration Bisection Task Revealed by Contingent Negative Variation

The duration bisection paradigm is a classic task used to examine how humans and other animals perceive time. Typically, participants first learn short and long anchor durations and are subsequently asked to classify probe durations as closer to the short or long anchor duration. However, the specific representations of time and the decision rules applied in this task remain the subject of debate. For example, researchers have questioned whether participants actually use representations of the short and long anchor durations in the decision process rather than merely a response threshold that is derived from those anchor durations. Electroencephalographic (EEG) measures, like the contingent negative variation (CNV), can provide information about the perceptual and cognitive processes that occur between the onset of the timing stimulus and the motor response. The CNV has been implicated as an electrophysiological marker of interval timing processes such as temporal accumulation, representation of the target duration, and the decision that the target duration has been attained. We used the CNV to investigate which durations are involved in the bisection categorization decision. The CNV increased in amplitude up to the value of the short anchor, remained at a constant level until about the geometric mean (GM) of the short and long anchors, and then began to resolve. These results suggest that the short anchor and the GM of the short and long anchors are critical target durations used in the bisection categorization decision process. In addition, larger mean N1P2 amplitude differences were associated with larger amplitude CNVs, which may reflect the participant’s precision in initiating timing on each trial across a test session. Overall, the results demonstrate the value of using scalp-recorded EEG to address basic questions about interval timing