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

description: Data, model runs and codes used in 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. 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 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 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 Fc time series with 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 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 model domain.<br> 3.1. <em>SUEWS_SourceCode</em> is a folder including SUEWS V2020b source Fortran codes. For detailed description, 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 model domain, covering circles with a 500 m, 750 m, 1000 m, 1500 m radius, respectively. Note that except the size of model domain, the parameters are the same as model run <em>VegetationSensitivity\case-gs_LAI</em>. <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>.