An ensemble approach to uncertainty estimation for satellite-based rainfall estimates

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.

Uncertainty in the estimation of potential evapotranspiration under climate change.

21st century climate change is projected to result in an intensification of the global hydrological cycle, but there is substantial uncertainty in how this will impact freshwater availability. A relatively overlooked aspect of this uncertainty pertains to how different methods of estimating potential evapotranspiration (PET) respond to changing climate. Here we investigate the global response of six different PET methods to a 2 °C rise in global mean temperature. All methods suggest an increase in PET associated with a warming climate. However, differences in PET climate change signal of over 100% are found between methods. Analysis of a precipitation/PET aridity index and regional water surplus indicates that for certain regions and GCMs, choice of PET method can actually determine the direction of projections of future water resources. As such, method dependence of the PET climate change signal is an important source of uncertainty in projections of future freshwater availability.

Diffusion of the synthetic pyrethroid permethrin into bed-sediments

Bed-sediments are a sink for many micro-organic contaminants in aquatic environments. The impact of toxic contaminants on benthic fauna often depends on their spatial distribution, and the fate of the parent compounds and their metabolites. The distribution of a synthetic pyrethroid, permethrin, a compound known to be toxic to aquatic invertebrates, was studied using river bed-sediments in lotic flume channels. trans/cis-Permethrin diagnostic ratios were used to quantify the photoisomerization of the trans isomer in water. Rates were affected by the presence of sediment particles and colloids when compared to distilled water alone. Two experiments in dark/light conditions with replicate channels were undertaken using natural sediment, previously contaminated with permethrin, to examine the effect of the growth of an algal biofilm at the sediment-water interface on diffusive fluxes of permethrin into the sediment. After 42 days, the bulk water was removed, allowing a fine sectioning of the sediment bed (i.e., every mm down to 5 mm and then 5-10 mm, then every 10 mm down to 50 mm). Permethrin was detected in all cases down to a depth of 5-10 mm, in agreement with estimates by the Millington and Quirk model, and measurements of concentrations in pore water produced a distribution coefficient (K-d) for each section, High K-d's were observed for the top layers, mainly as a result of high organic matter and specific surface area. Concentrations in the algal biofilm measured at the end of the experiment under light conditions, and increases in concentration in the top 1 mm of the sediment, demonstrated that algal/bacterial biofilm material was responsible for high K-d's at the sediment surface, and for the retardation of permethrin diffusion. This specific partition of permethrin to fine sediment particles and algae may enhance its threat to benthic invertebrates. In addition,the analysis of trans/cis-permethrin isomer ratios in sediment showed greater losses of trans-permethrin in the experiment under light conditions, which may have also resulted from enhanced biological activity at the sediment surface.

Uncertainties in annual riverine phosphorus load estimation: Impact of load estimation methodology, sampling frequency, baseflow index and catchment population density

Models developed to identify the rates and origins of nutrient export from land to stream require an accurate assessment of the nutrient load present in the water body in order to calibrate model parameters and structure. These data are rarely available at a representative scale and in an appropriate chemical form except in research catchments. Observational errors associated with nutrient load estimates based on these data lead to a high degree of uncertainty in modelling and nutrient budgeting studies. Here, daily paired instantaneous P and flow data for 17 UK research catchments covering a total of 39 water years (WY) have been used to explore the nature and extent of the observational error associated with nutrient flux estimates based on partial fractions and infrequent sampling. The daily records were artificially decimated to create 7 stratified sampling records, 7 weekly records, and 30 monthly records from each WY and catchment. These were used to evaluate the impact of sampling frequency on load estimate uncertainty. The analysis underlines the high uncertainty of load estimates based on monthly data and individual P fractions rather than total P. Catchments with a high baseflow index and/or low population density were found to return a lower RMSE on load estimates when sampled infrequently than those with a tow baseflow index and high population density. Catchment size was not shown to be important, though a limitation of this study is that daily records may fail to capture the full range of P export behaviour in smaller catchments with flashy hydrographs, leading to an underestimate of uncertainty in Load estimates for such catchments. Further analysis of sub-daily records is needed to investigate this fully. Here, recommendations are given on load estimation methodologies for different catchment types sampled at different frequencies, and the ways in which this analysis can be used to identify observational error and uncertainty for model calibration and nutrient budgeting studies. (c) 2006 Elsevier B.V. All rights reserved.