Scattered power from non-thermal, F-region plasma observed by EISCAT—evidence for coherent echoes?

Three rapid, poleward bursts of plasma flow, observed by the U.K.-POLAR EISCAT experiment, are studied in detail. In all three cases the large ion velocities (> 1 kms−1) are shown to drive the ion velocity distribution into a non-Maxwellian form, identified by the characteristic shape of the observed spectra and the fact that analysis of the spectra with the assumption of a Maxwellian distribution leads to excessive rises in apparent ion temperature, and an anticorrelation of apparent electron and ion temperatures. For all three periods the total scattered power is shown to rise with apparent ion temperature by up to 6 dB more than is expected for an isotropic Maxwellian plasma of constant density and by an even larger factor than that expected for non-thermal plasma. The anomalous increases in power are only observed at the lower altitudes (< 300 km). At greater altitudes the rise in power is roughly consistent with that simulated numerically for homogeneous, anisotropic, non-Maxwellian plasma of constant density, viewed using the U.K.-POLAR aspect angle. The spectra at times of anomalously high power are found to be asymmetric, showing an enhancement near the downward Doppler-shifted ion-acoustic frequency. Although it is not possible to eliminate completely rapid plasma density fluctuations as a cause of these power increases, such effects cannot explain the observed spectra and the correlation of power and apparent ion temperature without an unlikely set of coincidences. The observations are made along a beam direction which is as much as 16.5° from orthogonality with the geomagnetic field. Nevertheless, some form of coherent-like echo contamination of the incoherent scatter spectrum is the most satisfactory explanation of these data.

Model predictions of the occurrence of non-Maxwellian plasmas, and analysis of their effects on EISCAT data

The recent identification of non-thermal plasmas using EISCAT data has been made possible by their occurrence during large, short-lived flow bursts. For steady, yet rapid, ion convection the only available signature is the shape of the spectrum, which is unreliable because it is open to distortion by noise and sampling uncertainty and can be mimicked by other phenomena. Nevertheless, spectral shape does give an indication of the presence of non-thermal plasma, and the characteristic shape has been observed for long periods (of the order of an hour or more) in some experiments. To evaluate this type of event properly one needs to compare it to what would be expected theoretically. Predictions have been made using the coupled thermosphere-ionosphere model developed at University College London and the University of Sheffield to show where and when non-Maxwellian plasmas would be expected in the auroral zone. Geometrical and other factors then govern whether these are detectable by radar. The results are applicable to any incoherent scatter radar in this area, but the work presented here concentrates on predictions with regard to experiments on the EISCAT facility.

Non-Maxwellian ion velocity distributions observed using EISCAT

Recent observations from the EISCAT incoherent scatter radar have revealed bursts of poleward ion flow in the dayside auroral ionosphere which are consistent with the ionospheric signature of flux transfer events at the magnetopause. These bursts frequently contain ion drifts which exceed the neutral thermal speed and, because the neutral thermospheric wind is incapable of responding sufficiently rapidly, toroidal, non-Maxwellian ion velocity distributions are expected. The EISCAT observations are made with high time resolution (15 seconds) and at a large angle to the geomagnetic field (73.5°), allowing the non-Maxwellian nature of the distribution to be observed remotely for the first time. The observed features are also strongly suggestive of a toroidal distribution: characteristic spectral shape, increased scattered power (both consistent with reduced Landau damping and enhanced electric field fluctuations) and excessively high line-of-sight ion temperatures deduced if a Maxwellian distribution is assumed. These remote sensing observations allow the evolution of the distributions to be observed. They are found to be non-Maxwellian whenever the ion drift exceeds the neutral thermal speed, indicating that such distributions can exist over the time scale of the flow burst events (several minutes).

Characterising phosphorus and nitrate inputs to a rural river using high-frequency concentration-flow relationships

The total reactive phosphorus (TRP) and nitrate concentrations of the River Enborne, southern England, were monitored at hourly interval between January 2010 and December 2011. The relationships between these high-frequency nutrient concentration signals and flow were used to infer changes in nutrient source and dynamics through the annual cycle and each individual storm event, by studying hysteresis patterns. TRP concentrations exhibited strong dilution patterns with increasing flow, and predominantly clockwise hysteresis through storm events. Despite the Enborne catchment being relatively rural for southern England, TRP inputs were dominated by constant, non-rain-related inputs from sewage treatment works (STW) for the majority of the year, producing the highest phosphorus concentrations through the spring–summer growing season. At higher river flows, the majority of the TRP load was derived from within-channel remobilisation of phosphorus from the bed sediment, much of which was also derived from STW inputs. Therefore, future phosphorus mitigation measures should focus on STW improvements. Agricultural diffuse TRP inputs were only evident during storms in the May of each year, probably relating to manure application to land. The nitrate concentration–flow relationship produced a series of dilution curves, indicating major inputs from groundwater and to a lesser extent STW. Significant diffuse agricultural inputs with anticlockwise hysteresis trajectories were observed during the first major storms of the winter period. The simultaneous investigation of high-frequency time series data, concentration–flow relationships and hysteresis behaviour through multiple storms for both phosphorus and nitrate offers a simple and innovative approach for providing new insights into nutrient sources and dynamics.

Kraichnan–Leith–Batchelor similarity theory and two-dimensional inverse cascades

We study the scaling properties and Kraichnan–Leith–Batchelor (KLB) theory of forced inverse cascades in generalized two-dimensional (2D) fluids (α-turbulence models) simulated at resolution 8192x8192. We consider α=1 (surface quasigeostrophic flow), α=2 (2D Euler flow) and α=3. The forcing scale is well resolved, a direct cascade is present and there is no large-scale dissipation. Coherent vortices spanning a range of sizes, most larger than the forcing scale, are present for both α=1 and α=2. The active scalar field for α=3 contains comparatively few and small vortices. The energy spectral slopes in the inverse cascade are steeper than the KLB prediction −(7−α)/3 in all three systems. Since we stop the simulations well before the cascades have reached the domain scale, vortex formation and spectral steepening are not due to condensation effects; nor are they caused by large-scale dissipation, which is absent. One- and two-point p.d.f.s, hyperflatness factors and structure functions indicate that the inverse cascades are intermittent and non-Gaussian over much of the inertial range for α=1 and α=2, while the α=3 inverse cascade is much closer to Gaussian and non-intermittent. For α=3 the steep spectrum is close to that associated with enstrophy equipartition. Continuous wavelet analysis shows approximate KLB scaling ℰ(k)∝k−2 (α=1) and ℰ(k)∝k−5/3 (α=2) in the interstitial regions between the coherent vortices. Our results demonstrate that coherent vortex formation (α=1 and α=2) and non-realizability (α=3) cause 2D inverse cascades to deviate from the KLB predictions, but that the flow between the vortices exhibits KLB scaling and non-intermittent statistics for α=1 and α=2.

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.

Revisiting the ‘Venturi effect’ in passage ventilation between two non-parallel buildings

A recent study conducted by Blocken et al. (Numerical study on the existence of the Venturi effect in passages between perpendicular buildings. Journal of Engineering Mechanics, 2008,134: 1021-1028) challenged the popular view of the existence of the ‘Venturi effect’ in building passages as the wind is exposed to an open boundary. The present research extends the work of Blocken et al. (2008a) into a more general setup with the building orientation varying from 0° to 180° using CFD simulations. Our results reveal that the passage flow is mainly determined by the combination of corner streams. It is also shown that converging passages have a higher wind-blocking effect compared to diverging passages, explained by a lower wind speed and higher drag coefficient. Fluxes on the top plane of the passage volume reverse from outflow to inflow in the cases of α=135°, 150° and 165°. A simple mathematical expression to explain the relationship between the flux ratio and the geometric parameters has been developed to aid wind design in an urban neighborhood. In addition, a converging passage with α=15° is recommended for urban wind design in cold and temperate climates since the passage flow changes smoothly and a relatively lower wind speed is expected compared with that where there are no buildings. While for the high-density urban area in (sub)tropical climates such as Hong Kong where there is a desire for more wind, a diverging passage with α=150° is a better choice to promote ventilation at the pedestrian level.

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.

Comparison of terrain following and cut cell grids using a non-hydrostatic model

Terrain following coordinates are widely used in operational models but the cut cell method has been proposed as an alternative that can more accurately represent atmospheric dynamics over steep orography. Because the type of grid is usually chosen during model implementation, it becomes necessary to use different models to compare the accuracy of different grids. In contrast, here a C-grid finite volume model enables a like-for-like comparison of terrain following and cut cell grids. A series of standard two-dimensional tests using idealised terrain are performed: tracer advection in a prescribed horizontal velocity field, a test starting from resting initial conditions, and orographically induced gravity waves described by nonhydrostatic dynamics. In addition, three new tests are formulated: a more challenging resting atmosphere case, and two new advection tests having a velocity field that is everywhere tangential to the terrain following coordinate surfaces. These new tests present a challenge on cut cell grids. The results of the advection tests demonstrate that accuracy depends primarily upon alignment of the flow with the grid rather than grid orthogonality. A resting atmosphere is well-maintained on all grids. In the gravity waves test, results on all grids are in good agreement with existing results from the literature, although terrain following velocity fields lead to errors on cut cell grids. Due to semi-implicit timestepping and an upwind-biased, explicit advection scheme, there are no timestep restrictions associated with small cut cells. We do not find the significant advantages of cut cells or smoothed coordinates that other authors find.

The longitudinal variation of equatorial waves due to propagation on a varying zonal flow

The general 1-D theory of waves propagating on a zonally varying flow is developed from basic wave theory, and equations are derived for the variation of wavenumber and energy along ray paths. Different categories of behaviour are found, depending on the sign of the group velocity (cg) and a wave property, B. For B positive the wave energy and the wave number vary in the same sense, with maxima in relative easterlies or westerlies, depending on the sign of cg. Also the wave accumulation of Webster and Chang (1988) occurs where cg goes to zero. However for B negative they behave in opposite senses and wave accumulation does not occur. The zonal propagation of the gravest equatorial waves is analysed in detail using the theory. For non-dispersive Kelvin waves, B reduces to 2, and analytic solution is possible. B is positive for all the waves considered, except for the westward moving mixed Rossby-gravity (WMRG) wave which can have negative as well as positive B. Comparison is made between the observed climatologies of the individual equatorial waves and the result of pure propagation on the climatological upper tropospheric flow. The Kelvin wave distribution is in remarkable agreement, considering the approximations made. Some aspects of the WMRG and Rossby wave distributions are also in qualitative agreement. However the observed maxima in these waves in the winter westerlies in the eastern Pacific and Atlantic are not consistent with the theory. This is consistent with the importance of the sources of equatorial waves in these westerly duct regions due to higher latitude wave activity.

Shape-controlled continuous synthesis of metal nanostructures

A segmented flow-based microreactor is used for the continuous production of faceted nanocrystals. Flow segmentation is proposed as a versatile tool to manipulate the reduction kinetics and control the growth of faceted nanostructures; tuning the size and shape. Switching the gas from oxygen to carbon monoxide permits the adjustment in nanostructure growth from 1D (nanorods) to 2D (nanosheets). CO is a key factor in the formation of Pd nanosheets and Pt nanocubes; operating as a second phase, a reductant, and a capping agent. This combination confines the growth to specific structures. In addition, the segmented flow microfluidic reactor inherently has the ability to operate in a reproducible manner at elevated temperatures and pressures whilst confining potentially toxic reactants, such as CO, in nanoliter slugs. This continuous system successfully synthesised Pd nanorods with an aspect ratio of 6; thin palladium nanosheets with a thickness of 1.5 nm; and Pt nanocubes with a 5.6 nm edge length, all in a synthesis time as low as 150 s.