The residual circulation of the Southern Ocean: which spatio-temporal scales are needed?

The Southern Ocean circulation consists of a complicated mixture of processes and phenomena that arise at different time and spatial scales which need to be parametrized in the state-of-the-art climate models. The temporal and spatial scales that give rise to the present-day residual mean circulation are here investigated by calculating the Meridional Overturning Circulation (MOC) in density coordinates from an eddy-permitting global model. The region sensitive to the temporal decomposition is located between 38°S and 63°S, associated with the eddy-induced transport. The ‘‘Bolus’’ component of the residual circulation corresponds to the eddy-induced transport. It is dominated by timescales between 1 month and 1 year. The temporal behavior of the transient eddies is examined in splitting the ‘‘Bolus’’ component into a ‘‘Seasonal’’, an ‘‘Eddy’’ and an ‘‘Inter-monthly’’ component, respectively representing the correlation between density and velocity fluctuations due to the average seasonal cycle, due to mesoscale eddies and due to large-scale motion on timescales longer than one month that is not due to the seasonal cycle. The ‘‘Seasonal’’ bolus cell is important at all latitudes near the surface. The ‘‘Eddy’’ bolus cell is dominant in the thermocline between 50°S and 35°S and over the whole ocean depth at the latitude of the Drake Passage. The ‘‘Inter-monthly’’ bolus cell is important in all density classes and is maximal in the Brazil–Malvinas Confluence and the Agulhas Return Current. The spatial decomposition indicates that a large part of the Eulerian mean circulation is recovered for spatial scales larger than 11.25°, implying that small-scale meanders in the Antarctic Circumpolar Current (ACC), near the Subantarctic and Polar Fronts, and near the Subtropical Front are important in the compensation of the Eulerian mean flow.

Spatio-temporal analysis of an invasive plant pathogen (Phytophthora ramorum) in England and Wales

Phytophthora ramorum is a damaging invasive plant pathogen and was first discovered in the UK in 2002. Spatial point analyses were applied to the occurrence of this disease in England and Wales during the period of 2003-2006 in order to assess its spatio-temporal spread. Out of the 4301 garden centres and nurseries (GCN) surveyed, there were 164, 105, 123 and 41 sites with P. ramorum in 2003, 2004, 2005 and 2006, respectively. Spatial analysis of the observed point patterns of GCN outbreaks suggested that these sites were significantly clumped within a radius of ca 60 km in 2003, but not in later years. Further analyses were conducted to determine the relationship of GCN outbreak sites over two consecutive years and thus to infer possible disease spread over time. This analysis suggested that disease spread among GCN sites was most likely to have occurred within a distance of 60 km for 2003-2004, but not for the later years. There were 35, 63, 81 and 58 sites with P. ramorum in the semi-natural environment (SNE). Analyses were carried out to assess whether infected GCN sites could act as an inoculum source of infected SNE plants or vice versa. In all years, there was a significant spatial closeness among GCN and SNE outbreak sites within a distance of 1 km. But a significant relationship over a longer distance (within 60 km) was only observed between cases in 2003 and 2004. These analyses suggest that statutory actions taken so far appear to have reduced the extent of long-distance spread of P. ramorum among garden centres and nurseries, but not the disease spread at a shorter distance between GCN and SNE sites.

Modeling data with structural and temporal correlation using lower level and higher level multilevel models

Novel imaging techniques are playing an increasingly important role in drug development, providing insight into the mechanism of action of new chemical entities. The data sets obtained by these methods can be large with complex inter-relationships, but the most appropriate statistical analysis for handling this data is often uncertain - precisely because of the exploratory nature of the way the data are collected. We present an example from a clinical trial using magnetic resonance imaging to assess changes in atherosclerotic plaques following treatment with a tool compound with established clinical benefit. We compared two specific approaches to handle the correlations due to physical location and repeated measurements: two-level and four-level multilevel models. The two methods identified similar structural variables, but higher level multilevel models had the advantage of explaining a greater proportion of variation, and the modeling assumptions appeared to be better satisfied.

Spatio-temporal surface soil heat flux estimates from satellite data: results for the AMMA experiment at the Fakara (Niger) supersite

This paper describes a method that employs Earth Observation (EO) data to calculate spatiotemporal estimates of soil heat flux, G, using a physically-based method (the Analytical Method). The method involves a harmonic analysis of land surface temperature (LST) data. It also requires an estimate of near-surface soil thermal inertia; this property depends on soil textural composition and varies as a function of soil moisture content. The EO data needed to drive the model equations, and the ground-based data required to provide verification of the method, were obtained over the Fakara domain within the African Monsoon Multidisciplinary Analysis (AMMA) program. LST estimates (3 km × 3 km, one image 15 min−1) were derived from MSG-SEVIRI data. Soil moisture estimates were obtained from ENVISAT-ASAR data, while estimates of leaf area index, LAI, (to calculate the effect of the canopy on G, largely due to radiation extinction) were obtained from SPOT-HRV images. The variation of these variables over the Fakara domain, and implications for values of G derived from them, were discussed. Results showed that this method provides reliable large-scale spatiotemporal estimates of G. Variations in G could largely be explained by the variability in the model input variables. Furthermore, it was shown that this method is relatively insensitive to model parameters related to the vegetation or soil texture. However, the strong sensitivity of thermal inertia to soil moisture content at low values of relative saturation (<0.2) means that in arid or semi-arid climates accurate estimates of surface soil moisture content are of utmost importance, if reliable estimates of G are to be obtained. This method has the potential to improve large-scale evaporation estimates, to aid land surface model prediction and to advance research that aims to explain failure in energy balance closure of meteorological field studies.

Spatio-temporal dependencies in functional connectivity in rodent cortical cultures

Models of functional connectivity in cortical cultures on multi-electrodes arrays may aid in understanding how cognitive pathways form and improve techniques that aim to interface with neuronal systems. To enable research on such models, this study uses both data- and model-driven approaches to determine what dependencies are present in and between functional connectivity networks derived from bursts of extracellularly recorded activity. Properties of excitation in bursts were analysed using correlative techniques to assess the degree of linear dependence and then two parallel techniques were used to assess functional connectivity. Three models presenting increasing levels of spatio-temporal dependency were used to capture the dynamics of individual functional connections and their consistencies were verified using surrogate data. By comparing network-wide properties between model generated networks and functional networks from data, complex interdependencies were revealed. This indicates the persistent co-activation of neuronal pathways in spontaneous bursts, as can be found in whole brain structures.

Spatio-temporal variation in predation by urban domestic cats (Felis catus) and the acceptability of possible management actions in the UK

Urban domestic cat (Felis catus) populations can attain exceedingly high densities and are not limited by natural prey availability. This has generated concerns that they may negatively affect prey populations, leading to calls for management. We enlisted cat-owners to record prey returned home to estimate patterns of predation by free-roaming pets in different localities within the town of Reading, UK and questionnaire surveys were used to quantify attitudes to different possible management strategies. Prey return rates were highly variable: only 20% of cats returned ≥4 dead prey annually. Consequently, approximately 65% of owners received no prey in a given season, but this declined to 22% after eight seasons. The estimated mean predation rate was 18.3 prey cat−1 year−1 but this varied markedly both spatially and temporally: per capita predation rates declined with increasing cat density. Comparisons with estimates of the density of six common bird prey species indicated that cats killed numbers equivalent to adult density on c. 39% of occasions. Population modeling studies suggest that such predation rates could significantly reduce the size of local bird populations for common urban species. Conversely, most urban residents did not consider cat predation to be a significant problem. Collar-mounted anti-predation devices were the only management action acceptable to the majority of urban residents (65%), but were less acceptable to cat-owners because of perceived risks to their pets; only 24% of cats were fitted with such devices. Overall, cat predation did appear to be of sufficient magnitude to affect some prey populations, although further investigation of some key aspects of cat predation is warranted. Management of the predation behavior of urban cat populations in the UK is likely to be challenging and achieving this would require considerable engagement with cat owners.

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.

In-situ spatio-temporal measurements of the detailed azimuthal substructure of the substorm current wedge

The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100~km at altitudes of 4,000-7,000~km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240~s after Cluster 4 at 1,300-2,000~km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction, in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven ``wedgelets". Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

Spatio-temporal variations of ozone and nitrogen dioxide concentrations under urban trees and in a nearby open area

The evergreen Quercus ilex L. is one of the most common trees in Italian urban environments and is considered effective in the uptake of particulate and gaseous atmospheric pollutants. However, the few available estimates on O3 and NO2 removal by urban Q. ilex originate from model-based studies (which indicate NO2/O3 removal capacity of Q. ilex) and not from direct measurements of air pollutant concentrations. Thus, in the urban area of Siena (central Italy) we began long-term monitoring of O3/NO2 concentrations using passive samplers at a distance of 1, 5, 10 m from a busy road, under the canopies of Q. ilex and in a nearby open-field. Measurements performed in the period June 2011-October 2013 showed always a greater decrease of NO2 concentrations under the Q. ilex canopy than in the open-field transect. Conversely, a decrease of average O3 concentrations under the tree canopy was found only in autumn after the typical Mediterranean post-summer rainfalls. Our results indicate that interactions between O3/NO2 concentrations and trees in Mediterranean urban ecosystems are affected by temporal variations in climatic conditions. We argue therefore that the direct measurement of atmospheric pollutant concentrations should be chosen to describe local changes of aerial pollution.

Spatio-temporal variability of polynya dynamics and ice production in the Laptev Sea between the winters of 1979/80 and 2007/08

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km39 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

Large-scale and synoptic meteorology in the south-east Pacific during the observations campaign VOCALS-REx in austral Spring 2008

We present a descriptive overview of the meteorology in the south eastern subtropical Pacific (SEP) during the VOCALS-REx intensive observations campaign which was carried out between October and November 2008. Mainly based on data from operational analyses, forecasts, reanalysis, and satellite observations, we focus on spatio-temporal scales from synoptic to planetary. A climatological context is given within which the specific conditions observed during the campaign are placed, with particular reference to the relationships between the large-scale and the regional circulations. The mean circulations associated with the diurnal breeze systems are also discussed. We then provide a summary of the day-to-day synoptic-scale circulation, air-parcel trajectories, and cloud cover in the SEP during VOCALS-REx. Three meteorologically distinct periods of time are identified and the large-scale causes for their different character are discussed. The first period was characterised by significant variability associated with synoptic-scale systems interesting the SEP; while the two subsequent phases were affected by planetary-scale disturbances with a slower evolution. The changes between initial and later periods can be partly explained from the regular march of the annual cycle, but contributions from subseasonal variability and its teleconnections were important. Across the whole of the two months under consideration we find a significant correlation between the depth of the inversion-capped marine boundary layer (MBL) and the amount of low cloud in the area of study. We discuss this correlation and argue that at least as a crude approximation a typical scaling may be applied relating MBL and cloud properties with the large-scale parameters of SSTs and tropospheric temperatures. These results are consistent with previously found empirical relationships involving lower-tropospheric stability.

Effects of urban ventilation patterns on the carbon monoxide concentration in a high-rise mega city

Carbon monoxide (CO) concentration data from 1999–2006, monitored at 5 different pollution stations in a high-rise mega city (Hong Kong), were collected and investigated. The spatio-temporal characteristics of urban CO concentration profiles were obtained. A new approach was put forward to examine the relationship between urban CO concentration and different wind flow patterns. Rather than relying on the meteorological data from a single weather station, usually adopted in previous studies, four weather stations on the boundary of Hong Kong territory were used in the present study so as to identify 16 different wind flow patterns, among which a typical urban heat island circulation (UHIC) can be distinguished. Higher concentrations were observed to be associated with the flow pattern of an inflow from Lau Fau Shan (LFS) station which is located in the northwest of Hong Kong. This suggests that the ability of dilution for north-to-west wind is relatively weak due to the pollutants carried from outside Hong Kong. The effectiveness of wind speed on the alleviation of urban concentration is dependent on the initial concentration of the approaching wind. The increase of wind speed of north-to-west wind from 0 m/s to 6 m/s has little effect on the reduction of urban CO concentration, especially on the non-roadside stations. By contrast, for the southerly marine wind, pollution concentration decreases sharply with an increase in the wind speed. It was also found that urban heat island circulation (UHIC) is conducive of the accumulation of pollutants, especially at night. There exists a positive correlation between CO concentration and UHI intensity. This correlation is much stronger at night compared to during the day. Keywords: urban pollution monitoring, urban ventilation pattern, urban heat island circulation, mega city

Improved evaporative flux partitioning and carbon flux in the land surface model JULES: impact on the simulation of land surface processes in temperate Europe

The primary role of land surface models embedded in climate models is to partition surface available energy into upwards, radiative, sensible and latent heat fluxes. Partitioning of evapotranspiration, ET, is of fundamental importance: as a major component of the total surface latent heat flux, ET affects the simulated surface water balance, and related energy balance, and consequently the feedbacks with the atmosphere. In this context it is also crucial to credibly represent the CO2 exchange between ecosystems and their environment. In this study, JULES, the land surface model used in UK weather and climate models, has been evaluated for temperate Europe. Compared to eddy covariance flux measurements, the CO2 uptake by the ecosystem is underestimated and the ET overestimated. In addition, the contribution to ET from soil and intercepted water evaporation far outweighs the contribution of plant transpiration. To alleviate these biases, adaptations have been implemented in JULES, based on key literature references. These adaptations have improved the simulation of the spatio-temporal variability of the fluxes and the accuracy of the simulated GPP and ET, including its partitioning. This resulted in a shift of the seasonal soil moisture cycle. These adaptations are expected to increase the fidelity of climate simulations over Europe. Finally, the extreme summer of 2003 was used as evaluation benchmark for the use of the model in climate change studies. The improved model captures the impact of the 2003 drought on the carbon assimilation and the water use efficiency of the plants. It, however, underestimates the 2003 GPP anomalies. The simulations showed that a reduction of evaporation from the interception and soil reservoirs, albeit not of transpiration, largely explained the good correlation between the carbon and the water fluxes anomalies that was observed during 2003. This demonstrates the importance of being able to discriminate the response of individual component of the ET flux to environmental forcing.