The code and data to the article"Can carbon economy explain leaf dynamic seasonality in a tropical seasonal rain forest?"

This README file is associated with the code and data to the article"Can carbon economy explain leaf dynamic seasonality in a tropical seasonal rain forest?"<br>Description on the code and data which involved in the study. 1. The Forest-CEW modelThis model is a self-developed model by Jiye Zeng and Zheng-Hong Tan which first finished in 2018 and a updated version in 2022. The model is portable. Users double click the fcew.exe file, user-interface will be jump out. Users could easily use the model follow the user interfact. For advance user, the source code is provided which written C++. <br>2. The reproducing of figure 2In Figure 2, the leaf area index was estimated by couple of sensors which placed both above the canopy and below the canopy according to the Beer-Lambert Law. The matlab code for walsh seasonality index was also provided. The solar zenith angle is important in apply Beer-Lambert law in forest canopy. The funtion SolarZenith is used to accomplish this. there is a data save in WSIdata.csv. The meanning of each column are: column 1: time column 2: latent heat flux (w m-2) column 3: sensible heat flux (w m-2) column 4: friction velocity (m s-1) column 5: air density (kg m-3) column 6: downward shortwave radiation (w m-2) column 7: upward shortwave radiation (w m-2) column 8: downward longwave radiation (w m-2) column 9: upward longwave radiation (w m-2) column 10: canopy temperature measured by Apogee sensor (℃) column 11: surface temperature at forest floor measured by Apogee sensor (℃) column 12: soil temperature at 5 cm depth (℃) column 13: soil temperature at 20 cm depth (℃) column 14: soil temperature at 40 cm depth (℃) column 15: air temperature at level 4 (℃) column 16: air temperature at level 5 (℃) column 17: air temperature at level 6 (℃) column 18: air temperature at level 7 (℃) column 19: wind speed at level 4 (m s-1) column 20: wind speed at level 5 (m s-1) column 21: wind speed at level 6 (m s-1) column 22: wind speed at level 7 (m s-1) column 23: water vapor pressure at level 4 (kPa) column 24: water vapor pressure at level 5 (kPa) column 25: water vapor pressure at level 6 (kPa) column 26: water vapor pressure at level 7 (kPa) column 27: wind direction (degree) column 28: standard deviation of wind direction column 29: precipitation (mm) column 30: solar radiation at 41 m (w m-2) column 31: solar radiation at 70 m (w m-2) column 32: photosynthetically active radiation at 4 m (umol m-2 s-1) column 33: photosynthetically active radiation at 10 m (umol m-2 s-1) column 34: photosynthetically active radiation at 22 m (umol m-2 s-1) column 35: photosynthetically active radiation at 31 m (umol m-2 s-1) column 36: photosynthetically active radiation at 36 m (umol m-2 s-1) column 37: photosynthetically active radiation at 70 m (umol m-2 s-1)Litterfall data were saved in .csv format year-by-year. column 1: plot number column 5: leaf litter<br>3. The reproducing of figure 3(a) The leaf demography model could be run when P and L set as input as well as a and b are determined. The model could quickly stablized no matter how initail value of X,Y and Z were set (see main text for their meanings). (b) how optimal leaf longevity varied with constructive value and the date of leaf flush. The environmental driving factors was store in F55_2003.csv (meteorological data) column 1: time column 18: light intensity, umol PAR photons m-2 s-1 column 38: air temperature, ℃ column 40: relative humidity, %<br>4. The reproducing of figure 4 Open the agePhoto.m file. Figure 4 could be reproduced. Please modify the file direction in your case of using. <br>5. The reproducing of figure 5 Several matlab code involved in reproducing figure 5. CANOPYPHS_2LEAF.m is one the key main program which could calculate canopy photosynthesis using the sunlit-shade scaling model. The driving varibles as stored in ClimF0.csv (meteorological data) column 1: day of the year column 4: air temperature, ℃ column 5: vapour pressure deficit, kPa column 6: light intensity, umol PAR photons m-2 s-1 column 7: air CO2 concentration<br>6. The code for Section 2.10 statistical analysisThe code for the significance test and seasonality analysis is located in the Statistical_Analysis folder and is provided in MATLAB. The code is accompanied by data and can be used to reproduce the results of the manuscript.<br>Any question on using the code could email to [email protected]. Since this work is not a case study (data is avaiable for most of the tropical forest site), users who want to simulate with this method in their site are warmly welcome. <br>