A comparative analysis of projected impacts of climate change on river runoff from global and catchment-scale hydrological models

We present a comparative analysis of projected impacts of climate change on river runoff from two types of distributed hydrological model, a global hydrological model (GHM) and catchment-scale hydrological models (CHM). Analyses are conducted for six catchments that are global in coverage and feature strong contrasts in spatial scale as well as climatic and development conditions. These include the Liard (Canada), Mekong (SE Asia), Okavango (SW Africa), Rio Grande (Brazil), Xiangu (China) and Harper's Brook (UK). A single GHM (Mac-PDM.09) is applied to all catchments whilst different CHMs are applied for each catchment. The CHMs typically simulate water resources impacts based on a more explicit representation of catchment water resources than that available from the GHM, and the CHMs include river routing. Simulations of average annual runoff, mean monthly runoff and high (Q5) and low (Q95) monthly runoff under baseline (1961-1990) and climate change scenarios are presented. We compare the simulated runoff response of each hydrological model to (1) prescribed increases in global mean temperature from the HadCM3 climate model and (2)a prescribed increase in global-mean temperature of 2oC for seven GCMs to explore response to climate model and structural uncertainty. We find that differences in projected changes of mean annual runoff between the two types of hydrological model can be substantial for a given GCM, and they are generally larger for indicators of high and low flow. However, they are relatively small in comparison to the range of projections across the seven GCMs. Hence, for the six catchments and seven GCMs we considered, climate model structural uncertainty is greater than the uncertainty associated with the type of hydrological model applied. Moreover, shifts in the seasonal cycle of runoff with climate change are presented similarly by both hydrological models, although for some catchments the monthly timing of high and low flows differs.This implies that for studies that seek to quantify and assess the role of climate model uncertainty on catchment-scale runoff, it may be equally as feasible to apply a GHM as it is to apply a CHM, especially when climate modelling uncertainty across the range of available GCMs is as large as it currently is. Whilst the GHM is able to represent the broad climate change signal that is represented by the CHMs, we find, however, that for some catchments there are differences between GHMs and CHMs in mean annual runoff due to differences in potential evaporation estimation methods, in the representation of the seasonality of runoff, and in the magnitude of changes in extreme monthly runoff, all of which have implications for future water management issues.

The changing atmospheric water cycle in polar regions in a warmer climate

We have examined the atmospheric water cycle of both Polar Regions, pole wards of 60°N and 60°S, using the ERA-Interim re-analysis and high-resolution simulations with the ECHAM5 model for both the present and future climate based on the IPCC, A1B scenario, representative of the last three decades of the 21st century. The annual precipitation in ERA-Interim amounts to ~17000 km3 and is more or less the same in the Arctic and the Antarctic, but it is composed differently. In the Arctic the annual evaporation is some 8000 km3 but some 3000 km3 less in the Antarctica where the net horizontal transport is correspondingly larger. The net water transport of the model is more intense than in ERA-Interim, in the Arctic the difference is 2.5% and in the Antarctic it is 6.2%. Precipitation and net horizontal transport in the Arctic has a maximum in August and September. Evaporation peaks in June and July. The seasonal cycle is similar in Antarctica with the highest precipitation in the austral autumn. The largest net transport occurs at the end of the major extra-tropical storm tracks in the Northern Hemisphere such as the eastern Pacific and eastern north Atlantic. The variability of the model is virtually identical to that of the re-analysis and there are no changes in variability between the present climate and the climate at the end of the 21st century when normalized with the higher level of moisture. The changes from year to year are substantial with the 20 and 30-year records being generally too short to identify robust trends in the hydrological cycle. In the A1B climate scenario the strength of the water cycle increases by some 25% in the Arctic and by 19% in the Antarctica, as measured by annual precipitation. The increase in the net horizontal transport is 29% and 22% respectively, and the increase in evaporation correspondingly less. The net transport follows closely the Clausius-Clapeyron relation. There is 2 a minor change in the annual cycle of the Arctic atmospheric water cycle with the maximum transport and precipitation occurring later in the year. There is a small imbalance of some 4-6% between the net transport and precipitation minus evaporation. We suggest that this is mainly due to the fact the transport is calculated from instantaneous 6-hourly data while precipitation and evaporation is accumulated over a 6 hour period. The residual difference is proportionally similar for all experiments and hardly varies from year to year.

Mechanisms of burst release from pH-responsive polymeric microparticles.

Microencapsulation of drugs into preformed polymers is commonly achieved through solvent evaporation techniques or spray drying. We compared these encapsulation methods in terms of controlled drug release properties of the prepared microparticles and investigated the underlying mechanisms responsible for the “burst release” effect. Using two different pH-responsive polymers with a dissolution threshold of pH 6 (Eudragit L100 and AQOAT AS-MG), hydrocortisone, a model hydrophobic drug, was incorporated into microparticles below and above its solubility within the polymer matrix. Although, spray drying is an attractive approach due to rapid particle production and relatively low solvent waste, the oil-in-oil microencapsulation method is superior in terms of controlled drug release properties from the microparticles. Slow solvent evaporation during the oil-in-oil emulsification process allows adequate time for drug and polymer redistribution in the microparticles and reduces uncontrolled drug burst release. Electron microscopy showed that this slower manufacturing procedure generated non-porous particles whereas thermal analysis and X-ray diffractometry showed that drug loading above the solubility limit of the drug in the polymer generated excess crystalline drug on the surface of the particles. Raman spectral mapping illustrated that drug was homogeneously distributed as a solid solution in the particles when loaded below saturation in the polymer with consequently minimal burst release.

Remarkable dynamics of nanoparticles in the urban atmosphere

Nanoparticles emitted from road traffic are the largest source of respiratory exposure for the general public living in urban areas. It has been suggested that the adverse health effects of airborne particles may scale with the airborne particle number, which if correct, focuses attention on the nanoparticle (less than 100 nm) size range which dominates the number count in urban areas. Urban measurements of particle size distributions have tended to show a broadly similar pattern dominated by a mode centred on 20–30 nm diameter particles emitted by diesel engine exhaust. In this paper we report the results of measurements of particle number concentration and size distribution made in a major London park as well as on the BT Tower, 160 m high. These measurements taken during the REPARTEE project (Regents Park and BT Tower experiment) show a remarkable shift in particle size distributions with major losses of the smallest particle class as particles are advected away from the traffic source. In the Park, the traffic related mode at 20–30 nm diameter is much reduced with a new mode at <10 nm. Size distribution measurements also revealed higher number concentrations of sub-50 nm particles at the BT Tower during days affected by higher turbulence as determined by Doppler Lidar measurements and indicate a loss of nanoparticles from air aged during less turbulent conditions. These results suggest that nanoparticles are lost by evaporation, rather than coagulation processes. The results have major implications for understanding the impacts of traffic-generated particulate matter on human health.

Atmospheric conservation properties in ERA-Interim

We study the global atmospheric budgets of mass, moisture, energy and angular momentum in the latest reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA-Interim, for the period 1989–2008 and compare with ERA-40. Most of the measures we use indicate that the ERA-Interim reanalysis is superior in quality to ERA-40. In ERA-Interim the standard deviation of the monthly mean global dry mass of 0.7 kg m−2 (0.007%) is slightly worse than in ERA-40, and long time-scale variations in dry mass originate predominately in the surface pressure field. The divergent winds are improved in ERA-Interim: the global standard deviation of the time-averaged dry mass budget residual is 10 kg m−2 day−1 and the quality of the cross-equatorial mass fluxes is improved. The temporal variations in the global evaporation minus precipitation (E − P) are too large but the global moisture budget residual is 0.003 kg m−2 day−1 with a spatial standard deviation of 0.3 kg m−2 day−1. Both the E − P over ocean and P − E over land are about 15% larger than the 1.1 Tg s−1 transport of water from ocean to land. The top of atmosphere (TOA) net energy losses are improved, with a value of 1 W m−2, but the meridional gradient of the TOA net energy flux is smaller than that from the Clouds and the Earth's Radiant Energy System (CERES) data. At the surface the global energy losses are worse, with a value of 7 W m−2. Over land however, the energy loss is only 0.5 W m−2. The downwelling thermal radiation at the surface in ERA-Interim of 341 W m−2 is towards the higher end of previous estimates. The global mass-adjusted energy budget residual is 8 W m−2 with a spatial standard deviation of 11 W m−2, and the mass-adjusted atmospheric energy transport from low to high latitudes (the sum for the two hemispheres) is 9.5 PW

The role of water vapour in Earth’s energy flows

Water vapour modulates energy flows in Earth's climate system through transfer of latent heat by evaporation and condensation and by modifying the flows of radiative energy both in the longwave and shortwave portions of the electromagnetic spectrum. This article summarizes the role of water vapour in Earth's energy flows with particular emphasis on (1) the powerful thermodynamic constraint of the Clausius Clapeyron equation, (2) dynamical controls on humidity above the boundary layer (or free-troposphere), (3) uncertainty in continuum absorption in the relatively transparent "window" regions of the radiative spectrum and (4) implications for changes in the atmospheric hydrological cycle.

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.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

The impact of a seasonally ice free Arctic Ocean on the temperature, precipitation and surface mass balance of Svalbard

The observed decline in summer sea ice extent since the 1970s is predicted to continue until the Arctic Ocean is seasonally ice free during the 21st Century. This will lead to a much perturbed Arctic climate with large changes in ocean surface energy flux. Svalbard, located on the present day sea ice edge, contains many low lying ice caps and glaciers and is expected to experience rapid warming over the 21st Century. The total sea level rise if all the land ice on Svalbard were to melt completely is 0.02 m. The purpose of this study is to quantify the impact of climate change on Svalbard’s surface mass balance (SMB) and to determine, in particular, what proportion of the projected changes in precipitation and SMB are a result of changes to the Arctic sea ice cover. To investigate this a regional climate model was forced with monthly mean climatologies of sea surface temperature (SST) and sea ice concentration for the periods 1961–1990 and 2061–2090 under two emission scenarios. In a novel forcing experiment, 20th Century SSTs and 21st Century sea ice were used to force one simulation to investigate the role of sea ice forcing. This experiment results in a 3.5 m water equivalent increase in Svalbard’s SMB compared to the present day. This is because over 50 % of the projected increase in winter precipitation over Svalbard under the A1B emissions scenario is due to an increase in lower atmosphere moisture content associated with evaporation from the ice free ocean. These results indicate that increases in precipitation due to sea ice decline may act to moderate mass loss from Svalbard’s glaciers due to future Arctic warming.

Hydrological cycle in the Danube basin in present-day and XXII century simulations by IPCCAR4 global climate models

We present an intercomparison and verification analysis of 20 GCMs (Global Circulation Models) included in the 4th IPCC assessment report regarding their representation of the hydrological cycle on the Danube river basin for 1961–2000 and for the 2161–2200 SRESA1B scenario runs. The basin-scale properties of the hydrological cycle are computed by spatially integrating the precipitation, evaporation, and runoff fields using the Voronoi-Thiessen tessellation formalism. The span of the model- simulated mean annual water balances is of the same order of magnitude of the observed Danube discharge of the Delta; the true value is within the range simulated by the models. Some land components seem to have deficiencies since there are cases of violation of water conservation when annual means are considered. The overall performance and the degree of agreement of the GCMs are comparable to those of the RCMs (Regional Climate Models) analyzed in a previous work, in spite of the much higher resolution and common nesting of the RCMs. The reanalyses are shown to feature several inconsistencies and cannot be used as a verification benchmark for the hydrological cycle in the Danubian region. In the scenario runs, for basically all models the water balance decreases, whereas its interannual variability increases. Changes in the strength of the hydrological cycle are not consistent among models: it is confirmed that capturing the impact of climate change on the hydrological cycle is not an easy task over land areas. Moreover, in several cases we find that qualitatively different behaviors emerge among the models: the ensemble mean does not represent any sort of average model, and often it falls between the models’ clusters.

Incorporation of water vapour transfer in the JULES Land Surface Model: implications for key soil variables and land surface fluxes

This study focuses on the mechanisms underlying water and heat transfer in upper soil layers, and their effects on soil physical prognostic variables and the individual components of the energy balance. The skill of the JULES (Joint UK Land Environment Simulator) land surface model (LSM) to simulate key soil variables, such as soil moisture content and surface temperature, and fluxes such as evaporation, is investigated. The Richards equation for soil water transfer, as used in most LSMs, was updated by incorporating isothermal and thermal water vapour transfer. The model was tested for three sites representative of semi-arid and temperate arid climates: the Jornada site (New Mexico, USA), Griffith site (Australia) and Audubon site (Arizona, USA). Water vapour flux was found to contribute significantly to the water and heat transfer in the upper soil layers. This was mainly due to isothermal vapour diffusion; thermal vapour flux also played a role at the Jornada site just after rainfall events. Inclusion of water vapour flux had an effect on the diurnal evolution of evaporation, soil moisture content and surface temperature. The incorporation of additional processes, such as water vapour flux among others, into LSMs may improve the coupling between the upper soil layers and the atmosphere, which in turn could increase the reliability of weather and climate predictions.

The impact of 21st century Arctic sea ice decline on the climate and mass balance of the Greenland ice sheet and Svalbard’s glaciers and ice caps

The Arctic is a region particularly susceptible to rapid climate change. General circulation models (GCMs) suggest a polar amplification of any global warming signal by a factor of about 1.5 due, in part, to sea ice feedbacks. The dramatic recent decline in multi-year sea ice cover lies outside the standard deviation of the CMIP3 ensemble GCM predictions. Sea ice acts as a barrier between cold air and warmer oceans during winter, as well as inhibiting evaporation from the ocean surface water during the summer. An ice free Arctic would likely have an altered hydrological cycle with more evaporation from the ocean surface leading to changes in precipitation distribution and amount. Using the U.K. Met Office Regional Climate Model (RCM), HadRM3, the atmospheric effects of the observed and projected reduction in Arctic sea ice are investigated. The RCM is driven by the atmospheric GCM HadAM3. Both models are forced with sea surface temperature and sea ice for the period 2061-2090 from the CMIP3 HadGEM1 experiments. Here we use an RCM at 50km resolution over the Arctic and 25km over Svalbard, which captures well the present-day pattern of precipitation and provides a detailed picture of the projected changes in the behaviour of the oceanic-atmosphere moisture fluxes and how they affect precipitation. These experiments show that the projected 21stCentury sea ice decline alone causes large impacts to the surface mass balance (SMB) on Svalbard. However Greenland’s SMB is not significantly affected by sea ice decline alone, but responds with a strongly negative shift in SMB when changes to SST are incorporated into the experiments. This is the first study to characterise the impact of changes in future sea ice to Arctic terrestrial cryosphere mass balance.

Post-frying oil drainage from potato chips and french fries: a comparative study of atmospheric and vacuum drainage

This paper represents a study of the transient changes occurring in temperature, and moisture and oil contents during the so called “post-frying drainage”—which is the duration for which a product is held in the head space of the fryer after it is removed from the oil. Since most of the oil adhering to the product penetrates into the structure during this period, this paper examines the effects of applying vacuum during drainage (1.33 kPa) to maintain the product temperature consistently above the water saturation temperature corresponding to the prevailing pressure (11 °C), which potentially eliminates water condensation and prevents the occluded surface oil from penetrating into the product structure. Draining under vacuum significantly lowers the oil content of potato chips by 38% compared to atmospheric drainage. This phenomenon can be further confirmed by confocal laser scanning microscopy (CLSM) images, which show that the boundary between the core and the crust regions is clearly visible in the case of vacuum drainage, whereas in the case of atmospheric drainage, the oil is distributed throughout the structure. Unfortunately, the same approach did not reduce the oil content of French fries—the oil content of vacuum-drained product was found similar to the product obtained by draining under atmospheric pressure. This is because the reduction in oil content only occurs when there is net moisture evaporation from the product and the evaporation rate is sufficient to force out the oil from the product; this was clearly not the case with French fries. The CLSM images show that the oil distribution in the products drained under atmospheric pressure and vacuum was similar.

Under what conditions does embedded convection enhance orographic precipitation?

Idealised convection-permitting simulations are used to quantify the impact of embedded convection on the precipitation generated by moist flow over midlatitude mountain ridges. A broad range of mountain dimensions and moist stabilities are considered to encompass a spectrum of physically plausible flows. The simulations reveal that convection only enhances orographic precipitation in cap clouds that are otherwise unable to efficiently convert cloud condensate into precipitate. For tall and wide mountains (e.g. the Washington Cascades or the southern Andes), precipitate forms efficiently through vapour deposition and collection, even in the absence of embedded convection. When embedded convection develops in such clouds, it produces competing effects (enhanced condensation in updraughts and enhanced evaporation through turbulent mixing and compensating subsidence) that cancel to yield little net change in precipitation. By contrast, convection strongly enhances precipitation over short and narrow mountains (e.g. the UK Pennines or the Oregon Coastal Range) where precipitation formation is otherwise highly inefficient. Although cancellation between increased condensation and evaporation still occurs, the enhanced precipitation formation within the convective updraughts leads to a net increase in precipitation efficiency. The simulations are physically interpreted through non-dimensional diagnostics and relevant time-scales that govern advective, microphysical, and convective processes.

Faunal evidence of a Holocene pluvial phase in Southern Arabia with remarks on the morphological variability of Helenina anderseni

Although foraminifera have been found living in inland saline lakes isolated from the sea, this phenomenon has rarely been recognized in the fossil record. This study documents the occurrence of benthic foraminifera in Holocene lake sediments located nearly 500 km inland from the Red Sea, in the Al-Mundafan region of southern Saudi Arabia. The lake formed during a regional pluvial period, 10,500–6000 yr BP. The presence of foraminifera and brackish charophytes in the studied section represent an interval when the lake was slightly brackish due to high evaporation. The studied sediments yielded a bispecific benthic foraminiferal fauna comprised of Helenina anderseni and Trichohyalus aguayoi, as well as the brackish charophyte genus Lamprothamnium. The benthic foraminifera are species characteristic of mangrove swamps, salt marshes, and lagoons, which are environments currently widespread along the Red Sea coasts. Because the Al Mundafan area was never connected to the sea during the Quaternary, wading birds must have been the vector that transported the foraminifera to the paleolake.