Ed Gerber

Ed Gerber

Professor at the Courant Institute

How does the numerical formulation of an atmospheric model affect trace gas transport?

2 minute read

Published:

Please see our new paper exploring the role of model numerics on trace gas transport, just submitted to the Quarterly Journal of the Royal Meteorological Society.

This study, lead by Aman Gupta, develops new tests, or recipes, to assess the impact of model numerics on the coupled stratosphere-troposphere system. The transport of trace gases by the atmospheric circulation – in particular, radiatively active species such as water vapor and ozone, and their chemical precursors – plays an important role in the climate system and its response to external forcing. Transport presents a challenge for Atmospheric General Circulation Models (AGCMs), where errors in both the resolved circulation and the numerical representation of trace gases can bias transport.

We develope two benchmark tests to assess transport by the “dynamical core” of an AGCM. A dynamical core is is the numerical model which integrates the equations of motion, forming the “core” of a AGCM, which couples all the “physics” of the atmosphere to the dynamics. To separate transport from chemistry, the tests focus on the “age-of-air”, an estimate of the mean transport time by the circulation. The overturning circulation of the stratosphere, or Brewer-Dobson Circulation, provides a stiff test of model numerics, as the transport time scales are on the order of months to years. Small errors in the numerical representation of trace gases have plenty of time to add up into large biases.

To establish a benchmark, four dynamical cores employing different numerical schemes (finite-volume, pseudo- spectral, and spectral-element) and discretizations (cubed sphere vs. latitude-longitude) are compared across a range of resolutions. We found that the subtle momentum balance of the tropical stratosphere is very sensitive to model numerics, and the first intercomparison reveals stark differences in climatological circulation, particularly at high vertical resolution; some cores develop westerly jets and others, easterly jets. This leads to substantial spread in transport, biasing the age-of-air by up to 25% relative to its climatological mean, making it difficult to assess the impact of the numerical representation of trace gases on transport.

This uncertainty is removed by constraining the tropical winds in the second intercomparison test, akin to specifying the Quasi-Biennial Oscillation in an AGCM. The dynamical cores exhibit qualitative agreement on the structure of atmospheric transport in the second test, with evidence of convergence as the horizontal and vertical resolution is increased. Significant quantitative differences remain, however, particularly between models employing spectral vs. finite-volume numerics, even in state-of-the-art cores.

You May Also Enjoy

Stressed out jet streams!

1 minute read

Published:

Xingjian (Ken) Yan, a precocious undergraduate (now bound for a PhD at MIT) working with Lei Wang and I just submitted a paper exploring the utility of the traffic jam theory of blocking onset for perdiction to Geophysical Research Letters. Ken defined and explored “flux exceedance events”, meteorological situations where the jet stream gets overloaded with storm activity. Nakamura and Huang suggested that this overloaded jet situation creates a pile up storm activity – a traffic jam – leading to blocking events. Ken found that the climatological structure of exceedance events is remarkably similar to that of blocks, but that they appear to be distinct phenomenon: an overloaded jet stream is unfortunately not a reliable harbinger of an atmospheric block.

Protection without poison

1 minute read

Published:

Aaron Match developed a new theory for why the ozone layer is up in the stratosphere, reaching a maximum 26 km above the surface in paper submitted to the journal of Atmospheric Chemistry and Physics. In 1880, Walter Hartley deduced that ozone must be absorbing UV-B and UV-C radiation from the sun, but since it isn’t present at the surface (except in polluted air), this ozone must be somewhere “up there”. The stratosphere wasn’t discovered for a couple more decades, but it turns out that ozone that shields us from this harmful UV radiation is largely between 16 and 40 km above us. It’s a good thing its there: ozone is toxic. If it were uniformly distributed through the atmosphere, the concentration at the surface would be 8 times the EPA safety limit. Why is ozone safely up in the stratosphere, where it protects us without poisoning us? Advanced chemistry climate models can accurately predict the distribution of ozone, but in turns out our text book understanding of the ozone layer were incomplete and gave the wrong explanation for why it reaches a maximum in the stratosphere at 26 km.

Using Explainable AI and Transfer Learning to understand and predict the maintenance of Atlantic blocking with limited observational data

1 minute read

Published:

Huan Zhang submitted our paper on how to use machine learning to predict the persistence of blocking events to the new Journal of Geophysical Research Machine Learning and Computation, published by the AGU. Blocking events are persistent high pressure systems that “block” the flow of the jet stream. They are associated with extreme weather, as they shift the direction of storms, and, in summer, create heat domes that drive heat waves. A key element of a block is its persistence. Huan developed a convolutional neural network to predict whether an nascent blocking anomaly would persist, or fade away, and then interrogated the network to understand why it worked, and how it could be trained with the short observational record.

So, you have an imperfect data-driven parameterization for climate modeling. How can you make it better?

1 minute read

Published:

Minah Yang submitted a paper on how to combat data imbalance in regression problems to JAMES. The goal is to improve the performance of data-driven paramterizations, particularly for profiles that are rare, but important. This is a data imbalance problem: we need ensure that the parameterization works well on input-output pairs that are seldom seen in training. Minah proposes a technique based on histogram equalization, visualized below with help from Cece, Minah’s faithful companion! The idea is to oversample or reweight these rare cases during training, to ensure the method learns from them.