A spatio-temporal structure-based approach to drought characterisation

Drought characterisation is an intrinsically spatio-temporal problem. A limitation of previous approaches to characterisation is that they discard much of the spatio-temporal information by reducing events to a lower-order subspace. To address this, an explicit 3-dimensional (longitude, latitude, time) structure-based method is described in which drought events are defined by a spatially and temporarily coherent set of points displaying standardised precipitation below a given threshold. Geometric methods can then be used to measure similarity between individual drought structures. Groupings of these similarities provide an alternative to traditional methods for extracting recurrent space-time signals from geophysical data. The explicit consideration of structure encourages the construction of summary statistics which relate to the event geometry. Example measures considered are the event volume, centroid, and aspect ratio. The utility of a 3-dimensional approach is demonstrated by application to the analysis of European droughts (15 °W to 35°E, and 35 °N to 70°N) for the period 1901–2006. Large-scale structure is found to be abundant with 75 events identified lasting for more than 3 months and spanning at least 0.5 × 106 km2. Near-complete dissimilarity is seen between the individual drought structures, and little or no regularity is found in the time evolution of even the most spatially similar drought events. The spatial distribution of the event centroids and the time evolution of the geographic cross-sectional areas strongly suggest that large area, sustained droughts result from the combination of multiple small area (∼106 km2) short duration (∼3 months) events. The small events are not found to occur independently in space. This leads to the hypothesis that local water feedbacks play an important role in the aggregation process.

A robust mechanism for strengthening of the Brewer–Dobson circulation in response to climate change: critical-layer control of subtropical wave breaking

Climate models consistently predict a strengthened Brewer–Dobson circulation in response to greenhouse gas (GHG)-induced climate change. Although the predicted circulation changes are clearly the result of changes in stratospheric wave drag, the mechanism behind the wave-drag changes remains unclear. Here, simulations from a chemistry–climate model are analyzed to show that the changes in resolved wave drag are largely explainable in terms of a simple and robust dynamical mechanism, namely changes in the location of critical layers within the subtropical lower stratosphere, which are known from observations to control the spatial distribution of Rossby wave breaking. In particular, the strengthening of the upper flanks of the subtropical jets that is robustly expected from GHG-induced tropospheric warming pushes the critical layers (and the associated regions of wave drag) upward, allowing more wave activity to penetrate into the subtropical lower stratosphere. Because the subtropics represent the critical region for wave driving of the Brewer–Dobson circulation, the circulation is thereby strengthened. Transient planetary-scale waves and synoptic-scale waves generated by baroclinic instability are both found to play a crucial role in this process. Changes in stationary planetary wave drag are not so important because they largely occur away from subtropical latitudes.

The 20 February 2010 Madeira flash-floods: synoptic analysis and extreme rainfall assessment

This study aims to characterise the rainfall exceptionality and the meteorological context of the 20 February 2010 flash-floods in Madeira (Portugal). Daily and hourly precipitation records from the available rain-gauge station networks are evaluated in order to reconstitute the temporal evolution of the rainstorm, as its geographic incidence, contributing to understand the flash-flood dynamics and the type and spatial distribution of the associated impacts. The exceptionality of the rainstorm is further confirmed by the return period associated with the daily precipitation registered at the two long-term record stations, with 146.9 mm observed in the city of Funchal and 333.8 mm on the mountain top, corresponding to an estimated return period of approximately 290 yr and 90 yr, respectively. Furthermore, the synoptic associated situation responsible for the flash-floods is analysed using different sources of information, e.g., weather charts, reanalysis data, Meteosat images and radiosounding data, with the focus on two main issues: (1) the dynamical conditions that promoted such anomalous humidity availability over the Madeira region on 20 February 2010 and (2) the uplift mechanism that induced deep convection activity.

Cyclones causing wind storms in the Mediterranean: characteristics, trends and links to large-scale patterns

A climatology of cyclones with a focus on their relation to wind storm tracks in the Mediterranean region (MR) is presented. Trends in the frequency of cyclones and wind storms, as well as variations associated with the North Atlantic Oscillation (NAO), the East Atlantic/West Russian (EAWR) and the Scandinavian variability pattern (SCAND) are discussed. The study is based on the ERA40 reanalysis dataset. Wind storm tracks are identified by tracking clusters of adjacent grid boxes characterised by extremely high local wind speeds. The wind track is assigned to a cyclone track independently identified with an objective scheme. Areas with high wind activity – quantified by extreme wind tracks – are typically located south of the Golf of Genoa, south of Cyprus, southeast of Sicily and west of the Iberian Peninsula. About 69% of the wind storms are caused by cyclones located in the Mediterranean region, while the remaining 31% can be attributed to North Atlantic or Northern European cyclones. The North Atlantic Oscillation, the East Atlantic/West Russian pattern and the Scandinavian pattern all influence the amount and spatial distribution of wind inducing cyclones and wind events in the MR. The strongest signals exist for the NAO and the EAWR pattern, which are both associated with an increase in the number of organised strong wind events in the eastern MR during their positive phase. On the other hand, the storm numbers decrease over the western MR for the positive phase of the NAO and over the central MR during the positive phase of the EAWR pattern. The positive phase of the Scandinavian pattern is associated with a decrease in the number of winter wind storms over most of the MR. A third of the trends in the number of wind storms and wind producing cyclones during the winter season of the ERA40 period may be attributed to the variability of the North Atlantic Oscillation.

Parameterised diabatic processes in numerical simulations of an extratropical cyclone

The parameterisation of diabatic processes in numerical models is critical for the accuracy of weather forecasts and for climate projections. A novel approach to the evaluation of these processes in models is introduced in this contribution. The approach combines a suite of on-line tracer diagnostics with off-line trajectory calculations. Each tracer tracks accumulative changes in potential temperature associated with a particular parameterised diabatic process in the model. A comparison of tracers therefore allows the identification of the most active diabatic processes and their downstream impacts. The tracers are combined with trajectories computed using model-resolved winds, allowing the various diabatic contributions to be tracked back to their time and location of occurrence. We have used this approach to investigate diabatic processes within a simulated extratropical cyclone. We focus on the warm conveyor belt, in which the dominant diabatic contributions come from large-scale latent heating and parameterised convection. By contrasting two simulations, one with standard convection parameterisation settings and another with reduced parameterised convection, the effects of parameterised convection on the structure of the cyclone have been determined. Under reduced parameterised convection conditions, the large-scale latent heating is forced to release convective instability that would otherwise have been released by the convection parameterisation. Although the spatial distribution of precipitation depends on the details of the split between parameterised convection and large-scale latent heating, the total precipitation amount associated with the cyclone remains largely unchanged. For reduced parameterised convection, a more rapid and stronger latent heating episode takes place as air ascends within the warm conveyor belt.

A multithickness sea ice model accounting for sliding friction

A multithickness sea ice model explicitly accounting for the ridging and sliding friction contributions to sea ice stress is developed. Both ridging and sliding contributions depend on the deformation type through functions adopted from the Ukita and Moritz kinematic model of floe interaction. In contrast to most previous work, the ice strength of a uniform ice sheet of constant ice thickness is taken to be proportional to the ice thickness raised to the 3/2 power, as is revealed in discrete element simulations by Hopkins. The new multithickness sea ice model for sea ice stress has been implemented into the Los Alamos “CICE” sea ice model code and is shown to improve agreement between model predictions and observed spatial distribution of sea ice thickness in the Arctic.

Simulations of the impact of orbital forcing and ocean on the Asian summer monsoon during the Holocene

The importance of orbital forcing and ocean impact on the Asian summer monsoon in the Holocene is investigated by comparing simulations with a fully coupled ocean–atmosphere general circulation model (FOAM) and with the atmospheric component of this model (FSSTAM) forced with prescribed modern sea-surface temperatures (SSTs). The results show: (1) the ocean amplifies the orbitally-induced increase in African monsoon precipitation, makes somewhat increase in southern India and damps the increase over the southeastern China. (2) The ocean could change the spatial distribution and local intensity of the orbitally-induced latitudinal atmospheric oscillation over the southeastern China and the subtropical western Pacific Ocean. (3) The orbital forcing mostly enhances the Asian summer precipitation in the FOAM and FSSTAM simulations. However, the ocean reduces the orbitally-induced summer precipitation and postpones the time of summer monsoon onset over the Asian monsoon region. (4) The orbital forcing considerably enhances the intensity of upper divergence, which is amplified by ocean further, over the eastern hemisphere. But the divergence is weaker in the FOAM simulations than in the FSSTAM simulations when the orbital forcing is fixed. (5) The orbital forcing can enhance the amplitude of precipitation variability over the subtropical Africa, the southeastern China and northwestern China, inversely, reduce it over central India and North China in the FOAM and FSSTAM simulations. The ocean obviously reduces the amplitude of precipitation variability over most of the Asian monsoon regions in the fixed orbital forcing simulations. (6) The areas characterized by increased summer precipitation in the long-term mean are mostly characterized by increased amplitude of short-term variability, whereas regions characterized by decreased precipitation are primarily characterized by decreased amplitude of short-term variability. However, the influences of orbital forcing or dynamical ocean on regional climate depend on the model.

The water balance of northern Africa during the mid-Holocene: an evaluation of the 6ka BP PMIP simulations

Runoff fields over northern Africa (10–25°N, 20°W–30°E) derived from 17 atmospheric general circulation models driven by identical 6 ka BP orbital forcing, sea surface temperatures, and CO2 concentration have been analyzed using a hydrological routing scheme (HYDRA) to simulate changes in lake area. The AGCM-simulated runoff produced six-fold differences in simulated lake area between models, although even the largest simulated changes considerably underestimate the observed changes in lake area during the mid-Holocene. The inter-model differences in simulated lake area are largely due to differences in simulated runoff (the squared correlation coefficient, R2, is 0.84). Most of these differences can be attributed to differences in the simulated precipitation (R2=0.83). The higher correlation between runoff and simulated lake area (R2=0.92) implies that simulated differences in evaporation have a contributory effect. When runoff is calculated using an offline land-surface scheme (BIOME3), the correlation between runoff and simulated lake area is (R2=0.94). Finally, the spatial distribution of simulated precipitation can exert an important control on the overall response.

Formulation and implementation of a residual-mean ocean circulation model

A parameterization of mesoscale eddies in coarse-resolution ocean general circulation models (GCM) is formulated and implemented using a residual-mean formalism. In that framework, mean buoyancy is advected by the residual velocity (the sum of the Eulerian and eddy-induced velocities) and modified by a residual flux which accounts for the diabatic effects of mesoscale eddies. The residual velocity is obtained by stepping forward a residual-mean momentum equation in which eddy stresses appear as forcing terms. Study of the spatial distribution of eddy stresses, derived by using them as control parameters to ‘‘fit’’ the residual-mean model to observations, supports the idea that eddy stresses can be likened to a vertical down-gradient flux of momentum with a coefficient which is constant in the vertical. The residual eddy flux is set to zero in the ocean interior, where mesoscale eddies are assumed to be quasi-adiabatic, but is parameterized by a horizontal down-gradient diffusivity near the surface where eddies develop a diabatic component as they stir properties horizontally across steep isopycnals. The residual-mean model is implemented and tested in the MIT general circulation model. It is shown that the resulting model (1) has a climatology that is superior to that obtained using the Gent and McWilliams parameterization scheme with a spatially uniform diffusivity and (2) allows one to significantly reduce the (spurious) horizontal viscosity used in coarse resolution GCMs.

Urban microclimates and urban heat island in Chongqing, China

Chongqing is the largest directly-controlled municipality in China, which is now undergoing a rapid urbanization. The urbanization rate increased from 35.6% in 2000 to 48.3% in 2007, and it is estimated to reach at least 70% by 2020. The question remains open: What are the consequences of such rapid urbanization in Chongqing in terms of urban microclimate? Furthermore, Chongqing is located within the Three Gorges Reservoir (TGR) region and the upper Yangtze River, where the Three Gorges Reservoir (TGR) project started in 1993 and was completed in 2010. As one of the biggest construction projects in the world with a rising water level of 175m and water storage capacity of about 39.3 billion m3, it would be interesting to investigate how such a gigantic project impacts the surrounding micro-environment, especially in Chongqing. Different research approaches are adopted in the study. Our literature review indicates present studies on the urban climate in Chongqing are mainly confined within the historical trend analysis of several weather stations operated by the Chongqing government, little is known about the spatial distribution of urban air temperature and how the local land cover influences the air temperature, especially when there are rivers running through the Chongqing urban area. To contribute to the present knowledge, a series of field measurement campaigns and numerical simulations were carried out. Two complementary types of field measurements are included: fixed weather stations and mobile transverse measurement. Numerical simulations using a house-developed program are able to predict the urban air temperature in Chongqing.

Computational fluid dynamics modelling of the air movement in an environmental test chamber with a respiring manikin

In recent years, computational fluid dynamics (CFD) has been widely used as a method of simulating airflow and addressing indoor environment problems. The complexity of airflows within the indoor environment would make experimental investigation difficult to undertake and also imposes significant challenges on turbulence modelling for flow prediction. This research examines through CFD visualization how air is distributed within a room. Measurements of air temperature and air velocity have been performed at a number of points in an environmental test chamber with a human occupant. To complement the experimental results, CFD simulations were carried out and the results enabled detailed analysis and visualization of spatial distribution of airflow patterns and the effect of different parameters to be predicted. The results demonstrate the complexity of modelling human exhalation within a ventilated enclosure and shed some light into how to achieve more realistic predictions of the airflow within an occupied enclosure.

Effects of noise and by-catch on a Danish harbour porpoise population

Ships and wind turbines generate noise, which can have a negative impact on marine mammal populations by scaring animals away. Effective modelling of how this affects the populations has to take account of the location and timing of disturbances. Here we construct an individual-based model of harbour porpoises in the Inner Danish Waters. Individuals have their own energy budgets constructed using established principles of physiological ecology. Data are lacking on the spatial distribution of food which is instead inferred from knowledge of time-varying porpoise distributions. The model produces plausible patterns of population dynamics and matches well the age distribution of porpoises caught in by-catch. It estimates the effect of existing wind farms as a 10% reduction in population size when food recovers fast (after two days). Proposed new wind farms and ships do not result in further population declines. The population is however sensitive to variations in mortality resulting from by-catch and to the speed at which food recovers after being depleted. If food recovers slowly the effect of wind turbines becomes negligible, whereas ships are estimated to have a significant negative impact on the population. Annual by-catch rates ≥10% lead to monotonously decreasing populations and to extinction, and even the estimated by-catch rate from the adjacent area (approximately 4.1%) has a strong impact on the population. This suggests that conservation efforts should be more focused on reducing by-catch in commercial gillnet fisheries than on limiting the amount of anthropogenic noise. Individual-based models are unique in their ability to take account of the location and timing of disturbances and to show their likely effects on populations. The models also identify deficiencies in the existing database and can be used to set priorities for future field research.

Re-fitting Etton: space, time and material culture within a causewayed enclosure in Cambridgeshire

This paper considers the dynamics of deposition around and across the causewayed enclosure at Etton, Cambridgeshire. As a result of detailed re-analysis (particularly refitting) of the pottery and flint assemblages from the site, it proved possible to shed new light both on the temporality of occupation and the character of deposition there. Certain aspects of our work challenge previous interpretations of the site, and of causewayed enclosures in general; but, just as importantly, others confirm materially what has previously been suggested. The quantities of material deposited at Etton reveal that the enclosure was occupied only very intermittently and certainly less regularly than other contemporary sites in the region. The spatial distribution of material suggests that the enclosure ditch lay open for the entirety of the monument's life, but that acts of deposition generally focused on a specific part of the monument at any one time. As well as enhancing our knowledge of one particular causewayed enclosure, it is hoped that this paper – in combination with our earlier analysis of the pit site at Kilverstone – makes clear the potential that detailed material analysis has to offer in relation to our understanding of the temporality of occupation on prehistoric sites in general.

A numerical model of the ionospheric signatures of time-varying magneticreconnection: I. ionospheric convection

This paper presents a numerical model for predicting the evolution of the pattern of ionospheric convection in response to general time-dependent magnetic reconnection at the dayside magnetopause and in the cross-tail current sheet of the geomagnetic tail. The model quantifies the concepts of ionospheric flow excitation by Cowley and Lockwood (1992), assuming a uniform spatial distribution of ionospheric conductivity. The model is demonstrated using an example in which travelling reconnection pulses commence near noon and then move across the dayside magnetopause towards both dawn and dusk. Two such pulses, 8 min apart, are used and each causes the reconnection to be active for 1 min at every MLT that they pass over. This example demonstrates how the convection response to a given change in the interplanetary magnetic field (via the reconnection rate) depends on the previous reconnection history. The causes of this effect are explained. The inherent assumptions and the potential applications of the model are discussed.

The influence of anisotropic F region ion velocity distributions on ionospheric ion outflows into the magnetosphere

The contribution to the field-aligned ionospheric ion momentum equation, due to coupling between pressure anisotropy and the inhomogeneous geomagnetic field, is investigated. We term this contribution the “hydrodynamic mirror force” and investigate its dependence on the ion drift and the resulting deformations of the ion velocity distribution function from an isotropic form. It is shown that this extra upforce increases rapidly with ion drift relative to the neutral gas but is not highly dependent on the ion-neutral collision model employed. An example of a burst of flow observed by EISCAT, thought to be the ionospheric signature of a flux transfer event at the magnetopause, is studied in detail and it is shown that the nonthermal plasma which results is subject to a hydrodynamic mirror force which is roughly 10% of the gravitational downforce. In addition, predictions by the coupled University College London-Sheffield University model of the ionosphere and thermosphere show that the hydrodynamic mirror force in the auroral oval is up to 3% of the gravitational force for Kp of about 3, rising to 10% following a sudden increase in cross-cap potential. The spatial distribution of the upforce shows peaks in the cusp region and in the post-midnight auroral oval, similar to that of observed low-energy heavy ion flows from the ionosphere into the magnetosphere. We suggest the hydrodynamic mirror force may modulate these outflows by controlling the supply of heavy ions to regions of ion acceleration and that future simulations of the effects of Joule heating on ion outflows should make allowance for it.