The GlobMODEL demonstrator: assimilation of new satellite products in an operational meteorological center

As part of its Data User Element programme, the European Space Agency funded the GlobMODEL project which aimed at investigating the scientific, technical, and organizational issues associated with the use and exploitation of remotely-sensed observations, particularly from new sounders. A pilot study was performed as a "demonstrator" of the GlobMODEL idea, based on the use of new data, with a strong European heritage, not yet assimilated operationally. Two parallel assimilation experiments were performed, using either total column ozone or ozone profiles retrieved at the Royal Netherlands Meteorological Institute (KNMI) from the Ozone Monitoring Instrument (OMI). In both cases, the impact of assimilating OMI data in addition to the total ozone columns from the SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) on the European Centre for Medium Range Weather Forecasts (ECMWF) ozone analyses was assessed by means of independent measurements. We found that the impact of OMI total columns is mainly limited to the region between 20 and 80 hPa, and is particularly important at high latitudes in the Southern hemisphere where the stratospheric ozone transport and chemical depletion are generally difficult to model with accuracy. Furthermore, the assimilation experiments carried out in this work suggest that OMI DOAS (Differential Optical Absorption Spectroscopy) total ozone columns are on average larger than SCIAMACHY total columns by up to 3 DU, while OMI total columns derived from OMI ozone profiles are on average about 8 DU larger than SCIAMACHY total columns. At the same time, the demonstrator brought to light a number of issues related to the assimilation of atmospheric composition profiles, such as the shortcomings arising when the vertical resolution of the instrument is not properly accounted for in the assimilation. The GlobMODEL demonstrator accelerated scientific and operational utilization of new observations and its results - prompted ECMWF to start the operational assimilation of OMI total column ozone data.

Storm track sensitivity to sea surface temperature resolution in a regional atmosphere model

A high resolution regional atmosphere model is used to investigate the sensitivity of the North Atlantic storm track to the spatial and temporal resolution of the sea surface temperature (SST) data used as a lower boundary condition. The model is run over an unusually large domain covering all of the North Atlantic and Europe, and is shown to produce a very good simulation of the observed storm track structure. The model is forced at the lateral boundaries with 15–20 years of data from the ERA-40 reanalysis, and at the lower boundary by SST data of differing resolution. The impacts of increasing spatial and temporal resolution are assessed separately, and in both cases increasing the resolution leads to subtle, but significant changes in the storm track. In some, but not all cases these changes act to reduce the small storm track biases seen in the model when it is forced with low-resolution SSTs. In addition there are several clear mesoscale responses to increased spatial SST resolution, with surface heat fluxes and convective precipitation increasing by 10–20% along the Gulf Stream SST gradient.

A geomatics-based approach for the derivation of the spatial distribution of sediment transport processes in periglacial mountain environments

Within this paper modern techniques such as satellite image analysis and tools provided by geographic information systems (GIS.) are exploited in order to extend and improve existing techniques for mapping the spatial distribution of sediment transport processes. The processes of interest comprise mass movements such as solifluction, slope wash, dirty avalanches and rock- and boulder falls. They differ considerably in nature and therefore different approaches for the derivation of their spatial extent are required. A major challenge is addressing the differences between the comparably coarse resolution of the available satellite data (Landsat TM/ETM+, 30 in x 30 m) and the actual scale of sediment transport in this environment. A three-stepped approach has been developed which is based on the concept of Geomorphic Process Units (GPUs): parameterization, process area delineation and combination. Parameters include land cover from satellite data and digital elevation model derivatives. Process areas are identified using a hierarchical classification scheme utilizing thresholds and definition of topology. The approach has been developed for the Karkevagge in Sweden and could be successfully transferred to the Rabotsbekken catchment at Okstindan, Norway using similar input data. Copyright (C) 2008 John Wiley & Sons, Ltd.

Causes of the apparent scale independence of fractal indices associated with forest fragmentation in Bolivia

In a recent investigation, Landsat TM and ETM+ data were used to simulate different resolutions of remotely-sensed images (from 30 to 1100 m) and to analyze the effect of resolution on a range of landscape metrics associated with spatial patterns of forest fragmentation in Chapare, Bolivia since the mid-1980s. Whereas most metrics were found to be highly dependent on pixel size, several fractal metrics (DLFD, MPFD, and AWMPFD) were apparently independent of image resolution, in contradiction with a sizeable body of literature indicating that fractal dimensions of natural objects depend strongly on image characteristics. The present re-analysis of the Chapare images, using two alternative algorithms routinely used for the evaluation of fractal dimensions, shows that the values of the box-counting and information fractal dimensions are systematically larger, sometimes by as much as 85%, than the "fractal" indices DLFD, MPFD, and AWMFD for the same images. In addition, the geometrical fractal features of the forest and non-forest patches in the Chapare region strongly depend on the resolution of images used in the analysis. The largest dependency on resolution occurs for the box-counting fractal dimension in the case of the non-forest patches in 1993, where the difference between the 30 and I 100 m-resolution images corresponds to 24% of the full theoretical range (1.0 to 2.0) of the mass fractal dimension. The observation that the indices DLFD, MPFD, and AWMPFD, unlike the classical fractal dimensions, appear relatively unaffected by resolution in the case of the Chapare images seems due essentially to the fact that these indices are based on a heuristic, "non-geometric" approach to fractals. Because of their lack of a foundation in fractal geometry, nothing guarantees that these indices will be resolution-independent in general. (C) 2006 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.

Using the variogram to explore imagery of two different spatial resolutions

The resolution of remotely sensed data is becoming increasingly fine, and there are now many sources of data with a pixel size of 1 m x 1 m. This produces huge amounts of data that have to be stored, processed and transmitted. For environmental applications this resolution possibly provides far more data than are needed: data overload. This poses the question: how much is too much? We have explored two resolutions of data-20 in pixel SPOT data and I in pixel Computerized Airborne Multispectral Imaging System (CAMIS) data from Fort A. P. Hill (Virginia, USA), using the variogram of geostatistics. For both we used the normalized difference vegetation index (NDVI). Three scales of spatial variation were identified in both the SPOT and 1 in data: there was some overlap at the intermediate spatial scales of about 150 in and of 500 m-600 in. We subsampled the I in data and scales of variation of about 30 in and of 300 in were identified consistently until the separation between pixel centroids was 15 in (or 1 in 225pixels). At this stage, spatial scales of about 100m and 600m were described, which suggested that only now was there a real difference in the amount of spatial information available from an environmental perspective. These latter were similar spatial scales to those identified from the SPOT image. We have also analysed I in CAMIS data from Fort Story (Virginia, USA) for comparison and the outcome is similar.:From these analyses it seems that a pixel size of 20m is adequate for many environmental applications, and that if more detail is required the higher resolution data could be sub-sampled to a 10m separation between pixel centroids without any serious loss of information. This reduces significantly the amount of data that needs to be stored, transmitted and analysed and has important implications for data compression.

GC/multiple collector-ICPMS method for chlorine stable isotope analysis of chlorinated aliphatic hydrocarbons

Stable isotopic characterization of chlorine in chlorinated aliphatic pollution is potentially very valuable for risk assessment and monitoring remediation or natural attenuation. The approach has been underused because of the complexity of analysis and the time it takes. We have developed a new method that eliminates sample preparation. Gas chromatography produces individually eluted sample peaks for analysis. The He carrier gas is mixed with Ar and introduced directly into the torch of a multicollector ICPMS. The MC-ICPMS is run at a high mass resolution of >= 10 000 to eliminate interference of mass 37 ArH with Cl. The standardization approach is similar to that for continuous flow stable isotope analysis in which sample and reference materials are measured successively. We have measured PCE relative to a laboratory TCE standard mixed with the sample. Solvent samples of 200 nmol to 1.3 mu mol ( 24- 165 mu g of Cl) were measured. The PCE gave the same value relative to the TCE as measured by the conventional method with a precision of 0.12% ( 2 x standard error) but poorer precision for the smaller samples.

Diurnal cycle of land surface temperature in a desert encroachment zone as observed from satellites

Climate variability in the African Soudano-Sahel savanna zone has attracted much attention because of the persistence of anomalously low rainfall. Past efforts to monitor the climate of this region have focused on rainfall and vegetation conditions, while land surface temperature (LST) has received less attention. Remote sensing of LST is feasible and possible at global scale. Most remotely sensed estimates of LST are based on the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) that are limited in their ability to capture the full diurnal cycle. Although more frequent observations are available from past geostationary satellites, their spatial resolution is coarser than that of polar orbiting satellites. In this study, the improved capabilities of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the METEOSAT Second Generation (MSG) instrument are used to remotely sense the LST in the African Soudano-Sahel savanna zone at a resolution of 3 km and 15 minutes. In support of the Radiative Atmospheric Divergence using the ARM Mobile Facility (AMF), GERB and AMMA Stations (RADAGAST) project, African Monsoon Multidisciplinary Analyses (AMMA) project and the Department of Energy's Atmospheric Radiation Measurement (ARM) program, the ARM Mobile Facility was deployed during 2006 in this climatically sensitive region, thereby providing a unique opportunity to evaluate remotely sensed algorithms for deriving LST.

Elimination of the Greenland Ice Sheet in a High CO2 Climate

Projections of future global sea level depend on reliable estimates of changes in the size of polar ice sheets. Calculating this directly from global general circulation models (GCMs) is unreliable because the coarse resolution of 100 km or more is unable to capture narrow ablation zones, and ice dynamics is not usually taken into account in GCMs. To overcome these problems a high-resolution (20 km) dynamic ice sheet model has been coupled to the third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3). A novel feature is the use of two-way coupling, so that climate changes in the GCM drive ice mass changes in the ice sheet model that, in turn, can alter the future climate through changes in orography, surface albedo, and freshwater input to the model ocean. At the start of the main experiment the atmospheric carbon dioxide concentration was increased to 4 times the preindustrial level and held constant for 3000 yr. By the end of this period the Greenland ice sheet is almost completely ablated and has made a direct contribution of approximately 7 m to global average sea level, causing a peak rate of sea level rise of 5 mm yr-1 early in the simulation. The effect of ice sheet depletion on global and regional climate has been examined and it was found that apart from the sea level rise, the long-term effect on global climate is small. However, there are some significant regional climate changes that appear to have reduced the rate at which the ice sheet ablates.

Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

The expected imprint of flux rope geometry on suprathermal electrons in magnetic clouds

Magnetic clouds are a subset of interplanetary coronal mass ejections characterized by a smooth rotation in the magnetic field direction, which is interpreted as a signature of a magnetic flux rope. Suprathermal electron observations indicate that one or both ends of a magnetic cloud typically remain connected to the Sun as it moves out through the heliosphere. With distance from the axis of the flux rope, out toward its edge, the magnetic field winds more tightly about the axis and electrons must traverse longer magnetic field lines to reach the same heliocentric distance. This increased time of flight allows greater pitch-angle scattering to occur, meaning suprathermal electron pitch-angle distributions should be systematically broader at the edges of the flux rope than at the axis. We model this effect with an analytical magnetic flux rope model and a numerical scheme for suprathermal electron pitch-angle scattering and find that the signature of a magnetic flux rope should be observable with the typical pitch-angle resolution of suprathermal electron data provided ACE's SWEPAM instrument. Evidence of this signature in the observations, however, is weak, possibly because reconnection of magnetic fields within the flux rope acts to intermix flux tubes.

Fourier transform infrared spectra of H2CS and D2CS

Infrared spectra of thoformaldehyde, H2CS and D2CS, were observed in the gas phase at a resolution of better than 0.1 cm−1 from 4000 to 400 cm−1 using a Nicolet FTIR system. Vibrational band origins and rotational constants were determined for ν2, ν3, ν4, and ν6 of H2CS and for ν1, ν2, ν3, ν4, and ν6 of D2CS. The ν3, ν4, and ν6 bands of H2CS were analyzed as a set of three Coriolis interacting bands, and three Coriolis constants were determined; similarly the ν4 and ν6 bands of D2CS were analyzed as a pair of interacting bands and one Coriolis constant was determined. A general harmonic force field was determined, without constraints, to fit the vibrational wavenumbers, Coriolis constants, and centrifugal distortion constants. A zero-point (rz) structure was determined from the ground-state rotational constants, and the equilibrium (re) bond lengths were estimated.

The infrared vibration-rotation spectrum of trans and cis nitrous acid

Infrared spectra of the trans and the cis isomers of nitrous acid, both HONO and DONO, have been observed in the gas phase using a Fourier transform interferometer with a resolution of about 0.05 cm−1 from 4000 to 500 cm−1. Rotational analyses are reported on eleven of the fundamentals and some overtones.

The ν2, 2ν2, 3ν2, ν4 and ν2+ν4 bands of 15NH3

The ir absorption of gaseous 15NH3 between 510 and 3040 cm−1 was recorded with a resolution of 0.06 cm−1. The ν2, 2ν2, 3ν2, ν4, and ν2 + ν4 bands were measured and analyzed on the basis of the vibration-rotation Hamiltonian developed by V. Špirko, J. M. R. Stone, and D. Papoušek (J. Mol. Spectrosc. 60, 159–178 (1976)). A set of effective molecular parameters for the ν2 = 1, 2, 3 states was derived, which reproduced the transition frequencies within the accuracy of the experimental measurements. For ν4 and ν2 + ν4 bands the standard deviation of the calculated spectrum is about four times larger than the measurements accuracy: a similar result was found for ν4 in 14NH3 by Š. Urban et al. (J. Mol. Spectrosc. 79, 455–495 (1980)). This result suggests that the present treatment takes into account only the most significant part of the rovibration interaction in the doubly degenerate vibrational states of ammonia.

Vibration-rotation spectra and the harmonic force field of HOCl

Vibration-rotation spectra of HOCl have been measured at a resolution of 0.05 cm−1 to determine vibration rotation constants, and 35–37 Cl isotope shifts in the vibration frequencies. The spectrum of DOCl has also been recorded, and a preliminary analysis for the band origins has been made. The vibrational frequency data and centrifugal distortion constants have been used to determine the harmonic force field in a least-squares refinement; the force field obtained also gives a good fit to data on the vibrational contributions to the inertial defect. The equilibrium rotational constants of HOCl have been obtained, and an equilibrium structure has been estimated.

Second order Coriolis resonance in the v2, v5 Raman bands of CH3F

The vibration-rotation Raman spectrum of the ν2 and ν5 fundamentals of CH3F is reported, from 1320 to 1640 cm−1, with a resolution of about 0.3 cm−1. The Coriolis resonance between the two bands leads to many perturbation-allowed transitions. Where the resonance is still sufficiently weak that the quantum number K′ retains its meaning, perturbation-allowed transitions are observed for all values of ΔK from +4 to −4; in regions of strong resonance, however, we can only say that the observed transitions obey the selection rule Δ(k−l) = 0 or ±3. The spectrum has been analyzed by band contour simulation using a computer program based on exact diagonalization of the Hamiltonian within the ν2, ν5 vibrational levels, and improved vibration-rotation constants for these bands are reported. The relative magnitudes and relative sings of polarizability derivatives involved in these vibrations are also reported.