Efficient accommodation of local minima in watershed model calibration

The Gauss-Marquardt-Levenberg (GML) method of computer-based parameter estimation, in common with other gradient-based approaches, suffers from the drawback that it may become trapped in local objective function minima, and thus report optimized parameter values that are not, in fact, optimized at all. This can seriously degrade its utility in the calibration of watershed models where local optima abound. Nevertheless, the method also has advantages, chief among these being its model-run efficiency, and its ability to report useful information on parameter sensitivities and covariances as a by-product of its use. It is also easily adapted to maintain this efficiency in the face of potential numerical problems (that adversely affect all parameter estimation methodologies) caused by parameter insensitivity and/or parameter correlation. The present paper presents two algorithmic enhancements to the GML method that retain its strengths, but which overcome its weaknesses in the face of local optima. Using the first of these methods an intelligent search for better parameter sets is conducted in parameter subspaces of decreasing dimensionality when progress of the parameter estimation process is slowed either by numerical instability incurred through problem ill-posedness, or when a local objective function minimum is encountered. The second methodology minimizes the chance of successive GML parameter estimation runs finding the same objective function minimum by starting successive runs at points that are maximally removed from previous parameter trajectories. As well as enhancing the ability of a GML-based method to find the global objective function minimum, the latter technique can also be used to find the locations of many non-global optima (should they exist) in parameter space. This can provide a useful means of inquiring into the well-posedness of a parameter estimation problem, and for detecting the presence of bimodal parameter and predictive probability distributions. The new methodologies are demonstrated by calibrating a Hydrological Simulation Program-FORTRAN (HSPF) model against a time series of daily flows. Comparison with the SCE-UA method in this calibration context demonstrates a high level of comparative model run efficiency for the new method. (c) 2006 Elsevier B.V. All rights reserved.

Individual calibration for estimating free-living walking speed using the MTI monitor

Purpose: This study was conducted to devise a new individual calibration method to enhance MTI accelerometer estimation of free-living level walking speed. Method: Five female and five male middle-aged adults walked 400 m at 3.5, 4.5, and 5.5 km(.)h(-1), and 800 in at 6.5 km(.)h(-1) on an outdoor track, following a continuous protocol. Lap speed was controlled by a global positioning system (GPS) monitor. MTI counts-to-speed calibration equations were derived for each trial, for each subject for four such trials with each of four MTI, for each subject for the average MTI. and for the pooled data. Standard errors of the estimate (SEE) with and without individual calibration were compared. To assess accuracy of prediction of free-living walking speed, subjects also completed a self-paced, brisk 3-km walk wearing one of the four MTI, and differences between actual and predicted walking speed with and without individual calibration were examined. Results: Correlations between MTI counts and walking speed were 0.90 without individual calibration, 0.98 with individual calibration for the average MTI. and 0.99 with individual calibration for a specific MTI. The SEE (mean +/- SD) was 0.58 +/- 0.30 km(.)h(-1) without individual calibration, 0.19 +/- 0.09 km h(-1) with individual calibration for the average MTI monitor, and 0.16 +/- 0.08 km(.)h(-1) with individual calibration for a specific MTI monitor. The difference between actual and predicted walking speed on the brisk 3-km walk was 0.06 +/- 0.25 km(.)h(-1) using individual calibration and 0.28 +/- 0.63 km(.)h(-1) without individual calibration (for specific accelerometers). Conclusion: MTI accuracy in predicting walking speed without individual calibration might be sufficient for population-based studies but not for intervention trials. This individual calibration method will substantially increase precision of walking speed predicted from MTI counts.

Predictive validity of three ActiGraph energy expenditure equations for children

Purpose: This Study evaluated the predictive validity of three previously published ActiGraph energy expenditure (EE) prediction equations developed for children and adolescents. Methods: A total of 45 healthy children and adolescents (mean age: 13.7 +/- 2.6 yr) completed four 5-min activity trials (normal walking. brisk walking, easy running, and fast running) in ail indoor exercise facility. During each trial, participants were all ActiGraph accelerometer oil the right hip. EE was monitored breath by breath using the Cosmed K4b(2) portable indirect calorimetry system. Differences and associations between measured and predicted EE were assessed using dependent t-tests and Pearson correlations, respectively. Classification accuracy was assessed using percent agreement, sensitivity, specificity, and area under the receiver operating characteristic (ROC) curve, Results: None of the equations accurately predicted mean energy expenditure during each of the four activity trials. Each equation, however, accurately predicted mean EE in at least one activity trial. The Puyau equation accurately predicted EE during slow walking. The Trost equation accurately predicted EE during slow running. The Freedson equation accurately predicted EE during fast running. None of the three equations accurately predicted EE during brisk walking. The equations exhibited fair to excellent classification accuracy with respect to activity intensity. with the Trost equation exhibiting the highest classification accuracy and the Puyau equation exhibiting the lowest. Conclusions: These data suggest that the three accelerometer prediction equations do not accurately predict EE on a minute-by-minute basis in children and adolescents during overground walking and running. The equations maybe, however, for estimating participation in moderate and vigorous activity.

Accuracy of strain rate techniques for ientification of viability at dobutamine stress echo: A follow-up study after revascularization

Background. Myocardial viability (VM) assessment based on wall motion scoring (WMS) with dobutamine echo (DbE) is difficult and subjective. New quantitative techniques such as strain rate imaging (SRI) correspond with isotopic techniques but their ability to predict functional recovery (FR) after revascularization is unclear. Methods. Stable post-MI pts (n=43, age 63±9, EF 36±6%) underwent SRI during DbE. WMS evidence of VM was based on lowdose augmentation at DbE. SR, end-systolic strain (ESS), post-systolic strain (PSS) and timing were analyzed at rest and low dose in abnormal segts. Pts were followed for 9±12 months; FR was defined as segt improvement on post-revascularization images. Results: Of 180 segts with abnormal resting function, 83 showed FR and 97 did not. Resting parameters were not predictive of recovery; resting post-systolic shortening had a sensitivity and specificity

Major Issues in Business Process Management: An Expert Perspective

The results presented in this report form a part of a larger global study on the major issues in BPM. Only one part of the larger study is reported here, viz. interviews with BPM experts. Interviews of BPM tool vendors together with focus groups involving user organizations, are continuing in parallel and will set the groundwork for the identification of BPM issues on a global scale via a survey (including a Delphi study). Through this multi-method approach, we identify four distinct sets of outcomes. First, as is the focus of this report, we identify the BPM issues as perceived by BPM experts. Second, the research design allows us to gain insight into the opinions of organisations deploying BPM solutions. Third, an understanding of organizations’ misconceptions of BPM technologies, as confronted by BPM tool vendors is obtained. Last, we seek to gain an understanding of BPM issues on a global scale, together with knowledge of matters of concern. This final outcome is aimed to produce an industry driven research agenda which will inform practitioners and in particular, the research community world-wide on issues and challenges that are prevalent or emerging in BPM and related areas.

A Novel Nonparametric Density Estimator

We present a novel nonparametric density estimator and a new data-driven bandwidth selection method with excellent properties. The approach is in- spired by the principles of the generalized cross entropy method. The pro- posed density estimation procedure has numerous advantages over the tra- ditional kernel density estimator methods. Firstly, for the first time in the nonparametric literature, the proposed estimator allows for a genuine incor- poration of prior information in the density estimation procedure. Secondly, the approach provides the first data-driven bandwidth selection method that is guaranteed to provide a unique bandwidth for any data. Lastly, simulation examples suggest the proposed approach outperforms the current state of the art in nonparametric density estimation in terms of accuracy and reliability.

The use of transcriptional profiles to predict adult mosquito age under field conditions

Age is a critical determinant of an adult female mosquito's ability to transmit a range of human pathogens. Despite its central importance, relatively few methods exist with which to accurately determine chronological age of field-caught mosquitoes. This fact is a major constraint on our ability to fully understand the relative importance of vector longevity to disease transmission in different ecological contexts. It also limits our ability to evaluate novel disease control strategies that specifically target mosquito longevity. We report the development of a transcriptional profiling approach to determine age of adult female Aedes aegypti under field conditions. We demonstrate that this approach surpasses current cuticular hydrocarbon methods for both accuracy of predicted age as well as the upper limits at which age can be reliably predicted. The method is based on genes that display age-dependent expression in a range of dipteran insects and, as such, is likely to be broadly applicable to other disease vectors.

Constraints on the spatiotemporal accuracy of interceptive action: effects of target size on hitting a moving object

Results of two experiments are reported that examined how people respond to rectangular targets of different sizes in simple hitting tasks. If a target moves in a straight line and a person is constrained to move along a linear track oriented perpendicular to the targetrsquos motion, then the length of the target along its direction of motion constrains the temporal accuracy and precision required to make the interception. The dimensions of the target perpendicular to its direction of motion place no constraints on performance in such a task. In contrast, if the person is not constrained to move along a straight track, the targetrsquos dimensions may constrain the spatial as well as the temporal accuracy and precision. The experiments reported here examined how people responded to targets of different vertical extent (height): the task was to strike targets that moved along a straight, horizontal path. In experiment 1 participants were constrained to move along a horizontal linear track to strike targets and so target height did not constrain performance. Target height, length and speed were co-varied. Movement time (MT) was unaffected by target height but was systematically affected by length (briefer movements to smaller targets) and speed (briefer movements to faster targets). Peak movement speed (Vmax) was influenced by all three independent variables: participants struck shorter, narrower and faster targets harder. In experiment 2, participants were constrained to move in a vertical plane normal to the targetrsquos direction of motion. In this task target height constrains the spatial accuracy required to contact the target. Three groups of eight participants struck targets of different height but of constant length and speed, hence constant temporal accuracy demand (different for each group, one group struck stationary targets = no temporal accuracy demand). On average, participants showed little or no systematic response to changes in spatial accuracy demand on any dependent measure (MT, Vmax, spatial variable error). The results are interpreted in relation to previous results on movements aimed at stationary targets in the absence of visual feedback.

Experimental study of the quantum driven pendulum and its classical analog in atom optics

We present experimental results for the dynamics of cold atoms in a far detuned amplitude-modulated optical standing wave. Phase-space resonances constitute distinct peaks in the atomic momentum distribution containing up to 65% of all atoms resulting from a mixed quantum chaotic phase space. We characterize the atomic behavior in classical and quantum regimes and we present the applicable quantum and classical theory, which we have developed and refined. We show experimental proof that the size and the position of the resonances in phase space can be controlled by varying several parameters, such as the modulation frequency, the scaled well depth, the modulation amplitude, and the scaled Planck’s constant of the system. We have found a surprising stability against amplitude noise. We present methods to accurately control the momentum of an ensemble of atoms using these phase-space resonances which could be used for efficient phase-space state preparation.

Magnetic susceptibility of the normal-superconducting transition in high-purity single-crystal α-uranium

We report complex ac magnetic susceptibility measurements of a superconducting transition in very high-quality single-crystal alpha-uranium using microfabricated coplanar magnetometers. We identify an onset of superconductivity at Tapproximate to0.7 K in both the real and imaginary components of the susceptibility which is confirmed by resistivity data. A superconducting volume fraction argument, based on a comparison with a calibration YBa2Cu3O7-delta sample, indicates that superconductivity in these samples may be filamentary. Our data also demonstrate the sensitivity of the coplanar micro-magnetometers, which are ideally suited to measurements in pulsed magnetic fields exceeding 100 T.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR

Analytical expressions are developed for the time-dependent reactant concentration and catalyst activity in an isothermal CSTR with Langmuir-Hinshelwood kinetics of deactivation and reaction. Several parallel and series posioning mechanisms are considered for a reactor which, without poisoning, would operate at a unique steady state. The use of matched asymptotic expansions and abandonment of the usual initial-steady-state assumption give results, valid from startup to final loss of activity, whose accuracy can be improved systematically.

New calibration technique for multiple-component stress wave force balances

Force measurement in hypervelocity expansion tubes is not possible using conventional techniques. The stress wave force balance technique can be applied in expansion tubes to measure forces despite the short test times involved. This paper presents a new calibration technique for multiple-component stress wave force balances where an impulse response created using a load distribution is required and no orthogonal surfaces on the model exist.. This new technique relies on the tensorial superposition of single-component impulse responses analogous to the vectorial superposition of the calibration loads. The example presented here is that of a scale model of the Mars Pathfinder, but the technique is applicable to any geometry and may be useful for cases where orthogonal loads cannot be applied.

Examination into the accuracy of exchangeable cation measurement in saline soils

Despite the increasing prevalence of salinity world-wide, the measurement of exchangeable cation concentrations in saline soils remains problematic. Two soil types (Mollisol and Vertisol) were equilibrated with a range of sodium adsorption ratio (SAR) solutions at various ionic strengths. The concentrations of exchangeable cations were then determined using several different types of methods, and the measured exchangeable cation concentrations compared to reference values. At low ionic strength (low salinity), the concentration of exchangeable cations can be accurately estimated from the total soil extractable cations. In saline soils, however, the presence of soluble salts in the soil solution precludes the use of this method. Leaching of the soil with a pre-wash solution (such as alcohol) was found to effectively remove the soluble salts from the soil, thus allowing the accurate measurement of the effective cation exchange capacity (ECEC). However, the dilution associated with this pre-washing increased the exchangeable Ca concentrations while simultaneously decreasing exchangeable Na. In contrast, when calculated as the difference between the total extractable cations and the soil solution cations, good correlations were found between the calculated exchangeable cation concentrations and the reference values for both Na (Mollisol: y=0.873x and Vertisol: y=0.960x) and Ca (Mollisol: y=0.901x and Vertisol: y=1.05x). Therefore, for soils with a soil solution ionic strength greater than 50 mM (electrical conductivity of 4 dS/m) (in which exchangeable cation concentrations are overestimated by the assumption they can be estimated as the total extractable cations), concentrations can be calculated as the difference between total extractable cations and soluble cations.

Bifurcation in growth patterns for arrays of parallel Griffith, edge and sliding cracks

This paper presents the recent finding by Muhlhaus et al [1] that bifurcation of crack growth patterns exists for arrays of two-dimensional cracks. This bifurcation is a result of the nonlinear effect due to crack interaction, which is, in the present analysis, approximated by the dipole asymptotic or pseudo-traction method. The nonlinear parameter for the problem is the crack length/ spacing ratio lambda = a/h. For parallel and edge crack arrays under far field tension, uniform crack growth patterns (all cracks having same size) yield to nonuniform crack growth patterns (i.e. bifurcation) if lambda is larger than a critical value lambda(cr) (note that such bifurcation is not found for collinear crack arrays). For parallel and edge crack arrays respectively, the value of lambda(cr) decreases monotonically from (2/9)(1/2) and (2/15.096)(1/2) for arrays of 2 cracks, to (2/3)(1/2)/pi and (2/5.032)(1/2)/pi for infinite arrays of cracks. The critical parameter lambda(cr) is calculated numerically for arrays of up to 100 cracks, whilst discrete Fourier transform is used to obtain the exact solution of lambda(cr) for infinite crack arrays. For geomaterials, bifurcation can also occurs when array of sliding cracks are under compression.