Changes in area and geodetic mass balance of small glaciers, polar Urals, Russia, 1950-2008

Changes in area of 30 small glaciers (mostly <1 km2) in the northern Polar Urals (67.5-68.25 °N) between 1953 and 2000 were assessed using historic aerial photography from 1953 and 1960, ASTER and panchromatic Landsat ETM+ imagery from 2000, and data from 1981 and 2008 terrestrial surveys. Changes in volume and geodetic mass balance of IGAN and Obruchev glaciers were calculated using data from terrestrial surveys in 1963 and 2008. In total, glacier area declined by 22.3 ± 3.9% in the 1953/60-2000 period. The areas of individual glaciers decreased by 4-46%. Surfaces of Obruchev and IGAN glaciers lowered by 22.5 ± 1.7 m and 14.9 ± 2.1 m. Over 45 years, geodetic mass balances of Obruchev and IGAN glaciers were -20.66 ± 2.91 and -13.54 ± 2.57 m w.e. respectively. Glacier shrinkage in the Polar Urals is related to a summer warming of 1 °C between 1953-81 and 1981-2008 and its rates are consistent with other regions of northern Asia but are higher than in Scandinavia. While glacier shrinkage intensified in the 1981-2000 period relative to 1953-81, increasing winter precipitation and shading effects slowed glacier wastage in 2000-08.

An Earth Observation Land Data Assimilation System (EO-LDAS)

Current methods for estimating vegetation parameters are generally sub-optimal in the way they exploit information and do not generally consider uncertainties. We look forward to a future where operational dataassimilation schemes improve estimates by tracking land surface processes and exploiting multiple types of observations. Dataassimilation schemes seek to combine observations and models in a statistically optimal way taking into account uncertainty in both, but have not yet been much exploited in this area. The EO-LDAS scheme and prototype, developed under ESA funding, is designed to exploit the anticipated wealth of data that will be available under GMES missions, such as the Sentinel family of satellites, to provide improved mapping of land surface biophysical parameters. This paper describes the EO-LDAS implementation, and explores some of its core functionality. EO-LDAS is a weak constraint variational dataassimilationsystem. The prototype provides a mechanism for constraint based on a prior estimate of the state vector, a linear dynamic model, and EarthObservationdata (top-of-canopy reflectance here). The observation operator is a non-linear optical radiative transfer model for a vegetation canopy with a soil lower boundary, operating over the range 400 to 2500 nm. Adjoint codes for all model and operator components are provided in the prototype by automatic differentiation of the computer codes. In this paper, EO-LDAS is applied to the problem of daily estimation of six of the parameters controlling the radiative transfer operator over the course of a year (> 2000 state vector elements). Zero and first order process model constraints are implemented and explored as the dynamic model. The assimilation estimates all state vector elements simultaneously. This is performed in the context of a typical Sentinel-2 MSI operating scenario, using synthetic MSI observations simulated with the observation operator, with uncertainties typical of those achieved by optical sensors supposed for the data. The experiments consider a baseline state vector estimation case where dynamic constraints are applied, and assess the impact of dynamic constraints on the a posteriori uncertainties. The results demonstrate that reductions in uncertainty by a factor of up to two might be obtained by applying the sorts of dynamic constraints used here. The hyperparameter (dynamic model uncertainty) required to control the assimilation are estimated by a cross-validation exercise. The result of the assimilation is seen to be robust to missing observations with quite large data gaps.

A comparison of the performance of Kalpana and Meteosat-7 WV radiances in the ECMWF NWP model

The Kalpana Very High Resolution Radiometer (VHRR) water vapour (WV) channel is very similar to the WV channel of the Meteosat Visible and Infrared Radiation Imager (MVIRI) on Meteosat-7, and both satellites observe the Indian subcontinent. Thus it is possible to compare the performance of VHRR and MVIRI in numerical weather prediction (NWP) models. In order to do so, the impact of Kalpana- and Meteosat-7-measured WV radiances was evaluated using analyses and forecasts of moisture, temperature, geopotential and winds, using the European Centre for Medium-range Weather Forecasts (ECMWF) NWP model. Compared with experiments using Meteosat-7, the experiments using Kalpana WV radiances show a similar fit to all observations and produce very similar forecasts.

On instabilities in data assimilation algorithms

Data assimilation algorithms are a crucial part of operational systems in numerical weather prediction, hydrology and climate science, but are also important for dynamical reconstruction in medical applications and quality control for manufacturing processes. Usually, a variety of diverse measurement data are employed to determine the state of the atmosphere or to a wider system including land and oceans. Modern data assimilation systems use more and more remote sensing data, in particular radiances measured by satellites, radar data and integrated water vapor measurements via GPS/GNSS signals. The inversion of some of these measurements are ill-posed in the classical sense, i.e. the inverse of the operator H which maps the state onto the data is unbounded. In this case, the use of such data can lead to significant instabilities of data assimilation algorithms. The goal of this work is to provide a rigorous mathematical analysis of the instability of well-known data assimilation methods. Here, we will restrict our attention to particular linear systems, in which the instability can be explicitly analyzed. We investigate the three-dimensional variational assimilation and four-dimensional variational assimilation. A theory for the instability is developed using the classical theory of ill-posed problems in a Banach space framework. Further, we demonstrate by numerical examples that instabilities can and will occur, including an example from dynamic magnetic tomography.

‘NSC 68 and the Soviet Threat: A New Perspective on Western Threat Perception and Policy Making’

It has been argued that the US strategy paper NSC 68 of 1950 which ushered in the great Cold War rearmament process first in the US and then in NATO was a vast overreaction. This paper argues, by contrast, that given the intelligence about the Soviet Union's and its satellites' military buildup and the role of China in that period, this was a reasonable strategy to embrace.

Accelerated loss of alpine glaciers in the Kodar Mountains, south-eastern Siberia

The recession of mountain glaciers around the world has been linked to anthropogenic climate change and small glaciers (e.g. < 2 km2) are thought to be particularly vulnerable, with reports of their disappearance from several regions. However, the response of small glaciers to climate change can be modulated by non-climatic factors such as topography and debris cover and there remain a number of regions where their recent change has evaded scrutiny. This paper presents results of the first multi-year remote sensing survey of glaciers in the Kodar Mountains, the only glaciers in SE Siberia, which we compare to previous glacier inventories from this continental setting that reported total glacier areas of 18.8 km2 in ca. 1963 (12.6 km2 of exposed ice) and 15.5 km2 in 1974 (12 km2 of exposed ice). Mapping their debris-covered termini is difficult but delineation of debris-free ice on Landsat imagery reveals 34 glaciers with a total area of 11.72 ± 0.72 km2 in 1995, followed by a reduction to 9.53 ± 0.29 km2 in 2001 and 7.01 ± 0.23 km2 in 2010. This represents a ~ 44% decrease in exposed glacier ice between ca. 1963 and 2010, but with 40% lost since 1995 and with individual glaciers losing as much as 93% of their exposed ice. Thus, although continental glaciers are generally thought to be less sensitive than their maritime counterparts, a recent acceleration in shrinkage of exposed ice has taken place and we note its coincidence with a strong summer warming trend in the region initiated at the start of the 1980s. Whilst smaller and shorter glaciers have, proportionally, tended to shrink more rapidly, we find no statistically significant relationship between shrinkage and elevation characteristics, aspect or solar radiation. This is probably due to the small sample size, limited elevation range, and topographic setting of the glaciers in deep valleys-heads. Furthermore, many of the glaciers possess debris-covered termini and it is likely that the ablation of buried ice is lagging the shrinkage of exposed ice, such that a growth in the proportion of debris cover is occurring, as observed elsewhere. If recent trends continue, we hypothesise that glaciers could evolve into a type of rock glacier within the next few decades, introducing additional complexity in their response and delaying their potential demise.

Foreword: International Space Science Institute (ISSI) workshop on observing and predicting Earth’s energy flows

The Earth’s climate, as well as planetary climates in general, is broadly regulated by three fundamental parameters: the total solar irradiance, the planetary albedo and the planetary emissivity. Observations from series of different satellites during the last three decades indicate that these three quantities are generally very stable. The total solar irradiation of some 1,361 W/m2 at 1 A.U. varies within 1 W/m2 during the 11-year solar cycle (Fröhlich 2012). The albedo is close to 29 % with minute changes from year to year but with marked zonal differences (Stevens and Schwartz 2012). The only exception to the overall stability is a minor decrease in the planetary emissivity (the ratio between the radiation to space and the radiation from the surface of the Earth). This is a consequence of the increase in atmospheric greenhouse gas amounts making the atmosphere gradually more opaque to long-wave terrestrial radiation. As a consequence, radiation processes are slightly out of balance as less heat is leaving the Earth in the form of thermal radiation than the amount of heat from the incoming solar radiation. Present space-based systems cannot yet measure this imbalance, but the effect can be inferred from the increase in heat in the oceans where most of the heat accumulates. Minor amounts of heat are used to melt ice and to warm the atmosphere and the surface of the Earth.

Interpreting Differential Temperature Trends at the Surface and in the Lower Troposphere

Estimated global-scale temperature trends at Earth's surface (as recorded by thermometers) and in the lower troposphere (as monitored by satellites) diverge by up to 0.14°C per decade over the period 1979 to 1998. Accounting for differences in the spatial coverage of satellite and surface measurements reduces this differential, but still leaves a statistically significant residual of roughly 0.1°C per decade. Natural internal climate variability alone, as simulated in three state-of-the-art coupled atmosphere-ocean models, cannot completely explain this residual trend difference. A model forced by a combination of anthropogenic factors and volcanic aerosols yields surface-troposphere temperature trend differences closest to those observed.

Assessing levels of uncertainty in recent temperature time series

We examine to what degree we can expect to obtain accurate temperature trends for the last two decades near the surface and in the lower troposphere. We compare temperatures obtained from surface observations and radiosondes as well as satellite-based measurements from the Microwave Soundings Units (MSU), which have been adjusted for orbital decay and non-linear instrument-body effects, and reanalyses from the European Centre for Medium-Range Weather Forecasts (ERA) and the National Centre for Environmental Prediction (NCEP). In regions with abundant conventional data coverage, where the MSU has no major influence on the reanalysis, temperature anomalies obtained from microwave sounders, radiosondes and from both reanalyses agree reasonably. Where coverage is insufficient, in particular over the tropical oceans, large differences are found between the MSU and either reanalysis. These differences apparently relate to changes in the satellite data availability and to differing satellite retrieval methodologies, to which both reanalyses are quite sensitive over the oceans. For NCEP, this results from the use of raw radiances directly incorporated into the analysis, which make the reanalysis sensitive to changes in the underlying algorithms, e.g. those introduced in August 1992. For ERA, the bias-correction of the one-dimensional variational analysis may introduce an error when the satellite relative to which the correction is calculated is biased itself or when radiances change on a time scale longer than a couple of months, e.g. due to orbit decay. ERA inhomogeneities are apparent in April 1985, October/November 1986 and April 1989. These dates can be identified with the replacements of satellites. It is possible that a negative bias in the sea surface temperatures (SSTs) used in the reanalyses may have been introduced over the period of the satellite record. This could have resulted from a decrease in the number of ship measurements, a concomitant increase in the importance of satellite-derived SSTs, and a likely cold bias in the latter. Alternately, a warm bias in SSTs could have been caused by an increase in the percentage of buoy measurements (relative to deeper ship intake measurements) in the tropical Pacific. No indications for uncorrected inhomogeneities of land surface temperatures could be found. Near-surface temperatures have biases in the boundary layer in both reanalyses, presumably due to the incorrect treatment of snow cover. The increase of near-surface compared to lower tropospheric temperatures in the last two decades may be due to a combination of several factors, including high-latitude near-surface winter warming due to an enhanced NAO and upper-tropospheric cooling due to stratospheric ozone decrease.

Assimilation of non-synoptic observations

A system for continuous data assimilation described recently (Bengtsson & Gustavsson, 1971) has been further developed and tested under more realistic conditions. A balanced barotropic model is used and the integration is performed over an octagon covering the area to the north of 20° N. Comparisons have been made between using data from the actual aerological network and data from a satellite in a polar orbit. The result of the analyses has been studied in different subregions situated in data sparse as well as in data dense areas. The errors of the analysis have also been studied in the wave spectrum domain. Updating is performed using data generated by the model but also by model-independent data. Rather great differences are obtained between the two experiments especially with respect to the ultra-long waves. The more realistic approach gives much larger analysis error. In general the satellite updating yields somewhat better result than the updating from the conventional aerological network especially in the data sparse areas over the oceans. Most of the experiments are performed by a satellite making 200 observations/track, a sidescan capability of 40° and with a RMS-error of 20 m. It is found that the effect of increasing the number of satellite observations from 100 to 200 per orbit is almost negligible. Similarly the effect is small of improving the observations by diminishing the RMS-error below a certain value. An observing system using two satellites 90° out of phase has also been investigated. This is found to imply a substantial improvement. Finally an experiment has been performed using actual SIRS-soundings from NIMBUS IV. With respect to the very small number of soundings at 500 mb, 142 during 48 hours, the result can be regarded as quite satisfactory.

Numerical prediction of atmospheric blocking-A case study

A series of numerical models have been used to investigate the predictability of atmospheric blocking for an episode selected from FGGE Special Observing Period I. Level II-b FGGE data have been used in the experiment. The blocking took place over the North Atlantic region and is a very characteristic example of high winter blocking. It is found that the very high resolution models developed at ECMWF, in a remarkable way manage to predict the blocking event in great detail, even beyond 1 week. Although models with much less resolution manage to predict the blocking phenomenon as such, the actual evolution differs very much from the observed and consequently the practical value is substantially reduced. Wind observations from the geostationary satellites are shown to have a substantial impact on the forecast beyond 5 days, as well as an extension of the integration domain to the whole globe. Quasi-geostrophic baroclinic models and, even more, barotropic models, are totally inadequate to predict blocking except in its initial phase. The prediction experiment illustrates clearly that efforts which have gone into the improvement of numerical prediction models in the last decades have been worth while.

Scheduling satellite-based SAR acquisition for sequential assimilation of water level observations into flood modelling

Satellite-based Synthetic Aperture Radar (SAR) has proved useful for obtaining information on flood extent, which, when intersected with a Digital Elevation Model (DEM) of the floodplain, provides water level observations that can be assimilated into a hydrodynamic model to decrease forecast uncertainty. With an increasing number of operational satellites with SAR capability, information on the relationship between satellite first visit and revisit times and forecast performance is required to optimise the operational scheduling of satellite imagery. By using an Ensemble Transform Kalman Filter (ETKF) and a synthetic analysis with the 2D hydrodynamic model LISFLOOD-FP based on a real flooding case affecting an urban area (summer 2007,Tewkesbury, Southwest UK), we evaluate the sensitivity of the forecast performance to visit parameters. We emulate a generic hydrologic-hydrodynamic modelling cascade by imposing a bias and spatiotemporal correlations to the inflow error ensemble into the hydrodynamic domain. First, in agreement with previous research, estimation and correction for this bias leads to a clear improvement in keeping the forecast on track. Second, imagery obtained early in the flood is shown to have a large influence on forecast statistics. Revisit interval is most influential for early observations. The results are promising for the future of remote sensing-based water level observations for real-time flood forecasting in complex scenarios.

Intensification of winter transatlantic aviation turbulence in response to climate change

Atmospheric turbulence causes most weather-related aircraft incidents1. Commercial aircraft encounter moderate-or-greater turbulence tens of thousands of times each year worldwide, injuring probably hundreds of passengers (occasionally fatally), costing airlines tens of millions of dollars and causing structural damage to planes1, 2, 3. Clear-air turbulence is especially difficult to avoid, because it cannot be seen by pilots or detected by satellites or on-board radar4, 5. Clear-air turbulence is linked to atmospheric jet streams6, 7, which are projected to be strengthened by anthropogenic climate change8. However, the response of clear-air turbulence to projected climate change has not previously been studied. Here we show using climate model simulations that clear-air turbulence changes significantly within the transatlantic flight corridor when the concentration of carbon dioxide in the atmosphere is doubled. At cruise altitudes within 50–75° N and 10–60° W in winter, most clear-air turbulence measures show a 10–40% increase in the median strength of turbulence and a 40–170% increase in the frequency of occurrence of moderate-or-greater turbulence. Our results suggest that climate change will lead to bumpier transatlantic flights by the middle of this century. Journey times may lengthen and fuel consumption and emissions may increase. Aviation is partly responsible for changing the climate9, but our findings show for the first time how climate change could affect aviation.

Ultralow-frequency modulation of whistler-mode wave growth

Measurements from ground-based magnetometers and riometers at auroral latitudes have demonstrated that energetic (~30-300keV) electron precipitation can be modulated in the presence of magnetic field oscillations at ultra-low frequencies. It has previously been proposed that an ultra-low frequency (ULF) wave would modulate field and plasma properties near the equatorial plane, thus modifying the growth rates of whistler-mode waves. In turn, the resulting whistler-mode waves would mediate the pitch-angle scattering of electrons resulting in ionospheric precipitation. In this paper, we investigate this hypothesis by quantifying the changes to the linear growth rate expected due to a slow change in the local magnetic field strength for parameters typical of the equatorial region around 6.6RE radial distance. To constrain our study, we determine the largest possible ULF wave amplitudes from measurements of the magnetic field at geosynchronous orbit. Using nearly ten years of observations from two satellites, we demonstrate that the variation in magnetic field strength due to oscillations at 2mHz does not exceed ±10% of the background field. Modifications to the plasma density and temperature anisotropy are estimated using idealised models. For low temperature anisotropy, there is little change in the whistler-mode growth rates even for the largest ULF wave amplitude. Only for large temperature anisotropies can whistler-mode growth rates be modulated sufficiently to account for the changes in electron precipitation measured by riometers at auroral latitudes.