The use of land-sea warming contrast under climate change to improve impact metrics

A favoured method of assimilating information from state-of-the-art climate models into integrated assessment models of climate impacts is to use the transient climate response (TCR) of the climate models as an input, sometimes accompanied by a pattern matching approach to provide spatial information. More recent approaches to the problem use TCR with another independent piece of climate model output: the land-sea surface warming ratio (φ). In this paper we show why the use of φ in addition to TCR has such utility. Multiple linear regressions of surface temperature change onto TCR and φ in 22 climate models from the CMIP3 multi-model database show that the inclusion of φ explains a much greater fraction of the inter-model variance than using TCR alone. The improvement is particularly pronounced in North America and Eurasia in the boreal summer season, and in the Amazon all year round. The use of φ as the second metric is beneficial for three reasons: firstly it is uncorrelated with TCR in state-of-the-art climate models and can therefore be considered as an independent metric; secondly, because of its projected time-invariance, the magnitude of φ is better constrained than TCR in the immediate future; thirdly, the use of two variables is much simpler than approaches such as pattern scaling from climate models. Finally we show how using the latest estimates of φ from climate models with a mean value of 1.6—as opposed to previously reported values of 1.4—can significantly increase the mean time-integrated discounted damage projections in a state-of-the-art integrated assessment model by about 15 %. When compared to damages calculated without the inclusion of the land-sea warming ratio, this figure rises to 65 %, equivalent to almost 200 trillion dollars over 200 years.

Land-atmosphere energy exchange in Arctic tundra and boreal forest: available data and feedbacks to climate

This paper summarizes and analyses available data on the surface energy balance of Arctic tundra and boreal forest. The complex interactions between ecosystems and their surface energy balance are also examined, including climatically induced shifts in ecosystem type that might amplify or reduce the effects of potential climatic change. High latitudes are characterized by large annual changes in solar input. Albedo decreases strongly from winter, when the surface is snow-covered, to summer, especially in nonforested regions such as Arctic tundra and boreal wetlands. Evapotranspiration (QE) of high-latitude ecosystems is less than from a freely evaporating surface and decreases late in the season, when soil moisture declines, indicating stomatal control over QE, particularly in evergreen forests. Evergreen conifer forests have a canopy conductance half that of deciduous forests and consequently lower QE and higher sensible heat flux (QH). There is a broad overlap in energy partitioning between Arctic and boreal ecosystems, although Arctic ecosystems and light taiga generally have higher ground heat flux because there is less leaf and stem area to shade the ground surface, and the thermal gradient from the surface to permafrost is steeper. Permafrost creates a strong heat sink in summer that reduces surface temperature and therefore heat flux to the atmosphere. Loss of permafrost would therefore amplify climatic warming. If warming caused an increase in productivity and leaf area, or fire caused a shift from evergreen to deciduous forest, this would increase QE and reduce QH. Potential future shifts in vegetation would have varying climate feedbacks, with largest effects caused by shifts from boreal conifer to shrubland or deciduous forest (or vice versa) and from Arctic coastal to wet tundra. An increase of logging activity in the boreal forests appears to reduce QE by roughly 50% with little change in QH, while the ground heat flux is strongly enhanced.

The Greenland Ice Sheet's surface mass balance in a seasonally sea-ice free Arctic

General circulation models predict a rapid decrease in sea ice extent with concurrent increases in near surface air temperature and precipitation in the Arctic over the 21st century. This has led to suggestions that some Arctic land ice masses may experience an increase in accumulation due to enhanced evaporation from a seasonally sea ice free Arctic Ocean. To investigate the impact of this phenomenon on Greenland ice sheet climate and surface mass balance (SMB) a regional climate model, HadRM3, was used to force an insolation-temperature melt SMB model. A set of experiments designed to investigate the role of sea ice independently from sea surface temperature (SST) forcing are described. In the warmer and wetter SI + SST simulation Greenland experiences a 23% increase in winter SMB but 65% reduced summer SMB, resulting in a net decrease in the annual value. This study shows that sea ice decline contributes to the increased winter balance, causing 25% of the increase in winter accumulation; this is largest in eastern Greenland as the result of increased evaporation in the Greenland Sea. These results indicate that the seasonal cycle of Greenland's SMB will increase dramatically as global temperatures increase, with the largest changes in temperature and precipitation occurring in winter. This demonstrates that the accurate prediction of changes in sea ice cover is important for predicting Greenland SMB and ice sheet evolution.

Late-Holocene changes in character and behaviour of land-terminating glaciers on James Ross Island, Antarctica

Virtually no information is available on the response of land-terminating Antarctic Peninsula glaciers to climate change on a centennial timescale. This paper analyses the topography, geomorphology and sedimentology of prominent moraines on James Ross Island, Antarctica, to determine geometric changes and to interpret glacier behaviour. The moraines are very likely due to a late-Holocene phase of advance and featured (1) shearing and thrusting within the snout, (2) shearing and deformation of basal sediment, (3) more supraglacial debris than at present and (4) short distances of sediment transport. Retreat of ∼100 m and thinning of 15–20 m has produced a loss of 0.1 km3 of ice. The pattern of surface lowering is asymmetric. These geometrical changes are suggested most simply to be due to a net negative mass balance caused by a drier climate. Comparisons of the moraines with the current glaciological surface structure of the glaciers permits speculation of a transition from a polythermal to a cold-based thermal regime. Small land-terminating glaciers in the northern Antarctic Peninsula region could be cooling despite a warming climate.

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.

An investigation into the impact of using various techniques to estimate Arctic surface air temperature anomalies

Time series of global and regional mean Surface Air Temperature (SAT) anomalies are a common metric used to estimate recent climate change. Various techniques can be used to create these time series from meteorological station data. The degree of difference arising from using five different techniques, based on existing temperature anomaly dataset techniques, to estimate Arctic SAT anomalies over land and sea ice were investigated using reanalysis data as a testbed. Techniques which interpolated anomalies were found to result in smaller errors than non-interpolating techniques relative to the reanalysis reference. Kriging techniques provided the smallest errors in estimates of Arctic anomalies and Simple Kriging was often the best kriging method in this study, especially over sea ice. A linear interpolation technique had, on average, Root Mean Square Errors (RMSEs) up to 0.55 K larger than the two kriging techniques tested. Non-interpolating techniques provided the least representative anomaly estimates. Nonetheless, they serve as useful checks for confirming whether estimates from interpolating techniques are reasonable. The interaction of meteorological station coverage with estimation techniques between 1850 and 2011 was simulated using an ensemble dataset comprising repeated individual years (1979-2011). All techniques were found to have larger RMSEs for earlier station coverages. This supports calls for increased data sharing and data rescue, especially in sparsely observed regions such as the Arctic.

Land use and climate change impacts on soil organic carbon stocks in semi-arid Spain

Purpose The sensitivity of soil organic carbon to global change drivers, according to the depth profile, is receiving increasing attention because of its importance in the global carbon cycle and its potential feedback to climate change. A better knowledge of the vertical distribution of SOC and its controlling factors—the aim of this study—will help scientists predict the consequences of global change. Materials and methods The study area was the Murcia Province (S.E. Spain) under semiarid Mediterranean conditions. The database used consists of 312 soil profiles collected in a systematic grid, each 12 km2 covering a total area of 11,004 km2. Statistical analysis to study the relationships between SOC concentration and control factors in different soil use scenarios was conducted at fixed depths of 0–20, 20–40, 40–60, and 60–100 cm. Results and discussion SOC concentration in the top 40 cm ranged between 6.1 and 31.5 g kg−1, with significant differences according to land use, soil type and lithology, while below this depth, no differences were observed (SOC concentration 2.1–6.8 g kg−1). The ANOVA showed that land use was the most important factor controlling SOC concentration in the 0–40 cm depth. Significant differences were found in the relative importance of environmental and textural factors according to land use and soil depth. In forestland, mean annual precipitation and texture were the main predictors of SOC, while in cropland and shrubland, the main predictors were mean annual temperature and lithology. Total SOC stored in the top 1 m in the region was about 79 Tg with a low mean density of 7.18 kg Cm−3. The vertical distribution of SOC was shallower in forestland and deeper in cropland. A reduction in rainfall would lead to SOC decrease in forestland and shrubland, and an increase of mean annual temperature would adversely affect SOC in croplands and shrubland. With increasing depth, the relative importance of climatic factors decreases and texture becomes more important in controlling SOC in all land uses. Conclusions Due to climate change, impacts will be much greater in surface SOC, the strategies for C sequestration should be focused on subsoil sequestration, which was hindered in forestland due to bedrock limitations to soil depth. In these conditions, sequestration in cropland through appropriate management practices is recommended.

Effects of land management on different forms of soil carbon in olive groves in Mediterranean areas

This study analyses soil organic carbon (SOC) and hot-water extractable carbon (HWC), both measures of soil quality, under different land management: (1) conventional tillage (CT); (2) CT plus the addition of oil mill waste alperujo (A); (3) CT plus the addition of oil mill waste olive leaves (L); (4) no tillage with chipped pruned branches (NT1); and (5) no tillage with chipped pruned branches and weeds (NT2); in a typical Mediterranean agricultural area; the olive groves of Andalucía, southern Spain. SOC values in CT, A, NT1 and NT2 decreased with depth, but in NT2 the surface horizon (0-5 cm) had higher values than the other treatments, 47% more than the average values in the other three soils. In L, SOC also decreased with depth, although there was an increase of 88.5% from the first (0-10 cm) to the second horizon (10-16 cm). Total SOC stock values were very similar under A (101.9 Mg ha−1), CT (101.7 Mg ha−1), NT1 (105.8 Mg ha−1) and NT2 (111.3 Mg ha−1, if we consider the same depth of the others). However, SOC under L was significantly higher (p < 0.05) at 250.2 Mg ha−1. HWC decreased with depth in A, CT and NT1. NT2 and L followed the same pattern as the other management types but with a higher value in the surface horizon (2.3 and 4.9 mg g−1 respectively). Overall, our results indicate that application of oil mill waste olive leaves under CT (L) is a good management practice to improve SOC and reduce waste.

Evidence of a topographic signal in surface soil moisture derived from ENVISAT ASAR wide swath data

The susceptibility of a catchment to flooding is affected by its soil moisture prior to an extreme rainfall event. While soil moisture is routinely observed by satellite instruments, results from previous work on the assimilation of remotely sensed soil moisture into hydrologic models have been mixed. This may have been due in part to the low spatial resolution of the observations used. In this study, the remote sensing aspects of a project attempting to improve flow predictions from a distributed hydrologic model by assimilating soil moisture measurements are described. Advanced Synthetic Aperture Radar (ASAR) Wide Swath data were used to measure soil moisture as, unlike low resolution microwave data, they have sufficient resolution to allow soil moisture variations due to local topography to be detected, which may help to take into account the spatial heterogeneity of hydrological processes. Surface soil moisture content (SSMC) was measured over the catchments of the Severn and Avon rivers in the South West UK. To reduce the influence of vegetation, measurements were made only over homogeneous pixels of improved grassland determined from a land cover map. Radar backscatter was corrected for terrain variations and normalized to a common incidence angle. SSMC was calculated using change detection. To search for evidence of a topographic signal, the mean SSMC from improved grassland pixels on low slopes near rivers was compared to that on higher slopes. When the mean SSMC on low slopes was 30–90%, the higher slopes were slightly drier than the low slopes. The effect was reversed for lower SSMC values. It was also more pronounced during a drying event. These findings contribute to the scant information in the literature on the use of high resolution SAR soil moisture measurement to improve hydrologic models.

Eemian and Holocene sea-surface conditions in the southern Black Sea: organic-walled dinoflagellate cyst record from core 22-GC3/3

In order to compare the sea-surface conditions in the Black Sea during the Holocene and Eemian, sapropelic parts of marine core 22-GC3 (42°13.53′N/36°29.55′E, 838 m water depth) were studied for organic-walled dinoflagellate cyst content. The record shows a change from freshwater/brackish assemblages (Pyxidinopsis psilata, Spiniferites cruciformis, and Caspidinium rugosum) to more marine assemblages (Lingulodinium machaerophorum and Spiniferites ramosus complex) during each interglacial, due to the inflow of saline Mediterranean water. The lacustrine–marine transitions in 22-GC3 occurred at ~ 8.3 cal kyr BP during the early Holocene and ~ 128 kyr BP during the early Eemian, slightly later compared to the onset of interglacial conditions on the adjacent land. Dinoflagellate cyst assemblages reveal higher sea-surface salinity (~ 28–30) (e.g. Spiniferites pachydermus, Bitectatodinium tepikiense, and Spiniferites mirabilis) around ~ 126.5–121 kyr BP in comparison to the Holocene (~ 15–20) as well as relatively high sea-surface temperature (e.g. Tuberculodinium vancampoae, S. pachydermus, and S. mirabilis) especially at ~ 127.6–125.3 kyr BP. Establishment of high sea-surface salinity during the Eemian correlates very well with reconstructed relatively high global sea-level and is explained as a combined effect of increased Mediterranean supply and high temperatures at the beginning of the last interglacial. The observed changes in the dinocyst record highlight the importance of nutrients for the composition of the Eemian and Holocene dinocyst assemblages.

Combining satellite observations and reanalysis energy transports to estimate global net surface energy fluxes 1985-2012

Two methods are developed to estimate net surface energy fluxes based upon satellite-based reconstructions of radiative fluxes at the top of atmosphere and the atmospheric energy tendencies and transports from the ERA-Interim reanalysis. Method 1 applies the mass adjusted energy divergence from ERA-Interim while method 2 estimates energy divergence based upon the net energy difference at the top of atmosphere and the surface from ERA-Interim. To optimise the surface flux and its variability over ocean, the divergences over land are constrained to match the monthly area mean surface net energy flux variability derived from a simple relationship between the surface net energy flux and the surface temperature change. The energy divergences over the oceans are then adjusted to remove an unphysical residual global mean atmospheric energy divergence. The estimated net surface energy fluxes are compared with other data sets from reanalysis and atmospheric model simulations. The spatial correlation coefficients of multi-annual means between the estimations made here and other data sets are all around 0.9. There are good agreements in area mean anomaly variability over the global ocean, but discrepancies in the trend over the eastern Pacific are apparent.

Representing the Arctic in global surface temperature time series of recent climate change

The Arctic is an important region in the study of climate change, but monitoring surface temperatures in this region is challenging, particularly in areas covered by sea ice. Here in situ, satellite and reanalysis data were utilised to investigate whether global warming over recent decades could be better estimated by changing the way the Arctic is treated in calculating global mean temperature. The degree of difference arising from using five different techniques, based on existing temperature anomaly dataset techniques, to estimate Arctic SAT anomalies over land and sea ice were investigated using reanalysis data as a testbed. Techniques which interpolated anomalies were found to result in smaller errors than non-interpolating techniques. Kriging techniques provided the smallest errors in anomaly estimates. Similar accuracies were found for anomalies estimated from in situ meteorological station SAT records using a kriging technique. Whether additional data sources, which are not currently utilised in temperature anomaly datasets, would improve estimates of Arctic surface air temperature anomalies was investigated within the reanalysis testbed and using in situ data. For the reanalysis study, the additional input anomalies were reanalysis data sampled at certain supplementary data source locations over Arctic land and sea ice areas. For the in situ data study, the additional input anomalies over sea ice were surface temperature anomalies derived from the Advanced Very High Resolution Radiometer satellite instruments. The use of additional data sources, particularly those located in the Arctic Ocean over sea ice or on islands in sparsely observed regions, can lead to substantial improvements in the accuracy of estimated anomalies. Decreases in Root Mean Square Error can be up to 0.2K for Arctic-average anomalies and more than 1K for spatially resolved anomalies. Further improvements in accuracy may be accomplished through the use of other data sources.

On the parameterization of turbulent surface fluxes over heterogeneous sea ice surfaces

Turbulent surface fluxes of momentum and sensible and latent heat as well as surface temperature, air temperature, air humidity, and wind speed were measured by the German Falcon research aircraft over the marginal ice zone (MIZ) of the northern Baltic Sea and the Fram Strait. Applying the bulk formulas and the stability functions to the measurements, the roughness lengths for momentum z0, sensible heat zT, and latent heat zq were calculated. As mean values over a wide range of sea ice conditions, we obtain z0 = 5 � 10�4 m, zT = 1 � 10�8 m, and zq = 1 � 10�7 m. These correspond to the following mean values (± standard deviations) of neutral transfer coefficients reduced to 10 m height, CDN10 = (1.9 ± 0.8) � 10�3, CHN10 = (0.9 ± 0.3) � 10�3, and CEN10 = (1.0 ± 0.2) � 10�3. An average ratio of z0/zT � 104 was observed over the range of 10�6 m < z0 < 10�2 m and differs from previously published results over compact sea ice (10�1 < z0/zT < 103). Other observational results over heterogeneous sea ice do not exist. However, our z0/zT ratio approximately agrees with observations over heterogeneous land surfaces. Flux parameterizations based on commonly used roughness lengths ratios (z0 = zT = zq) overestimate the surface heat fluxes compared to our measurements by more than 100%.

Temporal and spatial variability of surface fluxes over the ice edge zone in the northern Baltic Sea

Three land-fast ice stations (one of them was the Finnish research ice breaker Aranda) and the German research aircraft Falcon were applied to measure the turbulent and radiation fluxes over the ice edge zone in the northern Baltic Sea during the Baltic Air-Sea-Ice Study (BASIS) field experiment from 16 February to 6 March 1998. The temporal and spatial variability of the surface fluxes is discussed. Synoptic weather systems passed the experimental area in a rapid sequence and dominated the conditions (wind speed, airsurface temperature difference, cloud field) for the variability of the turbulent and radiation fluxes. At the ice stations, the largest upward sensible heat fluxes of about 100 Wm�2 were measured during the passage of a cold front when the air cooled faster (�5 K per hour) than the surface. The largest downward flux of about �200 Wm�2 occurred during warm air advection when the air temperature reached +10�C but the surface temperature remained at 0�C. Spatial variability of fluxes was observed from the small scale (scale of ice floes and open water spots) to the mesoscale (width of the ice edge zone). The degree of spatial variability depends on the synoptic situation: during melting conditions downward heat fluxes were the same over ice and open water, whereas during strong cold-air advection upward heat fluxes differed by more than 100 Wm�2. A remarkable amount of grey ice with intermediate surface temperature was observed. The ice in the Baltic Sea cannot be described by one ice type only.

Influence of ground surface characteristics on the mean radiant temperature in urban areas

The effect of variations in land cover on mean radiant surface temperature (Tmrt) is explored through a simple scheme developed within the radiation model SOLWEIG. Outgoing longwave radiation is parameterised using surface temperature observations on a grass and an asphalt surface, whereas outgoing shortwave radiation is modelled through variations in albedo for the different surfaces. The influence of surface materials on Tmrt is small compared to the effects of shadowing. Nevertheless, altering ground surface materials could contribute to a reduction on Tmrt to reduce the radiant load during heat-wave episodes in locations where shadowing is not an option. Evaluation of the new scheme suggests that despite its simplicity it can simulate the outgoing fluxes well, especially during sunny conditions. However, it underestimates at night and in shadowed locations. One grass surface used to develop the parameterisation, with very different characteristics compared to an evaluation grass site, caused Tmrt to be underestimated. The implications of using high resolution (e.g. 15 minutes) temporal forcing data under partly cloudy conditions are demonstrated even for fairly proximal sites.