Predicting atmospheric optical properties for radiative transfer computations using neural networks.

The radiative transfer equations are well known, but radiation parametrizations in atmospheric models are computationally expensive. A promising tool for accelerating parametrizations is the use of machine learning techniques. In this study, we develop a machine learning-based parametrization for the gaseous optical properties by training neural networks to emulate a modern radiation parametrization (RRTMGP). To minimize computa- tional costs, we reduce the range of atmospheric conditions for which the neural networks are applicable and use machine-specific optimized BLAS functions to accelerate matrix computations. To generate training data, we use a set of randomly perturbed atmospheric profiles and calculate optical properties using RRTMGP. Predicted optical properties are highly accurate and the resulting radiative fluxes have average errors within 0.5 W m-2 compared to RRTMGP. Our neural network-based gas optics parametrization is up to four times faster than RRTMGP, depending on the size of the neural networks. We further test the trade-off between speed and accuracy by training neural networks for the narrow range of atmospheric conditions of a single large-eddy simulation, so smaller and therefore faster networks can achieve a desired accuracy. We conclude that our machine learning-based parametrization can speed-up radiative transfer computations while retaining high accuracy. This article is part of the theme issue 'Machine learning for weather and climate modelling'.

The Time Scale of Shallow Convective Self-Aggregation in Large-Eddy Simulations Is Sensitive to Numerics.

Numerical simulations of the tropical mesoscales often exhibit a self-reinforcing feedback between cumulus convection and shallow circulations, which leads to the self-aggregation of clouds into large clusters. We investigate whether this basic feedback can be adequately captured by large-eddy simulations (LESs). To do so, we simulate the non-precipitating, cumulus-topped boundary layer of the canonical "BOMEX" case over a range of numerical settings in two models. Since the energetic convective scales underpinning the self-aggregation are only slightly larger than typical LES grid spacings, aggregation timescales do not converge even at rather high resolutions (<100 m). Therefore, high resolutions or improved sub-filter scale models may be required to faithfully represent certain forms of trade-wind mesoscale cloud patterns and self-aggregating deep convection in large-eddy and cloud-resolving models, and to understand their significance relative to other processes that organize the tropical mesoscales.

Advancing understanding of land-atmosphere interactions by breaking discipline and scale barriers.

Vegetation and atmosphere processes are coupled through a myriad of interactions linking plant transpiration, carbon dioxide assimilation, turbulent transport of moisture, heat and atmospheric constituents, aerosol formation, moist convection, and precipitation. Advances in our understanding are hampered by discipline barriers and challenges in understanding the role of small spatiotemporal scales. In this perspective, we propose to study the atmosphere-ecosystem interaction as a continuum by integrating leaf to regional scales (multiscale) and integrating biochemical and physical processes (multiprocesses). The challenges ahead are (1) How do clouds and canopies affect the transferring and in-canopy penetration of radiation, thereby impacting photosynthesis and biogenic chemical transformations? (2) How is the radiative energy spatially distributed and converted into turbulent fluxes of heat, moisture, carbon, and reactive compounds? (3) How do local (leaf-canopy-clouds, 1 m to kilometers) biochemical and physical processes interact with regional meteorology and atmospheric composition (kilometers to 100 km)? (4) How can we integrate the feedbacks between cloud radiative effects and plant physiology to reduce uncertainties in our climate projections driven by regional warming and enhanced carbon dioxide levels? Our methodology integrates fine-scale explicit simulations with new observational techniques to determine the role of unresolved small-scale spatiotemporal processes in weather and climate models.

Soil drought can mitigate deadly heat stress thanks to a reduction of air humidity.

Global warming increases the number and severity of deadly heatwaves. Recent heatwaves often coincided with soil droughts that intensify air temperature but lower air humidity. Since lowering air humidity may reduce human heat stress, the net impact of soil desiccation on the morbidity and mortality of heatwaves remains unclear. Combining weather balloon and satellite observations, atmospheric modelling, and meta-analyses of heatwave mortality, we find that soil droughts­despite their warming effect­lead to a mild reduction in heatwave lethality. More specifically, morning dry soils attenuate afternoon heat stress anomaly by ~5%. This occurs because of reduced surface evaporation and increased entrainment of dry air aloft. The benefit appears more pronounced during specific events, such as the Chicago 1995 and Northern U.S. 2006 and 2012 heatwaves. Our findings suggest that irrigated agriculture may intensify lethal heat stress, and question recently proposed heatwave mitigation measures involving surface moistening to increase evaporative cooling.

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations.

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

Towards Adaptive Grids for Atmospheric Boundary-Layer Simulations.

We present a proof-of-concept for the adaptive mesh refinement method applied to atmospheric boundary-layer simulations. Such a method may form an attractive alternative to static grids for studies on atmospheric flows that have a high degree of scale separation in space and/or time. Examples include the diurnal cycle and a convective boundary layer capped by a strong inversion. For such cases, large-eddy simulations using regular grids often have to rely on a subgrid-scale closure for the most challenging regions in the spatial and/or temporal domain. Here we analyze a flow configuration that describes the growth and subsequent decay of a convective boundary layer using direct numerical simulation (DNS). We validate the obtained results and benchmark the performance of the adaptive solver against two runs using fixed regular grids. It appears that the adaptive-mesh algorithm is able to coarsen and refine the grid dynamically whilst maintaining an accurate solution. In particular, during the initial growth of the convective boundary layer a high resolution is required compared to the subsequent stage of decaying turbulence. More specifically, the number of grid cells varies by two orders of magnitude over the course of the simulation. For this specific DNS case, the adaptive solver was not yet more efficient than the more traditional solver that is dedicated to these types of flows. However, the overall analysis shows that the method has a clear potential for numerical investigations of the most challenging atmospheric cases.

Observational evidence for cloud cover enhancement over western European forests.

Forests impact regional hydrology and climate directly by regulating water and heat fluxes. Indirect effects through cloud formation and precipitation can be important in facilitating continental-scale moisture recycling but are poorly understood at regional scales. In particular, the impact of temperate forest on clouds is largely unknown. Here we provide observational evidence for a strong increase in cloud cover over large forest regions in western Europe based on analysis of 10 years of 15 min resolution data from geostationary satellites. In addition, we show that widespread windthrow by cyclone Klaus in the Landes forest led to a significant decrease in local cloud cover in subsequent years. Strong cloud development along the downwind edges of larger forest areas are consistent with a forest-breeze mesoscale circulation. Our results highlight the need to include impacts on cloud formation when evaluating the water and climate services of temperate forests, in particular around densely populated areas.