Turbulent structure and scaling of the inertial subrange in a stratocumulus-topped boundary layer observed by a Doppler lidar

The turbulent structure of a stratocumulus-topped marine boundary layer over a 2-day period is observed with a Doppler lidar at Mace Head in Ireland. Using profiles of vertical velocity statistics, the bulk of the mixing is identified as cloud driven. This is supported by the pertinent feature of negative vertical velocity skewness in the sub-cloud layer which extends, on occasion, almost to the surface. Both coupled and decoupled turbulence characteristics are observed. The length and timescales related to the cloud-driven mixing are investigated and shown to provide additional information about the structure and the source of the mixing inside the boundary layer. They are also shown to place constraints on the length of the sampling periods used to derive products, such as the turbulent dissipation rate, from lidar measurements. For this, the maximum wavelengths that belong to the inertial subrange are studied through spectral analysis of the vertical velocity. The maximum wavelength of the inertial subrange in the cloud-driven layer scales relatively well with the corresponding layer depth during pronounced decoupled structure identified from the vertical velocity skewness. However, on many occasions, combining the analysis of the inertial subrange and vertical velocity statistics suggests higher decoupling height than expected from the skewness profiles. Our results show that investigation of the length scales related to the inertial subrange significantly complements the analysis of the vertical velocity statistics and enables a more confident interpretation of complex boundary layer structures using measurements from a Doppler lidar.

Effects of agricultural management practices on earthworm populations and crop yield: validation and application of a mechanistic modelling approach

There is little consensus on how agriculture will meet future food demands sustainably. Soils and their biota play a crucial role by mediating ecosystem services that support agricultural productivity. However, a multitude of site-specific environmental factors and management practices interact to affect the ability of soil biota to perform vital functions, confounding the interpretation of results from experimental approaches. Insights can be gained through models, which integrate the physiological, biological and ecological mechanisms underpinning soil functions. We present a powerful modelling approach for predicting how agricultural management practices (pesticide applications and tillage) affect soil functioning through earthworm populations. By combining energy budgets and individual-based simulation models, and integrating key behavioural and ecological drivers, we accurately predict population responses to pesticide applications in different climatic conditions. We use the model to analyse the ecological consequences of different weed management practices. Our results demonstrate that an important link between agricultural management (herbicide applications and zero, reduced and conventional tillage) and earthworms is the maintenance of soil organic matter (SOM). We show how zero and reduced tillage practices can increase crop yields while preserving natural ecosystem functions. This demonstrates how management practices which aim to sustain agricultural productivity should account for their effects on earthworm populations, as their proliferation stimulates agricultural productivity. Synthesis and applications. Our results indicate that conventional tillage practices have longer term effects on soil biota than pesticide control, if the pesticide has a short dissipation time. The risk of earthworm populations becoming exposed to toxic pesticides will be reduced under dry soil conditions. Similarly, an increase in soil organic matter could increase the recovery rate of earthworm populations. However, effects are not necessarily additive and the impact of different management practices on earthworms depends on their timing and the prevailing environmental conditions. Our model can be used to determine which combinations of crop management practices and climatic conditions pose least overall risk to earthworm populations. Linking our model mechanistically to crop yield models would aid the optimization of crop management systems by exploring the trade-off between different ecosystem services.

Energy consistent discontinuous Galerkin methods for the Navier–Stokes–Korteweg system

We design consistent discontinuous Galerkin finite element schemes for the approximation of the Euler-Korteweg and the Navier-Stokes-Korteweg systems. We show that the scheme for the Euler-Korteweg system is energy and mass conservative and that the scheme for the Navier-Stokes-Korteweg system is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to viscous effects, that is, there is no numerical dissipation. In this sense the methods are consistent with the energy dissipation of the continuous PDE systems. - See more at: http://www.ams.org/journals/mcom/2014-83-289/S0025-5718-2014-02792-0/home.html#sthash.rwTIhNWi.dpuf

Energy consistent discontinuous Galerkin methods for a quasi-incompressible diffuse two phase flow model

We design consistent discontinuous Galerkin finite element schemes for the approximation of a quasi-incompressible two phase flow model of Allen–Cahn/Cahn–Hilliard/Navier–Stokes–Korteweg type which allows for phase transitions. We show that the scheme is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to discrete equivalents of those effects already causing dissipation on the continuous level, that is, there is no artificial numerical dissipation added into the scheme. In this sense the methods are consistent with the energy dissipation of the continuous PDE system.

Future changes in the Western North Pacific tropical cyclone activity projected by a multidecadal simulation with a 16-km global atmospheric GCM

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Niño–Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.

The solsticial pause on Mars: 2 modelling and investigation of causes

The martian solsticial pause, presented in a companion paper (Lewis et al., this issue), was investigated further through a series of model runs using the UK version of the LMD/UK Mars Global Climate Model. It was found that the pause could not be adequately reproduced if radiatively active water ice clouds were omitted from the model. When clouds were used, along with a realistic time-dependent dust opacity distribution, a substantial minimum in near-surface transient eddy activity formed around solstice in both hemispheres. The net effect of the clouds in the model is, by altering the thermal structure of the atmosphere, to decrease the vertical shear of the westerly jet near the surface around solstice, and thus reduce baroclinic growth rates. A similar effect was seen under conditions of large dust loading, implying that northern midlatitude eddy activity will tend to become suppressed after a period of intense flushing storm formation around the northern cap edge. Suppression of baroclinic eddy generation by the barotropic component of the flow and via diabatic eddy dissipation were also investigated as possible mechanisms leading to the formation of the solsticial pause but were found not to make major contributions. Zonal variations in topography were found to be important, as their presence results in weakened transient eddies around winter solstice in both hemispheres, through modification of the near-surface flow. The zonal topographic asymmetry appears to be the primary reason for the weakness of eddy activity in the southern hemisphere relative to the northern hemisphere, and the ultimate cause of the solsticial pause in both hemispheres. The meridional topographic gradient was found to exert a much weaker influence on near-surface transient eddies.

Thermal aspects of satellite downscaling

The strong trend toward nanosatellites creates new challenges in terms of thermal balance control. The thermal balance of a satellite is determined by the heat dissipation in its subsystems and by the thermal connections between them. As satellites become smaller, heat dissipation in their subsystems tends to decrease and thermal connectivity scales down with dimension. However, these two terms do not necessarily scale in the same way, and so the thermal balance may alter and the temperature of subsystems may reach undesired levels. This paper focuses on low-Earth-orbit satellites. We constructed a generalized lumped thermal model that combines a generalized low-Earth-orbit satellite configuration with scaling trends in subsystem heat dissipation and thermal connectivity. Using satellite mass as a scaling parameter, we show that subsystems do not become thermally critical by scaling mass alone.

The non-conservation of potential vorticity by a dynamical core compared with the effects of parametrized physical processes

Numerical models of the atmosphere combine a dynamical core, which approximates solutions to the adiabatic, frictionless governing equations for fluid dynamics, with tendencies arising from the parametrization of other physical processes. Since potential vorticity (PV) is conserved following fluid flow in adiabatic, frictionless circumstances, it is possible to isolate the effects of non-conservative processes by accumulating PV changes in an air-mass relative framework. This “PV tracer technique” is used to accumulate separately the effects on PV of each of the different non-conservative processes represented in a numerical model of the atmosphere. Dynamical cores are not exactly conservative because they introduce, explicitly or implicitly, some level of dissipation and adjustment of prognostic model variables which acts to modify PV. Here, the PV tracers technique is extended to diagnose the cumulative effect of the non-conservation of PV by a dynamical core and its characteristics relative to the PV modification by parametrized physical processes. Quantification using the Met Office Unified Model reveals that the magnitude of the non-conservation of PV by the dynamical core is comparable to those from physical processes. Moreover, the residual of the PV budget, when tracing the effects of the dynamical core and physical processes, is at least an order of magnitude smaller than the PV tracers associated with the most active physical processes. The implication of this work is that the non-conservation of PV by a dynamical core can be assessed in case studies with a full suite of physics parametrizations and directly compared with the PV modification by parametrized physical processes. The nonconservation of PV by the dynamical core is shown to move the position of the extratropical tropopause while the parametrized physical processes have a lesser effect at the tropopause level.

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.