Model type and unit issue 4

[jinshijian]

390

• Model has soil moisture component, add Q16 to it.

1384

• If use the reported Model_paramB values, the soil temperature has to be over 50 C to get the reported Rs_Annual;
• Model_paramB = Model_paramB*10? Because this give us a reasonable results (required soil temperature = 5 C to get the reported Rs_annual).

2182

• The model Parameters should have an issue, because Ts should be around 13C to get the reported Rs_annual according to the model, the required soil temperature is too higher than T_annual from the Delaware climate data (-2 C).
• Cannot figure it out what is the problem.

4756

• The model Parameters should have an issue, because Ts should be around 28C to get the reported Rs_annual according to the model, the required soil temperature is too higher than T_annual from the Delaware climate data (9 C).
• Cannot figure it out what is the problem.

2560

• Study_temp are available from 'Table 1', these information can be added:
• Study_temp = 26 WHERE Record_number = 1134
• Study_temp = 26.2 WHERE Record_number = 1135
• Study_temp = 26.6 WHERE Record_number = 1136
• Study_temp = 30 WHERE Record_number = 1137
• Study_temp = 18.4 WHERE Record_number = 1138
• Study_temp = 20 WHERE Record_number = 1139

3976

• Study_temp =13 according to Figure 1, this information can be added.

3649

• Study_temp = 6.2 according to Figure 1.

5162

• The model is a seasonal model, not annual, add Q15

7530

• Model_paramC = 10

7596

• The model in Figure 3 is a hourly model
• I re-simulated the model based on the data from figure 5&6
• Model type: "Polynomial, R=exp(a+b(T-d)+c(T-d)^2)"
• Model_paramA = -798.76, Model_paramB = 71.714, Model_paramC = -1.36 WHERE Record_number = 5822
• Model_paramA = -773.22, Model_paramB = 69.774, Model_paramC = -1.3233 WHERE Record_number = 5823

7704

• Model_type = 'Exponential, R=a exp(b(T-c))' according to Figure 5
• The unit for the Rs_annual should be reported wrong in the mamuscript, if use g c/m2/yr, change to umol CO2/m2/s, we get (780=2.06, 840=2.22, 1175=3.11,647=1.72,1448=3.83);
• But from Figure 5, we can see that the Rs_Annual is much smaller than above values.
• I guess the unit of Rs_annual should change to g co2/m2/yr, because the values perfectly mached the results from Figure 5.
• Rs_annual = Rs_annual(current values)*44/12

8007

• Model_type = 'Exponential, R=a exp(b(T-c))'.

8079

• According to Figure 8 upright panel, model only for summer and spring (record number 6621).
• I Re-simulated the model based on the data from Figure 8 (upright panel): Model_paramA = -0.0075, Model_paramB = 0.141, Model_paramC = -0.0031.
• Model_paramD = 0 Record_number = 6622.

8334

• Model_paramA values reported wrong in the manuscript, because if use these values, the temperature has to be over 60 C to get the annual_Rs.
• If Model_paramA devide by 10, the required temperature is about 17 C to get the reported annual_Rs.
• Model_paramA = Model_paramA(original)*10

8382

• The model reported in the paper has a soil water content component.
• I re-simulated the model (only for control, record number 5882) use the same model formate, but without SWC conponent,
• for control (Recordnumber = 5882) Model_paramA = 0.995, Model_paramB = 250.96, Model_type = 'R10 (L&T), R=a exp(b((1/c)-(1/(T-d))), T in K'

8866

• The unit may reported worong in the manuscript,
• If model_output_units is g co2/m2/hr, the temperature should be -0.6 C to get the reported Rs_Annual according to the model, the temperature is far way from the T_Annual (12 C) in this site,
• If the model_output_units change to g c/m2/d, the annual TS mean is 15 C, close to the T_Annual (12 C from delaware climate data).
• Update the unit?

8917

• Model_paramB = 0.0331 WHERE Record_number = 5973 (according to Table 2)
• Model_paramB = 0.0541 WHERE Record_number = 5979 (according toTable 2)

9563

• For record number 5962, 5963, 5964, according to the model parameters, the calculated temperature is around 5 C to get the reported Rs_annual, the required temperature is much lower than the annual mean temperature from Delaware climate data (21 C);
• But for record number 5959, 5960, and 5961, the required temperature is close to 20 C;
• I cannot find the model parameters in the manuscript, so I am not able to identify the problem, please let me know if you got those information from the manuscript.

[bpbond]

1384

If use the reported Model_paramB values, the soil temperature has to be over 50 C to get the reported Rs_Annual;
Model_paramB = Model_paramB*10? Because this give us a reasonable results (required soil temperature = 5 C to get the reported Rs_annual).

I looked at this and yes, must be. I'll fix and add a note.

[bpbond]

2182

The model Parameters should have an issue, because Ts should be around 13C to get the reported Rs_annual according to the model, the required soil temperature is too higher than T_annual from the Delaware climate data (-2 C).
Cannot figure it out what is the problem.

I checked and the parameters in Table 1 do produce roughly the curves shown in Figure 1. For the first equation, the parameters at -2 C = about 0.2 µmol/m2/s = ~90 gC/m2/yr via Bahn equation. But the paper reports much higher annual numbers. Right? Huh, I wonder what happened. 🤷‍♂️

[bpbond]

4756

The model Parameters should have an issue, because Ts should be around 28C to get the reported Rs_annual according to the model, the required soil temperature is too higher than T_annual from the Delaware climate data (9 C).
Cannot figure it out what is the problem.

This is listed in the database as first author "Brye", but we're talking about Cahill et al. (2009) right?

[bpbond]

7704

Model_type = 'Exponential, R=a exp(b(T-c))' according to Figure 5
The unit for the Rs_annual should be reported wrong in the mamuscript, if use g c/m2/yr, change to umol CO2/m2/s, we get (780=2.06, 840=2.22, 1175=3.11,647=1.72,1448=3.83);
But from Figure 5, we can see that the Rs_Annual is much smaller than above values.
I guess the unit of Rs_annual should change to g co2/m2/yr, because the values perfectly mached the results from Figure 5.
Rs_annual = Rs_annual(current values)*44/12

I'm a little confused here. I should (i) change the equation type and (ii) change the Rs_annual values to Rs_annual * 44/12? What about Model_output_units?

[bpbond]

8866

The unit may reported worong in the manuscript,
If model_output_units is g co2/m2/hr, the temperature should be -0.6 C to get the reported Rs_Annual according to the model, the temperature is far way from the T_Annual (12 C) in this site,
If the model_output_units change to g c/m2/d, the annual TS mean is 15 C, close to the T_Annual (12 C from delaware climate data).
Update the unit?

Agreed. Done.

[bpbond]

9563

For record number 5962, 5963, 5964, according to the model parameters, the calculated temperature is around 5 C to get the reported Rs_annual, the required temperature is much lower than the annual mean temperature from Delaware climate data (21 C);
But for record number 5959, 5960, and 5961, the required temperature is close to 20 C;
I cannot find the model parameters in the manuscript, so I am not able to identify the problem, please let me know if you got those information from the manuscript.

In the Notes field it says: "Q10 models from 8917" (8917 Han in Scientific Reports).

[jinshijian]

2182

Yes, the reported Rs_annual is too high if the model parameters are correct.

4756

Yes, Cahill et al. (2009).

7704

Sorry for the unclear description, if change the Rs_annual values, the Model_output_units should maintain the same. So, only need two changes: Model_type change to 'Exponential, R=a exp(b(T-c))', Rs_annual = Rs_annual(current values)/44*12

9563

Thanks, according table 2 in 8917, Model_paramB should change, and the estimates match well with the reported Rs_annual. SET Model_paramB = 0.0375 WHERE Record_number = 5959 SET Model_paramB = 0.0377 WHERE Record_number = 5960 SET Model_paramB = 0.0272 WHERE Record_number = 5961 SET Model_paramB = 0.0541 WHERE Record_number = 5962 SET Model_paramB = 0.045 WHERE Record_number = 5963 SET Model_paramB = 0.0311 WHERE Record_number = 5964

[bpbond]

Changed 7704 and 9563.

Re 2182 and 4756–let's go over later today.