Sensitivity analysis and parameter estimation for distributed hydrological modeling: potential of variational methods

Variational methods are widely used for the analysis and control of computationally intensive spatially distributed systems. In particular, the adjoint state method enables a very efficient calculation of the derivatives of an objective function (response function to be analysed or cost function to be optimised) with respect to model inputs. In this contribution, it is shown that the potential of variational methods for distributed catchment scale hydrology should be considered. A distributed flash flood model, coupling kinematic wave overland flow and Green Ampt infiltration, is applied to a small catchment of the Thoré basin and used as a relatively simple (synthetic observations) but didactic application case. It is shown that forward and adjoint sensitivity analysis provide a local but extensive insight on the relation between the assigned model parameters and the simulated hydrological response. Spatially distributed parameter sensitivities can be obtained for a very modest calculation effort (~6 times the computing time of a single model run) and the singular value decomposition (SVD) of the Jacobian matrix provides an interesting perspective for the analysis of the rainfall-runoff relation. For the estimation of model parameters, adjoint-based derivatives were found exceedingly efficient in driving a bound-constrained quasi-Newton algorithm. The reference parameter set is retrieved independently from the optimization initial condition when the very common dimension reduction strategy (i.e. scalar multipliers) is adopted. Furthermore, the sensitivity analysis results suggest that most of the variability in this high-dimensional parameter space can be captured with a few orthogonal directions. A parametrization based on the SVD leading singular vectors was found very promising but should be combined with another regularization strategy in order to prevent overfitting

Exploring parameter sensitivities of the land surface using a locally coupled land-atmosphere model

This paper presents a multicriteria analysis that explores the sensitivity of the land surface to changes in both land and atmospheric parameters, in terms of reproducing surface heat fluxes and ground temperature; for the land parameters, offline sensitivity analyses were also conducted for comparison to infer the influence of land-atmosphere interactions. A simple "one-at-a-time" sensitivity analysis was conducted first to filter out some insensitive parameters, followed by a multicriteria sensitivity analysis using the multiobjective generalized sensitivity analysis algorithm. The models used were the locally coupled National Center for Atmospheric Research (NCAR) single-column community climate model and the offline NCAR land surface model, driven and evaluated by a summer intensive operational periods (IOP) data set from the southern Great Plains. As expected, the results show that land-atmosphere interactions (with or without land-atmosphere parameter interactions) can have significant influences on the sensitivity of the land surface to changes in the land parameters, and the single-criterion sensitivities can be significantly different from the multicriteria sensitivity. These findings are mostly model and data independent and can be generally useful, regardless of the model/data dependence of the sensitivities of individual parameters. The exceptionally high sensitivities of the selected atmospheric parameters in a multicriteria sense (and in particular for latent heat) appeal for adequate attention to the specification of effective values of these parameters in an atmospheric model. Overall, this study proposes an effective framework of multicriteria sensitivity analysis beneficial to future studies in the development and parameter estimation of other complex (offline or coupled) land surface models. Copyright 2004 by the American Geophysical Union

Searching for sgluons in multitop events at a center-of-mass energy of 8 TeV

Large classes of new physics theories predict the existence of new scalar states, commonly dubbed sgluons, lying in the adjoint representation of the QCD gauge group. Since these new fields are expected to decay into colored Standard Model particles, and in particular into one or two top quarks, these theories predict a possible enhancement of the hadroproduction rate associated with multitop final states. We therefore investigate multitop events produced at the Large Hadron Collider, running at a center-of-mass energy of 8 TeV, and employ those events to probe the possible existence of color adjoint scalar particles. We first construct a simplified effective field theory motivated by R-symmetric supersymmetric models where sgluon fields decay dominantly into top quarks. We then use this model to analyze the sensitivity of the Large Hadron Collider in both a multilepton plus jets and a single lepton plus jets channel. After having based our event selection strategy on the possible presence of two, three and four top quarks in the final state, we find that sgluon-induced new physics contributions to multitop cross sections as low as 10-100 fb can be excluded at the 95% confidence level, assuming an integrated luminosity of 20 inverse fb. Equivalently, sgluon masses of about 500-700 GeV can be reached for several classes of benchmark scenarios.Comment: 26 pages; 8 figures, 6 tables; version accepted by JHE

Measurement network design including traveltime determinations to minimize model prediction uncertainty

Traveltime determinations have found increasing application in the characterization of groundwater systems. No algorithms are available, however, to optimally design sampling strategies including this information type. We propose a first-order methodology to include groundwater age or tracer arrival time determinations in measurement network design and apply the methodology in an illustrative example in which the network design is directed at contaminant breakthrough uncertainty minimization. We calculate linearized covariances between potential measurements and the goal variables of which we want to reduce the uncertainty: the groundwater age at the control plane and the breakthrough locations of the contaminant. We assume the traveltime to be lognormally distributed and therefore logtransform the age determinations in compliance with the adopted Bayesian framework. Accordingly, we derive expressions for the linearized covariances between the transformed age determinations and the parameters and states. In our synthetic numerical example, the derived expressions are shown to provide good first-order predictions of the variance of the natural logarithm of groundwater age if the variance of the natural logarithm of the conductivity is less than 3.0. The calculated covariances can be used to predict the posterior breakthrough variance belonging to a candidate network before samples are taken. A Genetic Algorithm is used to efficiently search, among all candidate networks, for a near-optimal one. We show that, in our numerical example, an age estimation network outperforms (in terms of breakthrough uncertainty reduction) equally sized head measurement networks and conductivity measurement networks even if the age estimations are highly uncertain