Parameter identifiability, Parameter identifiability, constraints, and constraints, and equifinality in data equifinality in data assimilation with ecosystem assimilation with ecosystem models models Dr. Yiqi Luo Dr. Yiqi Luo Botany and microbiology Botany and microbiology department department University of Oklahoma, University of Oklahoma, USA USA Land surface models and FluxNET data Edinburgh, 4-6 June 2008 (Luo et al. Ecol Appl. (Luo et al. Ecol Appl. In press In press ) )
34
Embed
Dr. Yiqi Luo Botany and microbiology department University of Oklahoma, USA
Parameter identifiability, constraints, and equifinality in data assimilation with ecosystem models. (Luo et al. Ecol Appl. In press ). Dr. Yiqi Luo Botany and microbiology department University of Oklahoma, USA. Land surface models and FluxNET data Edinburgh , 4-6 June 2008. - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Parameter identifiability, Parameter identifiability, constraints, and equifinality constraints, and equifinality
in data assimilation with in data assimilation with ecosystem modelsecosystem models
Dr. Yiqi LuoDr. Yiqi LuoBotany and microbiology Botany and microbiology
departmentdepartmentUniversity of Oklahoma, USAUniversity of Oklahoma, USA
Land surface models and FluxNET dataEdinburgh, 4-6 June 2008
(Luo et al. Ecol Appl. (Luo et al. Ecol Appl. In pressIn press))
Wang et al. (2001) ------ a maximum of Wang et al. (2001) ------ a maximum of 3 or 43 or 4 p parameters can be determined.arameters can be determined. Braswell et al. (2005) ------ Braswell et al. (2005) ------ 13 out of 2313 out of 23 parame parameters were well-constrained.ters were well-constrained. Xu et al. (2006) ------ Xu et al. (2006) ------ 4 or 3 out of 74 or 3 out of 7 parameters parameters can be constrained, respectively under ambiecan be constrained, respectively under ambient and elevated COnt and elevated CO22..
Identiable parameters
Three methods to exaThree methods to examine parameter identimine parameter identifiabilityfiability1.1. Search methodSearch method2.2. Model Model
structurestructure3.3. Data variabilityData variability
Harvard Forest EMS-Tower
Eddy flux data Eddy flux data
COCO22 flux flux HH22O fluxO flux Wind speedWind speed TemperatureTemperature PARPAR Relative humidityRelative humidityHourly or half-hourlyHourly or half-hourly
Eddy flux technology
Leaf-level Photosynthesis
Sub-model
Canopy-level Photosynthesis
Sub-model
System-level C balanceSub-model
ModelModel
Table 1 Parameters informationTable 1 Parameters information
Develop prior distributionDevelop prior distribution
a) generate candidate a) generate candidate pp from sample space from sample spaceb) input to model and calculate cost functionb) input to model and calculate cost functionc) select according to decision criterionc) select according to decision criteriond) repeatd) repeat
Construct posterior distributionConstruct posterior distribution
Fig. 2 Decrease of cost function with each step of conditional inversion
ConclusionsConclusions Conditional inversion can Conditional inversion can
substantially increase the number of substantially increase the number of constrained parameters.constrained parameters.
Cost function and information loss Cost function and information loss decrease with each step of conditional decrease with each step of conditional inversion.inversion.
Measurement errors Measurement errors and parameter identiand parameter identifiabilityfiability
Leaves X1 Woody X2 Fine Roots X3
Metabolic Litter X4 Structural Litter X5
Microbes X6
Slow SOM X7
Passive SOM X8
GPP
TECO – biogeochemical model
)(
100000010000
1000000101000010000001000000001000000001
8786
7675
68676564
5351
4341
cdiagC
ffff
ffffffff
A
0)0(
)()()(
XtX
tBPtACXtXdtd
TbbbB )00000( 321
No. of parameter
8
12
8
3
c1
010
2030
40
0.0E+00 1.0E-03 2.0E-03 3.0E-03 4.0E-03 5.0E-03
2.0-SD
1.0-SD
0.5-SD
c2
0
5
10
15
0 0.00005 0.0001 0.00015 0.0002 0.00025
2.0-SD
1.0-SD
0.5-SD
c3
05
10152025
0 0.002 0.004 0.006 0.008
2.0-SD
1.0-SD
0.5-SD
c4
0
1
2
3
4
0 0.01 0.02 0.03 0.04
c5
0
2
4
6
0 0.001 0.002 0.003
c6
0
20
40
60
80
0 0.1 0.2 0.3 0.4 0.5
c7
0
2
4
6
8
10
12
0 0.0005 0.001 0.0015
c8
0
0.5
1
1.5
2
2.5
0 2E-06 4E-06 6E-06 8E-06 0.00001
Exit rates
A1
0
1
2
3
0 0.5 1 1.5
2.0-SD
1.0-SD
0.5-SD
A4
012345
0 0.2 0.4 0.6 0.8 1
2.0-SD
1.0-SD
0.5-SD
A5
0
1
2
3
4
0 0.2 0.4 0.6
2.0-SD
1.0-SD
0.5-SD
A6
0
1
2
3
0 0.2 0.4 0.6
A7
0
1
2
3
0 0.2 0.4 0.6 0.8
A8
0
1
2
3
4
0 0.1 0.2 0.3 0.4
A9
0
1
2
3
0 0.2 0.4 0.6 0.8
A10
0
1
2
3
0 0.1 0.2 0.3 0.4
A11
0
1
2
3
0 0.5 1 1.5
Transfer coefficients
Xo1
02
46
8
0 100 200 300 400 500
2.0-SD
1.0-SD
0.5-SD
Xo2
0
510
15
20
0 2000 4000 6000
2.0-SD
1.0-SD
0.5-SD
Xo3
0
12
3
4
0 100 200 300 400
2.0-SD
1.0-SD
0.5-SD
Xo4
00.5
11.5
22.5
0 20 40 60 80
Xo5
0
1
2
3
4
0 100 200 300 400 500
Xo6
0
1
2
3
0 50 100 150
Xo7
0
5
10
15
20
0 1000 2000 3000 4000
Xo8
0
5
10
15
20
0 200 400 600 800
Initial values
X1
500 1000 1500 2000
Freq
uenc
y 10
2
05
10152025
X2
5000 6000 7000 8000 900001234
X3
500 1000 1500 2000
Freq
uenc
y 10
2
05
10152025 X4
500 1000 1500 20000
10203040
X5
500 1000 1500 2000
Freq
uenc
y 10
2
02468
10 X6
50 100 150 2000
20
40
X7
Carbon content (g C m-2)1000 2000 3000 4000
Freq
uenc
y 10
2
02468
10X8
Carbon content (g C m-2)
450 500 550 600 65002468
Pool sizes without data
X1
300 400 500 600 700 800
Freq
uenc
y 10
2
05
10152025
X2
5000 6000 7000 8000 900005
10152025
X3
100 200 300 400 500
Freq
uenc
y 10
2
05
10152025
X4
100 200 300 400 50005
10152025
X5
0 500 1000 1500 2000
Freq
uenc
y 10
2
05
10152025 X6
0 50 100 150 20005
10152025
X7
Carbon content (g C m-2)1500 2000 2500 3000 3500
Freq
uenc
y 10
2
05
10152025
halved SDambient SDdoubled SD
X8
Carbon content (g C m-2)
700 800 900 1000 1100 1200 130005
10152025
Pool sizes with data and different SD
ConclusionConclusion
Magnitudes of measurement errMagnitudes of measurement errors do not affect parameter idenors do not affect parameter identifiability but influence relative tifiability but influence relative constraints of parametersconstraints of parameters
Base modelBase modelGPP
Leaves X1 Stems X2 Roots X3
Metabolic L. X4 Struct. L. X5
Microbes X6
Slow SOM X7
Passive SOM X8
Simplified modelsSimplified models
Plant C
Litter C
GPP CO2
Soil C
Plant C
Litter C
GPP CO2
O Soil C
Miner. C
3P model 4P model
Simplified modelsSimplified models6P model 7P model
Differences in model structure Differences in model structure are corresponding to different are corresponding to different sets of parameters. The number sets of parameters. The number of constrained parameters varies of constrained parameters varies with data availability with data availability