An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Regulating vertical relations in the presence of retailer differentiation costs

We discuss public policy towards vertical relations, comparing different types of contracts between a manufacturer and a maximum of two retailers. Together with (potential) price competition between the retailers, we study the role of a (sunk) differentiation cost paid by them in order to relax competition in the retail market and broaden the market potential of the distributed product. This non-price competition element in the downstream market is responsible for our conclusion that, unlike in standard policy guidelines and previous theoretical analysis, restrictions in intra-brand competition may deserve a permissive treatment even in the absence of inter-brand competition, if retailer differentiation is costly.

Horizontal and vertical structure of the Eyjafjallajökull ash cloud over the UK: a comparison of airborne lidar observations and simulations

During April and May 2010 the ash cloud from the eruption of the Icelandic volcano Eyjafjallajökull caused widespread disruption to aviation over northern Europe. The location and impact of the eruption led to a wealth of observations of the ash cloud were being obtained which can be used to assess modelling of the long range transport of ash in the troposphere. The UK FAAM (Facility for Airborne Atmospheric Measurements) BAe-146-301 research aircraft overflew the ash cloud on a number of days during May. The aircraft carries a downward looking lidar which detected the ash layer through the backscatter of the laser light. In this study ash concentrations derived from the lidar are compared with simulations of the ash cloud made with NAME (Numerical Atmospheric-dispersion Modelling Environment), a general purpose atmospheric transport and dispersion model. The simulated ash clouds are compared to the lidar data to determine how well NAME simulates the horizontal and vertical structure of the ash clouds. Comparison between the ash concentrations derived from the lidar and those from NAME is used to define the fraction of ash emitted in the eruption that is transported over long distances compared to the total emission of tephra. In making these comparisons possible position errors in the simulated ash clouds are identified and accounted for. The ash layers seen by the lidar considered in this study were thin, with typical depths of 550–750 m. The vertical structure of the ash cloud simulated by NAME was generally consistent with the observed ash layers, although the layers in the simulated ash clouds that are identified with observed ash layers are about twice the depth of the observed layers. The structure of the simulated ash clouds were sensitive to the profile of ash emissions that was assumed. In terms of horizontal and vertical structure the best results were obtained by assuming that the emission occurred at the top of the eruption plume, consistent with the observed structure of eruption plumes. However, early in the period when the intensity of the eruption was low, assuming that the emission of ash was uniform with height gives better guidance on the horizontal and vertical structure of the ash cloud. Comparison of the lidar concentrations with those from NAME show that 2–5% of the total mass erupted by the volcano remained in the ash cloud over the United Kingdom.

Unsteady flow simulation of a vertical axis augmented wind turbine: a two-dimensional study

As the integration of vertical axis wind turbines in the built environment is a promising alternative to horizontal axis wind turbines, a 2D computational investigation of an augmented wind turbine is proposed and analysed. In the initial CFD analysis, three parameters are carefully investigated: mesh resolution; turbulence model; and time step size. It appears that the mesh resolution and the turbulence model affect result accuracy; while the time step size examined, for the unsteady nature of the flow, has small impact on the numerical results. In the CFD validation of the open rotor with secondary data, the numerical results are in good agreement in terms of shape. It is, however, observed a discrepancy factor of 2 between numerical and experimental data. Successively, the introduction of an omnidirectional stator around the wind turbine increases the power and torque coefficients by around 30–35% when compared to the open case; but attention needs to be given to the orientation of the stator blades for optimum performance. It is found that the power and torque coefficients of the augmented wind turbine are independent of the incident wind speed considered.

Vertical profile measurements of lower troposphere ionisation

Vertical soundings of the atmospheric ion production rate have been obtained from Geiger counters integrated with conventional meteorological radiosondes. In launches made from Reading (UK) during 2013-2014, the Regener-Pfotzer ionisation maximum was at an altitude equivalent to a pressure of (63.1±2.4) hPa, or, expressed in terms of the local air density, (0.101±0.005) kgm−3. The measured ionisation profiles have been evaluated against the Usoskin-Kovaltsov model and, separately, surface neutron monitor data from Oulu. Model ionisation rates agree well with the observed cosmic ray ionisation below 20 km altitude. Above 10 km, the measured ionisation rates also correlate well with simultaneous neutron monitor data, although, consistently with previous work, measured variability at the ionisation maximum is greater than that found by the neutron monitor. However, in the lower atmosphere (below 5 km altitude), agreement between the measurements and simultaneous neutron monitor data is poor. For studies of transient lower atmosphere phenomena associated with cosmic ray ionisation, this indicates the need for in situ ionisation measurements and improved lower atmosphere parameterisations.

Bio-inspired autonomous visual vertical control of a quadrotor UAV

Near ground maneuvers, such as hover, approach and landing, are key elements of autonomy in unmanned aerial vehicles. Such maneuvers have been tackled conventionally by measuring or estimating the velocity and the height above the ground often using ultrasonic or laser range finders. Near ground maneuvers are naturally mastered by flying birds and insects as objects below may be of interest for food or shelter. These animals perform such maneuvers efficiently using only the available vision and vestibular sensory information. In this paper, the time-to-contact (Tau) theory, which conceptualizes the visual strategy with which many species are believed to approach objects, is presented as a solution for Unmanned Aerial Vehicles (UAV) relative ground distance control. The paper shows how such an approach can be visually guided without knowledge of height and velocity relative to the ground. A control scheme that implements the Tau strategy is developed employing only visual information from a monocular camera and an inertial measurement unit. To achieve reliable visual information at a high rate, a novel filtering system is proposed to complement the control system. The proposed system is implemented on-board an experimental quadrotor UAV and shown not only to successfully land and approach ground, but also to enable the user to choose the dynamic characteristics of the approach. The methods presented in this paper are applicable to both aerial and space autonomous vehicles.

Sea ice–ocean coupling using a rescaled vertical coordinate z*

Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate ‘‘z*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of "z*" avoids the difficult issue of vanishing levels under thick ice. Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model.

Vertical structure of stratospheric water vapour trends derived from merged satellite data

Stratospheric water vapour is a powerful greenhouse gas. The longest available record from balloon observations over Boulder, Colorado, USA shows increases in stratospheric water vapour concentrations that cannot be fully explained by observed changes in the main drivers, tropical tropopause temperatures and methane. Satellite observations could help resolve the issue, but constructing a reliable long-term data record from individual short satellite records is challenging. Here we present an approach to merge satellite data sets with the help of a chemistry–climate model nudged to observed meteorology. We use the models’ water vapour as a transfer function between data sets that overcomes issues arising from instrument drift and short overlap periods. In the lower stratosphere, our water vapour record extends back to 1988 and water vapour concentrations largely follow tropical tropopause temperatures. Lower and mid-stratospheric long-term trends are negative, and the trends from Boulder are shown not to be globally representative. In the upper stratosphere, our record extends back to 1986 and shows positive long-term trends. The altitudinal differences in the trends are explained by methane oxidation together with a strengthened lower-stratospheric and a weakened upper stratospheric circulation inferred by this analysis. Our results call into question previous estimates of surface radiative forcing based on presumed global long-term increases in water vapour concentrations in the lower stratosphere.

Resonance response of electrorheological fluids in vertical oscillation squeeze flow

The resonance effect of microcrystalline cellulose/castor oil electrorheological (ER) suspensions was studied in a compressed oscillatory squeeze flow under external electric fields. The resonance frequency first increases linearly with increasing external held, and then shift to high-field plateau. The amplitudes of resonance peak increase sharply with the applied fields in the range of 0.17-1.67kV/mm. The phase difference of the reduced displacement relative to the excitation force inverses in the case of resonance. A viscoelasticity model of the ER suspensions, which offers both the equivalent stiffness and the viscous damping, should be responsible for the appearance of resonance. The influence of the electric field on the resonance frequency and the resonance hump is consistent qualitatively with the interpretation of our proposed model. Storage modulus G' was presented for the purpose of investigating this influence.

Ocean heat uptake processes: a model intercomparison

We compare the quasi-equilibrium heat balances, as well as their responses to 4×CO2 perturbation, among three global climate models with the aim to identify and explain inter-model differences in ocean heat uptake (OHU) processes. We find that, in quasi-equilibrium, convective and mixed layer processes, as well as eddy-related processes, cause cooling of the subsurface ocean. The cooling is balanced by warming caused by advective and diapycnally diffusive processes. We also find that in the CO2-perturbed climates the largest contribution to OHU comes from changes in vertical mixing processes and the mean circulation, particularly in the extra-tropics, caused both by changes in wind forcing, and by changes in high-latitude buoyancy forcing. There is a substantial warming in the tropics, a significant part of which occurs because of changes in horizontal advection in extra-tropics. Diapycnal diffusion makes only a weak contribution to the OHU, mainly in the tropics, due to increased stratification. There are important qualitative differences in the contribution of eddy-induced advection and isopycnal diffusion to the OHU among the models. The former is related to the different values of the coefficients used in the corresponding scheme. The latter is related to the different tapering formulations of the isopycnal diffusion scheme. These differences affect the OHU in the deep ocean, which is substantial in two of the models, the dominant region of deep warming being the Southern Ocean. However, most of the OHU takes place above 2000 m, and the three models are quantitatively similar in their global OHU efficiency and its breakdown among processes and as a function of latitude.

The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Vertical ion flow in the polar ionosphere

Outflowing ions from the polar ionosphere fall into two categories: the classical polar wind and the suprathermal ion flows. The flows in both these categories vary a great deal with altitude. The classical polar wind is supersonic at high altitude: at ∼3 RE geocentric, the observed polar wind is H+ dominated and has a Mach number of 2.5–5.1. At 400–600 km, thermal and suprathermal upward O+ ion fluxes frequently occur at the poleward edge of the nightside auroral oval during magnetically active times. Above 500 km, ions are accelerated transverse to the local geomagnetic field. At 1400 km, transversely accelerated ions are frequently observed in winter nights but rarely appear in the summer. In the dayside cleft above ∼2000 km, ions of all species are transversely heated and upwell with significant number and heat fluxes, forming a cleft ion fountain as they convect across the polar cap. Upwelling ions are observed most (least) frequently in the summer (winter). At yet higher altitudes, energetic (>10 eV to several kiloelectron volts) upflowing H+ and O+ ions are frequently observed, their active time occurrence frequency being as high as 0.7 at auroral latitudes and 0.3 in the polar cap. Their composition, intensity, and angular characteristics vary quantitatively with solar activity, being O+ dominant and more intense near solar maximum. Their resulting ion outflow is dominated by ions below 1 keV and reaches 3.5×10^26 O+ and 7×10^25 H+ ions s^{−1} at magnetically active times (Kp≥5) near solar maximum. In comparison, the estimated polar wind ion outflow at times of moderate solar activity is 7×10^25H+ and 4×10^24 He+ ions s^{−1}. The estimated <10-eV cleft ion fountain flow is 3.8×10^25 O+ and 8.6×10^23 H+ ions s^{−1} near solar maximum.

Bioinspired autonomous visual vertical control of a quadrotor unmanned aerial vehicle

Near-ground maneuvers, such as hover, approach, and landing, are key elements of autonomy in unmanned aerial vehicles. Such maneuvers have been tackled conventionally by measuring or estimating the velocity and the height above the ground, often using ultrasonic or laser range finders. Near-ground maneuvers are naturally mastered by flying birds and insects because objects below may be of interest for food or shelter. These animals perform such maneuvers efficiently using only the available vision and vestibular sensory information. In this paper, the time-tocontact (tau) theory, which conceptualizes the visual strategy with which many species are believed to approach objects, is presented as a solution for relative ground distance control for unmanned aerial vehicles. The paper shows how such an approach can be visually guided without knowledge of height and velocity relative to the ground. A control scheme that implements the tau strategy is developed employing only visual information from a monocular camera and an inertial measurement unit. To achieve reliable visual information at a high rate, a novel filtering system is proposed to complement the control system. The proposed system is implemented onboard an experimental quadrotor unmannedaerial vehicle and is shown to not only successfully land and approach ground, but also to enable the user to choose the dynamic characteristics of the approach. The methods presented in this paper are applicable to both aerial and space autonomous vehicles.

A traceable physical calibration of the vertical advection-diffusion equation for modelling ocean heat uptake

The classic vertical advection-diffusion (VAD) balance is a central concept in studying the ocean heat budget, in particular in simple climate models (SCMs). Here we present a new framework to calibrate the parameters of the VAD equation to the vertical ocean heat balance of two fully-coupled climate models that is traceable to the models’ circulation as well as to vertical mixing and diffusion processes. Based on temperature diagnostics, we derive an effective vertical velocity w∗ and turbulent diffusivity k∗ for each individual physical process. In steady-state, we find that the residual vertical velocity and diffusivity change sign in mid-depth, highlighting the different regional contributions of isopycnal and diapycnal diffusion in balancing the models’ residual advection and vertical mixing. We quantify the impacts of the time-evolution of the effective quantities under a transient 1%CO2 simulation and make the link to the parameters of currently employed SCMs.

Vertical structure and physical processes of the Madden-Julian oscillation: Linking hindcast fidelity to simulated diabatic heating and moistening

Many theories for the Madden-Julian oscillation (MJO) focus on diabatic processes, particularly the evolution of vertical heating and moistening. Poor MJO performance in weather and climate models is often blamed on biases in these processes and their interactions with the large-scale circulation. We introduce one of three components of a model-evaluation project, which aims to connect MJO fidelity in models to their representations of several physical processes, focusing on diabatic heating and moistening. This component consists of 20-day hindcasts, initialised daily during two MJO events in winter 2009-10. The 13 models exhibit a range of skill: several have accurate forecasts to 20 days' lead, while others perform similarly to statistical models (8-11 days). Models that maintain the observed MJO amplitude accurately predict propagation, but not vice versa. We find no link between hindcast fidelity and the precipitation-moisture relationship, in contrast to other recent studies. There is also no relationship between models' performance and the evolution of their diabatic-heating profiles with rain rate. A more robust association emerges between models' fidelity and net moistening: the highest-skill models show a clear transition from low-level moistening for light rainfall to mid-level moistening at moderate rainfall and upper-level moistening for heavy rainfall. The mid-level moistening, arising from both dynamics and physics, may be most important. Accurately representing many processes may be necessary, but not sufficient for capturing the MJO, which suggests that models fail to predict the MJO for a broad range of reasons and limits the possibility of finding a panacea.