DOI: 10.5281/zenodo.8012194
title: Datasets for simulating heat and CO2 fluxes in Beijing using SUEWS V2020b

description: Data, model runs and codes used in the manuscript "<em>Simulating heat and CO<sub>2</sub> fluxes in Beijing using SUEWS V2020b: sensitivity to vegetation phenology and maximum conductance</em>". The source code for Surface Urban Energy and Water balance Scheme (SUEWS) V2020b is also included. <strong>[Latest Version]</strong> June 9, 2023: Version 2 was created. (1) Modifications in <strong>LAI</strong>: the input data of LAI (<em>LAI/source_data/Beijing_LAI.csv</em>) was previously derived from a MODIS product MOD15A2H; it was replaced with the LAI time series at a higher spatial resolution derived from data provided by Landsat 7 ETM+; the code for data retrieval (<em>LAI/LAI_from_GEE.txt</em>) was added. (2) Modifications in <strong>ModelRuns</strong>: some of the external parameters given to SUEWS were adjusted, including the optimized LAI parameters (for case LAI and case gs_LAI), and Fr<sub>nonheat</sub> (from 0.2 to 0.5); a <em>read_me.txt</em> was added. (3) Modifications in <strong>Fig&Statistics</strong>: the shaded area for the line charts were changed from standard deviation to the interquartile range. <strong>[Data&Code Description]</strong> The folders are: <strong>1. Observations</strong><br> The units for CO2 flux (Fc), latent heat flux (QE), and sensible heat flux (QH) are μmol m<sup>-2</sup> s<sup>-1</sup>, W m<sup>-2</sup>, and W m<sup>-2</sup>, respectively, unless otherwise stated.<br> <strong>1.1. Observation datasets</strong><br> 1.1.1. <em>radiation_flux_2010-2011.csv</em> is radiation flux observations from May 2010 to July 2011 at 140 m on the IAP tower, Beijing.<br> 1.1.2. <em>turbulent_flux_before_QC_2016.csv</em> is turbulent flux observations before quality control (QC) obtained with the eddy covariance (EC) technique for the year 2016 at 47 m on the IAP tower.<br> 1.1.3. <em>turbulent_flux_after_QC_2016.csv</em> is turbulent flux observations after quality control (QC).<br> 1.1.4. <em>Fc_mean_diurnal_by_season.csv</em> is the mean diurnal seasonal cycle of Fc observations.<br> 1.1.5. <em>Fc_gapfilled_with_MDC.csv</em> is the Fc time series for the year 2016 gap-filled with Mean Diurnal Cycle (MDC) method.<br> 1.1.6. <em>meteorological_observations </em>is a folder including meteorological observations at different heights measured at the IAP tower for the years 2010-2012.<br> <strong>1.2 Codes for observations processing</strong><br> 1.2.1. <em>turbulent_flux_QC.py</em> is for turbulent flux QC covering wind direction filtering, stationarity test, friction velocity filtering, and nighttime filtering.<br> 1.2.2. <em>mean_diurnal_cycle_of_flux_observations.py</em> is for obtaining the mean diurnal cycle of flux observations.<br> 1.2.3. <em>Fc_gap_filling.py</em> is for filling the gaps in the Fc time series with the MDC method. <strong>2. LAI</strong><br> This folder includes a project aiming to obtain the parameters for the leaf area index (LAI) model adopted by SUEWS. In this LAI model, the LAI is related to air temperature only. Incorporating hourly air temperature observations and the 'real' LAI (e.g. observed LAI, remotely sensed LAI), the optimization method Covariance Matrix Adaptation Evolutionary Strategies (CMA-ES) will derive the optimized parameters for the LAI model. For more details, please read the description text file (<em>read_me.txt</em>) under this folder.<br> 2.1. <em>read_me.txt.</em><br> 2.2. <em>source_data</em> is a folder including the input for CMA-ES optimization. The input files of Beijing are included, allowing a quick run to reproduce the results demonstrated in the manuscript.<br> 2.3. <em>result_data</em> is a folder including the output for CMA-ES optimization.<br> 2.4. The code for data retrieval from the Google Earth Engine platform (<em>LAI/LAI_from_GEE.txt</em>) was added. The others are Python scripts regarding CMA-ES, including the LAI model, data de-spiking, interpolation, CMA-ES training, and result visualization. <strong>3. ModelRuns</strong><br> This folder includes SUEWS source code, four model runs conducted to test the sensitivity to vegetation-related parameters (i.e., maximum conductance, LAI model parameters), and four model runs to test sensitivity to the radius of the model domain.<br> 3.1. <em>SUEWS_SourceCode</em> is a folder including SUEWS V2020b source Fortran codes. For detailed descriptions, readers are referred to SUEWS webpage (https://suews.readthedocs.io/en/latest/).<br> 3.2. <em>Input</em> is a folder including the necessary input files shared by each of the model runs. Please copy all the files to a model run folder (e.g. <em>VegetationSensitivity\case-base\Input\</em>) before starting a SUEWS run.<br> 3.3. <em>VegetationSensitivity</em> is a folder including four SUEWS runs to test the model sensitivity to vegetation-related parameters: <em>case-base</em> is the control run, <em>case-gs</em> is a model run where maximum conductance of vegetation has been modified to more site-specific values, <em>case-LAI</em> is a model run where LAI parameters have been optimized with CMA-ES, and <em>case-gs_LAI</em> is a model run where both maximum conductance and LAI parameters have been modified. Note that all the model runs have selected 1000 m as the radius of the model domain. To conduct a quick model run to reproduce the output, enter the command "./SUEWS_V2020b" in the model run folder (e.g. <em>VegetationSensitivity\case-base\</em>) under a Linux environment (e.g., Ubuntu).<br> 3.4. <em>RadiusSensitivity</em> is a folder including four SUEWS runs to test the model sensitivity to the size of the model domain, covering circles with a 500 m, 750 m, 1000 m, 1500 m radius, respectively. Note that the vegetation-related parameters are the same as the model run <em>VegetationSensitivity\case-gs_LAI</em>. 3.5. read_me.txt is a note on how to reproduce the SUEWS cases shown in the manuscript. <strong>4. Fig&Statistics</strong><br> This folder includes Python scripts and data to reproduce a portion of the figures and statistics in the manuscript quickly. They are sorted into the following sub-folders:<br> 4.1 <em>validation_WFDE5</em>,<br> 4.2 <em>evaluation_radiation</em>,<br> 4.3 <em>evaluation_turbulent</em>,<br> 4.4 <em>accumulated_Fc</em>,<br> 4.5 <em>Fc_component_MDC</em>.