Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Search space pruning and global optimization of multiple gravity assist trajectories with deep space manoeuvres

This paper deals with the design of optimal multiple gravity assist trajectories with deep space manoeuvres. A pruning method which considers the sequential nature of the problem is presented. The method locates feasible vectors using local optimization and applies a clustering algorithm to find reduced bounding boxes which can be used in a subsequent optimization step. Since multiple local minima remain within the pruned search space, the use of a global optimization method, such as Differential Evolution, is suggested for finding solutions which are likely to be close to the global optimum. Two case studies are presented.

An improved algorithm for generating global window brightness temperatures from multiple satellite infrared imagery

An improved algorithm for the generation of gridded window brightness temperatures is presented. The primary data source is the International Satellite Cloud Climatology Project, level B3 data, covering the period from July 1983 to the present. The algorithm rakes window brightness, temperatures from multiple satellites, both geostationary and polar orbiting, which have already been navigated and normalized radiometrically to the National Oceanic and Atmospheric Administration's Advanced Very High Resolution Radiometer, and generates 3-hourly global images on a 0.5 degrees by 0.5 degrees latitude-longitude grid. The gridding uses a hierarchical scheme based on spherical kernel estimators. As part of the gridding procedure, the geostationary data are corrected for limb effects using a simple empirical correction to the radiances, from which the corrected temperatures are computed. This is in addition to the application of satellite zenith angle weighting to downweight limb pixels in preference to nearer-nadir pixels. The polar orbiter data are windowed on the target time with temporal weighting to account for the noncontemporaneous nature of the data. Large regions of missing data are interpolated from adjacent processed images using a form of motion compensated interpolation based on the estimation of motion vectors using an hierarchical block matching scheme. Examples are shown of the various stages in the process. Also shown are examples of the usefulness of this type of data in GCM validation.

The indoor localization system using vision for multiple mobile robots

This paper presents the development of an indoor localization system using camera vision. The localization system has a capability to determine 2D coordinate (x, y) for a team of mobile robots, Miabot. The experimental results show that the system outperforms our existing sonar localizer both in accuracy and a precision.

Multiple phases of North African humidity recorded in lacustrine sediments from the Fazzan Basin, Libyan Sahara

The Fazzan Basin of south-west Libya is at present arid with less than 20 mm of rainfall per annum. However, regionally extensive limestones, lacustrine sands and coquina (fossiliferous carbonate rock) deposits show that the Fazzan Basin previously contained a large palaeolake, indicating that the climate in the past was more humid. Optically stimulated luminescence (OSL) dating techniques have been applied to key lacustrine deposits within the basin in an attempt to provide an internally consistent chronology for this humidity record. Results indicate that palaeolake sediments within the Fazzan Basin record a very long history of palacohydrological change, ranging from present day and conditions to humidity capable of sustaining a lake with an approximate area of 76,250 km(2). The existence of humid periods in mid oxygen isotope stage 5 and the early Holocene is confirmed. An older lacustrine event, tentatively correlated to oxygen isotope stage 11, is also recognized. In addition, evidence is presented for at least two humid phases beyond the age range over which the conventional OSL dating technique is applicable. This study demonstrates that OSL dating of palaeolake sediments within the Fazzan Basin offers the potential to provide a detailed record of North African humidity spanning several glacial-interglacial cycles. (c) 2006 Elsevier Ltd. All rights reserved.

Single versus multiple occupancy - Effects on toxicity parameters measured on Eisenia fetida in lead nitrate-treated soil

The mortality (7 and 14 d), weight change (7 and 14 d), and metal uptake of Eisenia fetida (Savigny, 1826) kept in Pb(NO3)(2)-treated Kettering loam soil in single- and multiple-occupancy (10 earthworms) test containers were determined. The number of earthworms to dry mass (g) ratio of soil was 1:50 in both sets of test containers. Lead concentrations were in the nominal range of 0 to 10,000 mg Pb/kg soil (mg/kg hereafter). Levels of mortality at a given concentration were statistically identical between the single- and multiple-occupancy tests, except at 1,800 mg/kg, at which significantly (p less than or equal to 0.05) more mortality occurred in the multiple-occupancy tests. Death of individual earthworms in the multiple-occupancy tests did not trigger death of the other earthworms in that soil. The LC50 values (concentration statistically likely to kill 50% of the population) were identical between the multiple- and single-occupancy soils: 2,662 mg/kg (2,598-2,984, 7 d) and 2,589 mg/kg (2,251-3,013, 14 d) for the multiple-occupancy soils and 2,827 mg/kg (2,443-3,168, both 7 and 14 d) for the single-occupancy soils (values in brackets represent the 95% confidence intervals). Data were insufficient to calculate the concentration statistically likely to reduce individual earthworm mass by 50% (EC50), but after 14 d, the decrease in earthworm weight in the 1,800 and 3,000 mg/kg tests was significantly greater in the multiple- than in the single-occupancy soils. At 1,000, 1,800, and 3,000 mg/kg tests, earthworm Pb tissue concentration was significantly (p less than or equal to 0.05) greater in earthworms from the multiple-occupancy soils. The presence of earthworms increased the NH3 content of the soil; earthworm mortality increased NH3 concentrations further but not to toxic levels.

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.