“…Since data collection methods to support parameter estimation in two zone transient storage modeling are evolving (e.g., Briggs et al 2009;Neilson et al, 2010a,b), the need for flexibility when incorporating dynamic external information is underscored in model calibration particularly when dealing with both local and global scales. This type of flexibility is not available when optimization algorithms rely solely on the options encoded to solve the problem, which is the case for most single objective algorithms (e.g., nonlinear gradient-based search algorithms such as the Levenberg-Marquardt algorithm (Marquardt, 1963, used by Hil, 1998Doherty, 2005;Poeter et al, 2005), evolutionary algorithms (Duan et al, 1992;Deb, 2001) or Bayesian approaches (Metropolis et al, 1953;Hastings, 1970;Doherty, 2003).…”