Boundary layer dynamics over London, UK, as observed using Doppler lidar during REPARTEE-II

Urban boundary layers (UBLs) can be highly complex due to the heterogeneous roughness and heating of the surface, particularly at night. Due to a general lack of observations, it is not clear whether canonical models of boundary layer mixing are appropriate in modelling air quality in urban areas. This paper reports Doppler lidar observations of turbulence profiles in the centre of London, UK, as part of the second REPARTEE campaign in autumn 2007. Lidar-measured standard deviation of vertical velocity averaged over 30 min intervals generally compared well with in situ sonic anemometer measurements at 190 m on the BT telecommunications Tower. During calm, nocturnal periods, the lidar underestimated turbulent mixing due mainly to limited sampling rate. Mixing height derived from the turbulence, and aerosol layer height from the backscatter profiles, showed similar diurnal cycles ranging from c. 300 to 800 m, increasing to c. 200 to 850 m under clear skies. The aerosol layer height was sometimes significantly different to the mixing height, particularly at night under clear skies. For convective and neutral cases, the scaled turbulence profiles resembled canonical results; this was less clear for the stable case. Lidar observations clearly showed enhanced mixing beneath stratocumulus clouds reaching down on occasion to approximately half daytime boundary layer depth. On one occasion the nocturnal turbulent structure was consistent with a nocturnal jet, suggesting a stable layer. Given the general agreement between observations and canonical turbulence profiles, mixing timescales were calculated for passive scalars released at street level to reach the BT Tower using existing models of turbulent mixing. It was estimated to take c. 10 min to diffuse up to 190 m, rising to between 20 and 50 min at night, depending on stability. Determination of mixing timescales is important when comparing to physico-chemical processes acting on pollutant species measured simultaneously at both the ground and at the BT Tower during the campaign. From the 3 week autumnal data-set there is evidence for occasional stable layers in central London, effectively decoupling surface emissions from air aloft.

Remarkable dynamics of nanoparticles in the urban atmosphere

Nanoparticles emitted from road traffic are the largest source of respiratory exposure for the general public living in urban areas. It has been suggested that the adverse health effects of airborne particles may scale with the airborne particle number, which if correct, focuses attention on the nanoparticle (less than 100 nm) size range which dominates the number count in urban areas. Urban measurements of particle size distributions have tended to show a broadly similar pattern dominated by a mode centred on 20–30 nm diameter particles emitted by diesel engine exhaust. In this paper we report the results of measurements of particle number concentration and size distribution made in a major London park as well as on the BT Tower, 160 m high. These measurements taken during the REPARTEE project (Regents Park and BT Tower experiment) show a remarkable shift in particle size distributions with major losses of the smallest particle class as particles are advected away from the traffic source. In the Park, the traffic related mode at 20–30 nm diameter is much reduced with a new mode at <10 nm. Size distribution measurements also revealed higher number concentrations of sub-50 nm particles at the BT Tower during days affected by higher turbulence as determined by Doppler Lidar measurements and indicate a loss of nanoparticles from air aged during less turbulent conditions. These results suggest that nanoparticles are lost by evaporation, rather than coagulation processes. The results have major implications for understanding the impacts of traffic-generated particulate matter on human health.

System dynamics mapping of acute patient flows

Department of Health staff wished to use systems modelling to discuss acute patient flows with groups of NHS staff. The aim was to assess the usefulness of system dynamics (SD) in a healthcare context and to elicit proposals concerning ways of improving patient experience. Since time restrictions excluded simulation modelling, a hybrid approach using stock/flow symbols from SD was created. Initial interviews and hospital site visits generated a series of stock/flow maps. A ‘Conceptual Framework’ was then created to introduce the mapping symbols and to generate a series of questions about different patient paths and what might speed or slow patient flows. These materials formed the centre of three workshops for NHS staff. The participants were able to propose ideas for improving patient flows and the elicited data was subsequently employed to create a finalized suite of maps of a general acute hospital. The maps and ideas were communicated back to the Department of Health and subsequently assisted the work of the Modernization Agency.

Dynamics of the lower stratospheric circulation response to ENSO

A robust feature of the observed response to El Nin˜o–Southern Oscillation (ENSO) is an altered circulation in the lower stratosphere. When sea surface temperatures (SSTs) in the tropical Pacific are warmer there is enhanced upwelling and cooling in the tropical lower stratosphere and downwelling and warming in the midlatitudes, while the opposite is true of cooler SSTs. The midlatitude lower stratospheric response to ENSO is larger in the Southern Hemisphere (SH) than in the Northern Hemisphere (NH). In this study the dynamical version of the Canadian Middle Atmosphere Model (CMAM) is used to simulate 25 realizations of the atmospheric response to the 1982/83 El Nin˜o and the 1973/74 La Nin˜ a. This version ofCMAMis a comprehensive high-top general circulation model that does not include interactive chemistry. The observed lower stratospheric response to ENSO is well reproduced by the simulations, allowing them to be used to investigate the mechanisms involved. Both the observed and simulated responses maximize in December–March and so this study focuses on understanding the mechanisms involved in that season. The response in tropical upwelling is predominantly driven by anomalous transient synoptic-scale wave drag in the SH subtropical lower stratosphere, which is also responsible for the compensating SH midlatitude response. This altered wave drag stems from an altered upward flux of wave activity from the troposphere into the lower stratosphere between 208 and 408S. The altered flux of wave activity can be divided into two distinct components. In the Pacific, the acceleration of the zonal wind in the subtropics from the warmer tropical SSTs results in a region between the midlatitude and subtropical jets where there is an enhanced source of low phase speed eddies. At other longitudes, an equatorward shift of the midlatitude jet from the extratropical tropospheric response to El Nin˜o results in an enhanced source of waves of higher phase speeds in the subtropics. The altered resolved wave drag is only apparent in the SH and the difference between the two hemispheres can be related to the difference in the climatological jet structures in this season and the projection of the wind anomalies associated with ENSO onto those structures.

Stratospheric climate and variability from a general circulation model and observations

The climate and natural variability of the large-scale stratospheric circulation simulated by a newly developed general circulation model are evaluated against available global observations. The simulation consisted of a 30-year annual cycle integration performed with a comprehensive model of the troposphere and stratosphere. The observations consisted of a 15-year dataset from global operational analyses of the troposphere and stratosphere. The model evaluation concentrates on the simulation of the evolution of the extratropical stratospheric circulation in both hemispheres. The December–February climatology of the observed zonal mean winter circulation is found to be reasonably well captured by the model, although in the Northern Hemisphere upper stratosphere the simulated westerly winds are systematically stronger and a cold bias is apparent in the polar stratosphere. This Northern Hemisphere stratospheric cold bias virtually disappears during spring (March–May), consistent with a realistic simulation of the spring weakening of the mean westerly winds in the model. A considerable amount of monthly interannual variability is also found in the simulation in the Northern Hemisphere in late winter and early spring. The simulated interannual variability is predominantly caused by polar warmings of the stratosphere, in agreement with observations. The breakdown of the Northern Hemisphere stratospheric polar vortex appears therefore to occur in a realistic way in the model. However, in early winter the model severely underestimates the interannual variability, especially in the upper troposphere. The Southern Hemisphere winter (June–August) zonal mean temperature is systematically colder in the model, and the simulated winds are somewhat too strong in the upper stratosphere. Contrary to the results for the Northern Hemisphere spring, this model cold bias worsens during the Southern Hemisphere spring (September–November). Significant discrepancies between the model results and the observations are therefore found during the breakdown of the Southern Hemisphere polar vortex. For instance, the simulated Southern Hemisphere stratosphere westerly jet continuously decreases in intensity more or less in situ from June to November, while the observed stratospheric jet moves downward and poleward.

Climate predictability experiments with a general circulation model

The atmospheric response to the evolution of the global sea surface temperatures from 1979 to 1992 is studied using the Max-Planck-Institut 19 level atmospheric general circulation model, ECHAM3 at T 42 resolution. Five separate 14-year integrations are performed and results are presented for each individual realization and for the ensemble-averaged response. The results are compared to a 30-year control integration using a climate monthly mean state of the sea surface temperatures and to analysis data. It is found that the ECHAM3 model, by and large, does reproduce the observed response pattern to El Nin˜o and La Nin˜a. During the El Nin˜ o events, the subtropical jet streams in both hemispheres are intensified and displaced equatorward, and there is a tendency towards weak upper easterlies over the equator. The Southern Oscillation is a very stable feature of the integrations and is accurately reproduced in all experiments. The inter-annual variability at middle- and high-latitudes, on the other hand, is strongly dominated by chaotic dynamics, and the tropical SST forcing only modulates the atmospheric circulation. The potential predictability of the model is investigated for six different regions. Signal to noise ratio is large in most parts of the tropical belt, of medium strength in the western hemisphere and generally small over the European area. The ENSO signal is most pronounced during the boreal spring. A particularly strong signal in the precipitation field in the extratropics during spring can be found over the southern United States. Western Canada is normally warmer during the warm ENSO phase, while northern Europe is warmer than normal during the ENSO cold phase. The reason is advection of warm air due to a more intense Pacific low than normal during the warm ENSO phase and a more intense Icelandic low than normal during the cold ENSO phase, respectively.

Zonal-mean dynamics of extended recoveries from stratospheric sudden warmings

The recovery of the Arctic polar vortex following stratospheric sudden warmings is found to take upward of 3 months in a particular subset of cases, termed here polar-night jet oscillation (PJO) events. The anomalous zonal-mean circulation above the pole during this recovery is characterized by a persistently warm lower stratosphere, and above this a cold midstratosphere and anomalously high stratopause, which descends as the event unfolds. Composites of these events in the Canadian Middle Atmosphere Model show the persistence of the lower-stratospheric anomaly is a result of strongly suppressed wave driving and weak radiative cooling at these heights. The upper-stratospheric and lower-mesospheric anomalies are driven immediately following the warming by anomalous planetary-scale eddies, following which, anomalous parameterized nonorographic and orographic gravity waves play an important role. These details are found to be robust for PJO events (as opposed to sudden warmings in general) in that many details of individual PJO events match the composite mean. Azonal-mean quasigeostrophic model on the sphere is shown to reproduce the response to the thermal and mechanical forcings produced during a PJO event. The former is well approximated by Newtonian cooling. The response can thus be considered as a transient approach to the steady-state, downward control limit. In this context, the time scale of the lower-stratospheric anomaly is determined by the transient, radiative response to the extended absence of wave driving. The extent to which the dynamics of the wave-driven descent of the stratopause can be considered analogous to the descending phases of the quasi-biennial oscillation (QBO) is also discussed.

Lateral boundary contributions to wave-activity invariants and nonlinear stability theorems for balanced dynamics

The slow advective-timescale dynamics of the atmosphere and oceans is referred to as balanced dynamics. An extensive body of theory for disturbances to basic flows exists for the quasi-geostrophic (QG) model of balanced dynamics, based on wave-activity invariants and nonlinear stability theorems associated with exact symmetry-based conservation laws. In attempting to extend this theory to the semi-geostrophic (SG) model of balanced dynamics, Kushner & Shepherd discovered lateral boundary contributions to the SG wave-activity invariants which are not present in the QG theory, and which affect the stability theorems. However, because of technical difficulties associated with the SG model, the analysis of Kushner & Shepherd was not fully nonlinear. This paper examines the issue of lateral boundary contributions to wave-activity invariants for balanced dynamics in the context of Salmon's nearly geostrophic model of rotating shallow-water flow. Salmon's model has certain similarities with the SG model, but also has important differences that allow the present analysis to be carried to finite amplitude. In the process, the way in which constraints produce boundary contributions to wave-activity invariants, and additional conditions in the associated stability theorems, is clarified. It is shown that Salmon's model possesses two kinds of stability theorems: an analogue of Ripa's small-amplitude stability theorem for shallow-water flow, and a finite-amplitude analogue of Kushner & Shepherd's SG stability theorem in which the ‘subsonic’ condition of Ripa's theorem is replaced by a condition that the flow be cyclonic along lateral boundaries. As with the SG theorem, this last condition has a simple physical interpretation involving the coastal Kelvin waves that exist in both models. Salmon's model has recently emerged as an important prototype for constrained Hamiltonian balanced models. The extent to which the present analysis applies to this general class of models is discussed.

An examination of climate sensitivity for idealised climate change experiments in an intermediate general circulation model

Radiative forcing and climate sensitivity have been widely used as concepts to understand climate change. This work performs climate change experiments with an intermediate general circulation model (IGCM) to examine the robustness of the radiative forcing concept for carbon dioxide and solar constant changes. This IGCM has been specifically developed as a computationally fast model, but one that allows an interaction between physical processes and large-scale dynamics; the model allows many long integrations to be performed relatively quickly. It employs a fast and accurate radiative transfer scheme, as well as simple convection and surface schemes, and a slab ocean, to model the effects of climate change mechanisms on the atmospheric temperatures and dynamics with a reasonable degree of complexity. The climatology of the IGCM run at T-21 resolution with 22 levels is compared to European Centre for Medium Range Weather Forecasting Reanalysis data. The response of the model to changes in carbon dioxide and solar output are examined when these changes are applied globally and when constrained geographically (e.g. over land only). The CO2 experiments have a roughly 17% higher climate sensitivity than the solar experiments. It is also found that a forcing at high latitudes causes a 40% higher climate sensitivity than a forcing only applied at low latitudes. It is found that, despite differences in the model feedbacks, climate sensitivity is roughly constant over a range of distributions of CO2 and solar forcings. Hence, in the IGCM at least, the radiative forcing concept is capable of predicting global surface temperature changes to within 30%, for the perturbations described here. It is concluded that radiative forcing remains a useful tool for assessing the natural and anthropogenic impact of climate change mechanisms on surface temperature.

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988

As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

The effect of electric field direction on the-shear stress of ER fluids

The effect of the direction of external electric field on the shear stress of an ER fluid has been studied by molecular-dynamics simulation. Due to the formation of inclined chains, the shear stress strongly depends on the direction of the field, and it may be very large under some special field direction. And theoretical model of ideal microstructure of ER fluids has proved this result. Thus the ER effect may be greatly enhanced just by choosing an optimum direction for the field without any additional requirement, suggesting a promising way to the practical application of ER fluids.

Rerum cognoscere causas: part II - opportunities generated by the agency/structure debate and suggestions for clarifying the social theoretic position of system dynamics

This is the second half of a two-part paper dealing with the social theoretic assumptions underlying system dynamics. In the first half it was concluded that analysing system dynamics using traditional, paradigm-based social theories is highly problematic. An innovative and potentially fruitful resolution is now proposed to these problems. In the first section it is argued that in order to find an appropriate social theoretic home for system dynamics it is necessary to look to a key exchange in contemporary social science: the agency/structure debate. This debate aims to move beyond both the theories based only on the actions of individual human agents, and those theories that emphasise only structural influences. Emerging from this debate are various theories that instead aim to unite the human agent view of the social realm with views that concentrate solely on system structure. It is argued that system dynamics is best viewed as being implicitly grounded in such theories. The main conclusion is therefore that system dynamics can contribute to an important part of social thinking by providing a formal approach for explicating social mechanisms. This conclusion is of general significance for system dynamics. However, the over-arching aim of the two-part paper is to increase the understanding of system dynamics in related disciplines. Four suggestions are therefore offered for how the system dynamics method might be extended further into the social sciences. It is argued that, presented in the right way, the formal yet contingent feedback causality thinking of system dynamics should diffuse widely in the social sciences and make a distinctive and important contribution to them. Felix qui potuit rerum cognoscere causas Happy is he who comes to know the causes of things Virgil - Georgics, Book II, line 490. 29 BCE