Source author record

Tapio Schneider

Tapio Schneider appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.ao-ph physics.flu-dyn astro-ph.EP cond-mat.stat-mech physics.comp-ph physics.geo-ph astro-ph math.DS math.NA nlin.CD Numerical Analysis physics.data-an

Catalog footprint

What is connected

11works

12topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Iterated Kalman Methodology For Inverse Problems

This paper is focused on the optimization approach to the solution of inverse problems. We introduce a stochastic dynamical system in which the parameter-to-data map is embedded, with the goal of employing techniques from nonlinear Kalman filtering to estimate the parameter given the data. The extended Kalman filter (which we refer to as ExKI in the context of inverse problems) can be effective for some inverse problems approached this way, but is impractical when the forward map is not readily differentiable and is given as a black box, and also for high dimensional parameter spaces because of the need to propagate large covariance matrices. Application of ensemble Kalman filters, for example use of the ensemble Kalman inversion (EKI) algorithm, has emerged as a useful tool which overcomes both of these issues: it is derivative free and works with a low-rank covariance approximation formed from the ensemble. In this paper, we work with the ExKI, EKI, and a variant on EKI which we term unscented Kalman inversion (UKI). The paper contains two main contributions. Firstly, we identify a novel stochastic dynamical system in which the parameter-to-data map is embedded. We present theory in the linear case to show exponential convergence of the mean of the filtering distribution to the solution of a regularized least squares problem. This is in contrast to previous work in which the EKI has been employed where the dynamical system used leads to algebraic convergence to an unregularized problem. Secondly, we show that the application of the UKI to this novel stochastic dynamical system yields improved inversion results, in comparison with the application of EKI to the same novel stochastic dynamical system.

preprint2021arXiv

Large-eddy simulations with ClimateMachine: a new open-source code for atmospheric simulations on GPUs and CPUs

We introduce ClimateMachine, a new open-source atmosphere modeling framework using the Julia language to be performance portable on central processing units (CPUs) and graphics processing units (GPUs). ClimateMachine uses a common framework both for coarser-resolution global simulations and for high-resolution, limited-area large-eddy simulations (LES). Here, we demonstrate the LES configuration of the atmosphere model in canonical benchmark cases and atmospheric flows, using an energy-conserving nodal discontinuous-Galerkin (DG) discretization of the governing equations. Resolution dependence, conservation characteristics and scaling metrics are examined in comparison with existing LES codes. They demonstrate the utility of ClimateMachine as a modelling tool for limited-area LES flow configurations.

preprint2021arXiv

Parameter uncertainty quantification in an idealized GCM with a seasonal cycle

Climate models are generally calibrated manually by comparing selected climate statistics, such as the global top-of-atmosphere energy balance, to observations. The manual tuning only targets a limited subset of observational data and parameters. Bayesian calibration can estimate climate model parameters and their uncertainty using a larger fraction of the available data and automatically exploring the parameter space more broadly. In Bayesian learning, it is natural to exploit the seasonal cycle, which has large amplitude, compared with anthropogenic climate change, in many climate statistics. In this study, we develop methods for the calibration and uncertainty quantification (UQ) of model parameters exploiting the seasonal cycle, and we demonstrate a proof-of-concept with an idealized general circulation model (GCM). Uncertainty quantification is performed using the calibrate-emulate-sample approach, which combines stochastic optimization and machine learning emulation to speed up Bayesian learning. The methods are demonstrated in a perfect-model setting through the calibration and UQ of a convective parameterization in an idealized GCM with a seasonal cycle. Calibration and UQ based on seasonally averaged climate statistics, compared to annually averaged, reduces the calibration error by up to an order of magnitude and narrows the spread of posterior distributions by factors between two and five, depending on the variables used for UQ. The reduction in the size of the parameter posterior distributions leads to a reduction in the uncertainty of climate model predictions.

preprint2016arXiv

Cumulant expansions for atmospheric flows

The equations governing atmospheric flows are nonlinear. Consequently, the hierarchy of cumulant equations is not closed. But because atmospheric flows are inhomogeneous and anisotropic, the nonlinearity may manifest itself only weakly through interactions of mean fields with disturbances such as thermals or eddies. In such situations, truncations of the hierarchy of cumulant equations hold promise as a closure strategy. We review how truncations at second order can be used to model and elucidate the dynamics of atmospheric flows. Two examples are considered. First, we study the growth of a dry convective boundary layer, which is heated from below, leading to turbulent upward energy transport and growth of the boundary layer. We demonstrate that a quasilinear truncation of the equations of motion, in which interactions of disturbances among each other are neglected but interactions with mean fields are taken into account, can capture the growth of the convective boundary layer even if it does not capture important turbulent transport terms. Second, we study the evolution of two-dimensional large-scale waves representing waves in Earth's upper atmosphere. We demonstrate that a cumulant expansion truncated at second order (CE2) can capture the evolution of such waves and their nonlinear interaction with the mean flow in some circumstances, for example, when the wave amplitude is small enough or the planetary rotation rate is large enough. However, CE2 fails to capture the flow evolution when nonlinear eddy--eddy interactions in surf zones become important. Higher-order closures can capture these missing interactions. The results point to new ways in which the dynamics of turbulent boundary layers may be represented in climate models, and they illustrate different classes of nonlinear processes that can control wave dissipation and momentum fluxes in the troposphere.

preprint2015arXiv

Why Eddy Momentum Fluxes are Concentrated in the Upper Troposphere

The extratropical eddy momentum flux (EMF) is controlled by generation, propagation, and dissipation of large-scale eddies and is concentrated in Earth's upper troposphere. An idealized GCM is used to investigate how this EMF structure arises. In simulations in which the poles are heated more strongly than the equator, EMF is concentrated near the surface, demonstrating that surface drag generally is not responsible for the upper-tropospheric EMF concentration. Although Earth's upper troposphere favors linear wave propagation, quasi-linear simulations in which nonlinear eddy-eddy interactions are suppressed demonstrate that this is likewise not primarily responsible for the upper-tropospheric EMF concentration. The quasi-linear simulations reveal the essential role of nonlinear eddy-eddy interactions in the surf zone in the upper troposphere, where wave activity absorption away from the baroclinic generation regions occurs through the nonlinear generation of small scales. In Earth-like atmospheres, wave activity that is generated in the lower troposphere propagates upward, then turns meridionally, eventually to be absorbed nonlinearly in the upper troposphere. The level at which the wave activity begins to propagate meridionally appears to be set by the typical height reached by baroclinic eddies. This can coincide with the tropopause height but also can lie below it if convection controls the tropopause height. In the latter case, EMF is maximal well below the tropopause. The simulations suggest that EMF is concentrated in Earth's upper troposphere because typical baroclinic eddies reach the tropopause.

preprint2011arXiv

Changes in zonal surface temperature gradients and Walker circulations in a wide range of climates

Variations in zonal surface temperature gradients and zonally asymmetric tropical overturning circulations (Walker circulations) are examined over a wide range of climates simulated with an idealized atmospheric general circulation model (GCM). The asymmetry in the tropical climate is generated by an imposed ocean energy flux, which does not vary with climate. The range of climates is simulated by modifying the optical thickness of an idealized longwave absorber (representing greenhouse gases). The zonal surface temperature gradient in low latitudes generally decreases as the climate warms in the idealized GCM simulations. A scaling relationship based on a two-term balance in the surface energy budget accounts for the changes in the zonally asymmetric component of the GCM-simulated surface temperature gradients. The Walker circulation weakens as the climate warms in the idealized simulations, as it does in comprehensive simulations of climate change. The wide range of climates allows a systematic test of energetic arguments that have been proposed to account for these changes in the tropical circulation. The analysis shows that a scaling estimate based on changes in the hydrological cycle (precipitation rate and saturation specific humidity) accounts for the simulated changes in the Walker circulation. However, it must be evaluated locally, with local precipitation rates. If global-mean quantities are used, the scaling estimate does not generally account for changes in the Walker circulation, and the extent to which it does is the result of compensating errors in changes in precipitation and saturation specific humidity that enter the scaling estimate.

preprint2010arXiv

Atmospheric dynamics of Earth-like tidally locked aquaplanets

We present simulations of atmospheres of Earth-like aquaplanets that are tidally locked to their star, that is, planets whose orbital period is equal to the rotation period about their spin axis, so that one side always faces the star and the other side is always dark. As extreme cases illustrating the effects of slow and rapid rotation, we consider planets with rotation periods equal to one current Earth year and one current Earth day. The dynamics responsible for the surface climate (e.g., winds, temperature, precipitation) and the general circulation of the atmosphere are discussed in light of existing theories of atmospheric circulations. For example, as expected from the increasing importance of Coriolis accelerations relative to inertial accelerations as the rotation rate increases, the winds are approximately isotropic and divergent at leading order in the slowly rotating atmosphere but are predominantly zonal and rotational in the rapidly rotating atmosphere. Free-atmospheric horizontal temperature variations in the slowly rotating atmosphere are generally weaker than in the rapidly rotating atmosphere. Interestingly, the surface temperature on the night side of the planets does not fall below ~240 K in either the rapidly or slowly rotating atmosphere; that is, heat transport from the day side to the night side of the planets efficiently reduces temperature contrasts in either case. Rotational waves shape the distribution of winds, temperature, and precipitation in the rapidly rotating atmosphere; in the slowly rotating atmosphere, these distributions are controlled by simpler divergent circulations. The results are of interest in the study of tidally locked terrestrial exoplanets and as illustrations of how planetary rotation and the insolation distribution shape climate.

preprint2010arXiv

Mechanisms of jet formation on the giant planets

The giant planet atmospheres exhibit alternating prograde (eastward) and retrograde (westward) jets of different speeds and widths, with an equatorial jet that is prograde on Jupiter and Saturn and retrograde on Uranus and Neptune. The jets are variously thought to be driven by differential radiative heating of the upper atmosphere or by intrinsic heat fluxes emanating from the deep interior. But existing models cannot account for the different flow configurations on the giant planets in an energetically consistent manner. Here a three-dimensional general circulation model is used to show that the different flow configurations can be reproduced by mechanisms universal across the giant planets if differences in their radiative heating and intrinsic heat fluxes are taken into account. Whether the equatorial jet is prograde or retrograde depends on whether the deep intrinsic heat fluxes are strong enough that convection penetrates into the upper troposphere and generates strong equatorial Rossby waves there. Prograde equatorial jets result if convective Rossby wave generation is strong and low-latitude angular momentum flux divergence owing to baroclinic eddies generated off the equator is sufficiently weak (Jupiter and Saturn). Retrograde equatorial jets result if either convective Rossby wave generation is weak or absent (Uranus) or low-latitude angular momentum flux divergence owing to baroclinic eddies is sufficiently strong (Neptune). The different speeds and widths of the off-equatorial jets depend, among other factors, on the differential radiative heating of the atmosphere and the altitude of the jets, which are vertically sheared. The simulations have closed energy and angular momentum balances that are consistent with observations of the giant planets.

preprint2009arXiv

Water vapor and the dynamics of climate changes

Water vapor is not only Earth's dominant greenhouse gas. Through the release of latent heat when it condenses, it also plays an active role in dynamic processes that shape the global circulation of the atmosphere and thus climate. Here we present an overview of how latent heat release affects atmosphere dynamics in a broad range of climates, ranging from extremely cold to extremely warm. Contrary to widely held beliefs, atmospheric circulation statistics can change non-monotonically with global-mean surface temperature, in part because of dynamic effects of water vapor. For example, the strengths of the tropical Hadley circulation and of zonally asymmetric tropical circulations, as well as the kinetic energy of extratropical baroclinic eddies, can be lower than they presently are both in much warmer climates and in much colder climates. We discuss how latent heat release is implicated in such circulation changes, particularly through its effect on the atmospheric static stability, and we illustrate the circulation changes through simulations with an idealized general circulation model. This allows us to explore a continuum of climates, constrain macroscopic laws governing this climatic continuum, and place past and possible future climate changes in a broader context.

preprint2008arXiv

Formation of Jets and Equatorial Superrotation on Jupiter

The zonal flow in Jupiter's upper troposphere is organized into alternating retrograde and prograde jets, with a prograde (superrotating) jet at the equator. Existing models posit as the driver of the flow either differential radiative heating of the atmosphere or intrinsic heat fluxes emanating from the deep interior; however, they do not reproduce all large-scale features of Jupiter's jets and thermal structure. Here it is shown that the difficulties in accounting for Jupiter's jets and thermal structure resolve if the effects of differential radiative heating and intrinsic heat fluxes are considered together, and if upper-tropospheric dynamics are linked to a magnetohydrodynamic (MHD) drag that acts deep in the atmosphere. Baroclinic eddies generated by differential radiative heating can account for the off-equatorial jets; meridionally propagating equatorial Rossby waves generated by intrinsic convective heat fluxes can account for the equatorial superrotation. The zonal flow extends deeply into the atmosphere, with its speed changing with depth, up to depths at which the MHD drag acts. The theory is supported by simulations with an energetically consistent general circulation model of Jupiter's outer atmosphere. A simulation that incorporates differential radiative heating and intrinsic heat fluxes reproduces Jupiter's observed jets and thermal structure and makes testable predictions about as-yet unobserved aspects thereof. A control simulation that incorporates only differential radiative heating but not intrinsic heat fluxes produces off-equatorial jets but no equatorial superrotation; another control simulation that incorporates only intrinsic heat fluxes but not differential radiative heating produces equatorial superrotation but no off-equatorial jets.

preprint2007arXiv

Statistics of an Unstable Barotropic Jet from a Cumulant Expansion

Low-order equal-time statistics of a barotropic flow on a rotating sphere are investigated. The flow is driven by linear relaxation toward an unstable zonal jet. For relatively short relaxation times, the flow is dominated by critical-layer waves. For sufficiently long relaxation times, the flow is turbulent. Statistics obtained from a second-order cumulant expansion are compared to those accumulated in direct numerical simulations, revealing the strengths and limitations of the expansion for different relaxation times.