Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

A model of melt pond evolution on sea ice

A one-dimensional, thermodynamic, and radiative model of a melt pond on sea ice is presented that explicitly treats the melt pond as an extra phase. A two-stream radiation model, which allows albedo to be determined from bulk optical properties, and a parameterization of the summertime evolution of optical properties, is used. Heat transport within the sea ice is described using an equation describing heat transport in a mushy layer of a binary alloy (salt water). The model is tested by comparison of numerical simulations with SHEBA data and previous modeling. The presence of melt ponds on the sea ice surface is demonstrated to have a significant effect on the heat and mass balance. Sensitivity tests indicate that the maximum melt pond depth is highly sensitive to optical parameters and drainage. INDEX TERMS: 4207 Oceanography: General: Arctic and Antarctic oceanography; 4255 Oceanography: General: Numerical modeling; 4299 Oceanography: General: General or miscellaneous; KEYWORDS: sea ice, melt pond, albedo, Arctic Ocean, radiation model, thermodynamic

Does convective aggregation need to be represented in cumulus parameterizations?

Tropical deep convection exhibits a variety of levels of aggregation over a wide range of scales. Based on a multisatellite analysis, the present study shows at mesoscale that different levels of aggregation are statistically associated with differing large-scale atmospheric states, despite similar convective intensity and large-scale forcings. The more aggregated the convection, the dryer and less cloudy the atmosphere, the stronger the outgoing longwave radiation, and the lower the planetary albedo. This suggests that mesoscale convective aggregation has the potential to affect couplings between moisture and convection and between convection, radiation, and large-scale ascent. In so doing, aggregation may play a role in phenomena such as “hot spots” or the Madden-Julian Oscillation. These findings support the need for the representation of mesoscale organization in cumulus parameterizations; most parameterizations used in current climate models lack any such representation. The ability of a cloud system-resolving model to reproduce observed relationships suggests that such models may be useful to guide attempts at parameterizations of convective aggregation.

Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum, Part 2: feedbacks with emphasis on the location of the ITCZ and mid- and high latitudes heat budget

A set of coupled ocean-atmosphere(-vegetation) simulations using state of the art climate models is now available for the Last Glacial Maximum (LGM) and the Mid-Holocene (MH) through the second phase of the Paleoclimate Modeling Intercomparison Project (PMIP2). Here we quantify the latitudinal shift of the location of the Intertropical Convergence Zone (ITCZ) in the tropical regions during boreal summer and the change in precipitation in the northern part of the ITCZ. For both periods the shift is more pronounced over the continents and East Asia. The maritime continent is the region where the largest spread is found between models. We also clearly establish that the larger the increase in the meridional temperature gradient in the tropical Atlantic during summer at the MH, the larger the change in precipitation over West Africa. The vegetation feedback is however not as large as found in previous studies, probably due to model differences in the control simulation. Finally, we show that the feedback from snow and sea-ice at mid and high latitudes contributes for half of the cooling in the Northern Hemisphere for the LGM, with the remaining being achieved by the reduced CO2 and water vapour in the atmosphere. For the MH the snow and albedo feedbacks strengthen the spring cooling and enhance the boreal summer warming, whereas water vapour reinforces the late summer warming. These feedbacks are modest in the Southern Hemisphere. For the LGM most of the surface cooling is due to CO2 and water vapour.

Radiative forcing effect of ice-age dust

During glacial periods, dust deposition rates and inferred atmospheric concentrations were globally much higher than present. According to recent model results, the large enhancement of atmospheric dust content at the last glacial maximum (LGM) can be explained only if increases in the potential dust source areas are taken into account. Such increases are to be expected, due to effects of low precipitation and low atmospheric (CO2) on plant growth. Here the modelled three-dimensional dust fields from Mahowald et al. and modelled seasonally varying surface-albedo fields derived in a parallel manner, are used to quantify the mean radiative forcing due to modern (non-anthropogenic) and LGM dust. The effect of mineralogical provenance on the radiative properties of the dust is taken into account, as is the range of optical properties associated with uncertainties about the mixing state of the dust particles. The high-latitude (poleward of 45°) mean change in forcing (LGM minus modern) is estimated to be small (–0.9 to +0.2 W m–2), especially when compared to nearly –20 W m–2 due to reflection from the extended ice sheets. Although the net effect of dust over ice sheets is a positive forcing (warming), much of the simulated high-latitude dust was not over the ice sheets, but over unglaciated regions close to the expanded dust source region in central Asia. In the tropics the change in forcing is estimated to be overall negative, and of similarly large magnitude (–2.2 to –3.2 W m–2) to the radiative cooling effect of low atmospheric (CO2). Thus, the largest long-term climatic effect of the LGM dust is likely to have been a cooling of the tropics. Low tropical sea-surface temperatures, low atmospheric (CO2) and high atmospheric dust loading may be mutually reinforcing due to multiple positive feedbacks, including the negative radiative forcing effect of dust.

Energy exchange in a dense urban environment Part II: impact of spatial heterogeneity of the surface

The centre of cities, characterised by spatial and temporal complexity, are challenging environments for micrometeorological research. This paper considers the impact of sensor location and heterogeneity of the urban surface on flux observations in the dense city centre of London, UK. Data gathered at two sites in close vicinity, but with different measurement heights, were analysed to investigate the influence of source area characteristics on long-term radiation and turbulent heat fluxes. Combining consideration of diffuse radiation and effects of specular reflections, the non-Lambertian urban surface is found to impact the measurements of surface albedo. Comparisons of observations from the two sites reveal that turbulent heat fluxes are similar under some flow conditions. However, they mostly observe processes at different scales due to their differing measurement heights, highlighting the critical impact of siting sensors in urban areas. A detailed source area analysis is presented to investigate the surface controls influencing the energy exchanges at the different scales

Land surface feedbacks and palaeomonsoons in northern Africa

We ran a sequence of climate model experiments for 6000 years ago, with land-surface conditions based on a realistic map of palaeovegetation, lakes and wetlands, to quantify the effects of land-surface feedbacks in the Saharan region. Vegetation-induced albedo and moisture flux changes produced year-round warming, forced the monsoon to 17°–25°N two months earlier, and shifted the precipitation belt ≈300 km northwards compared to the effects of orbital forcing alone. The addition of lakes and wetlands produced localised changes in evaporation and precipitation, but caused no further extension of the monsoon belt. Diagnostic analyses with biome and continental hydrology models showed that the combined land-surface feedbacks, although substantial, could neither maintain grassland as far north as observed (≈26°N) nor maintain Lake “MegaChad” (330,000 km²).

September Arctic sea-ice minimum predicted by spring melt-pond fraction

The area of Arctic September sea ice has diminished from about 7 million km2 in the 1990s to less than 5 million km2 in five of the past seven years, with a record minimum of 3.6 million km2 in 2012 (ref. 1). The strength of this decrease is greater than expected by the scientific community, the reasons for this are not fully understood, and its simulation is an on-going challenge for existing climate models2, 3. With growing Arctic marine activity there is an urgent demand for forecasting Arctic summer sea ice4. Previous attempts at seasonal forecasts of ice extent were of limited skill5, 6, 7, 8, 9. However, here we show that the Arctic sea-ice minimum can be accurately forecasted from melt-pond area in spring. We find a strong correlation between the spring pond fraction and September sea-ice extent. This is explained by a positive feedback mechanism: more ponds reduce the albedo; a lower albedo causes more melting; more melting increases pond fraction. Our results help explain the acceleration of Arctic sea-ice decrease during the past decade. The inclusion of our new melt-pond model10 promises to improve the skill of future forecast and climate models in Arctic regions and beyond.

Synergistic angular and spectral estimation of aerosol properties using CHRIS/PROBA-1 and simulated Sentinel-3~data

A method has been developed to estimate Aerosol Optical Depth (AOD), Fine Mode Fraction (FMF) and Single Scattering Albedo (SSA) over land surfaces using simulated Sentinel-3 data. The method uses inversion of a coupled surface/atmosphere radiative transfer model, and includes a general physical model of angular surface reflectance. An iterative process is used to determine the optimum value of the aerosol properties providing the best fit of the corrected reflectance values for a number of view angles and wavelengths with those provided by the physical model. A method of estimating AOD using only angular retrieval has previously been demonstrated on data from the ENVISAT and PROBA-1 satellite instruments, and is extended here to the synergistic spectral and angular sampling of Sentinel-3 and the additional aerosol properties. The method is tested using hyperspectral, multi-angle Compact High Resolution Imaging Spectrometer (CHRIS) images. The values obtained from these CHRIS observations are validated using ground based sun-photometer measurements. Results from 22 image sets using the synergistic retrieval and improved aerosol models show an RMSE of 0.06 in AOD, reduced to 0.03 over vegetated targets.

Development of the Surface Urban Energy and Water balance Scheme (SUEWS) for cold climate cities

The Surface Urban Energy and Water Balance Scheme (SUEWS) is developed to include snow. The processes addressed include accumulation of snow on the different urban surface types: snow albedo and density aging, snow melting and re-freezing of meltwater. Individual model parameters are assessed and independently evaluated using long-term observations in the two cold climate cities of Helsinki and Montreal. Eddy covariance sensible and latent heat fluxes and snow depth observations are available for two sites in Montreal and one in Helsinki. Surface runoff from two catchments (24 and 45 ha) in Helsinki and snow properties (albedo and density) from two sites in Montreal are also analysed. As multiple observation sites with different land-cover characteristics are available in both cities, model development is conducted independent of evaluation. The developed model simulates snowmelt related runoff well (within 19% and 3% for the two catchments in Helsinki when there is snow on the ground), with the springtime peak estimated correctly. However, the observed runoff peaks tend to be smoother than the simulated ones, likely due to the water holding capacity of the catchments and the missing time lag between the catchment and the observation point in the model. For all three sites the model simulates the timing of the snow accumulation and melt events well, but underestimates the total snow depth by 18–20% in Helsinki and 29–33% in Montreal. The model is able to reproduce the diurnal pattern of net radiation and turbulent fluxes of sensible and latent heat during cold snow, melting snow and snow-free periods. The largest model uncertainties are related to the timing of the melting period and the parameterization of the snowmelt. The results show that the enhanced model can simulate correctly the exchange of energy and water in cold climate cities at sites with varying surface cover.

Key conclusions of the first international urban land surface model comparison project

he first international urban land surface model comparison was designed to identify three aspects of the urban surface-atmosphere interactions: (1) the dominant physical processes, (2) the level of complexity required to model these, and 3) the parameter requirements for such a model. Offline simulations from 32 land surface schemes, with varying complexity, contributed to the comparison. Model results were analysed within a framework of physical classifications and over four stages. The results show that the following are important urban processes; (i) multiple reflections of shortwave radiation within street canyons, (ii) reduction in the amount of visible sky from within the canyon, which impacts on the net long-wave radiation, iii) the contrast in surface temperatures between building roofs and street canyons, and (iv) evaporation from vegetation. Models that use an appropriate bulk albedo based on multiple solar reflections, represent building roof surfaces separately from street canyons and include a representation of vegetation demonstrate more skill, but require parameter information on the albedo, height of the buildings relative to the width of the streets (height to width ratio), the fraction of building roofs compared to street canyons from a plan view (plan area fraction) and the fraction of the surface that is vegetated. These results, whilst based on a single site and less than 18 months of data, have implications for the future design of urban land surface models, the data that need to be measured in urban observational campaigns, and what needs to be included in initiatives for regional and global parameter databases.

Contributions of different cloud types to feedbacks and rapid adjustments in CMIP5

Using five climate model simulations of the response to an abrupt quadrupling of CO2, the authors perform the first simultaneous model intercomparison of cloud feedbacks and rapid radiative adjustments with cloud masking effects removed, partitioned among changes in cloud types and gross cloud properties. Upon CO2 quadrupling, clouds exhibit a rapid reduction in fractional coverage, cloud-top pressure, and optical depth, with each contributing equally to a 1.1 W m−2 net cloud radiative adjustment, primarily from shortwave radiation. Rapid reductions in midlevel clouds and optically thick clouds are important in reducing planetary albedo in every model. As the planet warms, clouds become fewer, higher, and thicker, and global mean net cloud feedback is positive in all but one model and results primarily from increased trapping of longwave radiation. As was true for earlier models, high cloud changes are the largest contributor to intermodel spread in longwave and shortwave cloud feedbacks, but low cloud changes are the largest contributor to the mean and spread in net cloud feedback. The importance of the negative optical depth feedback relative to the amount feedback at high latitudes is even more marked than in earlier models. The authors show that the negative longwave cloud adjustment inferred in previous studies is primarily caused by a 1.3 W m−2 cloud masking of CO2 forcing. Properly accounting for cloud masking increases net cloud feedback by 0.3 W m−2 K−1, whereas accounting for rapid adjustments reduces by 0.14 W m−2 K−1 the ensemble mean net cloud feedback through a combination of smaller positive cloud amount and altitude feedbacks and larger negative optical depth feedbacks.

Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust

Saharan dust affects the climate by altering the radiation balance and by depositing minerals to the Atlantic Ocean. Both are dependent on particle size. We present aircraft measurements comprising 42 profiles of size distribution (0.1–300 µm), representing freshly uplifted dust, regional aged dust, and dust in the Saharan Air Layer (SAL) over the Canary Islands. The mean effective diameter of dust in SAL profiles is 4.5 µm smaller than that in freshly uplifted dust, while the vertical structure changes from a low shallow layer (0–1.5 km) to a well-mixed deep Saharan dust layer (0–5 km). Size distributions show a loss of 60 to 90% of particles larger than 30 µm 12 h after uplift. The single scattering albedo (SSA) increases from 0.92 to 0.94 to 0.95 between fresh, aged, and SAL profiles: this is enough to alter heating rates by 26%. Some fresh dust close to the surface shows SSA as low as 0.85

Regional and monthly and clear-sky aerosol direct radiative effect (and forcing) derived from the GlobAEROSOL-AATSR satellite aerosol product

Using the GlobAEROSOL-AATSR dataset, estimates of the instantaneous, clear-sky, direct aerosol radiative effect and radiative forcing have been produced for the year 2006. Aerosol Robotic Network sun-photometer measurements have been used to characterise the random and systematic error in the GlobAEROSOL product for 22 regions covering the globe. Representative aerosol properties for each region were derived from the results of a wide range of literature sources and, along with the de-biased GlobAEROSOL AODs, were used to drive an offline version of the Met Office unified model radiation scheme. In addition to the mean AOD, best-estimate run of the radiation scheme, a range of additional calculations were done to propagate uncertainty estimates in the AOD, optical properties, surface albedo and errors due to the temporal and spatial averaging of the AOD fields. This analysis produced monthly, regional estimates of the clear-sky aerosol radiative effect and its uncertainty, which were combined to produce annual, global mean values of (−6.7±3.9)Wm−2 at the top of atmosphere (TOA) and (−12±6)Wm−2 at the surface. These results were then used to give estimates of regional, clear-sky aerosol direct radiative forcing, using modelled pre-industrial AOD fields for the year 1750 calculated for the AEROCOM PRE experiment. However, as it was not possible to quantify the uncertainty in the pre-industrial aerosol loading, these figures can only be taken as indicative and their uncertainties as lower bounds on the likely errors. Although the uncertainty on aerosol radiative effect presented here is considerably larger than most previous estimates, the explicit inclusion of the major sources of error in the calculations suggest that they are closer to the true constraint on this figure from similar methodologies, and point to the need for more, improved estimates of both global aerosol loading and aerosol optical properties.

Aerosol scattering and absorption during the EUCAARI-LONGREX flights of the Facility for Airborne Atmospheric Measurements (FAAM) BAe-146: can measurements and models agree?

Scattering and absorption by aerosol in anthropogenically perturbed air masses over Europe has been measured using instrumentation flown on the UK’s BAe-146-301 large Atmospheric Research Aircraft (ARA) operated by the Facility for Airborne Atmospheric Measurements (FAAM) on 14 flights during the EUCAARI-LONGREX campaign in May 2008. The geographical and temporal variations of the derived shortwave optical properties of aerosol are presented. Values of single scattering albedo of dry aerosol at 550 nm varied considerably from 0.86 to near unity, with a campaign average of 0.93 ± 0.03. Dry aerosol optical depths ranged from 0.030 ± 0.009 to 0.24 ± 0.07. An optical properties closure study comparing calculations from composition data and Mie scattering code with the measured properties is presented. Agreement to within measurement uncertainties of 30% can be achieved for both scattering and absorption,but the latter is shown to be sensitive to the refractive indices chosen for organic aerosols, and to a lesser extent black carbon, as well as being highly dependent on the accuracy of the absorption measurements. Agreement with the measured absorption can be achieved either if organic carbon is assumed to be weakly absorbing, or if the organic aerosol is purely scattering and the absorption measurement is an overestimate due to the presence of large amounts of organic carbon. Refractive indices could not be inferred conclusively due to this uncertainty, despite the enhancement in methodology compared to previous studies that derived from the use of the black carbon measurements. Hygroscopic growth curves derived from the wet nephelometer indicate moderate water uptake by the aerosol with a campaign mean f (RH) value (ratio in scattering) of 1.5 (range from 1.23 to 1.63) at 80% relative humidity. This value is qualitatively consistent with the major chemical components of the aerosol measured by the aerosol mass spectrometer, which are primarily mixed organics and nitrate and some sulphate.