A review of urban roughness sublayer turbulence

It is becoming increasingly important that we can understand and model flow processes in urban areas. Applications such as weather forecasting, air quality and sustainable urban development rely on accurate modelling of the interface between an urban surface and the atmosphere above. This review gives an overview of current understanding of turbulence generated by an urban surface up to a few building heights, the layer called the roughness sublayer (RSL). High quality datasets are also identified which can be used in the development of suitable parameterisations of the urban RSL. Datasets derived from physical and numerical modelling, and full-scale observations in urban areas now exist across a range of urban-type morphologies (e.g. street canyons, cubes, idealised and realistic building layouts). Results show that the urban RSL depth falls within 2 – 5 times mean building height and is not easily related to morphology. Systematic perturbations away from uniform layouts (e.g. varying building heights) have a significant impact on RSL structure and depth. Considerable fetch is required to develop an overlying inertial sublayer, where turbulence is more homogeneous, and some authors have suggested that the “patchiness” of urban areas may prevent inertial sublayers from developing at all. Turbulence statistics suggest similarities between vegetation and urban canopies but key differences are emerging. There is no consensus as to suitable scaling variables, e.g. friction velocity above canopy vs. square root of maximum Reynolds stress, mean vs. maximum building height. The review includes a summary of existing modelling practices and highlights research priorities.

The first real-time worldwide ionospheric predictions network: an advance in support of spaceborne experimentation, on-line model validation, and space weather

We report on the first realtime ionospheric predictions network and its capabilities to ingest a global database and forecast F-layer characteristics and "in situ" electron densities along the track of an orbiting spacecraft. A global network of ionosonde stations reported around-the-clock observations of F-region heights and densities, and an on-line library of models provided forecasting capabilities. Each model was tested against the incoming data; relative accuracies were intercompared to determine the best overall fit to the prevailing conditions; and the best-fit model was used to predict ionospheric conditions on an orbit-to-orbit basis for the 12-hour period following a twice-daily model test and validation procedure. It was found that the best-fit model often provided averaged (i.e., climatologically-based) accuracies better than 5% in predicting the heights and critical frequencies of the F-region peaks in the latitudinal domain of the TSS-1R flight path. There was a sharp contrast however, in model-measurement comparisons involving predictions of actual, unaveraged, along-track densities at the 295 km orbital altitude of TSS-1R In this case, extrema in the first-principle models varied by as much as an order of magnitude in density predictions, and the best-fit models were found to disagree with the "in situ" observations of Ne by as much as 140%. The discrepancies are interpreted as a manifestation of difficulties in accurately and self-consistently modeling the external controls of solar and magnetospheric inputs and the spatial and temporal variabilities in electric fields, thermospheric winds, plasmaspheric fluxes, and chemistry.

Large plasma velocities along the magnetic field line in the auroral zone

In the auroral zone, ionospheric plasma often moves horizontally at more than 1 km s−1, driven by magnetospheric electric fields, but it has usually been assumed that vertical velocities are much smaller. On occasions, however, plasma has been seen to move upwards along the magnetic field line at several hundred m s−1. These upward velocities are associated with electric fields strong enough to heat the ion population and drive it into a non-thermal state1,2. Here we report observations of substantial upwards acceleration of plasma, to velocities as high as 500 m s−1. An initial upthrust was provided by a combined upwelling of the neutral atmosphere and ionosphere but the continued acceleration at greater heights is explained by a combination of enhanced plasma pressure and the 'hydrodynamic mirror force'3. This acceleration marks an important stage in the transport of plasma from the ionosphere into the magnetosphere.

Thermal ion flows in the topside auroral ionosphere and the effects of low-altitude, transverse acceleration

Topside ionospheric profiles are used to study the upward field-aligned flow of thermal O+ at high latitudes. On the majority of the field lines outside the plasmasphere, the mean flux is approximately equal to the mean polar wind measured by spacecraft at greater altitudes. This is consistent with the theory of thermal light ion escape supported, via charge exchange, by upward O+ flow at lower heights. Events of larger O+ flow are detected at auroral latitudes and their occurrence is found to agree with that of transversely accelerated ions within the topside ionosphere and the magnetosphere. The effects of low altitude heating of O+ by oxygen cyclotron waves, driven by downward field-aligned currents, are considered as a possible common cause of these two types of event.

Departures from diffusive equilibrium in the topside F-layer from satellite soundings

A method for the detection of O+ ion fluxes from topside soundings is described. The shape of the plasma scale-height profile is altered by such flows only at heights near the F2-peak, where ion-neutral drag is large. Model profiles are used to relate changes in scale height to the ratio (φ/φL) where φ is the field-aligned O+ flux (relative to the neutral air) and φL is the limiting value set by frictional drag. Values of (φ/φL) can then be determined to within a few per cent from experimental soundings, using the plasma temperature and its gradient (as deduced from the observed profile) and the MSIS model neutral temperature. It was found that 3700 topside profiles show departures from diffusive equilibrium, out of 10,000 used to obtain the global morphology of (φ/φL) near the sunspot minimum. Results reveal dynamic ion-flow effects such as the transequatorial breeze and the effects of the polar wind and protonospheric replenishment light-ion flows can be inferred.

Ionospheric origin of magnetospheric O+ ions

Flows of thermal atomic oxygen (O+) ions are deduced from topside ionospheric plasma density profiles. The mean flux within most of the polar cap is of the order of 10^12 m^{−2} s^{−1}, a figure which is consistent with both theoretical and experimental estimates of the light ion polar wind at greater altitudes. Larger flows (up to 6 × 10^12 m^{−2} s^{−1}) are observed near the poleward edge of the night-side statistical auroral oval, a feature not reproduced in the light ion flux. The implication is one of a low altitude acceleration mechanism, acting upon the O+ ions at these latitudes and at heights above that at which the fluxes are observed. Such a process would enable ions to escape from the ionosphere because they do not exchange charge with neutral hydrogen. The observations are in general agreement with energetic O+ ions as previously observed in various parts of the magnetosphere.

The plume also rises: trajectories of pheromone plumes issuing from point sources in an orchard canopy at night

The trajectories of pheromone plumes in canopied habitats, such as orchards, have been little studied. We documented the capture of male navel orangeworm moths, Amyelois transitella, in female-baited traps positioned at 5 levels, from ground level to the canopy top, at approximately 6 m above ground, in almond orchards. Males were captured in similar proportions at all levels, suggesting that they do not favor a particular height during ranging flight. A 3-D sonic anemometer was used to establish patterns of wind flow and temperature at 6 heights from 2.08 to 6.65 m in an almond orchard with a 5 m high canopy, every 3 h over 72 h. The horizontal velocity of wind flow was highest above the canopy, where its directionality also was the most consistent. During the time of A. transitella mating (0300–0600), there was a net vertical displacement upward. Vertical buoyancy combined with only minor reductions in the distance that plumes will travel in the lower compared to the upper canopy suggest that the optimal height for release of pheromone from high-release-rate sources, such as aerosol dispensers (“puffers”), that are deployed at low densities (e.g., 3 per ha.) would be at mid or low in the canopy, thereby facilitating dispersion of disruptant throughout the canopy. Optimal placement of aerosol dispensers will vary with the behavioral ecology of the target pest; however, our results suggest that current protocols, which generally propose dispenser placement in the upper third of the canopy, should be reevaluated.

An overview on the urban boundary-layer atmosphere network in Helsinki

The Helsinki Urban Boundary-Layer Atmosphere Network (UrBAN: http://urban.fmi.fi) is a dedicated research-grade observational network where the physical processes in the atmosphere above the city are studied. Helsinki UrBAN is the most poleward intensive urban research observation network in the world and thus will allow studying some unique features such as strong seasonality. The network's key purpose is for the understanding of the physical processes in the urban boundary layer and associated fluxes of heat, momentum, moisture, and other gases. A further purpose is to secure a research-grade database, which can be used internationally to validate and develop numerical models of air quality and weather prediction. Scintillometers, a scanning Doppler lidar, ceilometers, a sodar, eddy-covariance stations, and radiometers are used. This equipment is supplemented by auxiliary measurements, which were primarily set up for general weather and/or air-quality mandatory purposes, such as vertical soundings and the operational Doppler radar network. Examples are presented as a testimony to the potential of the network for urban studies, such as (i) evidence of a stable boundary layer possibly coupled to an urban surface, (ii) the comparison of scintillometer data with sonic anemometry above an urban surface, (iii) the application of scanning lidar over a city, and (iv) combination of sodar and lidar to give a fuller range of sampling heights for boundary layer profiling.

Wave height analysis from 10 years of observations in the Norwegian Sea

Large waves pose risks to ships, offshore structures, coastal infrastructure and ecosystems. This paper analyses 10 years of in-situ measurements of significant wave height (Hs) and maximum wave height (Hmax) from the ocean weather ship Polarfront in the Norwegian Sea. During the period 2000 to 2009, surface elevation was recorded every 0.59 s during sampling periods of 30 min. The Hmax observations scale linearly with Hs on average. A widely-used empirical Weibull distribution is found to estimate average values of Hmax/Hs and Hmax better than a Rayleigh distribution, but tends to underestimate both for all but the smallest waves. In this paper we propose a modified Rayleigh distribution which compensates for the heterogeneity of the observed dataset: the distribution is fitted to the whole dataset and improves the estimate of the largest waves. Over the 10-year period, the Weibull distribution approximates the observed Hs and Hmax well, and an exponential function can be used to predict the probability distribution function of the ratio Hmax/Hs. However, the Weibull distribution tends to underestimate the occurrence of extremely large values of Hs and Hmax. The persistence of Hs and Hmax in winter is also examined. Wave fields with Hs>12 m and Hmax>16 m do not last longer than 3 h. Low-to-moderate wave heights that persist for more than 12 h dominate the relationship of the wave field with the winter NAO index over 2000–2009. In contrast, the inter-annual variability of wave fields with Hs>5.5 m or Hmax>8.5 m and wave fields persisting over ~2.5 days is not associated with the winter NAO index.

Corporate responsibility research: past – present – future

The concept of corporate responsibility (CR) has moved a long way over the past six decades, since Bowen’s book titled Social Responsibilities of the Businessman marked a modern era of business and society research. In this article, we trace the development of business and society literature to provide contextual background to contemporary studies of CR. As well as looking backward, however, we also project forward arguing that fundamental changes in the current business environment are likely to propel CR to new heights in the future. In particular, we explore the mechanisms by which CR affects stakeholder behavior and suggest ways in which scholars and practitioners may want to apply advances from psychology to link CR strategy to human motivation and the notion of balance in organizations and society.

Earth hummocks in north-east Okstindan, northern Norway: morphology, distribution and environmental constraints

Earth hummocks (also termed pounus or thúfur) are a common form of periglacial non-sorted patterned ground. The study objectives were to determine the morphology, distribution and development on slopes of earth hummocks in north-east Okstindan, Norway, an area with many hummocks but few documented accounts. The methodology involved detailed geomorphological mapping and precise measurement with a profileometer. The internal structure of the hummocks was investigated through excavations and sediment sample analyses. Fourteen sites with well-developed earth hummocks (accounting for over 650 individual hummock forms) were investigated. The sites have an average altitude of 750 m and occur on slopes with an average gradient of 7°. The hummock heights are in the range 0.11–0.52 m and their diameters 0.7–1.5 m, although coalescent forms are up to 5 m in length. The hummock morphology is characterised by a variable plan form, asymmetry with respect to upslope and downslope forms, downslope elongation, coalescence, and superimposed microtopography. The hummocks’ distribution appeared to have been controlled by the existence of a frost-susceptible ‘host’ sediment, but moisture availability and topographic position played a role. The authors conclude that differential frost heave and vegetation cover stability are critical for the hummocks’ longevity in the studied landscape.

Improving the TanDEM-X Digital Elevation Model for flood modelling using flood extents from Synthetic Aperture Radar images

The topography of many floodplains in the developed world has now been surveyed with high resolution sensors such as airborne LiDAR (Light Detection and Ranging), giving accurate Digital Elevation Models (DEMs) that facilitate accurate flood inundation modelling. This is not always the case for remote rivers in developing countries. However, the accuracy of DEMs produced for modelling studies on such rivers should be enhanced in the near future by the high resolution TanDEM-X WorldDEM. In a parallel development, increasing use is now being made of flood extents derived from high resolution Synthetic Aperture Radar (SAR) images for calibrating, validating and assimilating observations into flood inundation models in order to improve these. This paper discusses an additional use of SAR flood extents, namely to improve the accuracy of the TanDEM-X DEM in the floodplain covered by the flood extents, thereby permanently improving this DEM for future flood modelling and other studies. The method is based on the fact that for larger rivers the water elevation generally changes only slowly along a reach, so that the boundary of the flood extent (the waterline) can be regarded locally as a quasi-contour. As a result, heights of adjacent pixels along a small section of waterline can be regarded as samples with a common population mean. The height of the central pixel in the section can be replaced with the average of these heights, leading to a more accurate estimate. While this will result in a reduction in the height errors along a waterline, the waterline is a linear feature in a two-dimensional space. However, improvements to the DEM heights between adjacent pairs of waterlines can also be made, because DEM heights enclosed by the higher waterline of a pair must be at least no higher than the corrected heights along the higher waterline, whereas DEM heights not enclosed by the lower waterline must in general be no lower than the corrected heights along the lower waterline. In addition, DEM heights between the higher and lower waterlines can also be assigned smaller errors because of the reduced errors on the corrected waterline heights. The method was tested on a section of the TanDEM-X Intermediate DEM (IDEM) covering an 11km reach of the Warwickshire Avon, England. Flood extents from four COSMO-SKyMed images were available at various stages of a flood in November 2012, and a LiDAR DEM was available for validation. In the area covered by the flood extents, the original IDEM heights had a mean difference from the corresponding LiDAR heights of 0.5 m with a standard deviation of 2.0 m, while the corrected heights had a mean difference of 0.3 m with standard deviation 1.2 m. These figures show that significant reductions in IDEM height bias and error can be made using the method, with the corrected error being only 60% of the original. Even if only a single SAR image obtained near the peak of the flood was used, the corrected error was only 66% of the original. The method should also be capable of improving the final TanDEM-X DEM and other DEMs, and may also be of use with data from the SWOT (Surface Water and Ocean Topography) satellite.

CO2 emission estimation in the urban environment: measurement of the CO2 storage term

Eddy covariance has been used in urban areas to evaluate the net exchange of CO2 between the surface and the atmosphere. Typically, only the vertical flux is measured at a height 2–3 times that of the local roughness elements; however, under conditions of relatively low instability, CO2 may accumulate in the airspace below the measurement height. This can result in inaccurate emissions estimates if the accumulated CO2 drains away or is flushed upwards during thermal expansion of the boundary layer. Some studies apply a single height storage correction; however, this requires the assumption that the response of the CO2 concentration profile to forcing is constant with height. Here a full seasonal cycle (7th June 2012 to 3rd June 2013) of single height CO2 storage data calculated from concentrations measured at 10 Hz by open path gas analyser are compared to a data set calculated from a concurrent switched vertical profile measured (2 Hz, closed path gas analyser) at 10 heights within and above a street canyon in central London. The assumption required for the former storage determination is shown to be invalid. For approximately regular street canyons at least one other measurement is required. Continuous measurements at fewer locations are shown to be preferable to a spatially dense, switched profile, as temporal interpolation is ineffective. The majority of the spectral energy of the CO2 storage time series was found to be between 0.001 and 0.2 Hz (500 and 5 s respectively); however, sampling frequencies of 2 Hz and below still result in significantly lower CO2 storage values. An empirical method of correcting CO2 storage values from under-sampled time series is proposed.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Shadow lines in the arithmetic of elliptic curves

Let E/Q be an elliptic curve and p a rational prime of good ordinary reduction. For every imaginary quadratic field K/Q satisfying the Heegner hypothesis for E we have a corresponding line in E(K)\otimes Q_p, known as a shadow line. When E/Q has analytic rank 2 and E/K has analytic rank 3, shadow lines are expected to lie in E(Q)\otimes Qp. If, in addition, p splits in K/Q, then shadow lines can be determined using the anticyclotomic p-adic height pairing. We develop an algorithm to compute anticyclotomic p-adic heights which we then use to provide an algorithm to compute shadow lines. We conclude by illustrating these algorithms in a collection of examples.