Aerosol impacts on drizzle properties in warm clouds from ARM Mobile Facility maritime and continental deployments

We have extensively evaluated the response of cloud-base drizzle rate (Rcb; mm day–1) in warm clouds to liquid water path (LWP; g m–2) and to cloud condensation nuclei (CCN) number concentration (NCCN; cm–3), an aerosol proxy. This evaluation is based on a 19-month long dataset of Doppler radar, lidar, microwave radiometers and aerosol observing systems from the Atmospheric Radiation Measurement (ARM) Mobile Facility deployments at the Azores and in Germany. Assuming 0.55% supersaturation to calculate NCCN, we found a power law , indicating that Rcb decreases by a factor of 2–3 as NCCN increases from 200 to 1000 cm–3 for fixed LWP. Additionally, the precipitation susceptibility to NCCN ranges between 0.5 and 0.9, in agreement with values from simulations and aircraft measurements. Surprisingly, the susceptibility of the probability of precipitation from our analysis is much higher than that from CloudSat estimates, but agrees well with simulations from a multi-scale high-resolution aerosol-climate model. Although scale issues are not completely resolved in the intercomparisons, our results are encouraging, suggesting that it is possible for multi-scale models to accurately simulate the response of LWP to aerosol perturbations.

Evaluation of seven European aerosol optical depth retrieval algorithms for climate analysis

Satellite data are increasingly used to provide observation-based estimates of the effects of aerosols on climate. The Aerosol-cci project, part of the European Space Agency's Climate Change Initiative (CCI), was designed to provide essential climate variables for aerosols from satellite data. Eight algorithms, developed for the retrieval of aerosol properties using data from AATSR (4), MERIS (3) and POLDER, were evaluated to determine their suitability for climate studies. The primary result from each of these algorithms is the aerosol optical depth (AOD) at several wavelengths, together with the Ångström exponent (AE) which describes the spectral variation of the AOD for a given wavelength pair. Other aerosol parameters which are possibly retrieved from satellite observations are not considered in this paper. The AOD and AE (AE only for Level 2) were evaluated against independent collocated observations from the ground-based AERONET sun photometer network and against “reference” satellite data provided by MODIS and MISR. Tools used for the evaluation were developed for daily products as produced by the retrieval with a spatial resolution of 10 × 10 km2 (Level 2) and daily or monthly aggregates (Level 3). These tools include statistics for L2 and L3 products compared with AERONET, as well as scoring based on spatial and temporal correlations. In this paper we describe their use in a round robin (RR) evaluation of four months of data, one month for each season in 2008. The amount of data was restricted to only four months because of the large effort made to improve the algorithms, and to evaluate the improvement and current status, before larger data sets will be processed. Evaluation criteria are discussed. Results presented show the current status of the European aerosol algorithms in comparison to both AERONET and MODIS and MISR data. The comparison leads to a preliminary conclusion that the scores are similar, including those for the references, but the coverage of AATSR needs to be enhanced and further improvements are possible for most algorithms. None of the algorithms, including the references, outperforms all others everywhere. AATSR data can be used for the retrieval of AOD and AE over land and ocean. PARASOL and one of the MERIS algorithms have been evaluated over ocean only and both algorithms provide good results.

The three-dimensional morphology of simulated and observed convective storms over southern England

A set of high-resolution radar observations of convective storms has been collected to evaluate such storms in the UK Met Office Unified Model during the DYMECS project (Dynamical and Microphysical Evolution of Convective Storms). The 3-GHz Chilbolton Advanced Meteorological Radar was set up with a scan-scheduling algorithm to automatically track convective storms identified in real-time from the operational rainfall radar network. More than 1,000 storm observations gathered over fifteen days in 2011 and 2012 are used to evaluate the model under various synoptic conditions supporting convection. In terms of the detailed three-dimensional morphology, storms in the 1500-m grid-length simulations are shown to produce horizontal structures a factor 1.5–2 wider compared to radar observations. A set of nested model runs at grid lengths down to 100m show that the models converge in terms of storm width, but the storm structures in the simulations with the smallest grid lengths are too narrow and too intense compared to the radar observations. The modelled storms were surrounded by a region of drizzle without ice reflectivities above 0 dBZ aloft, which was related to the dominance of ice crystals and was improved by allowing only aggregates as an ice particle habit. Simulations with graupel outperformed the standard configuration for heavy-rain profiles, but the storm structures were a factor 2 too wide and the convective cores 2 km too deep.

Instantaneous gelation in Smoluchowski’s coagulation equation revisited

We study the solutions of the Smoluchowski coagulation equation with a regularization term which removes clusters from the system when their mass exceeds a specified cutoff size, M. We focus primarily on collision kernels which would exhibit an instantaneous gelation transition in the absence of any regularization. Numerical simulations demonstrate that for such kernels with monodisperse initial data, the regularized gelation time decreasesas M increases, consistent with the expectation that the gelation time is zero in the unregularized system. This decrease appears to be a logarithmically slow function of M, indicating that instantaneously gelling kernels may still be justifiable as physical models despite the fact that they are highly singular in the absence of a cutoff. We also study the case when a source of monomers is introduced in the regularized system. In this case a stationary state is reached. We present a complete analytic description of this regularized stationary state for the model kernel, K(m1,m2)=max{m1,m2}ν, which gels instantaneously when M→∞ if ν>1. The stationary cluster size distribution decays as a stretched exponential for small cluster sizes and crosses over to a power law decay with exponent ν for large cluster sizes. The total particle density in the stationary state slowly vanishes as [(ν−1)logM]−1/2 when M→∞. The approach to the stationary state is nontrivial: Oscillations about the stationary state emerge from the interplay between the monomer injection and the cutoff, M, which decay very slowly when M is large. A quantitative analysis of these oscillations is provided for the addition model which describes the situation in which clusters can only grow by absorbing monomers.

The Holocene temperature conundrum

A recent temperature reconstruction of global annual temperature shows Early Holocene warmth followed by a cooling trend through the Middle to Late Holocene [Marcott SA, et al., 2013, Science 339(6124):1198–1201]. This global cooling is puzzling because it is opposite from the expected and simulated global warming trend due to the retreating ice sheets and rising atmospheric greenhouse gases. Our critical reexamination of this contradiction between the reconstructed cooling and the simulated warming points to potentially significant biases in both the seasonality of the proxy reconstruction and the climate sensitivity of current climate models.

A novel ensemble method for retrieving properties of warm cloud in 3-D using ground-based scanning radar and zenith radiances

We present a novel method for retrieving high-resolution, three-dimensional (3-D) nonprecipitating cloud fields in both overcast and broken-cloud situations. The method uses scanning cloud radar and multiwavelength zenith radiances to obtain gridded 3-D liquid water content (LWC) and effective radius (re) and 2-D column mean droplet number concentration (Nd). By using an adaption of the ensemble Kalman filter, radiances are used to constrain the optical properties of the clouds using a forward model that employs full 3-D radiative transfer while also providing full error statistics given the uncertainty in the observations. To evaluate the new method, we first perform retrievals using synthetic measurements from a challenging cumulus cloud field produced by a large-eddy simulation snapshot. Uncertainty due to measurement error in overhead clouds is estimated at 20% in LWC and 6% in re, but the true error can be greater due to uncertainties in the assumed droplet size distribution and radiative transfer. Over the entire domain, LWC and re are retrieved with average error 0.05–0.08 g m-3 and ~2 μm, respectively, depending on the number of radiance channels used. The method is then evaluated using real data from the Atmospheric Radiation Measurement program Mobile Facility at the Azores. Two case studies are considered, one stratocumulus and one cumulus. Where available, the liquid water path retrieved directly above the observation site was found to be in good agreement with independent values obtained from microwave radiometer measurements, with an error of 20 g m-2.

Challenges in quantifying changes in the global water cycle

Understanding observed changes to the global water cycle is key to predicting future climate changes and their impacts. While many datasets document crucial variables such as precipitation, ocean salinity, runoff, and humidity, most are uncertain for determining long-term changes. In situ networks provide long time-series over land but are sparse in many regions, particularly the tropics. Satellite and reanalysis datasets provide global coverage, but their long-term stability is lacking. However, comparisons of changes among related variables can give insights into the robustness of observed changes. For example, ocean salinity, interpreted with an understanding of ocean processes, can help cross-validate precipitation. Observational evidence for human influences on the water cycle is emerging, but uncertainties resulting from internal variability and observational errors are too large to determine whether the observed and simulated changes are consistent. Improvements to the in situ and satellite observing networks that monitor the changing water cycle are required, yet continued data coverage is threatened by funding reductions. Uncertainty both in the role of anthropogenic aerosols, and due to large climate variability presently limits confidence in attribution of observed changes.

Ocean heat uptake processes: a model intercomparison

We compare the quasi-equilibrium heat balances, as well as their responses to 4×CO2 perturbation, among three global climate models with the aim to identify and explain inter-model differences in ocean heat uptake (OHU) processes. We find that, in quasi-equilibrium, convective and mixed layer processes, as well as eddy-related processes, cause cooling of the subsurface ocean. The cooling is balanced by warming caused by advective and diapycnally diffusive processes. We also find that in the CO2-perturbed climates the largest contribution to OHU comes from changes in vertical mixing processes and the mean circulation, particularly in the extra-tropics, caused both by changes in wind forcing, and by changes in high-latitude buoyancy forcing. There is a substantial warming in the tropics, a significant part of which occurs because of changes in horizontal advection in extra-tropics. Diapycnal diffusion makes only a weak contribution to the OHU, mainly in the tropics, due to increased stratification. There are important qualitative differences in the contribution of eddy-induced advection and isopycnal diffusion to the OHU among the models. The former is related to the different values of the coefficients used in the corresponding scheme. The latter is related to the different tapering formulations of the isopycnal diffusion scheme. These differences affect the OHU in the deep ocean, which is substantial in two of the models, the dominant region of deep warming being the Southern Ocean. However, most of the OHU takes place above 2000 m, and the three models are quantitatively similar in their global OHU efficiency and its breakdown among processes and as a function of latitude.

The impact of salinity perturbations on the future uptake of heat by the Atlantic Ocean

Anthropogenic ocean heat uptake is a key factor in determining climate change and sea-level rise. There is considerable uncertainty in projections of freshwater forcing of the ocean, with the potential to influence ocean heat uptake. We investigatethis by adding either -0.1 Sv or +0.1 Sv freshwater to the Atlantic in global climate model simulations, simultaneously imposing an atmospheric CO2 increase. The resulting changes in the Atlantic meridional overturning circulation are roughly equal and opposite (±2Sv). The impact of the perturbation on ocean heat content is more complex, although it is relatively small (~5%) compared to the total anthropogenic heat uptake. Several competing processes either accelerate or retard warming at different depths. Whilst positive freshwater perturbations cause an overall heating of the Atlantic, negative perturbations produce insignificant net changes in heat content. The processes active in our model appear robust, although their net result is likely model- and experiment-dependent.

Extension of the TAMSAT satellite-based rainfall monitoring over Africa and from 1983 to present

Tropical Applications of Meteorology Using Satellite Data and Ground-Based Observations (TAMSAT) rainfall monitoring products have been extended to provide spatially contiguous rainfall estimates across Africa. This has been achieved through a new, climatology-based calibration, which varies in both space and time. As a result, cumulative estimates of rainfall are now issued at the end of each 10-day period (dekad) at 4-km spatial resolution with pan-African coverage. The utility of the products for decision making is improved by the routine provision of validation reports, for which the 10-day (dekadal) TAMSAT rainfall estimates are compared with independent gauge observations. This paper describes the methodology by which the TAMSAT method has been applied to generate the pan-African rainfall monitoring products. It is demonstrated through comparison with gauge measurements that the method provides skillful estimates, although with a systematic dry bias. This study illustrates TAMSAT’s value as a complementary method of estimating rainfall through examples of successful operational application.

Nonlinear regional warming with increasing CO₂ concentration

When considering adaptation measures and global climate mitigation goals, stakeholders need regional-scale climate projections, including the range of plausible warming rates. To assist these stakeholders, it is important to understand whether some locations may see disproportionately high or low warming from additional forcing above targets such as 2 K (ref. 1). There is a need to narrow uncertainty2 in this nonlinear warming, which requires understanding how climate changes as forcings increase from medium to high levels. However, quantifying and understanding regional nonlinear processes is challenging. Here we show that regional-scale warming can be strongly superlinear to successive CO2 doublings, using five different climate models. Ensemble-mean warming is superlinear over most land locations. Further, the inter-model spread tends to be amplified at higher forcing levels, as nonlinearities grow—especially when considering changes per kelvin of global warming. Regional nonlinearities in surface warming arise from nonlinearities in global-mean radiative balance, the Atlantic meridional overturning circulation, surface snow/ice cover and evapotranspiration. For robust adaptation and mitigation advice, therefore, potentially avoidable climate change (the difference between business-as-usual and mitigation scenarios) and unavoidable climate change (change under strong mitigation scenarios) may need different analysis methods.

The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

Care in the community? Interpretations of a fractured goat bone from Neolithic Jarmo, Iraq

A case study of a goat metatarsal exhibiting a complex diaphyseal fracture from Pottery Neolithic Jarmo in the Central Zagros region of the eastern Fertile Crescent is here described and analysed. The Central Zagros is one of the areas with the earliest evidence for goat domestication. The significance of the pathology may be viewed within the context of domestic goat ecology in the landscape of Jarmo, potentially impacting browsing behaviour (goats raise themselves on their hind limbs to browse) and movement with the herd in the landscape (the terrain around Jarmo is very steep in places, which would be difficult for an animal to navigate on three legs). In the light of this, possible levels of care that the Neolithic human community may have afforded this animal are discussed – from a situation where therapeutic intervention may have occurred, to one of stall confinement of the animal to allow the pathology to heal, to a position of simple awareness of the condition – and how this impacts on our understanding of changes in attitudes towards animals through the process of domestication.

Preliminary ethnoarchaeological research on modern animal husbandry in Bestansur, Iraqi Kurdistan: integrating animal, plant and environmental data

This paper presents preliminary results from an ethnoarchaeological study of animal husbandry in the modern village of Bestansur, situated in the lower Zagros Mountains of Iraqi Kurdistan. This research explores how modern families use and manage their livestock within the local landscape and identifies traces of this use. The aim is to provide the groundwork for future archaeological investigations focusing on the nearby Neolithic site of Bestansur. This is based on the premise that modern behaviours can suggest testable patterns for past practices within the same functional and ecological domains. Semi-structured interviews conducted with villagers from several households provided large amounts of information on modern behaviours that helped direct data collection, and which also illustrate notable shifts in practices and use of the local landscape over time. Strontium isotope analysis of modern plant material demonstrates that a measurable variation exists between the alluvial floodplain and the lower foothills, while analysis of modern dung samples shows clear variation between sheep/goat and cow dung, in terms of numbers of faecal spherulites. These results are specific to the local environment of Bestansur and can be used for evaluating and contextualising archaeological evidence as well as providing modern reference material for comparative purposes.

Terrestrial biosphere changes over the last 120 kyr and their impact on ocean δ 13C

A new global synthesis and biomization of long (>40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 ka BP, and between 60 and 65 ka BP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ13C based on modelled land carbon storage, and palaeo-archives of ocean δ13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ13C changes.