Variational data assimilation for parameter estimation: application to a simple morphodynamic model

Data assimilation is a sophisticated mathematical technique for combining observational data with model predictions to produce state and parameter estimates that most accurately approximate the current and future states of the true system. The technique is commonly used in atmospheric and oceanic modelling, combining empirical observations with model predictions to produce more accurate and well-calibrated forecasts. Here, we consider a novel application within a coastal environment and describe how the method can also be used to deliver improved estimates of uncertain morphodynamic model parameters. This is achieved using a technique known as state augmentation. Earlier applications of state augmentation have typically employed the 4D-Var, Kalman filter or ensemble Kalman filter assimilation schemes. Our new method is based on a computationally inexpensive 3D-Var scheme, where the specification of the error covariance matrices is crucial for success. A simple 1D model of bed-form propagation is used to demonstrate the method. The scheme is capable of recovering near-perfect parameter values and, therefore, improves the capability of our model to predict future bathymetry. Such positive results suggest the potential for application to more complex morphodynamic models.

Spatial calibration of an optical see-through head-mounted display

We present here a method for calibrating an optical see-through Head Mounted Display (HMD) using techniques usually applied to camera calibration (photogrammetry). Using a camera placed inside the HMD to take pictures simultaneously of a tracked object and features in the HMD display, we could exploit established camera calibration techniques to recover both the intrinsic and extrinsic properties of the~HMD (width, height, focal length, optic centre and principal ray of the display). Our method gives low re-projection errors and, unlike existing methods, involves no time-consuming and error-prone human measurements, nor any prior estimates about the HMD geometry.

Parameter sensitivity and predictive uncertainty in a new water quality model, Q(2)

A new dynamic model of water quality, Q(2), has recently been developed, capable of simulating large branched river systems. This paper describes the application of a generalized sensitivity analysis (GSA) to Q(2) for single reaches of the River Thames in southern England. Focusing on the simulation of dissolved oxygen (DO) (since this may be regarded as a proxy for the overall health of a river); the GSA is used to identify key parameters controlling model behavior and provide a probabilistic procedure for model calibration. It is shown that, in the River Thames at least, it is more important to obtain high quality forcing functions than to obtain improved parameter estimates once approximate values have been estimated. Furthermore, there is a need to ensure reasonable simulation of a range of water quality determinands, since a focus only on DO increases predictive uncertainty in the DO simulations. The Q(2) model has been applied here to the River Thames, but it has a broad utility for evaluating other systems in Europe and around the world.

Assessing the sources and magnitude of diurnal nitrate variability in the San Joaquin River (California) with an in-situ optical nitrate sensor and dual nitrate isotopes

We investigated diurnal nitrate (NO3-) concentration variability in the San Joaquin River using an in situ optical NO3- sensor and discrete sampling during a 5-day summer period characterized by high algal productivity. Dual NO3- isotopes (delta N-15(NO3) and delta O-18(NO3)) and dissolved oxygen isotopes (delta O-18(DO)) were measured over 2 days to assess NO3- sources and biogeochemical controls over diurnal time-scales. Concerted temporal patterns of dissolved oxygen (DO) concentrations and delta O-18(DO) were consistent with photosynthesis, respiration and atmospheric O-2 exchange, providing evidence of diurnal biological processes independent of river discharge. Surface water NO3- concentrations varied by up to 22% over a single diurnal cycle and up to 31% over the 5-day study, but did not reveal concerted diurnal patterns at a frequency comparable to DO concentrations. The decoupling of delta N-15(NO3) and delta O-18(NO3) isotopes suggests that algal assimilation and denitrification are not major processes controlling diurnal NO3- variability in the San Joaquin River during the study. The lack of a clear explanation for NO3- variability likely reflects a combination of riverine biological processes and time-varying physical transport of NO3- from upstream agricultural drains to the mainstem San Joaquin River. The application of an in situ optical NO3- sensor along with discrete samples provides a view into the fine temporal structure of hydrochemical data and may allow for greater accuracy in pollution assessment.

Trends in aerosol optical depth in the Russian Arctic and their links with synoptic climatology

Temporal and spatial variability of aerosol optical depth (AOD) are examined using observations of direct solar radiation in the Eurasian Arctic for 1940-1990. AOD is estimated using empirical methods for 14 stations located between 66.2 degrees N and 80.6 degrees N, from the Kara Sea to the Chukchi Sea. While AOD exhibits a well-known springtime maximum and summertime minimum at all stations, atmospheric turbidity is higher in spring in the western (Kara-Laptev) part of the Eurasian Arctic. Between June and August, the eastern (East Siberian-Chukchi) sector experiences higher transparency than the western part. A statistically significant positive trend in AOD was observed in the Kara-Laptev sector between the late 1950s and the early 1930s predominantly in spring when pollution-derived aerosol dominates the Arctic atmosphere but not in the eastern sector. Although all stations are remote, those with positive trends are located closer to the anthropogenic sources of air pollution. By contrast, a widespread decline in AOD was observed between 1982 and 1990 in the eastern Arctic in spring but was limited to two sites in the western Arctic. These results suggest that the post-1982 decline in anthropogenic emissions in Europe and the former Soviet Union has had a limited effect on aerosol load in the Arctic. The post-1982 negative trends in AOD in summer, when marine aerosol is present in the atmosphere, were more common in the west. The relationships between AOD and atmospheric circulation are examined using a synoptic climatology approach. In spring, AOD depends primarily on the strength and direction of air flow. Thus strong westerly and northerly flows result in low AOD values in the East Siberian-Chukchi sector. By contrast, strong southerly flow associated with the passage of depressions results in high A OD in the Kara-Laptev sector and trajectory analysis points to the contribution of industrial regions of the sub-Arctic. In summer, low pressure gradient or anticyclonic conditions result in high atmospheric turbidity. The frequency of this weather type has declined significantly since the early 1980s in the Kara-Laptev sector, which partly explains the decline in summer AOD values. (c) 2004 Elsevier B.V. All rights reserved.

Time-multiplexed Laguerre-Gaussian holographic optical tweezers for biological applications

A ferroelectric liquid crystal spatial light modulator is used to generate up to 24 independently controllable traps in a holographic optical tweezers system using time-multiplexed Fresnel zone plates. For use in biological applications, helical zone plates are used to generate Laguerre-Gaussian laser modes. The high speed switching of the ferroelectric device together with recent advances in computer technology enable fast, smooth movement of traps that can be independently controlled in real time. This is demonstrated by the trapping and manipulation of yeast cells and fungal spores. (c) 2006 Optical Society of America.

The significance of volcanic eruption strength and frequency for climate

A simple physical model of the atmospheric effects of large explosive volcanic eruptions is developed. Using only one input parameter - the initial amount of sulphur dioxide injected into the stratosphere - the global-average stratospheric optical-depth perturbation and surface temperature response are modelled. The simplicity of this model avoids issues of incomplete data (applicable to more comprehensive models), making it a powerful and useful tool for atmospheric diagnostics of this climate forcing mechanism. It may also provide a computationally inexpensive and accurate way of introducing volcanic activity into larger climate models. The modelled surface temperature response for an initial sulphur-dioxide injection, coupled with emission-history statistics, is used to demonstrate that the most climatically significant volcanic eruptions are those of sufficient explosivity to just reach into the stratosphere (and achieve longevity). This study also highlights the fact that this measure of significance is highly sensitive to the representation of the climatic response and the frequency data used, and that we are far from producing a definitive history of explosive volcanism for at least the past 1000 years. Given this high degree of uncertainty, these results suggest that eruptions that release around and above 0.1 Mt SO2 into the stratosphere have the maximum climatic impact.