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.

Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

Correlations between microphysical properties of large-scale semi-transparent cirrus and the state of the atmosphere

By making use of TOVS Path-B satellite retrievals and ECMWF reanalyses, correlations between bulk microphysical properties of large-scale semi-transparent cirrus (visible optical thickness between 0.7 and 3.8) and thermodynamic and dynamic properties of the surrounding atmosphere have been studied on a global scale. These clouds constitute about half of all high clouds. The global averages (from 60°N to 60°S) of mean ice crystal diameter, De, and ice water path (IWP) of these clouds are 55 μm and 30 g m−2, respectively. IWP of these cirrus is slightly increasing with cloud-top temperature, whereas De of cold cirrus does not depend on this parameter. Correlations between De and IWp of large-scale cirrus seem to be different in the midlatitudes and in the tropics. However, we observe in general stronger correlations between De and IWP and atmospheric humidity and winds deduced from the ECMWF reanalyses: De and IWP increase both with increasing atmospheric water vapour. There is also a good distinction between different dynamical situations: In humid situations, IWP is on average about 10 gm−2 larger in regions with strong large-scale vertical updraft only that in regions with strong large-scale horizontal winds only, whereas the mean De of cold large-scale cirrus decreases by about 10 μm if both strong large-scale updraft and horizontal winds are present.

The effects of environmental perturbation and measurement error on estimates of the shape parameter in the theta-logistic model of population regulation

The theta-logistic is a widely used generalisation of the logistic model of regulated biological processes which is used in particular to model population regulation. Then the parameter theta gives the shape of the relationship between per-capita population growth rate and population size. Estimation of theta from population counts is however subject to bias, particularly when there are measurement errors. Here we identify factors disposing towards accurate estimation of theta by simulation of populations regulated according to the theta-logistic model. Factors investigated were measurement error, environmental perturbation and length of time series. Large measurement errors bias estimates of theta towards zero. Where estimated theta is close to zero, the estimated annual return rate may help resolve whether this is due to bias. Environmental perturbations help yield unbiased estimates of theta. Where environmental perturbations are large, estimates of theta are likely to be reliable even when measurement errors are also large. By contrast where the environment is relatively constant, unbiased estimates of theta can only be obtained if populations are counted precisely Our results have practical conclusions for the design of long-term population surveys. Estimation of the precision of population counts would be valuable, and could be achieved in practice by repeating counts in at least some years. Increasing the length of time series beyond ten or 20 years yields only small benefits. if populations are measured with appropriate accuracy, given the level of environmental perturbation, unbiased estimates can be obtained from relatively short censuses. These conclusions are optimistic for estimation of theta. (C) 2008 Elsevier B.V All rights reserved.

Thermodynamic analysis of ferrous ion binding to Escherichia coli ferritin EcFtnA

Iron oxidation in the bacterial ferritin EcFtnA from Escherichia coli shows marked differences from its homologue human H-chain ferritin (HuHF). While the amino acid residues that constitute the dinuclear center in these proteins are highly conserved, EcFtnA has a third iron-binding site (C site) in close proximity to the dinuclear center that is seemingly responsible for these differences. Here, we describe the first thermodynamic study of Fe2+ binding to EcFtnA and its variants to determine the location of the primary ferrous ion-binding sites on the protein and to better understand the role of the third C site in iron binding. Isothermal titration calorimetric analyses of the wild-type protein reveal the presence of two main classes of binding sites in the pH range of 6.5-7.5, ascribed to Fe2+ binding, first at the A and then the B sites. Site-directed mutagenesis of ligands in the A, B, or C sites affects the apparent Fe2+-binding stoichiometries at the unaltered sites. The data imply some degree of inter- and intrasubunit negative cooperative interaction between sites. Unlike HuHF where only the A site initially binds Fe2+, both A and B sites in EcFtnA bind Fe2+, implying a role for the C site in influencing the binding of Fe2+ at the B site of the di-iron center of EcFtnA. The ITC equations describing a binding model for three classes of independent binding sites are reported here for the first time.