Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConstraining stochastic parametrisation schemes using high-resolution simulations
54
0
0.0
(
0
)
Authors
Hannah M. Christensen
Ask ChatGPT about the research
No Arabic abstract
Stochastic parametrisations are used in weather and climate models to improve the representation of unpredictable unresolved processes. When compared to a deterministic model, a stochastic model represents `model uncertainty, i.e., sources of error in the forecast due to the limitations of the forecast model. We present a technique for systematically deriving new stochastic parametrisations or for constraining existing stochastic approaches. A high-resolution model simulation is coarse-grained to the desired forecast model resolution. This provides the initial conditions and forcing data needed to drive a Single Column Model (SCM). By comparing the SCM parametrised tendencies with the evolution of the high resolution model, we can estimate the error in the SCM tendencies that a stochastic parametrisation seeks to represent. We use this approach to assess the physical basis of the widely used Stochastically Perturbed Parametrisation Tendencies (SPPT) scheme. We find justification for the multiplicative nature of SPPT, and for the use of spatio-temporally correlated stochastic perturbations. We find evidence that the stochastic perturbation should be positively skewed, indicating that occasional large-magnitude positive perturbations are physically realistic. However other key assumptions of SPPT are less well justified, including coherency of the stochastic perturbations with height, coherency of the perturbations for different physical parametrisation schemes, and coherency for different prognostic variables. Relaxing these SPPT assumptions allows for an error model that explains a larger fractional variance than traditional SPPT. In particular, we suggest that independently perturbing the tendencies associated with different parametrisation schemes is justifiable, and would improve the realism of the SPPT approach.
rate research
Read More
Observations of tropical convection from precipitation radar and the concurring large-scale atmospheric state at two locations (Darwin and Kwajalein) are used to establish effective stochastic models to parameterise subgrid-scale tropical convective activity. Two approaches are presented which rely on the assumption that tropical convection induces a stationary equilibrium distribution. In the first approach we parameterise convection variables such as convective area fraction as an instantaneous random realisation conditioned on the large-scale vertical velocities according to a probability density function estimated from the observations. In the second approach convection variables are generated in a Markov process conditioned on the large-scale vertical velocity, allowing for non-trivial temporal correlations. Despite the different prevalent atmospheric and oceanic regimes at the two locations, with Kwajalein being exposed to a purely oceanic weather regime and Darwin exhibiting land-sea interaction, we establish that the empirical measure for the convective variables conditioned on large-scale mid-level vertical velocities for the two locations are close. This allows us to train the stochastic models at one location and then generate time series of convective activity at the other location. The proposed stochastic subgrid-scale models adequately reproduce the statistics of the observed convective variables and we discuss how they may be used in future scale-independent mass-flux convection parameterisations.
The daily cycle of heating and cooling of the near-surface ocean may be quite different in a shallow lagoon with a few meters deep seafloor that can be heated directly by the sun. If important, the solar radiation will affect the local benthic communities. To study the physical processes associated with the daily cycle of south-Pacific lagoon Bora Bora, a vertical string of five high-resolution temperature sensors was moored at a 2-m deep site for 3 weeks. Besides the standard ocean warming (approximately during daytime) and cooling (approximately nighttime), the sensors show relatively highest temperature near the lagoon-floor during the warming phase and a weakly stable stratification towards the end of the cooling phase. During the warming phase, highly variable stratification is observed extending into the water column under calm weather and turbid waters, otherwise not. Under trade wind and clear waters, the lowest sensor(s) show(s) consistently higher temperature variability than sensors higher-up with spectral slopes indicative of shear- and/or convective turbulence. During the cooling phase, the lower sensor shows consistently very low variance (non-turbulent), while other sensors show a spectral slope around the buoyancy frequency evidencing weakly stratified waters supporting internal waves. These observations contrast with open-ocean near-surface observations of stable stratification during the warming phase and of turbulent free convection during the cooling phase. Thus, lagoons seem to more resemble the atmosphere than the ocean in daytime thermodynamics and possibly act as a natural solar pond with bottom conductive heating (when salinity compensates for unstable temperature variations).
With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two secondary initial leaders that developed out of stepped leaders.
We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge-Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov--Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov--Poisson simulations on $128^3 times 128^3$ mesh grids have been successfully performed on a massively parallel computer.
We study the density structures of dark matter subhalos for both cold dark matter and self-interacting dark matter models using high-resolution cosmological $N$-body simulations. We quantify subhalos central density at 150 pc from the center of each subhalo at the classical dwarf spheroidal and ultrafaint dwarf scales. By comparing them with observations, we find that the self-interacting scattering cross-section of $sigma/m<3 rm{cm^{2}g^{-1}}$ is favored. Due to the combination of hosts tide and self-interactions, the central density of subhalos with small pericenter shows a noticeable difference between the cold and the self-interacting models, indicating that dwarf satellites with small pericenter are ideal sites to further constrain the nature of dark matter by future large spectroscopic surveys.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
