Serial clustering of extratropical cyclones in a multi-model ensemble of historical and future simulations

This study has investigated serial (temporal) clustering of extra-tropical cyclones simulated by 17 climate models that participated in CMIP5. Clustering was estimated by calculating the dispersion (ratio of variance to mean) of 30 December-February counts of Atlantic storm tracks passing nearby each grid point. Results from single historical simulations of 1975-2005 were compared to those from historical ERA40 reanalyses from 1958-2001 ERA40 and single future model projections of 2069-2099 under the RCP4.5 climate change scenario. Models were generally able to capture the broad features in reanalyses reported previously: underdispersion/regularity (i.e. variance less than mean) in the western core of the Atlantic storm track surrounded by overdispersion/clustering (i.e. variance greater than mean) to the north and south and over western Europe. Regression of counts onto North Atlantic Oscillation (NAO) indices revealed that much of the overdispersion in the historical reanalyses and model simulations can be accounted for by NAO variability. Future changes in dispersion were generally found to be small and not consistent across models. The overdispersion statistic, for any 30 year sample, is prone to large amounts of sampling uncertainty that obscures the climate change signal. For example, the projected increase in dispersion for storm counts near London in the CNRMCM5 model is 0.1 compared to a standard deviation of 0.25. Projected changes in the mean and variance of NAO are insufficient to create changes in overdispersion that are discernible above natural sampling variations.

The potential of flood forecasting using a variable-resolution global Digital Terrain Model and flood extents from Synthetic Aperture Radar images

A basic data requirement of a river flood inundation model is a Digital Terrain Model (DTM) of the reach being studied. The scale at which modeling is required determines the accuracy required of the DTM. For modeling floods in urban areas, a high resolution DTM such as that produced by airborne LiDAR (Light Detection And Ranging) is most useful, and large parts of many developed countries have now been mapped using LiDAR. In remoter areas, it is possible to model flooding on a larger scale using a lower resolution DTM, and in the near future the DTM of choice is likely to be that derived from the TanDEM-X Digital Elevation Model (DEM). A variable-resolution global DTM obtained by combining existing high and low resolution data sets would be useful for modeling flood water dynamics globally, at high resolution wherever possible and at lower resolution over larger rivers in remote areas. A further important data resource used in flood modeling is the flood extent, commonly derived from Synthetic Aperture Radar (SAR) images. Flood extents become more useful if they are intersected with the DTM, when water level observations (WLOs) at the flood boundary can be estimated at various points along the river reach. To illustrate the utility of such a global DTM, two examples of recent research involving WLOs at opposite ends of the spatial scale are discussed. The first requires high resolution spatial data, and involves the assimilation of WLOs from a real sequence of high resolution SAR images into a flood model to update the model state with observations over time, and to estimate river discharge and model parameters, including river bathymetry and friction. The results indicate the feasibility of such an Earth Observation-based flood forecasting system. The second example is at a larger scale, and uses SAR-derived WLOs to improve the lower-resolution TanDEM-X DEM in the area covered by the flood extents. The resulting reduction in random height error is significant.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Uncertainties in tidally adjusted estimates of sea level rise flooding (bathtub model) for the Greater London

Sea-level rise (SLR) from global warming may have severe consequences for coastal cities, particularly when combined with predicted increases in the strength of tidal surges. Predicting the regional impact of SLR flooding is strongly dependent on the modelling approach and accuracy of topographic data. Here, the areas under risk of sea water flooding for London boroughs were quantified based on the projected SLR scenarios reported in Intergovernmental Panel on Climate Change (IPCC) fifth assessment report (AR5) and UK climatic projections 2009 (UKCP09) using a tidally-adjusted bathtub modelling approach. Medium- to very high-resolution digital elevation models (DEMs) are used to evaluate inundation extents as well as uncertainties. Depending on the SLR scenario and DEMs used, it is estimated that 3%–8% of the area of Greater London could be inundated by 2100. The boroughs with the largest areas at risk of flooding are Newham, Southwark, and Greenwich. The differences in inundation areas estimated from a digital terrain model and a digital surface model are much greater than the root mean square error differences observed between the two data types, which may be attributed to processing levels. Flood models from SRTM data underestimate the inundation extent, so their results may not be reliable for constructing flood risk maps. This analysis provides a broad-scale estimate of the potential consequences of SLR and uncertainties in the DEM-based bathtub type flood inundation modelling for London boroughs.

Atmospheric warming and the amplification of precipitation extremes

Climate models suggest that extreme precipitation events will become more common in an anthropogenically warmed climate. However, observational limitations have hindered a direct evaluation of model-projected changes in extreme precipitation. We used satellite observations and model simulations to examine the response of tropical precipitation events to naturally driven changes in surface temperature and atmospheric moisture content. These observations reveal a distinct link between rainfall extremes and temperature, with heavy rain events increasing during warm periods and decreasing during cold periods. Furthermore, the observed amplification of rainfall extremes is found to be larger than that predicted by models, implying that projections of future changes in rainfall extremes in response to anthropogenic global warming may be underestimated.

Modelled soil organic carbon stocks and changes in the Indo-Gangetic Plains, India from 1980 to 2030

The Global Environment Facility co-financed Soil Organic Carbon (GEFSOC) Project developed a comprehensive modelling system for predicting soil organic carbon (SOC) stocks and changes over time. This research is an effort to predict SOC stocks and changes for the Indian, Indo-Gangetic Plains (IGP), an area with a predominantly rice (Oryza sativa) - wheat (Triticum aestivum) cropping system, using the GEFSOC Modelling System and to compare output with stocks generated using mapping approaches based on soil survey data. The GEFSOC Modelling System predicts an estimated SOC stock for the IGP, India of 1.27, 1.32 and 1.27 Pg for 1990, 2000 and 2030, respectively, in the top 20 cm of soil. The SOC stock using a mapping approach based on soil survey data was 0.66 and 0.88 Pg for 1980 and 2000, respectively. The SOC stock estimated using the GEFSOC Modelling System is higher than the stock estimated using the mapping approach. This is due to the fact that while the GEFSOC System accounts for variation in crop input data (crop management), the soil mapping approach only considers regional variation in soil texture and wetness. The trend of overall change in the modelled SOC stock estimates shows that the IGP, India may have reached an equilibrium following 30-40 years of the Green Revolution. This can be seen in the SOC stock change rates. Various different estimation methods show SOC stocks of 0.57-1.44 Pg C for the study area. The trend of overall change in C stock assessed from the soil survey data indicates that the soils of the IGP, India may store a projected 1.1 Pg of C in 2030. (C) 2007 Elsevier B.V. All rights reserved.

Reflections on practical approaches to involving children and young people in the data analysis process

This article reflects on key methodological issues emerging from children and young people's involvement in data analysis processes. We outline a pragmatic framework illustrating different approaches to engaging children, using two case studies of children's experiences of participating in data analysis. The article highlights methods of engagement and important issues such as the balance of power between adults and children, training, support, ethical considerations, time and resources. We argue that involving children in data analysis processes can have several benefits, including enabling a greater understanding of children's perspectives and helping to prioritise children's agendas in policy and practice. (C) 2007 The Author(s). Journal compilation (C) 2007 National Children's Bureau.

Mesozoic climates: General circulation models and the rock record

General circulation models (GCMs) use the laws of physics and an understanding of past geography to simulate climatic responses. They are objective in character. However, they tend to require powerful computers to handle vast numbers of calculations. Nevertheless, it is now possible to compare results from different GCMs for a range of times and over a wide range of parameterisations for the past, present and future (e.g. in terms of predictions of surface air temperature, surface moisture, precipitation, etc.). GCMs are currently producing simulated climate predictions for the Mesozoic, which compare favourably with the distributions of climatically sensitive facies (e.g. coals, evaporites and palaeosols). They can be used effectively in the prediction of oceanic upwelling sites and the distribution of petroleum source rocks and phosphorites. Models also produce evaluations of other parameters that do not leave a geological record (e.g. cloud cover, snow cover) and equivocal phenomena such as storminess. Parameterisation of sub-grid scale processes is the main weakness in GCMs (e.g. land surfaces, convection, cloud behaviour) and model output for continental interiors is still too cold in winter by comparison with palaeontological data. The sedimentary and palaeontological record provides an important way that GCMs may themselves be evaluated and this is important because the same GCMs are being used currently to predict possible changes in future climate. The Mesozoic Earth was, by comparison with the present, an alien world, as we illustrate here by reference to late Triassic, late Jurassic and late Cretaceous simulations. Dense forests grew close to both poles but experienced months-long daylight in warm summers and months-long darkness in cold snowy winters. Ocean depths were warm (8 degrees C or more to the ocean floor) and reefs, with corals, grew 10 degrees of latitude further north and south than at the present time. The whole Earth was warmer than now by 6 degrees C or more, giving more atmospheric humidity and a greatly enhanced hydrological cycle. Much of the rainfall was predominantly convective in character, often focused over the oceans and leaving major desert expanses on the continental areas. Polar ice sheets are unlikely to have been present because of the high summer temperatures achieved. The model indicates extensive sea ice in the nearly enclosed Arctic seaway through a large portion of the year during the late Cretaceous, and the possibility of sea ice in adjacent parts of the Midwest Seaway over North America. The Triassic world was a predominantly warm world, the model output for evaporation and precipitation conforming well with the known distributions of evaporites, calcretes and other climatically sensitive facies for that time. The message from the geological record is clear. Through the Phanerozoic, Earth's climate has changed significantly, both on a variety of time scales and over a range of climatic states, usually baldly referred to as "greenhouse" and "icehouse", although these terms disguise more subtle states between these extremes. Any notion that the climate can remain constant for the convenience of one species of anthropoid is a delusion (although the recent rate of climatic change is exceptional). (c) 2006 Elsevier B.V. All rights reserved.

Utilizing chromophoric dissolved organic matter measurements to derive export and reactivity of dissolved organic carbon exported to the Arctic Ocean: A case study of the Yukon River, Alaska

The quality and quantity of dissolved organic matter (DOM) exported by Arctic rivers is known to vary with hydrology and this exported material plays a fundamental role in the biogeochemical cycling of carbon at high latitudes. We highlight the potential of optical measurements to examine DOM quality across the hydrograph in Arctic rivers. Furthermore, we establish chromophoric DOM (CDOM) relationships to dissolved organic carbon (DOC) and lignin phenols in the Yukon River and model DOC and lignin loads from CDOM measurements, the former in excellent agreement with long-term DOC monitoring data. Intensive sampling across the historically under-sampled spring flush period highlights the importance of this time for total export of DOC and particularly lignin. Calculated riverine DOC loads to the Arctic Ocean show an increase from previous estimates, especially when new higher discharge data are incorporated. Increased DOC loads indicate decreased residence times for terrigenous DOM in the Arctic Ocean with important implications for the reactivity and export of this material to the Atlantic Ocean.

Atmospheric processes and observations

The central role of the atmosphere in abrupt climate change is proposed and discussed. This discussion is given in the context of the poleward transport of energy in the climate system and of climate variability and change. A number of examples based on observational and model data are used to illustrate the ideas.

A kinematically distorted flux rope model for magnetic clouds

Constant-α force-free magnetic flux rope models have proven to be a valuable first step toward understanding the global context of in situ observations of magnetic clouds. However, cylindrical symmetry is necessarily assumed when using such models, and it is apparent from both observations and modeling that magnetic clouds have highly noncircular cross sections. A number of approaches have been adopted to relax the circular cross section approximation: frequently, the cross-sectional shape is allowed to take an arbitrarily chosen shape (usually elliptical), increasing the number of free parameters that are fit between data and model. While a better “fit” may be achieved in terms of reducing the mean square error between the model and observed magnetic field time series, it is not always clear that this translates to a more accurate reconstruction of the global structure of the magnetic cloud. We develop a new, noncircular cross section flux rope model that is constrained by observations of CMEs/ICMEs and knowledge of the physical processes acting on the magnetic cloud: The magnetic cloud is assumed to initially take the form of a force-free flux rope in the low corona but to be subsequently deformed by a combination of axis-centered self-expansion and heliocentric radial expansion. The resulting analytical solution is validated by fitting to artificial time series produced by numerical MHD simulations of magnetic clouds and shown to accurately reproduce the global structure.

An event-based approach to validating solar wind speed predictions: high-speed enhancements in the Wang-Sheeley-Arge model

One of the primary goals of the Center for Integrated Space Weather Modeling (CISM) effort is to assess and improve prediction of the solar wind conditions in near‐Earth space, arising from both quasi‐steady and transient structures. We compare 8 years of L1 in situ observations to predictions of the solar wind speed made by the Wang‐Sheeley‐Arge (WSA) empirical model. The mean‐square error (MSE) between the observed and model predictions is used to reach a number of useful conclusions: there is no systematic lag in the WSA predictions, the MSE is found to be highest at solar minimum and lowest during the rise to solar maximum, and the optimal lead time for 1 AU solar wind speed predictions is found to be 3 days. However, MSE is shown to frequently be an inadequate “figure of merit” for assessing solar wind speed predictions. A complementary, event‐based analysis technique is developed in which high‐speed enhancements (HSEs) are systematically selected and associated from observed and model time series. WSA model is validated using comparisons of the number of hit, missed, and false HSEs, along with the timing and speed magnitude errors between the forecasted and observed events. Morphological differences between the different HSE populations are investigated to aid interpretation of the results and improvements to the model. Finally, by defining discrete events in the time series, model predictions from above and below the ecliptic plane can be used to estimate an uncertainty in the predicted HSE arrival times.

Parcel theory in three dimensions and the calculation of SCAPE

Slantwise convective available potential energy (SCAPE) is a measure of the degree to which the atmosphere is unstable to conditional symmetric instability (CSI). It has, until now, been defined by parcel theory in which the atmosphere is assumed to be nonevolving and balanced, that is, two-dimensional. When applying this two-dimensional theory to three-dimensional evolving flows, these assumptions can be interpreted as an implicit assumption that a timescale separation exists between a relatively rapid timescale for slantwise ascent and a slower timescale for the development of the system. An approximate extension of parcel theory to three dimensions is derived and it is shown that calculations of SCAPE based on the assumption of relatively rapid slantwise ascent can be qualitatively in error. For a case study example of a developing extratropical cyclone, SCAPE calculated along trajectories determined without assuming the existence of the timescale separation show large SCAPE values for parcels ascending from the warm sector and along the warm front. These parcels ascend into the cloud head within which there is some evidence consistent with the release of CSI from observational and model cross sections. This region of high SCAPE was not found for calculations along the relatively rapidly ascending trajectories determined by assuming the existence of the timescale separation.

Glucosinolate-accumulating S-cells in Arabidopsis leaves and flower stalks undergo programmed cell death at early stages of differentiation

Plant secondary metabolites glucosinolates (GSL) have important functions in plant resistance to herbivores and pathogens. We identified all major GSL that are accumulated in S-cells in Arabidopsis by MALDI-TOF MS, and estimated by LC-MS that the total GSL concentration in these cells is above 130 mM. The precise locations of the S-cells outside phloem bundles in rosette and cauline leaves and in flower stalks were visualised using sulphur mapping by cryo-SEM/EDX. S-cells contain up to 40% of total sulphur in flower stalk tissues. S-cells in emerging flower stalks and developing leaf tissues show typical signs of Programmed Cell Death (PCD) or apoptosis, such as chromatin condensation in the nucleus and blebbing of the membranes. TUNEL staining for DNA double strand breaks confirmed PCD in S-cells in postmeristematic tissues in the flower stalk as well as in the leaf. Our results show that S-cells in postmeristematic tissues proceed to an extreme degree of metabolic specialisation besides PCD. Accumulation and maintenance of a high concentration of GSL in these cells are accompanied by degradation of a number of cell organelles. The substantial changes in the cell composition during S-cell differentiation indicate the importance of this particular GSL-based phloem defence system. The specific anatomy of the S-cells and ability to accumulate specialised secondary metabolites is similar to that of the non-articulated laticifer cells in latex plants and thus indicates a common evolutionary origin.

The extrapolation problem and how population modeling can help

We argue that population modeling can add value to ecological risk assessment by reducing uncertainty when extrapolating from ecotoxicological observations to relevant ecological effects. We review other methods of extrapolation, ranging from application factors to species sensitivity distributions to suborganismal (biomarker and "-omics'') responses to quantitative structure activity relationships and model ecosystems, drawing attention to the limitations of each. We suggest a simple classification of population models and critically examine each model in an extrapolation context. We conclude that population models have the potential for adding value to ecological risk assessment by incorporating better understanding of the links between individual responses and population size and structure and by incorporating greater levels of ecological complexity. A number of issues, however, need to be addressed before such models are likely to become more widely used. In a science context, these involve challenges in parameterization, questions about appropriate levels of complexity, issues concerning how specific or general the models need to be, and the extent to which interactions through competition and trophic relationships can be easily incorporated.