Department of Hydrology and Water Resources, University of Arizona, Tucson, Arizona ...... Walnut Gulch Experimental Watershed, Tombstone, Arizona,.
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 106, NO. D24, PAGES 33,421-33,433,DECEMBER 27, 2001
Multiobjective calibration and sensitivityof a distributed land surface water and energy balance model Paul R. Houser,1HoshinV. Gupta,andW. JamesShuttleworth Departmentof HydrologyandWaterResources, Universityof Arizona,Tucson,Arizona
JamesS. Famiglietti Departmentof Earth SystemScience,Universityof California,Irvine, California
Abstract. The feasibilityof usingspatiallydistributedinformationto improvethe predictive ability of a spatiallydistributedlandsurfacewater andenergybalancemodel(LSM) was exploredat the U.S. Departmentof AgricultureAgriculturalResearchService(USDA-ARS) WalnutGulchExperimentalWatershedin southeastern Arizona. The inclusionof spatially variablesoil andvegetationinformationproducedunrealisticsimulationsthat were inconsistent with observations, whichwaslikely an artifactof bothdiscretelyassigninga singlesetof parameters to a givenareaandinadequateknowledgeof spatiallyvarying parametervalues. Becausesomeof the modelparameterswere not measuredor are abstract quantitiesa multiobjectiveleastsquaresstrategywas usedto find catchmentaveraged parametervaluesthat minimizethe predictionerrorof latentheatflux, soil heatflux, and surfacesoil moisture. This resultedin a substantialimprovementin the model'sspatially distributedperformanceandyieldedvaluableinsightsinto the interactionandoptimal selectionof modelparameters. 1. Introduction
Many aspectsof LSM calibrationhave been exploredin other investigations. For example, Sellers et al. [1989] employedmanualcalibrationof nineparameters in the simple biosphere(SiB) modelto improveits comparison with field observations. Franks and Beven [1997] use a generalized likelihooduncertainty estimationapproach to reducetheoutput uncertaintyfor the TOPUP model using First International SatelliteLand SurfaceClimatologyProjectField Experiment (FIFE) and Anglo-Brazilian Climate Observation Study (ABRACOS) observations.Lettenmaieret al. [1996] reported that in an intercomparisonof various LSMs, those whose
A land surfacewater and energy-balancemodel (LSM) [Famigliettiand Wood, 1994] was used to simulatethe spatiallydistributedbehaviorof the U.S. Departmentof AgricultureAgricultural ResearchService (USDA-ARS) WalnutGulchExperimental Watershed(WGEW). This type of modelgenerallyhas a large numberof parameters that describesurfacevegetationand soilsthat mustbe correctly specifiedto produceaccuratewater and energy balance simulations.Many of theseparameters arephysicallyrealistic and observable;that is, they are dimensionsor capacities parameterswere calibratedperformedbetterthanthosewithout which can be measuredreliably. For these parameters, parametercalibration.For a reviewof conceptualhydrologic observedvaluescan be usedwhen available. Other parameters modelcalibrationto observeddata,seeGuptaet al. [ 1998]. are not observableor are not physicallyrealistic. Such Automatic spatially distributedLSM calibrationhas been parameters maywell be conceptual representations of abstract severelylimited by the absenceof adequatespatiallyvariable watershedcharacteristics or empiricalconstantsthat can be hydrologiccalibrationdata. Automaticcalibrationmethods usedto optimizethe simulationresultson a trial-and-error adjustLSM parameters to obtaina bestfit of LSM predictions basis to match the simulations to observations. with observations. Therefore,to automatically obtainspatially Two aspectsadd complexityto the problemof parameter distributedLSM parameters,spatiallydistributedobservations specificationin this study. First, becausethe surface of LSM predictions,suchas soil moisture,surfaceheatfluxes, conditions are spatiallyvariable,the parameters mightalsobe and runoff, must be available. Currently, distributed expectedto vary spatially. Se•nd, because the performance observationsof LSM predictions are scarce, and when of the model is to be judged in terms of its ability to available, have relatively high uncertainty [Franks et al., simultaneously predictthe soil moisturestatesand the latent, 1998]. Therefore spatially distributedLSM calibrationhas sensible, andsoilheatfluxes,the specification of parameters is generallybeenlimitedto (1) specification of spatiallyvariable a multiobjective calibrationandevaluationprocess. parameters from soil, vegetation,or topographic mapsderived from field surveysor remotesensingor (2) lumpedor area1NowatHydrological Sciences Branch, NASA/GSFC, Greenbelt,averagedcalibration,where parametersare assumedto be Maryland.
Copyright 2001bytheAmerican Geophysical Union. Papernumber2000JD900803. 0148-0227/01/2 000JD900803 $09.00
constant over an area and are adjusted coincidentallyto minimize error in streamflow[e.g. Storcket al. 1998]. An exceptionis a studyby Frankset al. [1998], wheredistributed model parameterizations were first conditionedon lumped discharges and then further conditioned on distributed 33,421
33,422
HOUSER ET AL.: MULTIOBJECTIVE
CALIBRATION
estimatesof saturatedareasestimatedfrom remotesensingand energy fluxesatthelandsurface andlocalverticalrecharge to topographyusingthe weightingof estimatesin a fuzzy set the watertable. It incorporates simplerepresentations of framework.
This paperconsistsof six sections.Section2 describesthe studysite andthe dataset. Section3 providesan overviewof the LSM. The fourthsectionexploresthe hypothesis that the useof spatiallydistributedsoilsandvegetationinformationcan help to improvethe predictiveability of the LSM. Section5 explores the hypothesisthat a multiobjectiveparameter calibrationmethodcan be used to improvethe estimatesof nonobserved modelparameters.The last sectionpresentsthe
atmosphericforcing, vertical soil moisturetransport, plant-controlledtranspiration,interception,evaporation, infiltration,surfacerunoff, and sensibleand groundheat fluxes. The LSM incorporates a diurnalcycleandis driven withstandard meteorological datawith an hourlytimestep, thisbeingconsidered sufficient to resolve thedynamics of the land-atmosphere interaction.
conclusions and some final remarks.
4. Use of Spatially Distributed Soilsand Vegetation Information
2. Walnut Gulch and Monsoon '90 Observations
This sectionexaminesthe hypothesisthat the use of spatiallydistributed soilsandvegetation information canhelp
The WalnutGulchExperimental Watershed (WGEW) is to improvethe predictiveability of the LSM. The LSM was thelandsurface waterandenergydynamics in operated by the Southwest Watershed Research Center usedto simulate (SWRC) of the Agricultural Research Service, U.S. a spatiallydistributedmannerfor the entire WGEW at a 40 m
resolution. Predictions were madehourly from July 22 to August15, or day of year (DOY) 204 to 228, of 1990. LSM parameter estimates were first specified by using catchment-averaged parametervalues based largely on groundwater variesfrom 45 m at the lower endto 145 m in the observations. Next, spatially distributedvalues for the centerof thewatershed.Soiltypesrangefromclaysandsilts parameterswere derived from geographic information to well-cemented boulderconglomerates [Kustas et al., 1992], system's(GIS) maps. withthe surfacesoiltextures beinggravellyandsandyloams containing, on average,30% rock and little organicmatter. 4.1. Spatially ConstantParameters Themixedgrass-brush rangeland vegetation ranges from20 to Spatiallyconstantparametersfor the LSM were specified 60% in coverage.Thisrangeland regionreceives250 to 500 primarily on the basisof observations made duringMonsoon mm of precipitationannually,typicallytwo thirdsof it as '90 [Daughtryet al., 1991; Kustaset al., 1992, 1994a, 1994b, convectiveprecipitation duringthe summermonsoonseason 1994c; Stannard et al., 1994; Moran et al., 1994a, 1994b;
Department of Agriculture.Thestudycatchment is a heavily instrumented areacomprising theupper148km: of theWalnut Gulchdrainagebasinin an alluvialfan portionof the San PedroRiver watershed in southeastern Arizona. Depthto
[Renardet al., 1993].
Humes et al., 1994a, 1994b; Weltz et al., 1994; Menenti and
Eighty-fiverecordingrain gauges,11 primarywatershed Ritchie,1994;Hippset al., 1994;Spangler,1969]. The albedo runoff-measuring flumes, and micrometeorologicalof dry soilwas specifiedby adoptingthe parameters suggested
observations make the WGEW
a valuable research location.
by Dicla'nsonet al. [1993]. Soil profile temperatures measured
Duringthe Monsoon '90 experiment (July23 throughAugust in three trenchesat the Lucky Hills and Kendall siteswere 10, 1990),extensiveremote-sensing observations weremade, analyzedand modeledto determinethe temperaturediurnal while eight micrometeorological energy flux (Metflux) dampingdepth. The saturationsoil moisturewas assumedto
instruments provided continuous measurementof local be close to the porosity calculated from bulk density
meteorological conditionsand the surfaceenergybalance measurements.The surfacesaturatedhydraulicconductivity [Kustasand Goodrich,1994].
wasbasedon porosityusingthe Kozeny-Carmen equationwith Precipitation is the mostimportantspatialforcingvariable a coarsefraction correction,following RawIs et al. [1993]. in semiaridregionsbecause of its highlyvariable,convective The depth of the surfacezone was specifiedequal to the nature;thusmucheffort was devotedto derivingspatially maximum reported L band microwave penetrationdepth distributedprecipitationdata sets for the Monsoon '90 [Jackson,1993] to aid compatibilitywith remotesensingdata. experiment. A multiquadric-biharmonic interpolation Becauseof the deep water table at Walnut Gulch, surface algorithm[Syed, 1994] was used to producespatially processesshow no sensitivity to the LSM water table distributedprecipitationvaluesfor the entire model domain parameters,so thesewere assignedto arbitraryvalues. The fromthe availableraingaugedata. All othermeteorologicalvegetationstresssoil moistureparameterwas assumedto be forcingwasassumed to be spatiallyconstant andderivedfrom halfway between the saturatedand residual soil moisture averaging observations at the eightMetfluxstationsin place values,while the wiltingpointwasassumed to be equalto the during theexperiment. A spatial averaged soilmoisture profile valueof residualsoilmoisturebecausedesertvegetationrarely was derived from several in situ profile observationswilts. The soil moistureat which evapotranspiration reaches [Schmugge et al., 1994]andusedto initializetheLSM onJuly thepotentialratewasassumed to be 10% belowsaturation, this 22 (Table 1). beingthe approximate ratio suggested by Sellerset al. [ 1986]. Values for albedoof dry vegetationfor semiaridareaswere obtainedfrom Dicla'nsonet al. [1993], and the albedoof wet vegetationwas assumedto be 5% lower than that for dry The land surfacewater and energybalancemodel(LSM) vegetation.The stomatalresistance, root activity factor,root [Famiglietti and Wood, 1994] simulatesland surfacerunoff, density,root resistivity,and criticalleaf water potentialwere energyfluxes,andsoilmoisturedynamicsin threelayers.This not measured; indeed,becausemanyof theseparameters were LSM is designedto predictdiurnaldynamicsof the water and notobservable quantities, estimates wereused. 3. Land
Surface
Model
HOUSER ET AL.: MULTIOBJECTIVE CALIBRATION
33,423
Table1. Catchment-Averaged LSMParameters forWalnut Gulch, Based onObservations, GISCoverage, and MultiobjectiveCalibration a Parameter
Source
Vegetationparameters vegetationheight,m leaf area index
minimumstomatal resistance, s m-• initial waterstoragein canopy,m unstressed transpiration soilmoisture,% wiltingpointsoilmoisture,% vegetationfraction albedo,wet vegetation albedo,dry vegetation albedo,bare soil
root activityfactor
rootdensity, m m'3 rootresistivity, s m4 criticalleaf waterpotential,m Soil parameters surfacezonedepth,m initial surfacesoil moisture, %
root depth,m initial root zone soil moisture, %
maximum rateof capillary rise,m s'• initial transmissionzone soil moisture, %
percentsand,% percentclay, %
bulkdensity, gcm'3 residual soil moisture, % saturated soil moisture, %
saturated hydraulic conductivity, m s'• baresoilroughness length,m baresoil zeroplanedisplacement, m Water tableparameters averagetopographic index K sexponentialdecayparameter initial watertabledepth,m Energybalanceparameters soil moisturefor calculationof PET, %
diurnalheatpenetration,m temperature of deepsoillayer,øK
Observed
GIS
Calibrated
Humeset al. [ 1994a] Daughtryet al. [1991] Chowet al. [ 1988] assumed dry at startof simulation
(0.22) 1.275 40 (0)
(1.31) 267
(574)
assumed halfwaybetween0sand0r assumed 0r (cactirarelywilt) Kustaset al. [1994a] Dickinsonet al. [ 1993] Dickinsonet al. [1993] Dickinsonet al. [1993] Famiglietti[ 1992] Famiglietti[1992] Fatniglietti[ 1992] Fatniglietti[ 1992]
(20) (1) (0.42) (0.2) (0.25) (0.33) 10,000 1 1e9 -210
Basedon PBMR sensitivitydepth Schmugge et al. [ 1994] Schrnugge et al. [1994] Schrnugge et al. [ 1994] Famiglietti[1992] Schmugge et al. [ 1994] Kustasand Goodrich[1994] Kustasand Goodrich[ 1994] Kustasand Goodrich[ 1994] Rawlset al. [1993] Kustasand Goodrich[1994] Rawlset al., [ 1993] assumption assumption
(0.1) (10.0) 1.5 (17.0) (0.1) (17.0) 70.9 8.5 (1.6) 1 38.0 6.9e-5 (0.001) (0.0)
calculationfromA(tan(B))image Famiglietti[ 1992] Gilbert[ 1996]
(8.314) (7.0) (100.0)
Sellerset al. [1986] from observations [Houser,1996] fromobservations [Houser,1996]
28
0.5 297.0
(348025) (86.5) (4.8e 11) (-500)
(0.8)
63.1 22.1
(14.4) (8.•)
5.1 33.6 4.7 e-6
(2.2) (30) (8.7e-6)
(47) (0.33) (287.6)
Values inparentheses produce the"best" Monsoon '90 simulation. Readl e9as1fro.
4.2. Spatially Variable Parameters
for eachsoil seriesin the soil survey. Bulk density,B, was
estimated usingan empiricalformula[Rawlset al., 1993]: Spatialdistributions of the soil andvegetation parameters wereestimated usingGIS mapsof WalnutGulchvegetation B = 1.51+0.0025(S)- 0.0013(S)(O)- 0.0006(C)(O) andsoils. The potentialrootingdepthsfor eachsoil series were taken from a recent Walnut Gulch soil survey (D. J.
Breckenfeld, unpublished document,1993). Valuesof the effectiveporosity,•be, residualsoil moisture,0r, and the saturatedhydraulicconductivity, Ks, were determinedby Mayeux[1995]usingvaluessuggested by Rawlset al. [1993] andBouwer[1966]. J. Berglund,unpublished report,(1995) calculatedsaturatedsoil moisture,0e, for each soil classas follows:
Oe--•e-{-Or.
- 0.0048(C)(E).
(2)
whereS is percent sand,C is percent clay,O is percentorganic matter, and E is the cation exchangecapacity of clay normalizedby percentclay. Measurements at the eight Monsoon'90 Metflux sites [Daughtryet al., 1991] and
parameter rangesgiven by Jones[1983] were usedin conjunction withtheWalnutGulchvegetation GIS coverage to estimate the spatialdistribution of leaf areaindex(LAI) and
minimum stomatalresistance. The topographicindex was fromdigitalelevation models(DEMs)derivedfrom Percentsandandclay werecalculated by J. Berglund(1995) calculated 23 1:5000 stereo aerial photos [Matthews,1992]. DEM fromvaluesof percentmaterialpassing throughcertainsieves
33,424
HOUSER ET AL.' MULTIOBJECTIVE CALIBRATION a) Topographic Index
b)1August 1990 Rainfall
ß '••
20
.•
..?
'•'• '•'•' •--' '" Omm
c) Leaf Area Index
d) Saturated Soil Moisture ...,--
.
.::.. :
:•::
o
0%
e) 7 August 1990 LSM Surface Soil Moisture(%)
f) 7 August 1990 LSM Surface Soil Moisture(%)
using Spatlally Constant Parameters ...•.• .•..• using Spatially Variable Parameters 't
• ....
g) 5 August 1990 PBMR Soil Moisture
½;½• '• •:•5•?"• .....
•:½:::,½;•...;•t•½ •
•"
h) 9 August 1990 PBMR Soil Moisture
4O
Scale(km •'• • •Scale(kr•• • 0%
o
5
0%
Figure 1. SpatiallydistributedWalnutGulchExperimentalWatershed:(a) topography,(b) precipitation, (c) vegetation, and(d) soils. LSM spatialpredictions of surfacesoilmoistureat 1200LT on August7 (DOY 219) using(e) spatiallyconstantand (f) spatiallyvariablesoilsand vegetationparameters (all simulations using spatiallyvariabletopographyand precipitation). Pushbroommicrowaveradiometer(PBMR)'derived soil moisturefor (g) August5, 1990 (DOY 217) and (h) August9, 1990 (DOY 221) [Houseret al., 1998]. The additionof spatiallyvariablesoilsandvegetation produces unrealisticpolygonartifactsin the simulation.
overlapswere averaged,then a smoothingalgorithm was appliedto reducenoise,and isolatedrunoffsinkswere filled. The watershedwas delineatedand the topographicindex was derived using the standardmethodologydetailedby Beven [ 1995]. A selectionof representative spatiallyvariableWalnut Gulchparameters is shownin Figure1a-1d, andthe watershed averages of theseparameters arereportedin Table 1. 4.3. Results and Discussion
1e'1fi It isevident thatthespatially variable parameters havea largeimpacton the spatialpatternsof the simulation.Further, the soil moisturepatternsderivedusingthe spatiallyvariable parametersappearto be unrealisticand do not comparewell with observedpush broom microwaveradiometer(PBMR) surfacesoil moisture(Figures1g- 1h). A seriesof simulations, whereonly one spatiallyvariable parametersetwasusedat a time (theotherparameters were set to theirestimated catchment-averaged values),wereperformed to determinewhich subsetof spatialparameterscontribute
Thesimulated spatial patterns of surface soilmoisture for mostto thesepatterns(the topographic index and the August 7, 1990(DOY219),using boththespatially constantprecipitation fieldswere spatiallyvariablein all these andthespatially variable parameters, areshown in Figuressimulations).The simulations with spatiallyvariable
HOUSER ET AL.: MULTIOBJECTIVE
vegetationparametersperform similarly to the control run, indicatingthat spatialvariationin theseparametershas little effect on predictions. All of the simulationsusing spatially variablesoil parametersshowdistinctpolygonpatterns. The parameterspecifyingsaturatedsoil moisturehas the most influence on simulatedspatial patterns,while those which specify the percentagesof sand and clay, the saturated hydraulicconductivity,and the residualsoil moisturehave a more
moderate
influence.
The
soil
characteristics
are
extremelyimportant,and the vegetationcharacteristics are less important in regulating soil moisture variability and other surfacewaterand energycomponents in this watershed.The resultsof this sensitivityalso showthat the spatiallyvariable topographicindex hasvery little influenceon the simulations; thisresultis a consequence of the largedepthto the watertable
CALIBRATION
33,425
Gupta et al. [1999], this situationinvolvesthe problemof finding parameterestimateswhich can simultaneously minimize several noncommensurable criteria.
This section
examinesthe hypothesisthat a multiobjectiveparameter calibrationmethodcanbe usedto improvethe estimates of the nonobserved modelparameters of the LSM model. 5.1. Data Available
for Calibration
The primary observationsused to calibrate the LSM parametersin this study were hourly time series of catchment-averaged, water balancecorrected,latent heat flux during unstableatmosphericconditions,soil heat flux, and gravimetric calibrated near-surfacesoil moisture measured
with resistance sensors (Table 1). Duringthe Monsoon'90 study period, three replicate gravimetricsurfacesoil moisture at Walnut Gulch. samples were collected each day at the eight Metflux sites The pattern of enhanced spatial soil and vegetation [Schrnuggeet al., 1994]. Resistancesensorscollected polygonsapparentin the simulationsis an artifact of both continuous time series of soil moisture at 2.5 cm and 5 cm assigninglumped values of the parametersto discretely below thesurfaceat all eightMetflux sites[Arneret al., 1994]. partitionedareasas well as the misspecification of parameter Becausesuchsensorsare generallydifficult to calibrateand values. A moreappropriatespecification of spatialparameters would reflecttheir continuous variabilityin space,as obtained tend to drift, they were recalibratedeach day against measurements for the purposeof this study. The with remotesensing. It is also clearthat the specificationof gravimetric
energy balance was determined from measurementsof net the spatialparametervaluesusingscatteredobservations, table radiation, soil heat flux, and estimatesof sensibleand latent lookup,andphysicalrelationships (as describedin section4.2) variance result in inadequateparameterspecification. These results heat flux by either eddy covarianceor temperature methods [Kustas et al., 1994a]. A detailed water balance study indicatethat the available WGEW parameterobservations simply have insufficientaccuracyand spatial frequencyto showedthat during Monsoon'90 the sensibleheat flux was becauseof eddycorrelation propellerstalls, realistically constrainthe LSM. Becausethe simulations underestimated source area mismatch, or horizontal flux divergence in hilly using spatially constant vegetation and soil parameters terrain [Houser et al., 1997; William& 1996; Keefer et al., comparewell to the PBMR patterns,the subsequent parameter calibration studies reported in section 5 assume spatially 1997]. Therefore,for this study the observedlatent and with the constantsoil and vegetationparametersacrossthe catchment, sensibleheat fluxeswere adjustedto be consistent water balance, and unreliable observations during stable leaving only topographicindex and precipitationas spatially varyingentities. The use of catchment-averaged parametersin this study is Table2: Calibration Parameter Ranges andObjective Functions a supported by Whiteet al. [ 1997],who foundthatthe difference betweenthe calculatedarea-averagedsurfaceenergy fluxes Calibration Parameter Possible Range givenfor the WGEW by a singlepointlandsurfacemodelwith 1 to 1e6 Root activityfactor aggregateparametersand that given by a distributedarray of 1 to 100 Rootdensity,m m-3 land surface models was small.
This conclusion
was tested
here; althoughthere was more variability than reportedby Whiteet al. [1997] becauseof the finer spatialresolution,there was overall agreement. The relevanceof this finding to the presentstudyis that it supportsthe useof catchment-averaged parametersand forcingfor the WGEW if catchment-averaged surfacefluxes are required. This conclusionis only valid when catchment-averaged fluxes are of interest; for cases wheredownslopeflows or water availabilityin valley bottoms is of interest,thenspatiallyvariableparameters and forcingare required. In this study, the spatially distributedLSM is requiredto allow predictionof spatialpatternsof soil moisture at the resolutionof the availableremotelysensedobservations.
Rootresistivity,s m4 Criticalleafwaterpotential,m Potentialevaporation soil moisture,% Penetration of diurnalheating,m Temperature of deepsoillayer,K Minimumstomatalresistance, s m'l
le6 to le12 - 1 to - 1000 10 to 50
0.1 toO.7
285 to 310 0 to 700
Percentsand,%
lto
Percentclay, % Residualsoil moisture,% Saturatedsoil moisture,%
1 to 100
100
1 to20
20 to 50
Surface saturated hydraulic conductivity, m s4 5e-5to 5e-7 ObjectiveFunction
5. Multiobjective Calibration of Nonobserved Model
Parameters
The purposeof this study was to obtain estimatesof the nonobserved modelparametersthat would result in improved simulations of the overall behavioral responsesof the watershed.In particular,the desirewasto selectvaluesfor the parameters that minimizedthe error in predictedsoil moisture and sensible,latent,and groundheat fluxes. As discussed by
sensible heatflux(eightsites, waterbalance corrected, daytime average fromtemperature variance)
latentheatflux(eightsites, waterbalance corrected, daytime averagefromtemperature variance)
soilheatflux(eightsites average fromsoilheatfluxplates) surface soilmoisture (eightsites, gravimetric corrected average fromresistance sensors) aThe sensible heatflux objectivefunctionwasincludedasa diagnostic;it was not used in the determinationof the Pareto set. Read 1e6 as 10'6.
33,426
HOUSERET AL.: MULTIOBJECTIVECALIBRATION
atmospheric conditions wereexcluded fromobjective functionFigure2a illustrates howeachof the meansquare error calculations.
functions f•(0) andf2(0) mightvarywith 19.Notethat the parameter value19= 02minimizes the function f• buthasa relatively poorvalueforf2, whilea differentvalue19= 02 minimizes thefunction f2 buthasa relatively poorvalueforf•.
5.2. Parameters Calibrated
The 13 parameters listedin Table2 wereidentified for In otherwords,the individualsolution19= 0• that givesthe calibration.Table2 alsoliststhetypicalrangeof theirvalues; best match to the first model flux might give rather poor
theserangeswerespecified largerthanmightnormallybe performance in termsof matching thesecond modelfluxand expected in orderto allow the calibration procedure some vice versa.However,in general,we wish to havea model leewayin accounting for modeland dataerror. Moreover, solutionthatsimultaneously givesacceptable "good"matching becausesome of the parametersare nonphysical, the of both fluxes. specification of theirranges is somewhat arbitrary.Notethat Tothisendwenotethatthereexists a region ofvaluesOv= thesoilheatfluxparameters (temperature of thedeepsoillayer { 01< 0< 02} consisting of thoseparameter values forwhichit andpenetration of diumalheating) wereincluded in thislist is possible to obtainimprovement in f2(0 ) (by varying0 ) eventhoughtheyweremeasured with confidence.This is whilesimultaneously allowingsomedeterioration inf2(19);this because the soil heatflux was poorlysimulated by the LSM regionis calledthe "Pareto,""non-inferior," or "trade-off' owingto the simplicity of its soilheatflux submodel; this solution.Similarly, all valuesof 19q•O•arecalled"inferior" suggests thatnonphysical valuesfortheseparameters maybe solutions. Figure2b shows theproblem plottedin thefunction necessary forthemodelto performadequately. Notealsothat space withf• (19)andf2(0) ontheaxes,andthecurverepresents the five soil water retentionparameters,which can be
the trajectoryobtainedby varying19overits feasiblerange. Clearly,theobjective isto find 0 suchthatwe comeascloseto calibration, because this greatlyimprovedthe simulation of theorigin• ,f2}={0,0}aspossible.Notethatthesegment of soilmoisture andyieldedinsightintotheuseof observed soil the curve betweenB and C represents the trade-offregion, parameters. whileotherportions of the curverepresent inferiorsolutions. Any one of the trade-off solutions can be selected by theuser 5.3. TheoreticalBackgroundto the Multiobjective as "best," depending on which trade-off in matching of model Calibration Approach performance is considered mostacceptable by theuser. A theoreticaland practicalbasisfor the calibration of models withmultipleandnoncommensurable outputfluxeshas 5.4. ObjectiveFunction
estimated from observations,were also included in the
beenpresented by Guptaet al. [1998,1999],Bastidas et al. Threeobjectivefunctions wereusedto calibrate the model [1999],andMeixneret al. [1999]. Because of errorsin the parameters: thehourlymeansquare errorin matching of latent model(structure) andin the datait is impossible to find a heat flux duringunstableatmospheric conditions, soil heat single"best" parameter setthatsimultaneously minimizes the flux, andgravimetric calibrated near-surface soilmoisture.In errorsin matching all of theobservable modeloutputs.This eachcasethe meansquareerror objectivefunction,F, was leadsto a multiobjective optimization problemfor whichthe computed usingthenormalized sumof thesquared difference "solution" is a regionof the parameter spacethat reflects betweenthe catchment-averaged, water balancecorrected differenttrade-offsin the matchingof the outputs. An observations, Oi, and the predictions,Pi, wheneverthe
acceptable solution can beselected bytheuser from this region observation was available; thus of trade-offsolutions usingsubjective judgment.
To clarify,a simpleone-parameter example is presented in
Figure 2;we assume that data are available for two model output fluxes thatmust bematched, andweletf•(O)and f2(0)
F=n1y,• (O•P•) i=•
(3)
represent the meansquareerror(MSE) in matching these fluxesfor selected valuesof a specificmodelparameter 19. wherei isthetimestepandn isthenumber of observations. (b)
(a)
)
f•
f1(02 )
B
f•(O•)
01 '- 02
f2(02 ) f2(01 )
Figure2. A simpleone-parameter, two-objective modelcalibration example: (a) howeachof themean square errorfunctions, f•(0) andf2(0), mightvarywithmodel parameter 19,and(b) problem plotted in the function space withthecurverepresenting thetrajectory obtained byvarying 19overitsfeasible range.
HOUSER ET AL.: MULTIOBJECTIVE
5.5.MultiobjectiveCalibration Methodology
CALIBRATION
33,427
shows thecorresponding trade-off among theparameters.
Inspeclion of theseresults shows thatthefive solutions giving
Theprocedure adopted forthemultiobjective calibration of bestmatching to soil moisture alsogive relativelygood the LSM parameters at WalnutGulchwasto conduct a matching of theotherfourfluxes.In addition, it wasdecided Monte-Carlo search of thefeasible parameter space indicatedthatthe soilmoisture objective wasmostimportant for the inTable2. A totalof 200,000 parameter setrealizations were purposes of thisstudy. Therefore the parameter setthat randomly generated, andthethreeobjective function valuesproduced the lowestsoil moisture objective was chosen werecomputed at eachparameter set. Fromthese,the (bolded linein Figure4). Theselected parameter values are "noninferior"trade-off parameterregion was estimatedby shown in Table 1. findingthosepointsfor which theredid not existat leastone It is interesting to notethatthecalibrated parameter values othermemberof the original200,000realizationsfor whichall arereasonably closeto the defaultvalues,with the exception three objective function values were better. Only 125 of the rootparameters, andare in generalmorerealistic. A parameter set realizationswere found to belong to this rootdensity of 86m m'3ismuchmorerealistic than1 m m'3, estimatednoninferiorregion, indicatinga rather substantial for instance, while the high minimumstomatalresistance is reductionin overallparameteruncertainty.In addition,at each realisticfor water-conserving semiaridplants. The default noninferiorpoint,the meansquareerror in matchingof hourly percentsandand clay parameters that were derivedfrom a sensibleheatflux wasalsocomputedfor diagnosticpurposes. numberof soil analyses(J. Berglund,unpublishedreport, The estimatednoninferiorregion is displayedin Figure 3 1995)were likely too highbecause the largecoarsefraction usingthe formatsintroducedby Sorooshianet al. [ 1993] and (D. J. Breckenfeld, unpublished document, 1993)wasremoved
Gupta etal. [1998,1999].Figure 3 (le•) shows thetrade-offpriorto analysis.Basically, if the sand,silt, andclay amongthe objectives, whileFigure3 (right)showsthe percentages (whichin thepublished soilanalysis addupto corresponding trade-off among theparameters. Notethatthe 100%)areadjusted sothattheyaddupto 100%minusthe objectivefunctionvalueshave been normalized onto a coarsefraction, thentheyget verycloseto the parameters zero-to-one scale, where zerocorresponds tothesmallest error found using calibration. These insights leadtotheinclusion of andonecorresponds to thelargest errorin thefeasible region. a coarsefractioncorrection in the modelphysics[Houser, Similarly,the parameter setsare alsonormalized ontoa 1996].Thelargedifferences in soilheatfluxparameters are zero-to-one scale, withzerocorresponding to theminimumprobably anartifact of theLSMhaving a simplified soilheat valueof theparameter andonecorresponding tothemaximumfluxsubmodel. Finally,thehighvalues of criticalleafwater valuein theirdefined ranges.Theseplotsenable easein potential andsoilmoisture for thecalculation of potential identifying whichparameters arepoorly orwellidentified and evapotranspiration (PET)reflect thefactthattheplants donot alsohow variousparameter valuesaffectprediction errors. wilt and that the environment rarelyexperiences potential
Theestimated noninferior solution region shows relatively low evaporation. values formostoftheobjectives, withthelargest spread being Thereis a clearcompromise in the calibration between in the matching of soilmoisture.Someof the calibratedlatentheatfluxprediction andsoilmoisture prediction (Figure (noninferior) parameters ranges arequitetight(e.g.,deepsoil 5). Therange of prediction shows thatlatentheatfluxcanbe temperature andpercent clay),showing thattheseparameters calculated wellbytheLSMbutattheexpense of soilmoisture. areidentified well. Thenoninferior parameter values alsoare By choosing to predict soilmoisture mostaccurately, rather generally confined to moderate ranges of theparameter space. thanlatentheatflux, we obtain"flat tops"or "hats"in the However,quitea lot of parameterinteractionis observed,with simulatedlatentheat flux time series. Accordingto Carlson two clustersof parametersevident(i.e., seeminimumstomatal [ 1991,p. 353],"thereis considerable observational evidence of resistance, potentialevaporationsoil moisture,diurnalheating [such]behaviorin the field." This LSM prediction of latent penetration, residualsoilmoisture,etc.). heatflux may not be entirelyincorrect, but it is certainlynot Despite the relatively low objective function values supported by theMonsoon '90 observations. It is difficultto producedby the noninferiorparameterrangethe corresponding drawanydefinitiveconclusions herebecause the accuracy of total rangeof LSM flux predictionremainshigh. The eight these latent heat flux observationshas been questioned Metflux site averagerangesof predictionobtainedfrom the [Williams,1996;Keeferet al., 1997;Houseret al., 1997]. Paretoparametersetsare shownin Figure 5. There is still a large amountof predictionvariability even within the Pareto 6. Conclusions set,andhencethe final parametersetmustbe chosencarefully. High-quality land surface water and energy balance All of the noninferiorparametersets are good but with different compromises as to which objective is more fully simulationsare essentialfor the predictionof drought,floods, minimized. To selecta single"best"parameterset from the agriculturalproduction,land surfaceinputsto the atmosphere, of landsurfacesto climatechange. However, 125 noninferiorsolutions,a necessarilysubjectiveprocedure andthe response was followed. For each of the three objectives,the five the land surfacewater and energybalancemodelsthat make generallyutilize largenumbersof parameters, noninferior parameter sets giving the best value of that thesepredictions objectivefunctionwere selected(a total of 15 solutions);in manyof which are not regularlymeasuredand someof which addition,the five membersof the noninferiorregiongivingthe describe abstract land surface characteristics that cannot be bestvalueof hourlymeansquareerrorin matchingof sensible measured. The inherent heterogeneityof land surfaces heat flux also were selectedas shownin Figure 4. Note that requires spatially distributed LSM predictions, which the sensibleheat flux objectivefunctionwas not usedin the complicatesthe opli•num selectionof parameters. Further, estimationof the noninferior region. Each row of plots LSMs are expectedto predictmultipleland surfacestatesand corresponds to a differentoptimizedobjectivefunction,while fluxes reasonably well, which requires the use of a Figure4b showsthe trade-offamongobjectivesand Figure4a multiobjective calibration procedure to achieve optimal
33,428
HOUSER ET AL.: MULTIOBJECTIVE
a) Best5 Sensible HeatFluxObjective Parameter Sets
CALIBRATION
1.0
•_1.0 b Best 5Sensible Heat Flux Objective Parameter Sets
0.5
a' 05•
0.0
z 0.0
13
ß
0
1.0 C Best 5Latent Heat Flux Objective Parameter Sets ,_
0.5 [
•
d) Best5 Latent HeatFluxObjective Parameter Sets
1.0
0.5
.__N_ o
z 0.0
o.o
e)Best 5Soil Heat Flux Objective Parameter Sets
1.0
,_ 1.
o.s
,•
..•
O f)Best 5Soil Heat Flux Objective Parameter Sets
o.s
0.0
z 0.0
1.0 'q)Best 5Soil Moisture Objective Parameter Sets
._ 1.0
h) Best 5SoilMoisture Objective Parameter Sets
0.5
0.5
o
z z
Figure 3. Paretosetnormalizedobjectivefunctionandparametervaluesfor the 13-parameter, threeobjective function Monte-Carlo calibration. The best five parametersets and objective function values for each objectivefunctionare shown,with the subjectivelychosen"best"parametersetshownin bold.Parameterand objectivefunctiondefinitionsare shownin Table 2.
performance. Therefore, thisstudyhasexplored thefeasibility distributed LSM should improve localpredictions, thisstudy of usingspatiallydistributed soilsandvegetation information showsthat the polygonnatureof thesedata setsresultsin anda multiobjective parameter calibration methodto improve unrealistic simulations which were inconsistent with the predictiveability of a spatiallydistributedLSM at the observations. A parameter sensitivity studyrevealedthat WGEW. spatialvariations in vegetation parameters havelittle effecton Although it seems reasonable thattheinclusion of spatially predictions, butthatsoilparameters havea largeeffect. This variablesoil and vegetationinformationin a spatially problem is bothanartifactof discretely assigning a singleset
HOUSER ET AL.' MULTIOBJECTIVE
CALIBRATION
33,429
1.0a) 0.8
0.6
0ø4 0.2
0.0
b)
0.8
0.6
0.4
0.2
Figure4. Pareto setnormalized (a) objective function and(b) parameter values forthe13-parameter, three objective function calibration. Parameter andobjective function definitions areshown inTable2. of parametersto large areas of the catchmentand a However, the unrealisticvery strong surface soil spatial misspecification of spatialparameter fieldsdueto inadequatepatternspredictedon August7 (Figure 1) would likely be enhanced givena spatiallyvariablesoilmoistureinitialization measurement of the requiredparameters. or longer termmodelspinup.It is possible thatthe spatially ß It is alsorecognized thatthechoiceof spatiallyconstant soil moistureinitializationon July 22, 1990 and the shortmodel distributedinitial conditionscouldbe manipulatedto producea prediction.However,we had no spinupinfluencesthe spatiallydistributedpredictionsand betterspatiallydistributed influences thus the calibration of the distributed model. spatiallydistributed observations on July 22, for validating
33,430
HOUSER ET AL.: MULTIOBJECTIVE Latent
a)
-100
' 204
• 206
•
I 208
•
I 210
•
I 212
•
I
,
214
CALIBRATION
Heat Flux
I
216
I
I
218
Day of Year
•
I
,
220
,
I
i
224
1
226
228
........ LSMPareto SetRange --- •-.
Sensible
I
222
Observed Calibrated LSM
Heat Flux
2OO
-4OO
204
206
208
210
212
214
216
218
Day of Year
220
222
224
226
228
...................... LSMPareto SetRange Observed :-- ;-' Calibrated LSM
Figure5. Catchment-averaged timeseriesof (a) latent,(b) sensible, (c) soilheatflux, and(d) near-surface soilmoisture andtheirassociated prediction rangesfor the 13-parameter, threeobjective function Paretoset. Thesubjectively chosen calibrated parameter setprediction andarea-averaged observations arealsoindicated.
such manipulations,and we seriouslydoubt that such only eight surfaceflux observationpoints. Given that the manipulation wouldchangeour conclusions giventhe very availablespatial informationis insufficientto accurately strongspatialpatterns observed in thesubsequent predictions. specifythe LSM spatialparameterdistributions, we choseto There is doubtthat the processof derivingthe spatially calibratethecatchment average, parameters to providethebest variable parameterfields from a few scatteredobservations possiblesimulation. Becauseprecipitationis the dominant usingestablished physicalrelationships andtablelookupis driverof spatialvariabilityin soilmoistureat theWGEW, this sufficientfor realistically constraining a complexdistributed approach resultsin reasonable LSM spatialpredictions. land surfacemodel. The spatiallyvariableparameter fields It is acknowledged that the simulationshouldimproveif were derived from a variety of sources,which leads to therewere sufficientinformationto truly calibratethe LSM inconsistencies in the resultingparameterset and degraded spatially. It may have beenpossibleto use the spatially LSM performance. Ideally,the modelshouldbe calibrated distributed PBMR soil moistureremote-sensing estimates to for eachpointin thedomainto produce consistent parameter calibratethe spatiallydistributed parameters.However,this sets that accountfor the spatial variability of surface was not done becauseour intentionin this studywas to characteristics.For the WGEW this is a poorly posed performa multiobjective calibration that wouldimprovethe problem,beingthat thereare 90,000modelgrid pointsand predictionof a numberof quantities.If the LSM had been
HOUSER ET AL.' MULTIOBJECTIVE CALIBRATION
33,431
Soil Heat Flux
• t
....
210
212
•
214
•.•
216
•
•
218
•
Day of Yoar
•.
•
i•: •.. :!•::'::• ::•:::
•
•4
d)
•8
......... LSMPareto •t Range Obsew• - Calibrat•
35
•6 LSM
SoilMoisture
25
20 o. 15 lO
o
204
206
208
21o
212
214
216
218
Day of Year
220
222
224
226
228
........ LSM Pareto SetRange
• Observed (gravimetric) Observed (resistance)
.......
Calibrated LSM
Figure 5. (continued)
calibrated spatially using only remotely sensed soilmoisture, werederived by optimization. A multiobjective LSM thenthesurface energy fluxprediction skillmayhavecalibration technique was used tocalibrate these parameters by degraded. estimating thenoninferior parameter region based ona large Theuseofcatchment-averaged parameters inthisstudy is number of Monte-Carlo parameter setrealizations, then also supported bythesmall difference between thecalculated selecting the"best" parameter setin a semisubjective way. surface energy fluxes given bya single LSMwithaggregate Thechosen parameter setcontained physically reasonable parameters andthatgiven bya distributed array ofLSMs.values which, inthecase ofthesoilwater retention parameters, to measured values (because of theremoval of This findingalso suggests that a spatiallydistributed LSM weresuperior appliedat theWGEWwiththeresolution usedin thisstudyis thecoarsefractionpriorto measurement). Severaladditionalinsightsinto multiobjective calibration probably notneededif catchment-averaged surfacefluxesare asa resultof additional experiments notreported required.However,a spatially distributed LSM isrequired for weregained andforwhichadditional investigation isrequired. prediction of spatialpatterns of soilmoisture at theresolutioninthispaper Theseincludethefollowing:(1) Someparameters tendto have of the available observations. LSMs have many parameters which must be carefully similarvaluesfor minimumobjectivesfrom severalsites, averageparameters may be specified fortheiroptimalimplementation. Severalparametersindicatingthat watershed However, thespatial patterns present inthePBMR were directlyobservable at Walnut Gulch,and their values acceptable. thatimproved distributed representations maybe couldbe readilyspecified,but otherswere nonphysical and datasuggest
33,432
HOUSER ET AL.: MULTIOBJECTIVE CALIBRATION
possible,and thereforea lumpedrepresentation may not be Chow, V. T., D. R. Maidment,andL. W. Mays, AppliedHydrology, acceptable for all the parameters. (2) Paretooptimalparameter McGraw-Hill, New York, 1988. Daughtry,C. S. T., M. A. Weltz, E. M. Perry, and W. P. Dulaney, sets often produce objective functionswith large error Direct andindirectestimatesof leaf areaindex,paperpresentedat indicating that some model componentsmay require Tenth Conferenceon Biometeorologyand Aerobiology,Special Sessionon Hydrometeorology,Am. Meteorol. Soc., Salt Lake improvement (the soil heatflux component for example).(3) City, Utah, 1991. Someobjectives thatarelinkedthroughmodelphysics,suchas Dickinson, R. E., A. Henderson-Sellers,and P. J. Kennedy, Biosphere-atmosphere transfer scheme(BATS) Version l e as make the multiobjectivecalibrationmore identifiable,only coupledto the NCAR CommunityClimateModel, NCAR Tech.
sensibleand latent heat fluxes, containsimilar information. To
objective functions with significantly different model performance informationshouldbe usedin the calibration
Note 387+STR, Natl. Cent. For Atmos. Res., Boulder, Colo.,
1993.
process. (4)Theinclusion ofsoilparameters inthecalibration Famiglietti, J. S.,Aggregation andScaling of spatially-variable hydrologicalprocesses:local, catchment-scale and macroscale results in significantly betterLSMperformance becausemodels ofwater and energy balance, Ph.D. dissertation, 207pp. measurements of these parameters werelikelybiased owingto
Dep.ofCiv.-Eng. andOper. Res., Princeton Univ., Princeton, N.
the removalof the coarsefraction(i.e., rocksand gravel), the
J., 1992.
presence ofmacropores, etc.(5)A large number ofobjectiveFamiglietti, J.S.,andE.F.Wood, Application ofmultiscale water
functions leads toalarge Pareto setsize, making identification and energy balance models onatallgrass prairie, Water Resour. Res.,30(1 i), 3061-3078, 1994.
of a preferred parameter setimpossible. Therefore, when Franks, S.W.,mad K.J.Beven, Bayesian estimation ofuncertainty in implementing a multiobjective calibration, it is criticalto select a minimumnumberof mostrelevantobjectives.(6) Because
landsurface-atmosphere fluxpredictions, J. Geophys. Res.,102, 23,991-23, 999,1997.
multiple paran•eter setscan be Pareto, multiobjective Franks, S.W.,P.Gineste, K.J.Beven, and P.Merot, Onconstraining the predictionsof a distributedmodel:The incorporation of fuzzy calibrationproblemsare inherentlysubjective.Procedures for estimatesof saturatedareasinto the calibrationprocess,[Vater findingthe Paretoparametersetsare relativelyautomaticand Resour.Res.,34(4), 787-797, 1998. will help to identify parametersets that minimize all the Gupta, H. ¾., S. Sorooshian,and P.O. Yapo, Toward improved calibration of hydrologic models: Multiple and objectivesto some extent. However, to choosea single noncommensurablemeasuresof information, Water Resour. Res., parameterset from thosethat are Pareto,a setof rulesabout 34(4), 751-763, 1998. whichobjective(or trade-offamongobjectives) is mostcritical Gupta,H. V., L. Bastidas,S. Sorooshian,W. J. Shuttleworthand Z. L. Yang, Parameterestimationof a land surfaceschemeusing to minimizemustbe adopted. Alternatively,in this context, multicrilefiamethods,d. Geophys.Res., 104, 19, 491-19, 503, the concept of a preferred parameter set may not be 1999. meaningful. It might be betterto understand the rangeof Hipps, L. E., E. Swiatek, and W. P. Kustas, Interactionsbetween model behaviorsgiven by the Pareto set. The range and regionalsurfacefluxesandthe atmospheric boundarylayerover a characterof the Pareto set may be an indicationof model heterogeneous watershed,Water Resour.Res., 30(5), 1387-1392, physicaldeficiencies, aswhenthe observations lie outsidethe 1994.
pareto setrange. However, a large Pareto setparameter rangeHouser, P.R.,Remote-sensing soilmoisture using four-dimensional dataassimilation,Ph.D. dissertation, 391pp.,Dep. of Hydrol. and may also result from a poorly posedparameterestimation problem.
Water Res., Univ. of Ariz., Tucson,Ariz., 1996. l-louser,P. R., C. Harlow, W. J. Shuttleworth, T. O. Keefer, W. E.
Emmerich, and D.C. Goodrich, Evaluation of multiple flux measurementtechniquesusing water balance informationat a Acknowledgments. This work was supported,fostered,and semiaridsite,paperpresentedat 13th Conferenceon Hydrology, carriedout to completionat the Universityof Arizona,Departmentof Am. Meteorol.Soc.,LongBeach,Calif., 1997. Hydrologyand Water Resources,with financialsupportprovided Houser,P. R., W. J. Shuttleworth, J. S. Famiglietti,H. V. Gupta,K. throughNASA grantsNAGW-4165andNAGW-4087andthroughthe H. Syed,and D.C. Goodrich,1998: Integrationof soil moisture awardof a NASA GlobalChangeFellowship.The studywouldnot remotesensingand hydrologicmodelingusingdata assimilation, havebeenpossiblewithoutthe dataprovidedby theUSDA-Southwest WaterResour.Res.,34(12), 3405-3420, 1998. WatershedResearchCenterandwithouttheir assistance in interpreting Humes,K. S., W. P. Kustas,andM. S. Moran, Useof remotesensing and reference site measurements to estimate instantaneous surface thesedata from the Walnut Gulch ExperimentalWatershed. We are energybalancecomponents over a semiaridrangelandwatershed, pleasedto thankCorrieThiesfor typographic suggestions andseveral WaterResour.Res.,30(5), 1363-1373, 1994a. anonymous reviewersfor very constructive suggestions. Humes, K. S., W. P. Kustas,M. S. Moran, W. D. Nichols, and M. A.
References Amer, S. A., T. O. Keefer, M. A. Weltz, D.C. Goodrich,and L. B. Bach, Soil moisturesensorsfor continuousmonitoring, Water Resour.Bull., 30(1), 69-83, 1994. Bastidas,L., H. V. Gupta, S. Sorooshian,W. J. Shuttleworth,and Z. L. Yang, Sensitivityanalysisof a land surfaceschemeusing multicriteriamethods,d. Geophys.Res., 104, 19, 481-19, 490,
Weitz, Variability of emissivityand surfacetemperatureover a sparselyvegetatedsurface,WaterResour.Res.,30(5), 1299-1310, 1994b.
Jackson,T. J., Measuring surface soil moisture using passive microwaveremotesensing,Hydrol.Processes, 7, 139-152,1993. Jones,H. G., Plantsand Microclimate:A QuantitativeApproachto Environmental Plant Physiology,CambridgeUniv. Press,New York, 1983.
Keefer,T., W. Kustas,D. Goodrich,andK. Quast,Estimatingwater 1999. use by semiaridrangelandvegetationduringthe summerin southernArizona, paper presentedat 13th Conferenceon Beven,K., The TOPMoDEL demonstration programs,Cent.for Res. on Environ. Syst. and Stat., Inst. of Environ. Syst.and Biol. Sci., Hydrology, Am.Meteorol.Soc.,LongBeach,Calif.,1997. LancasterUniv., U.K., 1995. Kustas,W. P. andD.C. Goodrich,Prefaceto the specialsectionon MONSOON 90, WaterResour.Res.,30(5), 1211-1225, 1994. Bouwer, H., Rapid field measurementof air entry value and Kustas,W. P., et al., An interdisciplinary field studyof the energy hydraulicconductivityof soil as significantparametersin flow and water fluxes in the atmosphere-biosphere system over systemanalysis, WaterResour.Res.,2, 729-738, 1966. Carlson,T. N., Recentadvancesin modelingthe infraredtemperature semiaridrangelands: Descriptionand somepreliminaryresults, Bull. Am.Meteorol.Soc.,72(11), 1682-1705,1992. of vegetation canopies, in Land Surface Evaporation; Measurementand Parameterization,editedby T. Schmuggeand Kustas,W. P, J. H. Bianford,D. I. Stannard,C. S. J. Daughtry,W. D. Nichols, and M. A. Weltz, Local energy flux estimatesfor J. C. Andre,pp. 349-358, Springer-Verlag, New York, 1991.
HOUSER ET AL.: MULTIOBJECTIVE
unstableconditionsusingvariancedata in semiaridrangelands. WaterResour.Res.,30(5), 1351-1361, 1994a.
CALIBRATION
33,433
(SiB) for usewithin generalcirculationmodels, JAtrnos.Sci., 43, 505-531, 1986.
Kustas,W. P., M. S. Moran, K. S. Humes,D. I. Stannard,P. J. Pinter,
Sellers,P. J., W. J. Shuttleworth,J. L. Dorman, A. Dalcher, and J. M. L. E. Hipps, E. Swiatek,and D.C. Goodrich,Surfaceenergy Roberts, Calibrating the simple biosphere model (SiB) for balanceestimatesat local andregionalscalesusingopticalremote amazoniantropicalforestusingfield and remotesensingdata, 1, sensingfrom an aircraftplatformand atmospheric datacollected Average calibrationwith field data, J. Appl. Meteorol., 28, 728756, 1989. over semiaridrangelands,WaterResour.Res.,30(5), 1241-1260, 1994b.
Kustas,W. P., T. J. Schmugge, K. S. Humes,T. J. Jackson,R. Parry, M. A. Weltz, and S. Moran, Relationshipsbetweenevaporative fractionand remotelysensedvegetationindex and microwave brightnesstemperature for semiarid rangelands, J. Appl. Meteorol.,32(12), 1781-1790, 1994c. Lettenmaier,D., D. Lohmann,E. F. Wood, and X. Liang, PILPS-2c workshopreport,PrincetonUniv., PrincetonN.J., October,1996. Matthews,N., Interpolatinglarge-scalecontourinformationfrom digital terrain data, U.S. Bureauof Land ManagementInternal Report,U.S. Dep. Of theInter.,Washington, D.C., 1992. Mayeux,B.C., The relativeimportance of rainfall intensityversus saturated hydraulicconductivity for runoffmodelingof a semiarid watershed,Master'sthesis,Dep. of Hydrol. andWater Res.,Univ. of Ariz., Tucson, Ariz., 1995.
Sorooshian, S., Q. Duan, and V. K. Gupta, Calibrationof rainfallrunoff models: application of global optimization to the sacramento soil moistureaccountingmodel, WaterResour.Res., 29(4), 1185-1194, 1993.
Spangler,D. P., A Geophysical studyof the hydrogeology of the Walnut Gulch ExperimentalWatershed,Tombstone,Arizona, Ph.D. dissertation, Dep. of Geol.,Univ. of Ariz., 1969. Stannard,D. I., J. H. Blanford,W. P. Kustas,W. D. Nichols, S. A. Amer, T. J. Schmugge, andM. A. Weltz, Interpretation of surface
flux measurements in heterogeneous terrainduringMonsoon'90, WaterResour.Res.,30(5), 1227-1239, 1994. Storck, P., L.C. Bowling, P. Wetherbee,and D.P. Lettenmaier, Applicationof a GIS-baseddistributedhydrologymodel for predictionof forestharvesteffectson peak streamflow in the PacificNorthwest,Hydrol. Processes,12, 889-904, 1998. Syed,K. H., Spatialstormcharacteristics and basinresponse,M.S. thesis,261 pp., Dep. of Hydrol.andWaterResour.,Univ. of Ariz.,
Meixner,T., H. V. Gupta,L. A. Bastidas,and R. C. Bales,The Use of massflux or concentration for sensitivityanalysisof the alpine hydrochemical model,in press,Hydrol.Processes,1999. Tucson, Ariz., 1994. Menenti, M. and J. C. Ritchie, Estimationof effective aerodynamic Weltz, M. A., J. C. Ritchie,and H. D. Fox, Comparison of laserand roughnessof Walnut Gulch watershedwith laser altimeter fieldmeasurements of vegetation heightandcanopycover,Water measurements,Water Resour.Res.,30(5), 1329-1337, 1994. Moran, M. S., T. R. Clarke, W. P. Kustas,M. Weltz, and S. A. Amer, Evaluationofhydrologicparameters in a semiaridrangelandusing
remotelysensedspectraldata, WaterResour.Res., 30(5), 12871297. 1994a.
Moran, M. S., W. P. Kustas,A. Vidal, De I. Stannard,J. H. Blanford,
andW. D. Nichols,Use of ground-based remotelysenseddatafor surfaceenergybalanceevaluationof a semiaridrangeland, Water Resour.Res.,30(5), 1339-1349, 1994b. RawIs,W. J., L. R. Ahuja,D. L. Brakensiek,andA. Shirmohammadi, Infiltrationand soil water movement,in Handbookof Hydrology, editedby D. R. Maidment,McGraw-Hill, pp. 5.1-5.51, New York,
Tucson, Ariz., ! 996.
J. S. Famiglietti,Departmentof EarthSystemScience,University
1993.
Renard, K. G., L. J. Lane, J. R. Simanton, W. B. Emmerich, J. J. Stone, M. A. Weltz, D.C. Goodrich, and S. D. Yakowitz,
Agricultural impacts in an arid environment:Walnut Gulch studies,Hydrol. $ci. Technol.,9(1-4), 145-190, 1993. Schmugge,T., T. J. Jackson,W. P. Kustas,R. Roberts,R. Parry,D. C. Goodrich, S. A. Amer, and M. A. Weltz, Push broom microwave
Resour.Res.,30(5), 1311-1319, 1994.
White,C. B., P. R. Houser,A.M. Arain,Z.-L. Yang, K. Syed,and W. J. Shuttleworth,The aggregatedescriptionof semiarid vegetation with precipitation-generated soil moisture heterogeneity,Hydrol Earth SystSci, 1,205-212, 1997. Williams, M.D., Evapotranspirationin southeasternArizona semiaridwatersheds: Walnut Gulch and CienegaCreek, M.S. thesis,178pp., Dep. of Hydrol.andWaterResour.,Univ. of Ariz.,
radiometer
of California, Irvine, CA 92697-3100.
H. Guptaand W. J. Shuttleworth,Departmentof Hydrologyand Water Resources,Universityof Arizona, Tucson,AZ 85721. P. R. Houser,HydrologicalSciencesBranch,NASA/GSFC, Code 974, Bldg. 33, Rm. A322, Greenbelt, MD 20771. (Paul.Houser•gsfc.nasa.gov)
observations of surface soil moisture in
Monsoon'90, WaterResour.Res.,30(5), 1321-1328, 1994. Sellers,P. J., Y. Mintz, and A. Dalcher, A simple biospheremodel
(ReceivedApril 22, 1999;revisedSeptember22, 2000; acceptedOctober5, 2000.)