From ensemble forecasts to predictive distribution functions

The translation of an ensemble of model runs into a probability distribution is a common task in model-based prediction. Common methods for such ensemble interpretations proceed as if verification and ensemble were draws from the same underlying distribution, an assumption not viable for most, if any, real world ensembles. An alternative is to consider an ensemble as merely a source of information rather than the possible scenarios of reality. This approach, which looks for maps between ensembles and probabilistic distributions, is investigated and extended. Common methods are revisited, and an improvement to standard kernel dressing, called ‘affine kernel dressing’ (AKD), is introduced. AKD assumes an affine mapping between ensemble and verification, typically not acting on individual ensemble members but on the entire ensemble as a whole, the parameters of this mapping are determined in parallel with the other dressing parameters, including a weight assigned to the unconditioned (climatological) distribution. These amendments to standard kernel dressing, albeit simple, can improve performance significantly and are shown to be appropriate for both overdispersive and underdispersive ensembles, unlike standard kernel dressing which exacerbates over dispersion. Studies are presented using operational numerical weather predictions for two locations and data from the Lorenz63 system, demonstrating both effectiveness given operational constraints and statistical significance given a large sample.

A ‘Boscastle-type’ quasi-stationary convective system over the UK Southwest Peninsula

An investigation is presented of a quasi-stationary convective system (QSCS) which occurred over the UK Southwest Peninsula on 21 July 2010. This system was remarkably similar in its location and structure to one which caused devastating flash flooding in the coastal village of Boscastle, Cornwall on 16 August 2004. However, in the 2010 case rainfall accumulations were around four times smaller and no flooding was recorded. The more extreme nature of the Boscastle case is shown to be related to three factors: (1) higher rain rates, associated with a warmer and moister tropospheric column and deeper convective clouds; (2) a more stationary system, due to slower evolution of the large-scale flow; and (3) distribution of the heaviest precipitation over fewer river catchments. Overall, however, the synoptic setting of the two events was broadly similar, suggesting that such conditions favour the development of QSCSs over the Southwest Peninsula. A numerical simulation of the July 2010 event was performed using a 1.5-km grid length configuration of the Met Office Unified Model. This reveals that convection was repeatedly initiated through lifting of low-level air parcels along a quasi-stationary coastal convergence line. Sensitivity tests are used to show that this convergence line was a sea breeze front which temporarily stalled along the coastline due to the retarding influence of an offshore-directed background wind component. Several deficiencies are noted in the 1.5-km model’s representation of the storm system, including delayed convective initiation; however, significant improvements are observed when the grid length is reduced to 500 m. These result in part from an improved representation of the convergence line, which enhances the associated low-level ascent allowing air parcels to more readily reach their level of free convection. The implications of this finding for forecasting convective precipitation are discussed.

What is the climate system able to do ‘on its own’?

The climate of the Earth, like planetary climates in general, is broadly controlled by solar irradiation, planetary albedo and emissivity as well as its rotation rate and distribution of land (with its orography) and oceans. However, the majority of climate fluctuations that affect mankind are internal modes of the general circulation of the atmosphere and the oceans. Some of these modes, such as El Nino-Southern Oscillation (ENSO), are quasi-regular and have some longer-term predictive skill; others like the Arctic and Antarctic Oscillation are chaotic and generally unpredictable beyond a few weeks. Studies using general circulation models indicate that internal processes dominate the regional climate and that some like ENSO events have even distinct global signatures. This is one of the reasons why it is so difficult to separate internal climate processes from external ones caused, for example, by changes in greenhouse gases and solar irradiation. However, the accumulation of the warmest seasons during the latest two decades is lending strong support to the forcing of the greenhouse gases. As models are getting more comprehensive, they show a gradually broader range of internal processes including those on longer time scales, challenging the interpretation of the causes of past and present climate events further.

Technical Note: The normal quantile transformation and its application in a flood forecasting system.

The Normal Quantile Transform (NQT) has been used in many hydrological and meteorological applications in order to make the Cumulated Distribution Function (CDF) of the observed, simulated and forecast river discharge, water level or precipitation data Gaussian. It is also the heart of the meta-Gaussian model for assessing the total predictive uncertainty of the Hydrological Uncertainty Processor (HUP) developed by Krzysztofowicz. In the field of geo-statistics this transformation is better known as the Normal-Score Transform. In this paper some possible problems caused by small sample sizes when applying the NQT in flood forecasting systems will be discussed and a novel way to solve the problem will be outlined by combining extreme value analysis and non-parametric regression methods. The method will be illustrated by examples of hydrological stream-flow forecasts.

Non-uniform data distribution for communication-efficient parallel clustering

Global communicationrequirements andloadimbalanceof someparalleldataminingalgorithms arethe major obstacles to exploitthe computational power of large-scale systems. This work investigates how non-uniform data distributions can be exploited to remove the global communication requirement and to reduce the communication costin parallel data mining algorithms and, in particular, in the k-means algorithm for cluster analysis. In the straightforward parallel formulation of the k-means algorithm, data and computation loads are uniformly distributed over the processing nodes. This approach has excellent load balancing characteristics that may suggest it could scale up to large and extreme-scale parallel computing systems. However, at each iteration step the algorithm requires a global reduction operationwhichhinders thescalabilityoftheapproach.Thisworkstudiesadifferentparallelformulation of the algorithm where the requirement of global communication is removed, while maintaining the same deterministic nature ofthe centralised algorithm. The proposed approach exploits a non-uniform data distribution which can be either found in real-world distributed applications or can be induced by means ofmulti-dimensional binary searchtrees. The approachcanalso be extended to accommodate an approximation error which allows a further reduction ofthe communication costs. The effectiveness of the exact and approximate methods has been tested in a parallel computing system with 64 processors and in simulations with 1024 processing element

Inhalation exposure to particulate matter in rooms with underfloor air distribution

This work investigated the personal exposure to indoor particulate matters using the intake fraction metric and provided a possible way to trace the particle inhaled from an indoor particle source. A turbulence model validated by the particle measurements in a room with underfloor air distribution (UFAD) system was used to predict the indoor particle concentrations. Inhalation intake fraction of indoor particles was defined and evaluated in two rooms equipped with the UFAD, i.e., the experimental room and a small office. According to the exposure characteristics and a typical respiratory rate, the intake fraction was determined in two rooms with a continuous and episodic (human cough) source of particles, respectively. The findings showed that the well-mixing assumption of indoor air failed to give an accurate estimation of inhalation exposure and the average concentration at return outlet or within the overall room could not relate well the intake fraction to the amount of particle emitted from an indoor source.

Tracing troposphere-to-stratosphere transport above a mid-latitude deep convective system

Within the project SPURT (trace gas measurements in the tropopause region) a variety of trace gases have been measured in situ in order to investigate the role of dynamical and chemical processes in the extra-tropical tropopause region. In this paper we report on a flight on 10 November 2001 leading from Hohn, Germany (52�N) to Faro, Portugal (37�N) through a strongly developed deep stratospheric intrusion. This streamer was associated with a large convective system over the western Mediterranean with potentially significant troposphere-to-stratosphere transport. Along major parts of the flight we measured unexpectedly high NOy mixing ratios. Also H2O mixing ratios were significantly higher than stratospheric background levels confirming the extraordinary chemical signature of the probed air masses in the interior of the streamer. Backward trajectories encompassing the streamer enable to analyze the origin and physical characteristics of the air masses and to trace troposphere-to-stratosphere transport. Near the western flank of the streamer features caused by long range transport, such as tropospheric filaments characterized by sudden drops in the O3 and NOy mixing ratios and enhanced CO and H2O can be reconstructed in great detail using the reverse domain filling technique. These filaments indicate a high potential for subsequent mixing with the stratospheric air. At the south-western edge of the streamer a strong gradient in the NOy and the O3 mixing ratios coincides very well with a sharp gradient in potential vorticity in the ECMWF fields. In contrast, in the interior of the streamer the observed highly elevated NOy and H2O mixing ratios up to a potential temperature level of 365K and potential vorticity values of maximum 10 PVU cannot be explained in terms of resolved troposphere-to-stratosphere transport along the backward trajectories. Also mesoscale simulations with a High Resolution Model reveal no direct evidence for convective H2O injection up to this level. Elevated H2O mixing ratios in the ECMWF and HRM are seen only up to about tropopause height at 340 hPa and 270 hPa, respectively, well below flight altitude of about 200 hPa. However, forward tracing of the convective influence as identified by satellite brightness temperature measurements and counts of lightning strokes shows that during this part of the flight the aircraft was closely following the border of an air mass which was heavily impacted by convective activity over Spain and Algeria. This is evidence that deep convection at mid-latitudes may have a large impact on the tracer distribution of the lowermost stratosphere reaching well above the thunderstorms anvils as claimed by recent studies using cloud-resolving models.

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop-soil system: validation and application

Models for water transfer in the crop-soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop-soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop-soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration. (c) 2009 Elsevier B.V. All rights reserved.

Fire in the Earth system

Fire is a worldwide phenomenon that appears in the geological record soon after the appearance of terrestrial plants. Fire influences global ecosystem patterns and processes, including vegetation distribution and structure, the carbon cycle, and climate. Although humans and fire have always coexisted, our capacity to manage fire remains imperfect and may become more difficult in the future as climate change alters fire regimes. This risk is difficult to assess, however, because fires are still poorly represented in global models. Here, we discuss some of the most important issues involved in developing a better understanding of the role of fire in the Earth system.

Fire in the Earth System

Fire is an important component of the Earth System that is tightly coupled with climate, vegetation, biogeochemical cycles, and human activities. Observations of how fire regimes change on seasonal to millennial timescales are providing an improved understanding of the hierarchy of controls on fire regimes. Climate is the principal control on fire regimes, although human activities have had an increasing influence on the distribution and incidence of fire in recent centuries. Understanding of the controls and variability of fire also underpins the development of models, both conceptual and numerical, that allow us to predict how future climate and land-use changes might influence fire regimes. Although fires in fire-adapted ecosystems can be important for biodiversity and ecosystem function, positive effects are being increasingly outweighed by losses of ecosystem services. As humans encroach further into the natural habitat of fire, social and economic costs are also escalating. The prospect of near-term rapid and large climate changes, and the escalating costs of large wildfires, necessitates a radical re-thinking and the development of approaches to fire management that promote the more harmonious co-existence of fire and people.

A review of the performance of different ventilation and airflow distribution systems in buildings

The objective of this article is to review the scientific literature on airflow distribution systems and ventilation effectiveness to identify and assess the most suitable room air distribution methods for various spaces. In this study, different ventilation systems are classified according to specific requirements and assessment procedures. This study shows that eight ventilation methods have been employed in the built environment for different purposes and tasks. The investigation shows that numerous studies have been carried out on ventilation effectiveness but few studies have been done regarding other aspects of air distribution. Amongst existing types of ventilation systems, the performance of each ventilation methods varies from one case to another due to different usages of the ventilation system in a room and the different assessment indices used. This review shows that the assessment of ventilation effectiveness or efficiency should be determined according to each task of the ventilation system, such as removal of heat, removal of pollutant, supply fresh air to the breathing zone or protecting the occupant from cross infection. The analysis results form a basic framework regarding the application of airflow distribution for the benefit of designers, architects, engineers, installers and building owners.

Trust in authorities monitoring the distribution of genetically modified foods: dimensionality, measurement issues, and determinants

Based on a combined internet and mail survey in Germany the independence of indica-tors of trust in public authorities from indicators of attitudes toward genetically modified food is tested. Despite evidence of a link between trust indicators on the one hand and evaluation of benefits and perceived likelihoods of risks, correlation with other factors is found to be moderate on average. But the trust indicators exhibit only a moderate relation with the re-spondents’ preference for either sole public control or a cooperation of public and private bodies in the monitoring of GM food distribution. Instead, age and location in either the New or the Old Lander are found to be significantly related with such preferences.

Freshwater transport in the coupled ocean-atmosphere system: a passive ocean

Conservation of water demands that meridional ocean and atmosphere freshwater transports (FWT) are of equal magnitude but opposite in direction. This suggests that the atmospheric FWT and its associated latent heat (LH) transport could be thought of as a \textquotedblleft coupled ocean/atmosphere mode\textquotedblright. But what is the true nature of this coupling? Is the ocean passive or active? Here we analyze a series of simulations with a coupled ocean-atmosphere-sea ice model employing highly idealized geometries but with markedly different coupled climates and patterns of ocean circulation. Exploiting streamfunctions in specific humidity coordinates for the atmosphere and salt coordinates for the ocean to represent FWT in their respective medium, we find that atmospheric FWT/LH transport is essentially independent of the ocean state. Ocean circulation and salinity distribution adjust to achieve a return freshwater pathway demanded of them by the atmosphere. So, although ocean and atmosphere FWTs are indeed coupled by mass conservation, the ocean is a passive component acting as a reservoir of freshwater.

Trends in Morbid Obesity and in Bariatric Surgeries Covered by the Brazilian Public Health System

Obesity is an increasingly serious public health problem on a global level. Morbid obesity, defined as a body mass index greater than 40 kg/m(2), is associated with increased mortality and a high burden of obesity-related morbidities. To study the prevalence of morbid obesity in Brazil, three national anthropometric surveys were reanalyzed. Data about bariatric surgeries were obtained from the Ministry of Health Hospital Information System, which is available online. A 255% rise in the prevalence of morbid obesity was observed, starting at 0.18% in 1975-1976 and growing to 0.33% in 1989 and 0.64% in 2002-2003. There was a higher rate in the South in the first two surveys, but the prevalence in the Southeast rose steadily, reaching 0.77% in 2002-2003 and overtaking the South. Since 1999, the Brazilian Unified Health System has covered surgical treatment for morbid obesity. From 2000 to 2006, there was a sixfold increase in the number of surgeries, which topped the 2,500 mark in 2006. The geographic distribution of these surgeries is heavily concentrated in the Southeast, the most developed region of Brazil, where there is also the highest prevalence of morbid obesity. This was followed by the Southern region. The figures for the rise in morbid obesity in Brazil are startling, especially the increase among men. This is a situation that calls for further study, alongside measures to encourage the adoption of healthy lifestyles. Preventive measures aimed at slowing down or reversing the obesity epidemic are urgently required.

Moisture transport and intraseasonal variability in the South America monsoon system

This paper examines moisture transport on intraseasonal timescales over the continent and over the South Atlantic convergence zone (SACZ) during the South America (SA) summer monsoon. Combined Empirical Orthogonal Function analysis (EOFc) of Global Precipitation Climatology Project pentad precipitation, specific humidity, air temperature, zonal and meridional winds at 850 hPa (NCEP/NCAR reanalysis) are performed to identify the large-scale variability of the South America monsoon system and the SACZ. The first EOFc was used as a large-scale index for the South American monsoon (LISAM), whereas the second EOFc characterized the SACZ. LISAM (SACZ) index showed spectral variance on 30-90 (15-20) days and were both band filtered (10-100 days). Intraseasonal wet anomalies were defined when LISAM and SACZ anomalies were above the 75th percentile of their respective distribution. LISAM and SACZ wet events were examined independently of each other and when they occur simultaneously. LISAM wet events were observed with the amplification of wave activity in the Northern Hemisphere and the enhancement of northwesterly cross-equatorial moisture transport over tropical continental SA. Enhanced SACZ was observed with moisture transport from the extratropics of the Southern Hemisphere. Simultaneous LISAM and SACZ wet events are associated with cross-equatorial moisture transport along with moisture transport from Subtropical Southwestern Atlantic.