Diagnostic analysis of atmospheric moisture and clear-sky radiative feedback in the Hadley Centre and Geophysical Fluid Dynamics Laboratory (GFDL) climate models

The interannual variability of the hydrological cycle is diagnosed from the Hadley Centre and Geophysical Fluid Dynamics Laboratory (GFDL) climate models, both of which are forced by observed sea surface temperatures. The models produce a similar sensitivity of clear-sky outgoing longwave radiation to surface temperature of ∼2 W m−2 K−1, indicating a consistent and positive clear-sky radiative feedback. However, differences between changes in the temperature lapse-rate and the height dependence of moisture fluctuations suggest that contrasting mechanisms bring about this result. The GFDL model appears to give a weaker water vapor feedback (i.e., changes in specific humidity). This is counteracted by a smaller upper tropospheric temperature response to surface warming, which implies a compensating positive lapse-rate feedback.

Exploitation of Geostationary Earth Radiation Budget data using simulations from a numerical weather prediction model: Methodology and data validation

We describe a new methodology for comparing satellite radiation budget data with a numerical weather prediction (NWP) model. This is applied to data from the Geostationary Earth Radiation Budget (GERB) instrument on Meteosat-8. The methodology brings together, in near-real time, GERB broadband shortwave and longwave fluxes with simulations based on analyses produced by the Met Office global NWP model. Results for the period May 2003 to February 2005 illustrate the progressive improvements in the data products as various initial problems were resolved. In most areas the comparisons reveal systematic errors in the model's representation of surface properties and clouds, which are discussed elsewhere. However, for clear-sky regions over the oceans the model simulations are believed to be sufficiently accurate to allow the quality of the GERB fluxes themselves to be assessed and any changes in time of the performance of the instrument to be identified. Using model and radiosonde profiles of temperature and humidity as input to a single-column version of the model's radiation code, we conduct sensitivity experiments which provide estimates of the expected model errors over the ocean of about ±5–10 W m−2 in clear-sky outgoing longwave radiation (OLR) and ±0.01 in clear-sky albedo. For the more recent data the differences between the observed and modeled OLR and albedo are well within these error estimates. The close agreement between the observed and modeled values, particularly for the most recent period, illustrates the value of the methodology. It also contributes to the validation of the GERB products and increases confidence in the quality of the data, prior to their release.

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.

Orographic precipitation in the tropics: experiments in Dominica

The “natural laboratory” of mountainous Dominica (15°N) in the trade wind belt is used to study the physics of tropical orographic precipitation in its purest form, unforced by weather disturbances or by the diurnal cycle of solar heating. A cross-island line of rain gauges and 5-min radar scans from Guadeloupe reveal a large annual precipitation at high elevation (7 m yr^{−1}) and a large orographic enhancement factor (2 to 8) caused primarily by repetitive convective triggering over the windward slope. The triggering is caused by terrain-forced lifting of the conditionally unstable trade wind cloud layer. Ambient humidity fluctuations associated with open-ocean convection may play a key role. The convection transports moisture upward and causes frequent brief showers on the hilltops. The drying ratio of the full air column from precipitation is less than 1% whereas the surface air dries by about 17% from the east coast to the mountain top. On the lee side, a plunging trade wind inversion and reduced instability destroys convective clouds and creates an oceanic rain shadow.

Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

The diurnal cycle of convection and atmospheric tides in an aquaplanet

The diurnal cycle of tropical convection and its relationship to the atmospheric tides is investigated using an aquaplanet GCM. The diurnal and semidiurnal harmonics of precipitation are both found to contribute significantly to the total diurnal variability of precipitation in the model, which is broadly consistent with observations of the diurnal cycle of convection over the open ocean. The semidiurnal tide is found to be the dominant forcing for the semidiurnal harmonic of precipitation. In contrast the diurnal tide plays only a small role in forcing the diurnal harmonic of precipitation, which is dominated by the variations in shortwave and longwave heating. In both the diurnal and semidiurnal harmonics, the feedback onto the convection by the humidity tendencies due to the convection is found to be important in determining the phase of the harmonics. Further experiments show that the diurnal cycle of precipitation is sensitive to the choice of closure in the convection scheme. While the surface pressure signal of the simulated atmospheric tides in the model agree well with both theory and observations in their magnitude and phase, sensitivity experiments suggest that the role of the stratospheric ozone in forcing the semidiurnal tide is much reduced compared to theoretical predictions. Furthermore, the influence of the cloud radiative effects seems small. It is suggested that the radiative heating profile in the troposphere, associated primarily with the water vapor distribution, is more important than previously thought for driving the semidiurnal tide. However, this result may be sensitive to the vertical resolution and extent of the model.

The cumulus-capped boundary layer. II: Interface fluxes

This paper considers the relationship between the mean temperature and humidity profiles and the fluxes of heat and moisture at cloud base and the base of the inversion in the cumulus-capped boundary layer. The relationships derived are based on an approximate form of the scalar-flux budget and the scaling properties of the turbulent kinetic energy (TKE) budget. The scalar-flux budget gives a relationship between the change in the virtual potential temperature across either the cloud base transition zone or the inversion and the flux at the base of the layer. The scaling properties of the TKE budget lead to a relationship between the heat and moisture fluxes and the mean subsaturation through the liquid-water flux. The 'jump relation' for the virtual potential temperature at cloud base shows the close connection between the cumulus mass flux in the cumulus-capped boundary layer and the entrainment velocity in the dry-convective boundary layer. Gravity waves are shown to be an important feature of the inversion.

Observations of the depth of ice particle evaporation beneath frontal cloud to improve NWP modelling

The evaporation (sublimation) of ice particles beneath frontal ice cloud can provide a significant source of diabatic cooling which can lead to enhanced slantwise descent below the frontal surface. The strength and vertical extent of the cooling play a role in determining the dynamic response of the atmosphere, and an adequate representation is required in numerical weather-prediction (NWP) models for accurate forecasts of frontal dynamics. In this paper, data from a vertically pointing 94 GHz radar are used to determine the characteristic depth-scale of ice particle sublimation beneath frontal ice cloud. A statistical comparison is made with equivalent data extracted from the NWP mesoscale model operational at the Met Office, defining the evaporation depth-scale as the distance for the ice water content to fall to 10% of its peak value in the cloud. The results show that the depth of the ice evaporation zone derived from observations is less than 1 km for 90% of the time. The model significantly overestimates the sublimation depth-scales by a factor of between two and three, and underestimates the local ice water content by a factor of between two and four. Consequently the results suggest the model significantly underestimates the strength of the evaporative cooling, with implications for the prediction of frontal dynamics. A number of reasons for the model discrepancy are suggested. A comparison with radiosonde relative humidity data suggests part of the overestimation in evaporation depth may be due to a high RH bias in the dry slot beneath the frontal cloud, but other possible reasons include poor vertical resolution and deficiencies in the evaporation rate or ice particle fall-speed parametrizations.

Can climate models capture the structure of extratropical cyclones?

Composites of wind speeds, equivalent potential temperature, mean sea level pressure, vertical velocity, and relative humidity have been produced for the 100 most intense extratropical cyclones in the Northern Hemisphere winter for the 40-yr ECMWF Re-Analysis (ERA-40) and the high resolution global environment model (HiGEM). Features of conceptual models of cyclone structure—the warm conveyor belt, cold conveyor belt, and dry intrusion—have been identified in the composites from ERA-40 and compared to HiGEM. Such features can be identified in the composite fields despite the smoothing that occurs in the compositing process. The surface features and the three-dimensional structure of the cyclones in HiGEM compare very well with those from ERA-40. The warm conveyor belt is identified in the temperature and wind fields as a mass of warm air undergoing moist isentropic uplift and is very similar in ERA-40 and HiGEM. The rate of ascent is lower in HiGEM, associated with a shallower slope of the moist isentropes in the warm sector. There are also differences in the relative humidity fields in the warm conveyor belt. In ERA-40, the high values of relative humidity are strongly associated with the moist isentropic uplift, whereas in HiGEM these are not so strongly associated. The cold conveyor belt is identified as rearward flowing air that undercuts the warm conveyor belt and produces a low-level jet, and is very similar in HiGEM and ERA-40. The dry intrusion is identified in the 500-hPa vertical velocity and relative humidity. The structure of the dry intrusion compares well between HiGEM and ERA-40 but the descent is weaker in HiGEM because of weaker along-isentrope flow behind the composite cyclone. HiGEM’s ability to represent the key features of extratropical cyclone structure can give confidence in future predictions from this model.

Quality of medicines stored together in multi-compartment compliance aids

Background and Objective: Dispensing medicines into compliance aids is a common practice in pharmacy contrary to manufacturers’ advice and studies have shown the appearance of light-sensitive tablets is compromised by such storage; we previously found evidence of reduced bioavailability at elevated temperature and humidity. Our objective was to examine the physicochemical stability of two generic atenolol tablets in different compliance aids and with aspirin co-storage at room temperature and at 40 °C/75% relative humidity. Methods: The physicochemical stability of atenolol tablets was evaluated after 28 days of storage and compared with controls by examining visual appearance, weight, disintegration, dissolution, friability and hardness to accepted standards and using a previously validated HPLC method for chemical assay. Results and Discussion: The response to storage was brand-dependent and not straightforward. With one make of atenolol (Alpharma), storage in compliance aids even at room temperature impacted on physical stability, reducing tablet hardness, with storage in Dosett® exerting a greater impact than storage in Medidos® (t-test P < 0·001). Co-storage at elevated temperature and humidity also impacted on the appearance of non-coated aspirin tablets (Angette™). The chemical stability of atenolol was not affected and we did not find evidence of changes to bioavailability with either make. Certainly data for one atenolol make (CP Pharmaceuticals) co-stored with aspirin (Angette™ and Nu-Seals) in both compliance aids at room temperature provided evidence of short-term stability. But medicines are dispensed into compliance aids in multi-factorial ways so our study highlights not only the lack of evidence but also a realization that evidence to support real practice may not be accomplished through research. Conclusion: Reassuring practitioners of the continued stability of medicines in compliance aids under the countless condition in which they are dispensed in practice may requires a different approach involving medical device regulators and more definitive professional guidance.

Pilot study of the short-term physico-chemical stability of atenolol tablets stored in a multi-compartment compliance aid

Study objectives: There is a possibility that lower air, moisture and light protection could impact on physico-chemical stability of medicines inside multi-compartment compliance aids (MCCAs), although this has not yet been proved. The objectives of the study were to examine the physico-chemical stability of atenolol tablets stored in a compliance aid at room temperature, and at elevated temperature and humidity to simulate practice conditions. Methods: Atenolol 100 mg tablets in 28-chamber, plastic compliance aids with transparent lids were stored for four weeks at room temperature and at 40°C with 75% relative humidity. Tablets were also stored at room temperature in original packaging and Petri dishes. Physical tests were conducted to standards as laid down in the British Pharmacopoeia 2005, and dissolution to those of the United States Pharmacopoeia volume 24. Chemical stability was assessed by a validated high-performance liquid chromatography (HPLC) method. Results: Tablets at room temperature in original packaging, in compliance aids and Petri dishes remained the same in appearance and passed physico-chemical tests. Tablets exposed to 40°C with 75% relative humidity in compliance aids passed tests for uniformity of weight, friability and chemical stability but became pale and moist, softer (82 newtons ± 4; p< 0.0001) than tablets in the original packaging (118 newtons ± 6), more friable (0.14% loss of mass) compared with other tablets (0.005%), and failed the tests for disintegration (>15 minutes) and dissolution (only 15% atenolol released at 30 minutes). Conclusion: Although chemical stability was unaffected, storage in compliance aids at 40°C with 75% relative humidity softened atenolol tablets, prolonged disintegration time and hindered dissolution which could significantly reduce bioavailability. This formulation could be suitable for storage in compliance aids at 25°C, but not in hotter, humid weather.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Hydrolysis of N2O5 on sub-micron mineral salt aerosols

The kinetics of reactive uptake of gaseous N2O5 on sub-micron aerosol particles composed of aqueous ammonium sulfate, ammonium hydrogensulfate and sodium nitrate has been investigated. Uptake was measured in a laminar flow reactor, coupled with a differential mobility analyser (DMA) to obtain the aerosol size distribution, with N2O5 detection using NO chemiluminescence. FTIR spectroscopy was used to obtain information about the composition and water content of the aerosol particles under the conditions used in the kinetic measurements. The aerosols were generated by the nebulisation of aqueous salt solutions. The uptake coefficient on the sulfate salts was in the range [gamma]=0.0015 to 0.033 depending on temperature, humidity and phase of the aerosol. On sodium nitrate aerosols the values were much lower, [gamma]<0.001, confirming the inhibition of N2O5 hydrolysis by nitrate ions. At high humidity (>50% r.h.) the uptake coefficient on liquid sulfate aerosols is independent of water content, but at lower humidity, especially below the efflorescence point, the reactivity of the aerosol declines, correlating with the lower water content. The lower uptake rate on solid aerosols may be due to limitations imposed by the liquid volume in the particles. Uptake on sulfate aerosols showed a negative temperature dependence at T>290 K but no significant temperature dependence at lower temperatures. The results are generally consistent with previous models of N2O5 hydrolysis where the reactive intermediate is NO2+ produced by autoionisation of nitrogen pentoxide in the condensed phase.

Reactive uptake coefficients for heterogeneous reaction of N2O5 with submicron aerosols of NaCl and natural sea salt

The kinetics of uptake of gaseous N2O5 on submicron aerosols containing NaCl and natural sea salt have been investigated in a flow reactor as a function of relative humidity (RH) in the range 30-80% at 295±2K and a total pressure of 1bar. The measured uptake coefficients, γ, were larger on the aerosols containing sea salt compared to those of pure NaCl, and in both cases increased with increasing RH. These observations are explained in terms of the variation in the size of the salt droplets, which leads to a limitation in the uptake rate into small particles. After correction for this effect the uptake coefficients are independent of relative humidity, and agree with those measured previously on larger droplets. A value of γ=0.025 is recommended for the reactive uptake coefficient for N2O5 on deliquesced sea salt droplets at 298K and RH>40%.

Stratospheric water vapour and high climate sensitivity in a version of the HadSM3 climate model

It has been shown previously that one member of the Met Office Hadley Centre single-parameter perturbed physics ensemble – the so-called "low entrainment parameter" member – has a much higher climate sensitivity than other individual parameter perturbations. Here we show that the concentration of stratospheric water vapour in this member is over three times higher than observations, and, more importantly for climate sensitivity, increases significantly when climate warms. The large surface temperature response of this ensemble member is more consistent with stratospheric humidity change, rather than upper tropospheric clouds as has been previously suggested. The direct relationship between the bias in the control state (elevated stratospheric humidity) and the cause of the high climate sensitivity (a further increase in stratospheric humidity) lends further doubt as to the realism of this particular integration. This, together with other evidence, lowers the likelihood that the climate system's physical sensitivity is significantly higher than the likely upper range quoted in the Intergovernmental Panel on Climate Change's Fourth Assessment Report.