A parametric model for the interval censored survival times of Acacia mangium plantation in a spacing trial

Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.

The impact of traffic emissions on atmospheric ozone and OH: results from QUANTIFY

To estimate the impact of emissions by road, aircraft and ship traffic on ozone and OH in the present-day atmosphere six different atmospheric chemistry models have been used. Based on newly developed global emission inventories for road, ship and aircraft emission data sets each model performed sensitivity simulations reducing the emissions of each transport sector by 5%. The model results indicate that on global annual average lower tropospheric ozone responds most sensitive to ship emissions (50.6%±10.9% of the total traffic induced perturbation), followed by road (36.7%±9.3%) and aircraft exhausts (12.7%±2.9%), respectively. In the northern upper troposphere between 200–300 hPa at 30–60° N the maximum impact from road and ship are 93% and 73% of the maximum effect of aircraft, respectively. The latter is 0.185 ppbv for ozone (for the 5% case) or 3.69 ppbv when scaling to 100%. On the global average the impact of road even dominates in the UTLS-region. The sensitivity of ozone formation per NOx molecule emitted is highest for aircraft exhausts. The local maximum effect of the summed traffic emissions on the ozone column predicted by the models is 0.2 DU and occurs over the northern subtropical Atlantic extending to central Europe. Below 800 hPa both ozone and OH respond most sensitively to ship emissions in the marine lower troposphere over the Atlantic. Based on the 5% perturbation the effect on ozone can exceed 0.6% close to the marine surface (global zonal mean) which is 80% of the total traffic induced ozone perturbation. In the southern hemisphere ship emissions contribute relatively strongly to the total ozone perturbation by 60%–80% throughout the year. Methane lifetime changes against OH are affected strongest by ship emissions up to 0.21 (± 0.05)%, followed by road (0.08 (±0.01)%) and air traffic (0.05 (± 0.02)%). Based on the full scale ozone and methane perturbations positive radiative forcings were calculated for road emissions (7.3±6.2 mWm−2) and for aviation (2.9±2.3 mWm−2). Ship induced methane lifetime changes dominate over the ozone forcing and therefore lead to a net negative forcing (−25.5±13.2 mWm−2).

Long-lived halocarbon trends and budgets from atmospheric chemistry modelling constrained with measurements in polar firn

The budgets of seven halogenated gases (CFC-11, CFC-12, CFC-113, CFC-114, CFC-115, CCl4 and SF6) are studied by comparing measurements in polar firn air from two Arctic and three Antarctic sites, and simulation results of two numerical models: a 2-D atmospheric chemistry model and a 1-D firn diffusion model. The first one is used to calculate atmospheric concentrations from emission trends based on industrial inventories; the calculated concentration trends are used by the second one to produce depth concentration profiles in the firn. The 2-D atmospheric model is validated in the boundary layer by comparison with atmospheric station measurements, and vertically for CFC-12 by comparison with balloon and FTIR measurements. Firn air measurements provide constraints on historical atmospheric concentrations over the last century. Age distributions in the firn are discussed using a Green function approach. Finally, our results are used as input to a radiative model in order to evaluate the radiative forcing of our target gases. Multi-species and multi-site firn air studies allow to better constrain atmospheric trends. The low concentrations of all studied gases at the bottom of the firn, and their consistency with our model results confirm that their natural sources are small. Our results indicate that the emissions, sinks and trends of CFC-11, CFC-12, CFC-113, CFC-115 and SF6 are well constrained, whereas it is not the case for CFC-114 and CCl4. Significant emission-dependent changes in the lifetimes of halocarbons destroyed in the stratosphere were obtained. Those result from the time needed for their transport from the surface where they are emitted to the stratosphere where they are destroyed. Efforts should be made to update and reduce the large uncertainties on CFC lifetimes.

Radiative forcing due to aviation water vapour emissions

Three emissions inventories have been used with a fully Lagrangian trajectory model to calculate the stratospheric accumulation of water vapour emissions from aircraft, and the resulting radiative forcing. The annual and global mean radiative forcing due to present-day aviation water vapour emissions has been found to be 0.9 [0.3 to 1.4] mW m^2. This is around a factor of three smaller than the value given in recent assessments, and the upper bound is much lower than a recently suggested 20 mW m^2 upper bound. This forcing is sensitive to the vertical distribution of emissions, and, to a lesser extent, interannual variability in meteorology. Large differences in the vertical distribution of emissions within the inventories have been identified, which result in the choice of inventory being the largest source of differences in the calculation of the radiative forcing due to the emissions. Analysis of Northern Hemisphere trajectories demonstrates that the assumption of an e-folding time is not always appropriate for stratospheric emissions. A linear model is more representative for emissions that enter the stratosphere far above the tropopause.

Accelerated loss of alpine glaciers in the Kodar Mountains, south-eastern Siberia

The recession of mountain glaciers around the world has been linked to anthropogenic climate change and small glaciers (e.g. < 2 km2) are thought to be particularly vulnerable, with reports of their disappearance from several regions. However, the response of small glaciers to climate change can be modulated by non-climatic factors such as topography and debris cover and there remain a number of regions where their recent change has evaded scrutiny. This paper presents results of the first multi-year remote sensing survey of glaciers in the Kodar Mountains, the only glaciers in SE Siberia, which we compare to previous glacier inventories from this continental setting that reported total glacier areas of 18.8 km2 in ca. 1963 (12.6 km2 of exposed ice) and 15.5 km2 in 1974 (12 km2 of exposed ice). Mapping their debris-covered termini is difficult but delineation of debris-free ice on Landsat imagery reveals 34 glaciers with a total area of 11.72 ± 0.72 km2 in 1995, followed by a reduction to 9.53 ± 0.29 km2 in 2001 and 7.01 ± 0.23 km2 in 2010. This represents a ~ 44% decrease in exposed glacier ice between ca. 1963 and 2010, but with 40% lost since 1995 and with individual glaciers losing as much as 93% of their exposed ice. Thus, although continental glaciers are generally thought to be less sensitive than their maritime counterparts, a recent acceleration in shrinkage of exposed ice has taken place and we note its coincidence with a strong summer warming trend in the region initiated at the start of the 1980s. Whilst smaller and shorter glaciers have, proportionally, tended to shrink more rapidly, we find no statistically significant relationship between shrinkage and elevation characteristics, aspect or solar radiation. This is probably due to the small sample size, limited elevation range, and topographic setting of the glaciers in deep valleys-heads. Furthermore, many of the glaciers possess debris-covered termini and it is likely that the ablation of buried ice is lagging the shrinkage of exposed ice, such that a growth in the proportion of debris cover is occurring, as observed elsewhere. If recent trends continue, we hypothesise that glaciers could evolve into a type of rock glacier within the next few decades, introducing additional complexity in their response and delaying their potential demise.

Characterisation of Bolivian savanna ecosystems by their modern pollen rain and implications for fossil pollen records

The majority of vegetation reconstructions from the Neotropics are derived from fossil pollen records extracted from lake sediments. However, the interpretation of these records is restricted by limited knowledge of the contemporary relationships between the vegetation and pollen rain of Neotropical ecosystems, especially for more open vegetation such as savannas. This research aims to improve the interpretation of these records by investigating the vegetation and modern pollen rain of different savanna ecosystems in Bolivia using vegetation inventories, artificial pollen traps and surface lake sediments. Two types of savanna were studied, upland savannas (cerrado), occurring on well drained soils, and seasonally-inundated savannas occurring on seasonally water-logged soils. Quantitative vegetation data are used to identify taxa that are floristically important in the different savanna types and to allow modern pollen/vegetation ratios to be calculated. Artificial pollen traps from the upland savanna site are dominated by Moraceae (35%), Poaceae (30%), Alchornea (6%) and Cecropia (4%). The two seasonally-inundated savanna sites are dominated by Moraceae (37%), Poaceae (20%), Alchornea (8%) and Cecropia (7%), and Moraceae (25%), Cyperaceae (22%), Poaceae (19%) and Cecropia (9%), respectively. The modern pollen rain of seasonally-inundated savannas from surface lake sediments is dominated by Cyperaceae (35%), Poaceae (33%), Moraceae (9%) and Asteraceae (5%). Upland and seasonally-flooded savannas were found to be only subtly distinct from each other palynologically. All sites have a high proportion of Moraceae pollen due to effective wind dispersal of this pollen type from areas of evergreen forest close to the study sites. Modern pollen/vegetation ratios show that many key woody plant taxa are absent/under-represented in the modern pollen rain (e.g., Caryocar and Tabebuia). The lower-than-expected percentages of Poaceae pollen, and the scarcity of savanna indicators, in the modern pollen rain of these ecosystems mean that savannas could potentially be overlooked in fossil pollen records without consideration of the full pollen spectrum available.

Differentiation between Neotropical rainforest, dry forest and savannah ecosystems by their modern pollen spectra and implications for the fossil pollen record

Accurate differentiation between tropical forest and savannah ecosystems in the fossil pollen record is hampered by the combination of: i) poor taxonomic resolution in pollen identification, and ii) the high species diversity of many lowland tropical families, i.e. with many different growth forms living in numerous environmental settings. These barriers to interpreting the fossil record hinder our understanding of the past distributions of different Neotropical ecosystems and consequently cloud our knowledge of past climatic, biodiversity and carbon storage patterns. Modern pollen studies facilitate an improved understanding of how ecosystems are represented by the pollen their plants produce and therefore aid interpretation of fossil pollen records. To understand how to differentiate ecosystems palynologically, it is essential that a consistent sampling method is used across ecosystems. However, to date, modern pollen studies from tropical South America have employed a variety of methodologies (e.g. pollen traps, moss polsters, soil samples). In this paper, we present the first modern pollen study from the Neotropics to examine the modern pollen rain from moist evergreen tropical forest (METF), semi-deciduous dry tropical forest (SDTF) and wooded savannah (cerradão) using a consistent sampling methodology (pollen traps). Pollen rain was sampled annually in September for the years 1999–2001 from within permanent vegetation study plots in, or near, the Noel Kempff Mercado National Park (NKMNP), Bolivia. Comparison of the modern pollen rain within these plots with detailed floristic inventories allowed estimates of the relative pollen productivity and dispersal for individual taxa to be made (% pollen/% vegetation or ‘p/v’). The applicability of these data to interpreting fossil records from lake sediments was then explored by comparison with pollen assemblages obtained from five lake surface samples.

The sources of phosphorus in the waters of Great Britain

Total phosphorus (TP) and soluble reactive phosphorus (SRP) loads to watercourses of the River Basin Districts (RBDs) of Great Britain (GB) were estimated using inventories of industrial P loads and estimates of P loads from sewage treatment works and diffuse P loads calculated using region-specific export coefficients for particular land cover classes combined with census data for agricultural stocking densities and human populations. The TP load to GB waters was estimated to be 60 kt yr(-1), of which households contributed 73, agriculture contributed 20, industry contributed 3, and 4 came from background sources. The SRP load to GB waters was estimated to be 47 kt yr(-1), of which households contributed 78, agriculture contributed 13, industry contributed 4, and 6 came from background Sources. The 'average' area-normalized TP and SRP loads to GB waters approximated 2.4 kg ha(-1) yr(-1) and 1.8 kg ha(-1) yr(-1), respectively. A consideration of uncertainties in the data contributing to these estimates suggested that the TP load to GB waters might lie between 33 and 68 kt yr(-1), with agriculture contributing between 10 and 28 of the TP load. These estimates are consistent with recent appraisals of annual TP and SRP loads to GB coastal waters and area-normalized TP loads from their catchments. Estimates of the contributions of RBDs to these P loads were consistent with the geographical distribution of P concentrations in GB rivers and recent assessments of surface waters at risk from P Pollution.

Atmospheric chemistry and physics in the atmosphere of a developed megacity (London): an overview of the REPARTEE experiment and its conclusions

The REgents PARk and Tower Environmental Experiment (REPARTEE) comprised two campaigns in London in October 2006 and October/November 2007. The experiment design involved measurements at a heavily trafficked roadside site, two urban background sites and an elevated site at 160–190 m above ground on the BT Tower, supplemented in the second campaign by Doppler lidar measurements of atmospheric vertical structure. A wide range of measurements of airborne particle physical metrics and chemical composition were made as well as measurements of a considerable range of gas phase species and the fluxes of both particulate and gas phase substances. Significant findings include (a) demonstration of the evaporation of traffic-generated nanoparticles during both horizontal and vertical atmospheric transport; (b) generation of a large base of information on the fluxes of nanoparticles, accumulation mode particles and specific chemical components of the aerosol and a range of gas phase species, as well as the elucidation of key processes and comparison with emissions inventories; (c) quantification of vertical gradients in selected aerosol and trace gas species which has demonstrated the important role of regional transport in influencing concentrations of sulphate, nitrate and secondary organic compounds within the atmosphere of London; (d) generation of new data on the atmospheric structure and turbulence above London, including the estimation of mixed layer depths; (e) provision of new data on trace gas dispersion in the urban atmosphere through the release of purposeful tracers; (f) the determination of spatial differences in aerosol particle size distributions and their interpretation in terms of sources and physico-chemical transformations; (g) studies of the nocturnal oxidation of nitrogen oxides and of the diurnal behaviour of nitrate aerosol in the urban atmosphere, and (h) new information on the chemical composition and source apportionment of particulate matter size fractions in the atmosphere of London derived both from bulk chemical analysis and aerosol mass spectrometry with two instrument types.

Effects of urban density on carbon dioxide exchanges: observations of dense urban, suburban and woodland areas of southern England

Anthropogenic and biogenic controls on the surface–atmosphere exchange of CO2 are explored for three different environments. Similarities are seen between suburban and woodland sites during summer, when photosynthesis and respiration determine the diurnal pattern of the CO2 flux. In winter, emissions from human activities dominate urban and suburban fluxes; building emissions increase during cold weather, while traffic is a major component of CO2 emissions all year round. Observed CO2 fluxes reflect diurnal traffic patterns (busy throughout the day (urban); rush-hour peaks (suburban)) and vary between working days and non-working days, except at the woodland site. Suburban vegetation offsets some anthropogenic emissions, but 24-h CO2 fluxes are usually positive even during summer. Observations are compared to estimated emissions from simple models and inventories. Annual CO2 exchanges are significantly different between sites, demonstrating the impacts of increasing urban density (and decreasing vegetation fraction) on the CO2 flux to the atmosphere.

A multi-model evaluation of aerosols over South Asia: common problems and possible causes

Atmospheric pollution over South Asia attracts special attention due to its effects on regional climate, water cycle and human health. These effects are potentially growing owing to rising trends of anthropogenic aerosol emissions. In this study, the spatio-temporal aerosol distributions over South Asia from seven global aerosol models are evaluated against aerosol retrievals from NASA satellite sensors and ground-based measurements for the period of 2000–2007. Overall, substantial underestimations of aerosol loading over South Asia are found systematically in most model simulations. Averaged over the entire South Asia, the annual mean aerosol optical depth (AOD) is underestimated by a range 15 to 44% across models compared to MISR (Multi-angle Imaging SpectroRadiometer), which is the lowest bound among various satellite AOD retrievals (from MISR, SeaWiFS (Sea-Viewing Wide Field-of-View Sensor), MODIS (Moderate Resolution Imaging Spectroradiometer) Aqua and Terra). In particular during the post-monsoon and wintertime periods (i.e., October–January), when agricultural waste burning and anthropogenic emissions dominate, models fail to capture AOD and aerosol absorption optical depth (AAOD) over the Indo–Gangetic Plain (IGP) compared to ground-based Aerosol Robotic Network (AERONET) sunphotometer measurements. The underestimations of aerosol loading in models generally occur in the lower troposphere (below 2 km) based on the comparisons of aerosol extinction profiles calculated by the models with those from Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) data. Furthermore, surface concentrations of all aerosol components (sulfate, nitrate, organic aerosol (OA) and black carbon (BC)) from the models are found much lower than in situ measurements in winter. Several possible causes for these common problems of underestimating aerosols in models during the post-monsoon and wintertime periods are identified: the aerosol hygroscopic growth and formation of secondary inorganic aerosol are suppressed in the models because relative humidity (RH) is biased far too low in the boundary layer and thus foggy conditions are poorly represented in current models, the nitrate aerosol is either missing or inadequately accounted for, and emissions from agricultural waste burning and biofuel usage are too low in the emission inventories. These common problems and possible causes found in multiple models point out directions for future model improvements in this important region.

Cities reducing their greenhouse gas emissions

The study asks how well are cities doing in reducing their greenhouse gas emissions. Data from six cities with repeat GHG emission inventories for the period 2004–2009 is examined: Berlin, Boston, Greater Toronto, London, New York City and Seattle. All of the cities are reducing their per capita GHG emissions, primarily through changes to stationary combustion. On average the cities are reducing per capita emissions by 0.27 t CO2e/capita per year; this is about the same average rate as the cities nation states, although the cities are reducing emissions faster in percentage terms.

Greenhouse gas emissions from waste management—assessment of quantification methods

Of the many sources of urban greenhouse gas (GHG) emissions, solid waste is the only one for which management decisions are undertaken primarily by municipal governments themselves and is hence often the largest component of cities’ corporate inventories. It is essential that decision-makers select an appropriate quantification methodology and have an appreciation of methodological strengths and shortcomings. This work compares four different waste emissions quantification methods, including Intergovernmental Panel on Climate Change (IPCC) 1996 guidelines, IPCC 2006 guidelines, U.S. Environmental Protection Agency (EPA) Waste Reduction Model (WARM), and the Federation of Canadian Municipalities- Partners for Climate Protection (FCM-PCP) quantification tool. Waste disposal data for the greater Toronto area (GTA) in 2005 are used for all methodologies; treatment options (including landfill, incineration, compost, and anaerobic digestion) are examined where available in methodologies. Landfill was shown to be the greatest source of GHG emissions, contributing more than three-quarters of total emissions associated with waste management. Results from the different landfill gas (LFG) quantification approaches ranged from an emissions source of 557 kt carbon dioxide equivalents (CO2e) (FCM-PCP) to a carbon sink of −53 kt CO2e (EPA WARM). Similar values were obtained between IPCC approaches. The IPCC 2006 method was found to be more appropriate for inventorying applications because it uses a waste-in-place (WIP) approach, rather than a methane commitment (MC) approach, despite perceived onerous data requirements for WIP. MC approaches were found to be useful from a planning standpoint; however, uncertainty associated with their projections of future parameter values limits their applicability for GHG inventorying. MC and WIP methods provided similar results in this case study; however, this is case specific because of similarity in assumptions of present and future landfill parameters and quantities of annual waste deposited in recent years being relatively consistent.

The Earls of Kildare and their books at the end of the Middle Ages

The two surviving inventories of the library of the Fitzgerald Earls of Kildare bear witness to a particularly large and diverse collection of books in the Earls' castle at Maynooth, Co. Kildare. Between them, the lists record well over one hundred separate items in four languages: Latin, French, English and Irish. This paper traces the history of the library and analyses the Fitzgeralds' particular interests as book collectors and as readers. It provides the first full published set of suggested identifications and bibliographical details for the books at Maynooth. It also includes a fresh transcription of the library lists and a discussion of the manuscript context in which they are preserved. Sources like the Kildare library lists provide valuable evidence for the potential circulation of a wide range of non-native manuscripts and prints in late-medieval Ireland.

Ratings of age of acquisition of 299 words across 25 languages: is there a cross-linguistic order of words?

We present a new set of subjective age-of-acquisition (AoA) ratings for 299 words (158 nouns, 141 verbs) in 25 languages from five language families (Afro-Asiatic: Semitic languages; Altaic: one Turkic language: Indo-European: Baltic, Celtic, Germanic, Hellenic, Slavic, and Romance languages; Niger-Congo: one Bantu language; Uralic: Finnic and Ugric languages). Adult native speakers reported the age at which they had learned each word. We present a comparison of the AoA ratings across all languages by contrasting them in pairs. This comparison shows a consistency in the orders of ratings across the 25 languages. The data were then analyzed (1) to ascertain how the demographic characteristics of the participants influenced AoA estimations and (2) to assess differences caused by the exact form of the target question (when did you learn vs. when do children learn this word); (3) to compare the ratings obtained in our study to those of previous studies; and (4) to assess the validity of our study by comparison with quasi-objective AoA norms derived from the MacArthur–Bates Communicative Development Inventories (MB-CDI). All 299 words were judged as being acquired early (mostly before the age of 6 years). AoA ratings were associated with the raters’ social or language status, but not with the raters’ age or education. Parents reported words as being learned earlier, and bilinguals reported learning them later. Estimations of the age at which children learn the words revealed significantly lower ratings of AoA. Finally, comparisons with previous AoA and MB-CDI norms support the validity of the present estimations. Our AoA ratings are available for research or other purposes.