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.

Modelling nitrogen fluxes from the land surface to the coastal zone in European systems: the perspective of the INCA project

This paper describes the results and conclusions of the INCA (Integrated Nitrogen Model for European CAtchments) project and sets the findings in the context of the ELOISE (European Land-Ocean Interaction Studies) programme. The INCA project was concerned with the development of a generic model of the major factors and processes controlling nitrogen dynamics in European river systems, thereby providing a tool (a) to aid the scientific understanding of nitrogen transport and retention in catchments and (b) for river-basin management and policy-making. The findings of the study highlight the heterogeneity of the factors and processes controlling nitrogen dynamics in freshwater systems. Nonetheless, the INCA model was able to simulate the in-stream nitrogen concentrations and fluxes observed at annual and seasonal timescales in Arctic, Continental and Maritime-Temperate regimes. This result suggests that the data requirements and structural complexity of the INCA model are appropriate to simulate nitrogen fluxes across a wide range of European freshwater environments. This is a major requirement for the production of coupled fiver-estuary-coastal shelf models for the management of our aquatic environment. With regard to river-basin management, to achieve an efficient reduction in nutrient fluxes from the land to the estuarine and coastal zone, the model simulations suggest that management options must be adaptable to the prevailing environmental and socio-economic factors in individual catchments: 'Blanket approaches' to environmental policy appear too simple. (c) 2004 Elsevier B.V. All rights reserved.

Uses of airborne scanning laser altimetry for river flood modelling and coastal geomorphology

Two ongoing projects at ESSC that involve the development of new techniques for extracting information from airborne LiDAR data and combining this information with environmental models will be discussed. The first project in conjunction with Bristol University is aiming to improve 2-D river flood flow models by using remote sensing to provide distributed data for model calibration and validation. Airborne LiDAR can provide such models with a dense and accurate floodplain topography together with vegetation heights for parameterisation of model friction. The vegetation height data can be used to specify a friction factor at each node of a model’s finite element mesh. A LiDAR range image segmenter has been developed which converts a LiDAR image into separate raster maps of surface topography and vegetation height for use in the model. Satellite and airborne SAR data have been used to measure flood extent remotely in order to validate the modelled flood extent. Methods have also been developed for improving the models by decomposing the model’s finite element mesh to reflect floodplain features such as hedges and trees having different frictional properties to their surroundings. Originally developed for rural floodplains, the segmenter is currently being extended to provide DEMs and friction parameter maps for urban floods, by fusing the LiDAR data with digital map data. The second project is concerned with the extraction of tidal channel networks from LiDAR. These networks are important features of the inter-tidal zone, and play a key role in tidal propagation and in the evolution of salt-marshes and tidal flats. The study of their morphology is currently an active area of research, and a number of theories related to networks have been developed which require validation using dense and extensive observations of network forms and cross-sections. The conventional method of measuring networks is cumbersome and subjective, involving manual digitisation of aerial photographs in conjunction with field measurement of channel depths and widths for selected parts of the network. A semi-automatic technique has been developed to extract networks from LiDAR data of the inter-tidal zone. A multi-level knowledge-based approach has been implemented, whereby low level algorithms first extract channel fragments based mainly on image properties then a high level processing stage improves the network using domain knowledge. The approach adopted at low level uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels. The higher level processing includes a channel repair mechanism.

Anthropogenic nutrient sources and loads from a Mediterranean catchment into a coastal lagoon: Mar Menor, Spain

The Mar Menor is a coastal lagoon increasingly threatened by urban and agricultural pressures. The main watercourse draining into the lagoon is the Rambla del Albujón. A fortnightly campaign carried out over one annual cycle enabled us to characterize the treated urban sewage effluents and agricultural sources which contribute to the nutrient fluxes in the watercourse. Multivariate analysis provided information for establishing chemical signatures and for assessing the relative influence of the various sources on the water quality at the outlet. Mass balances were used to examine net gains and losses, and cross-correlations with rainfall to analyze climatic influence and control factors in the trends of the nutrient flux. The rainfall pattern was significantly cross-correlated with nitrate and phosphorus fluxes from agricultural sources, while fluctuations in the resident population explained the phosphorus flux trend in urban sources. 50% of dissolved inorganic nitrogen was from agricultural sources, while 70% of total phosphate and 91% of total organic carbon were from urban point sources. The net amounts of all the nutrients fell as a result of plant uptake and/or denitrification in the channel. The control of urban point sources (phosphorus-enriched) is suggested as a promptly action for improving the health of the coastal lagoon.

Understanding 2H/1H systematics of leaf wax n-alkanes in coastal plants at Stiffkey saltmarsh, Norfolk, UK

Interpretation of sedimentary n-alkyl lipid d2H data is complicated by a limited understanding of factors controlling interspecies variation in biomarker 2H/1H composition. To distinguish between the effects of interrelated environmental, physical and biochemical controls on the hydrogen isotope composition of n-alkyl lipids, we conducted linked d2H analyses of soil water, xylem water, leaf water and n-alkanes from a range of C3 and C4 plants growing at a UK saltmarsh (i) across multiple sampling sites, (ii) throughout the 2012 growing season, and (iii) at different times of the day. Soil waters varied isotopically by up to 35& depending on marsh sub-environment, and exhibited site-specific seasonal shifts in d2H up to a maximum of 31 per mil. Maximum interspecies variation in xylem water was 38 per mil, while leaf waters differed seasonally by a maximum of 29 per mil. Leaf wax n-alkane 2H/1H, however, consistently varied by over 100 per mil throughout the 2012 growing season, resulting in an interspecies range in the ewax/leaf water values of -79 per mil to –227 per mil. From the discrepancy in the magnitude of these isotopic differences, we conclude that mechanisms driving variation in the 2H/1H composition of leaf water, including (i) spatial changes in soil water 2H/1H, (ii) temporal changes in soil water 2H/1H, (iii) differences in xylem water 2H/1H, and (iv) differences in leaf water evaporative 2H-enrichment due to varied plant life forms, cannot explain the range of n-alkane d2H values we observed. Results from this study suggests that accurate reconstructions of palaeoclimate regimes from sedimentary n-alkane d2H require further research to constrain those biological mechanisms influencing species-specific differences in 2H/1H fractionation during lipid biosynthesis, in particular where plants have developed biochemical adaptations to water-stressed conditions. Understanding how these mechanisms interact with environmental conditions will be crucial to ensure accurate interpretation of hydrogen isotope signals from the geological record.

Sequential hydrocarbon biodegradation of a soil from arid coastal Australia oil-treated under laboratory controlled conditions

A general consistency in the sequential order of petroleum hydrocarbon reduction in previous biodegradation studies has led to the proposal of several molecularly based biodegradation scales. Few studies have investigated the biodegradation susceptibility of petroleum hydrocarbon products in soil media, however, and metabolic preferences can change with habitat type. A laboratory based study comprising gas chromatography–mass spectrometry (GC–MS) analysis of extracts of oil-treated soil samples incubated for up to 161 days was conducted to investigate the biodegradation of crude oil exposed to sandy soils of Barrow Island, home to both a Class ‘‘A” nature reserve and Australia’s largest on-shore oil field. Biodegradation trends of the hydrocarbon-treated soils were largely consistent with previous reports but some unusual behaviour was recognised both between and within hydrocarbon classes. For example, the n-alkanes persisted at trace levels from day 86 to 161 following the removal of typically more stable dimethyl naphthalenes and methyl phenanthrenes. The relative susceptibility to biodegradation of different di- tri- and tetramethylnaphthalene isomers also showed several features distinct from previous reports. The unique biodegradation behaviour of Barrow Is. soil likely reflects difference in microbial functioning with physiochemical variation in the environment. Correlation of molecular parameters, reduction rates of selected alkyl naphthalene isomers and CO2 respiration values with a delayed (61 d) oil-treated soil identified a slowing of biodegradation with microcosm incubation; a reduced function or population of incubated soil flora might also influence the biodegradation patterns observed.