Modeled and observed atmospheric radiation balance during the West African dry season: role of mineral dust, biomass burning aerosol, and surface albedo

The global radiation balance of the atmosphere is still poorly observed, particularly at the surface. We investigate the observed radiation balance at (1) the surface using the ARM Mobile Facility in Niamey, Niger, and (2) the top of the atmosphere (TOA) over West Africa using data from the Geostationary Earth Radiation Budget (GERB) instrument on board Meteosat-8. Observed radiative fluxes are compared with predictions from the global numerical weather prediction (NWP) version of the Met Office Unified Model (MetUM). The evaluation points to major shortcomings in the NWP model's radiative fluxes during the dry season (December 2005 to April 2006) arising from (1) a lack of absorbing aerosol in the model (mineral dust and biomass burning aerosol) and (2) a poor specification of the surface albedo. A case study of the major Saharan dust outbreak of 6–12 March 2006 is used to evaluate a parameterization of mineral dust for use in the NWP models. The model shows good predictability of the large-scale flow out to 4–5 days with the dust parameterization providing reasonable dust uplift, spatial distribution, and temporal evolution for this strongly forced dust event. The direct radiative impact of the dust reduces net downward shortwave (SW) flux at the surface (TOA) by a maximum of 200 W m−2 (150 W m−2), with a SW heating of the atmospheric column. The impacts of dust on terrestrial radiation are smaller. Comparisons of TOA (surface) radiation balance with GERB (ARM) show the “dusty” forecasts reduce biases in the radiative fluxes and improve surface temperatures and vertical thermodynamic structure.

On spectral invariance of single scattering albedo for water droplets and ice crystals at weakly absorbing wavelengths

The single scattering albedo w_0l in atmospheric radiative transfer is the ratio of the scattering coefficient to the extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength l and droplet size r. This note shows that for water droplets at weakly absorbing wavelengths, the ratio w_0l(r)/w_0l(r0) of two single scattering albedo spectra is a linear function of w_0l(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo spectrum w_0l(r) via one known spectrum w_0l(r0). We provide a simple physical explanation of the discovered relationship. Similar linear relationships were found for the single scattering albedo spectra of non-spherical ice crystals.

Comparison of MODIS broadband albedo over an agricultural site with ground measurements and values derived from Earth observation data at a range of spatial scales

Land surface albedo is dependent on atmospheric state and hence is difficult to validate. Over the UK persistent cloud cover and land cover heterogeneity at moderate (km-scale) spatial resolution can also complicate comparison of field-measured albedo with that derived from instruments such as the Moderate Resolution Imaging Spectrometer (MODIS). A practical method of comparing moderate resolution satellite-derived albedo with ground-based measurements over an agricultural site in the UK is presented. Point measurements of albedo made on the ground are scaled up to the MODIS resolution (1 km) through reflectance data obtained at a range of spatial scales. The point measurements of albedo agreed in magnitude with MODIS values over the test site to within a few per cent, despite problems such as persistent cloud cover and the difficulties of comparing measurements made during different years. Albedo values derived from airborne and field-measured data were generally lower than the corresponding satellite-derived values. This is thought to be due to assumptions made regarding the ratio of direct to diffuse illumination used when calculating albedo from reflectance. Measurements of albedo calculated for specific times fitted closely to the trajectories of temporal albedo derived from both Systeme pour l'Observation de la Terre (SPOT) Vegetation (VGT) and MODIS instruments.

Tackling regional climate change by leaf albedo bio-geoengineering

The likelihood that continuing greenhouse-gas emissions will lead to an unmanageable degree of climate change [1] has stimulated the search for planetary-scale technological solutions for reducing global warming [2] (“geoengineering”), typically characterized by the necessity for costly new infrastructures and industries [3]. We suggest that the existing global infrastructure associated with arable agriculture can help, given that crop plants exert an important influence over the climatic energy budget 4 and 5 because of differences in their albedo (solar reflectivity) compared to soils and to natural vegetation [6]. Specifically, we propose a “bio-geoengineering” approach to mitigate surface warming, in which crop varieties having specific leaf glossiness and/or canopy morphological traits are specifically chosen to maximize solar reflectivity. We quantify this by modifying the canopy albedo of vegetation in prescribed cropland areas in a global-climate model, and thereby estimate the near-term potential for bio-geoengineering to be a summertime cooling of more than 1°C throughout much of central North America and midlatitude Eurasia, equivalent to seasonally offsetting approximately one-fifth of regional warming due to doubling of atmospheric CO2[7]. Ultimately, genetic modification of plant leaf waxes or canopy structure could achieve greater temperature reductions, although better characterization of existing intraspecies variability is needed first.

Assessing the benefits of crop albedo bio-geoengineering

It has been proposed that growing crop varieties with higher canopy albedo would lower summer-time temperatures over North America and Eurasia and provide a partial mitigation of global warming ('bio-geoengineering') (Ridgwell et al 2009 Curr. Biol. 19 1–5). Here, we use a coupled ocean–atmosphere–vegetation model (HadCM3) with prescribed agricultural regions, to investigate to what extent the regional effectiveness of crop albedo bio-geoengineering might be influenced by a progressively warming climate as well as assessing the impacts on regional hydrological cycling and primary productivity. Consistent with previous analysis, we find that the averted warming due to increasing crop canopy albedo by 0.04 is regionally and seasonally specific, with the largest cooling of ~1 °C for Europe in summer whereas in the low latitude monsoonal SE Asian regions of high density cropland, the greatest cooling is experienced in winter. In this study we identify potentially important positive impacts of increasing crop canopy albedo on soil moisture and primary productivity in European cropland regions, due to seasonal increases in precipitation. We also find that the background climate state has an important influence on the predicted regional effectiveness of bio-geoengineering on societally-relevant timescales (ca 100 years). The degree of natural climate variability and its dependence on greenhouse forcing that are evident in our simulations highlights the difficulties faced in the detection and verification of climate mitigation in geoengineering schemes. However, despite the small global impact, regionally focused schemes such as crop albedo bio-geoengineering have detection advantages.

Global analysis of radiative forcing from fire-induced shortwave albedo change

Land surface albedo, a key parameter to derive Earth's surface energy balance, is used in the parameterization of numerical weather prediction, climate monitoring and climate change impact assessments. Changes in albedo due to fire have not been fully investigated on a continental and global scale. The main goal of this study, therefore, is to quantify the changes in instantaneous shortwave albedo produced by biomass burning activities and their associated radiative forcing. The study relies on the MODerate-resolution Imaging Spectroradiometer (MODIS) MCD64A1 burned-area product to create an annual composite of areas affected by fire and the MCD43C2 bidirectional reflectance distribution function (BRDF) albedo snow-free product to compute a bihemispherical reflectance time series. The approximate day of burning is used to calculate the instantaneous change in shortwave albedo. Using the corresponding National Centers for Environmental Prediction (NCEP) monthly mean downward solar radiation flux at the surface, the global radiative forcing associated with fire was computed. The analysis reveals a mean decrease in shortwave albedo of −0.014 (1σ = 0.017), causing a mean positive radiative forcing of 3.99 Wm−2 (1σ = 4.89) over the 2002–20012 time period in areas affected by fire. The greatest drop in mean shortwave albedo change occurs in 2002, which corresponds to the highest total area burned (378 Mha) observed in the same year and produces the highest mean radiative forcing (4.5 Wm−2). Africa is the main contributor in terms of burned area, but forests globally give the highest radiative forcing per unit area and thus give detectable changes in shortwave albedo. The global mean radiative forcing for the whole period studied (~0.0275 Wm−2) shows that the contribution of fires to the Earth system is not insignificant.

Quantifying sources of inter-model diversity in the cloud albedo effect

There is large diversity in simulated aerosol forcing among models that participated in the fifth Coupled Model Intercomparison Project (CMIP5), particularly related to aerosol interactions with clouds. Here we use the reported model data and fitted aerosol-cloud relations to separate the main sources of inter-model diversity in the magnitude of the cloud albedo effect. There is large diversity in the global load and spatial distribution of sulfate aerosol, as well as in global-mean cloud-top effective radius. The use of different parameterizations of aerosol-cloud interactions makes the largest contribution to diversity in modeled radiative forcing (up to -39%, +48% about the mean estimate). Uncertainty in pre-industrial sulfate load also makes a substantial contribution (-15%, +61% about the mean estimate), with smaller contributions from inter-model differences in the historical change in sulfate load and in mean cloud fraction.

Snow-albedo feedback and Swiss spring temperature trends

We quantify the effect of the snow-albedo feedback on Swiss spring temperature trends using daily temperature and snow depth measurements from six station pairs for the period 1961–2011. We show that the daily mean 2-m temperature of a spring day without snow cover is on average 0.4 °C warmer than one with snow cover at the same location. This estimate is comparable with estimates from climate modelling studies. Caused by the decreases in snow pack, the snow-albedo feedback amplifies observed temperature trends in spring. The influence is small and confined to areas around the upward-moving snow line in spring and early summer. For the 1961–2011 period, the related temperature trend increases are in the order of 3–7 % of the total observed trend.

Seasonal variations of the physical and optical characteristics of Saharan Dust: results from the Dust outflow and Deposition to the Ocean (DODO) experiment

North African dust is important for climate through its direct radiative effect on solar and terrestrial radiation and its role in the biogeochemical system. The Dust Outflow and Deposition to the Ocean project (DODO) aimed to characterize the physical and optical properties of airborne North African dust in two seasons and to use these observations to constrain model simulations, with the ultimate aim of being able to quantify the deposition of iron to the North Atlantic Ocean. The in situ properties of dust from airborne campaigns measured during February and August 2006, based at Dakar, Senegal, are presented here. Average values of the single scattering albedo (0.99, 0.98), mass specific extinction (0.85 m^2 g^-1 , 1.14 m^2 g^-1 ), asymmetry parameter (0.68, 0.68), and refractive index (1.53--0.0005i,1.53--0.0014i) for the accumulation mode were found to differ by varying degrees between the dry and wet season, respectively. It is hypothesized that these differences are due to different source regions and transport processes which also differ between the DODO campaigns. Elemental ratios of Ca/Al were found to differ between the dry and wet season (1.1 and 0.5, respectively). Differences in vertical profiles are found between seasons and between land and ocean locations and reflect the different dynamics of the seasons. Using measurements of the coarse mode size distribution and illustrative Mie calculations, the optical properties are found to be very sensitive to the presence and amount of coarse mode of mineral dust, and the importance of accurate measurements of the coarse mode of dust is highlighted.

Tripleclouds: an efficient method for representing horizontal cloud inhomogeneity in 1D radiation schemes by using three regions at each height

Radiation schemes in general circulation models currently make a number of simplifications when accounting for clouds, one of the most important being the removal of horizontal inhomogeneity. A new scheme is presented that attempts to account for the neglected inhomogeneity by using two regions of cloud in each vertical level of the model as opposed to one. One of these regions is used to represent the optically thinner cloud in the level, and the other represents the optically thicker cloud. So, along with the clear-sky region, the scheme has three regions in each model level and is referred to as “Tripleclouds.” In addition, the scheme has the capability to represent arbitrary vertical overlap between the three regions in pairs of adjacent levels. This scheme is implemented in the Edwards–Slingo radiation code and tested on 250 h of data from 12 different days. The data are derived from cloud retrievals using radar, lidar, and a microwave radiometer at Chilbolton, southern United Kingdom. When the data are grouped into periods equivalent in size to general circulation model grid boxes, the shortwave plane-parallel albedo bias is found to be 8%, while the corresponding bias is found to be less than 1% using Tripleclouds. Similar results are found for the longwave biases. Tripleclouds is then compared to a more conventional method of accounting for inhomogeneity that multiplies optical depths by a constant scaling factor, and Tripleclouds is seen to improve on this method both in terms of top-of-atmosphere radiative flux biases and internal heating rates.

Sampling uncertainties in surface radiation budget calculations in RADAGAST

In the Radiative Atmospheric Divergence Using ARM Mobile Facility GERB and AMMA Stations (RADAGAST) project we calculate the divergence of radiative flux across the atmosphere by comparing fluxes measured at each end of an atmospheric column above Niamey, in the African Sahel region. The combination of broadband flux measurements from geostationary orbit and the deployment for over 12 months of a comprehensive suite of active and passive instrumentation at the surface eliminates a number of sampling issues that could otherwise affect divergence calculations of this sort. However, one sampling issue that challenges the project is the fact that the surface flux data are essentially measurements made at a point, while the top-of-atmosphere values are taken over a solid angle that corresponds to an area at the surface of some 2500 km2. Variability of cloud cover and aerosol loading in the atmosphere mean that the downwelling fluxes, even when averaged over a day, will not be an exact match to the area-averaged value over that larger area, although we might expect that it is an unbiased estimate thereof. The heterogeneity of the surface, for example, fixed variations in albedo, further means that there is a likely systematic difference in the corresponding upwelling fluxes. In this paper we characterize and quantify this spatial sampling problem. We bound the root-mean-square error in the downwelling fluxes by exploiting a second set of surface flux measurements from a site that was run in parallel with the main deployment. The differences in the two sets of fluxes lead us to an upper bound to the sampling uncertainty, and their correlation leads to another which is probably optimistic as it requires certain other conditions to be met. For the upwelling fluxes we use data products from a number of satellite instruments to characterize the relevant heterogeneities and so estimate the systematic effects that arise from the flux measurements having to be taken at a single point. The sampling uncertainties vary with the season, being higher during the monsoon period. We find that the sampling errors for the daily average flux are small for the shortwave irradiance, generally less than 5 W m−2, under relatively clear skies, but these increase to about 10 W m−2 during the monsoon. For the upwelling fluxes, again taking daily averages, systematic errors are of order 10 W m−2 as a result of albedo variability. The uncertainty on the longwave component of the surface radiation budget is smaller than that on the shortwave component, in all conditions, but a bias of 4 W m−2 is calculated to exist in the surface leaving longwave flux.

Some reports of snowfall from fog during the UK winter of 2008/09

Snowfall during anticyclonic, non-frontal, and foggy conditions is surprising. Because it is often not forecast, it can present a hazard to transport and modify the surface albedo. In this report, we present some observations of snowfall during conditions of freezing fog in the UK during the winter of 2008/09.

How many facets are needed to represent the surface energy balance of an urban area?

We investigate the question of how many facets are needed to represent the energy balance of an urban area by developing simplified 3-, 2- and 1-facet versions of a 4-facet energy balance model of two-dimensional streets and buildings. The 3-facet model simplifies the 4-facet model by averaging over the canyon orientation, which results in similar net shortwave and longwave balances for both wall facets, but maintains the asymmetry in the heat fluxes within the street canyon. For the 2-facet model, on the assumption that the wall and road temperatures are equal, the road and wall facets can be combined mathematically into a single street-canyon facet with effective values of the heat transfer coefficient, albedo, emissivity and thermodynamic properties, without further approximation. The 1-facet model requires the additional assumption that the roof temperature is also equal to the road and wall temperatures. Idealised simulations show that the geometry and material properties of the walls and road lead to a large heat capacity of the combined street canyon, whereas the roof behaves like a flat surface with low heat capacity. This means that the magnitude of the diurnal temperature variation of the street-canyon facets are broadly similar and much smaller than the diurnal temperature variation of the roof facets. Consequently, the approximation that the street-canyon facets have similar temperatures is sound, and the road and walls can be combined into a single facet. The roof behaves very differently and a separate roof facet is required. Consequently, the 2-facet model performs similarly to the 4-facet model, while the 1-facet model does not. The models are compared with previously published observations collected in Mexico City. Although the 3- and 2-facet models perform better than the 1-facet model, the present models are unable to represent the phase of the sensible heat flux. This result is consistent with previous model comparisons, and we argue that this feature of the data cannot be produced by a single column model. We conclude that a 2-facet model is necessary, and for numerical weather prediction sufficient, to model an urban surface, and that this conclusion is robust and therefore applicable to more general geometries.

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.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.