This dataset contains the flux measurements from the large aperture scintillometer (LAS) at Sidaoqiao Superstation in the Heihe integrated observatory network from January 1 to December 31 in 2018. There were one German BLS900 at Sidaoqiao Superstation. The north tower was set up with the BLS900 receiver and the south tower was equipped with the BLS900 transmitter. The site (north: 101.137° E, 42.008° N; south: 101.131° E, 41.987 N) was located in Ejinaqi, Inner Mongolia. The underlying surfaces between the two towers were tamarisk, populus, bare land and farmland. The elevation is 873 m. The effective height of the LAS was 25.5 m, and the path length was 2350 m. The data were sampled 1 minute. The raw data acquired at 1 min intervals were processed and quality controlled. The data were subsequently averaged over 30 min periods, in which sensible heat flux was iteratively calculated by combining Cn2 with meteorological data according to the Monin-Obukhov similarity theory. The main quality control steps were as follows: (1) The data were rejected when Cn2 exceeded the saturated criterion (Cn2>7.58E-14). (2) The data were rejected when the demodulation signal was small (Average X Intensity<1000). (3) The data were rejected when collected during precipitation. (4) The data were rejected if collected at night when weak turbulence occurred (u* was less than 0.1 m/s). In the iteration process, the universal functions of Thiermann and Grassl, 1992 was selected. Detailed can refer to Liu et al. (2011, 2013). Several instructions were included with the released data. (1) The missing data from the BLS900 instrument were denoted by -6999. (2) The dataset contained the following variables: Date/time (yyyy/m/d h:mm), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). In this dataset, a time of 0:30 corresponds to the average data for the period between 0:00 and 0:30, and the data were stored in *.xlsx format. Moreover, suspicious data were marked in red. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) (for data processing) in the Citation section.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This dataset contains the flux measurements from the large aperture scintillometer (LAS) at Arou Superstation in the Heihe integrated observatory network from January 1 to December 31 in 2018. There were two types of LASs at Arou Superstation: BLS450 and zzlas, produced by Germany and China, respectively. The north tower was set up with the zzlas receiver and the BLS450 transmitter, and the south tower was equipped with the zzlas transmitter and the BLS450 receiver. The site (north: 100.471° E, 38.057° N; south: 100.457° E, 38.038° N) was located in Caodaban village of A’rou town in Qilian county, Qinghai Province. The underlying surface between the two towers was alpine meadow. The elevation is 3033 m. The effective height of the LASs was 9.5 m, and the path length was 2390 m. The data were sampled 1 minute at both BLS450 and zzlas. The raw data acquired at 1 min intervals were processed and quality controlled. The data were subsequently averaged over 30 min periods, in which sensible heat flux was iteratively calculated by combining Cn2 with meteorological data according to the Monin-Obukhov similarity theory. The main quality control steps were as follows: (1) The data were rejected when Cn2 exceeded the saturated criterion (BLS450: Cn2>7.25E-14, zzlas: Cn2>7.84E-14). (2) The data were rejected when the demodulation signal was small (BLS450: Mininum X Intensity<50; zzlas: Demod>-20mv). (3) The data were rejected when collected during precipitation. (4) The data were rejected if collected at night when weak turbulence occurred (u* was less than 0.1 m/s). In the iteration process, the universal functions of Thiermann and Grassl, 1992 and Andreas, 1988 were selected for BLS450 and zzlas, respectively. Detailed can refer to Liu et al. (2011, 2013). Several instructions were included with the released data. (1) The data were primarily obtained from BLS450 measurements, and missing flux measurements from the BLS450 instrument were substituted with measurements from the zzlas instrument. The missing data were denoted by -6999. Due to the problems of storing and wireless transmission, data from 5 July to 24 August, were not collected. (2) The dataset contained the following variables: Date/time (yyyy/m/d h:mm), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). In this dataset, a time of 0:30 corresponds to the average data for the period between 0:00 and 0:30, and the data were stored in *.xlsx format. Moreover, suspicious data were marked in red. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) (for data processing) in the Citation section.
LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This dataset contains the flux measurements from the Sidaoqiao superstation eddy covariance system (EC) in the downstream reaches of the Heihe integrated observatory network from January 1 to December 31 in 2018. The site (101.1374° E, 42.0012° N) was located in the Ejina Banner in Inner Mongolia Autonomous Region . The elevation is 873 m. The EC was installed at a height of 3.2 m, and the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (CSAT3&Li7500A) was 0.15 m. The raw data acquired at 10 Hz were processed using the Eddypro post-processing software, including the spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. The observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC): class 1-3 (high quality), class 4-6 (good), class 7-8 (poor, better than gap filling data), class9 (rejected). In addition to the above processing steps, the half-hourly flux data were screened in a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data collected before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 10% of the 30 min raw record. There were 48 records per day, and the missing data were replaced with -6999. Suspicious data were marked in red. Latent heat flux during November 9 to 21, 2018 were missing due to the sensor malfunction of CO2/H2O gas analyzer. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m3), CO2 mass density (CO2, mg/m3), friction velocity (ustar, m/s), stability (L), sensible heat flux (Hs, W/m2), latent heat flux (LE, W/m2), carbon dioxide flux (Fc, mg/ (m2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xls format. Detailed information can be found in the suggested references. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) for data processing) in the Citation section.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This dataset contains the flux measurements from the Jingyangling station eddy covariance system (EC) in the upperstream reaches of the Heihe integrated observatory network from August 28 to December 31 in 2018. The site (101.1160E, 37.8384N) was located in the Jingyangling, near Qilian County in Qinghai Province. The elevation is 3750 m. The EC was installed at a height of 4.5 m, and the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (CSAT3&Li7500) was 0.15 m. The raw data acquired at 10 Hz were processed using the Eddypro post-processing software, including the spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. The observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC): class 1-3 (high quality), class 4-6 (good), class 7-8 (poor, better than gap filling data), class9 (rejected). In addition to the above processing steps, the half-hourly flux data were screened in a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data collected before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 10% of the 30 min raw record. There were 48 records per day, and the missing data were replaced with -6999. Suspicious data were marked in red. Data during insufficient power supply, data were missing occasionally. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m3), CO2 mass density (CO2, mg/m3), friction velocity (ustar, m/s), stability (L), sensible heat flux (Hs, W/m2), latent heat flux (LE, W/m2), carbon dioxide flux (Fc, mg/ (m2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xls format. Detailed information can be found in the suggested references. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) for data processing) in the Citation section.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This dataset contains the flux measurements from the large aperture scintillometer (LAS) at Daman Superstation in the Heihe integrated observatory network from January 1 to December 31 in 2018. There were two types of LASs at Daman Superstation: BLS450 and BLS900, produced by Germany. The north tower was set up with the BLS450 receiver and the BLS900 transmitter, and the south tower was equipped with the BLS450 transmitter and the BLS900 receiver. The site (north: 100.379° E, 38.861° N; south: 100.369° E, 38.847° N) was located in Daman irrigation district, which is near Zhangye, Gansu Province. The underlying surfaces between the two towers were corn, orchard, and greenhouse. The elevation is 1556 m. The effective height of the LASs was 22.45 m, and the path length was 1854 m. The data were sampled 1 minute at both BLS450 and BLS900. The raw data acquired at 1 min intervals were processed and quality controlled. The data were subsequently averaged over 30 min periods, in which sensible heat flux was iteratively calculated by combining Cn2 with meteorological data according to the Monin-Obukhov similarity theory. The main quality control steps were as follows: (1) The data were rejected when Cn2 exceeded the saturated criterion (Cn2>1.43E-13). (2) The data were rejected when the demodulation signal was small (Average X Intensity<1000). (3) The data were rejected when collected during precipitation. (4) The data were rejected if collected at night when weak turbulence occurred (u* was less than 0.1 m/s). In the iteration process, the universal functions of Thiermann and Grassl, 1992 was selected. Detailed can refer to Liu et al. (2011, 2013). Several instructions were included with the released data. (1) The data were primarily obtained from BLS900 measurements, and missing flux measurements from the BLS900 instrument were substituted with measurements from the BLS450 instrument. The missing data were denoted by -6999. (2) The dataset contained the following variables: Date/time (yyyy/m/d h:mm), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). In this dataset, a time of 0:30 corresponds to the average data for the period between 0:00 and 0:30, and the data were stored in *.xlsx format. Moreover, suspicious data were marked in red. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) (for data processing) in the Citation section.
LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This dataset contains the flux measurements from the Huazhaizi station eddy covariance system (EC) in the midstream reaches of the Heihe integrated observatory network from January 1 to December 31 in 2018. The site (100.3201° E, 38.7659° N) was located in the Zhangye City in Gansu Province. The elevation is 1731 m. The EC was installed at a height of 4.5 m, and the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (CSAT3&Li7500A) was 0.15 m. The raw data acquired at 10 Hz were processed using the Eddypro post-processing software, including the spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. The observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC): class 1-3 (high quality), class 4-6 (good), class 7-8 (poor, better than gap filling data), class9 (rejected). In addition to the above processing steps, the half-hourly flux data were screened in a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data collected before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 10% of the 30 min raw record. There were 48 records per day, and the missing data were replaced with -6999. Suspicious data were marked in red. Data during May 13 to July 12 and July 16 to August 21, 2018 were missing due to the malfunction of CO2/H2O gas analyzer. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m3), CO2 mass density (CO2, mg/m3), friction velocity (ustar, m/s), stability (z/L), sensible heat flux (Hs, W/m2), latent heat flux (LE, W/m2), carbon dioxide flux (Fc, mg/ (m2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xls format. Detailed information can be found in the suggested references. For more information, please refer to Liu et al. (2018) (for sites information), Liu et al. (2011) for data processing) in the Citation section.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
The data set contains eddy covariance System observation data of Daman super station which is located in the middle reaches of the Heihe Hydro-meteorological Observation Network from January 1, 2017 to December 31, 2017. The site is located in Daman Irrigation District, Zhangye, Gansu Province, and the underlying surface is corn. The latitude and longitude of the observation point is 100.37223E, 38.85551N, and the altitude is 1556.06m. The mount height of the Eddy Covariance System is 4.5 m, the sampling frequency is 10 Hz, the ultrasonic orientation is positive North, and the distance between the ultrasonic wind speed temperature meter (CSAT3) and the CO2/H2O analyzer (Li7500) is 17 cm. The original observation data of the Eddy Covariance System is 10 Hz, and the released data is a 30-minute data processed by Eddypro software. The main steps of the processing include: outlier eliminating, delay time correction, coordinates rotation (secondary coordinates rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction, etc. Meanwhile, the quality evaluation of each flux value was performed,mainly includes atmospheric stability (Δst) test and turbulence similarity (ITC) test. The 30-min flux value output of Eddypro software was also screened: (1) Data from the instrument error was eliminated; (2) Data obtained with one hour before and after precipitation was removed; (3) Data with a deletion rate greater than 10% of the 10 Hz raw data every 30 minutes was eliminated; (4) Observation data of weak turbulence at night (u* less than 0.1 m/s) was excluded. The average period of observation data is 30 minutes, 48 data per day, and the missing data is marked as -6999. The data of April 3 and 4 was missing due to Li7500 calibration of the eddy system; data from August 29 to September 5 was missing due to collector problem. Published observation data include: Date/Time, wind direction(°), horizontal wind speed(m/s), lateral wind speed standard deviation(m/s), ultrasonic virtual temperature (°C), water vapor density (g/m3), carbon dioxide concentration(mg/m3), friction velocity (m/s), length (m), sensible heat flux(W/m2), latent heat flux (W/m2), carbon dioxide flux (mg/(m2s)), sensible heat flux quality identification QA_Hs, latent heat flux quality identification QA_LE, carbon dioxide flux quality identification QA_Fc. The quality identification of sensible heat, latent heat, and carbon dioxide flux is divided into three levels (quality mark 0: (Δst <30, ITC<30); 1: (Δst <100, ITC<100); the rest is 2). The meaning of the data time, such as 0:30 represents an average data of 0:00-0:30; the data is stored in *.xls format. For hydro-meteorological network or station information, please refer to Liu et al. (2018). For observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
The data set contains the observation data of the eddy covariance system of Sidaoqiao superstation which is located along the lower reaches of the Heihe Hydrometeorological observation network, and the data set covers data from January 1, 2017 to December 31, 2017. The station is located in Sidao Bridge, Ejina Banner, Inner Mongolia, and the underlying surface is Tamarix. The latitude and longitude of the observation station is 101.1374E, 42.0012N, and the altitude is 873 m. The height of the eddy covariance system is 8 meters, the sampling frequency is 10Hz, the ultrasonic orientation is positive north, and the distance between the ultrasonic wind speed and temperature monitor (CSAT3) and the CO2/H2O analyzer (Li7500) is 15cm. The original observation data of the eddy covariance system is 10 Hz, and the released data is a 30-minute data processed by Eddypro software. The main steps of the processing include: outlier eliminating, delay time correction, coordinates rotation (secondary coordinates rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction, etc. Meanwhile, the quality evaluation of each flux value was performed,mainly includes atmospheric stability (Δst) test and turbulence similarity (ITC) test. The 30-min flux value output of Eddypro software was also screened: (1) Data from the instrument error was eliminated; (2) Data obtained with one hour before and after precipitation was removed; (3) Data with a deletion rate greater than 10% of the 10 Hz raw data every 30 minutes was eliminated; (4) Observation data of weak turbulence at night (u* less than 0.1 m/s) was excluded. The average period of observation data is 30 minutes, 48 data per day, and the missing data is marked as -6999. The data was missing due to Li7500 calibration of the eddy system on April 7 and 8; the suspicious data caused by instrument drift and other reasons was marked by red fonts. Published observation data include: date/time Date/Time, wind direction(°), horizontal wind speed(m/s), lateral wind speed standard deviation(m/s), ultrasonic virtual temperature (°C), water vapor density (g/m3), carbon dioxide concentration(mg/m3), friction velocity (m/s), length (m), sensible heat flux(W/m2), latent heat flux (W/m2), carbon dioxide flux (mg/(m2s)), sensible heat flux quality identification QA_Hs, latent heat flux quality identification QA_LE, carbon dioxide flux quality identification QA_Fc. The quality identification of sensible heat, latent heat, and carbon dioxide flux is divided into three levels (quality mark 0: (Δst <30, ITC<30); 1: (Δst <100, ITC<100); the rest is 2). The meaning of the data time, such as 0:30 represents an average data of 0:00-0:30; the data is stored in *.xls format. For hydrometeorological network or station information, please refer to Li et al. (2013). For observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This data set contains the observation data of vortex correlativity instrument at yakou station upstream of heihe hydrometeorological observation network from January 1, 2017 to December 31, 2017.The station is located in qilian county, qinghai province.The latitude and longitude of the observation point is 100.2421, 38.0142N, and the altitude is 4148 m.The height of the vortex correlation instrument is 3.2m, the sampling frequency is 10Hz, the ultrasonic direction is due to the north, and the distance between the ultrasonic wind speed and temperature instrument (CSAT3) and the CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of vorticity correlativity is 10Hz, and the released data is the data of 30 minutes processed by Eddypro software. The main steps of its processing include: outfield value elimination, delay time correction, coordinate rotation (secondary coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened.(2) data of 1h before and after precipitation were excluded;(3) the missing rate of 10Hz original data is more than 10% every 30min;(4) the observed data of weak turbulence at night were excluded (u* less than 0.1m/s).The average period of observation data was 30 minutes, 48 data a day, and the missing data was marked as -6999.April 13-14, instrument calibration, data missing.Suspicious data caused by instrument drift shall be identified in red.The eddy current correlator will be short of electricity at night in winter, resulting in the loss of data. Observations published include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Carbon dioxide flux mass identification QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest are 2).The meaning of data time, such as 0:30 represents the average between 0:00 and 0:30;The data is stored in *.xls format. Please refer to Liu et al. (2018) for hydrometeorological network or site information, and Liu et al. (2011) for observation data processing.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This data set contains the eddy correlativity observation data of the downstream desert station of heihe hydrometeorological observation network from January 1, 2017 to December 31, 2017.The station is located in ejin banner, Inner Mongolia.The longitude and latitude of the observation point are 100.9872e, 42.1135n and 1054m above sea level.The frame height of the vortex correlativity instrument is 4.7m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.The vortex correlator's Li7500 was calibrated on April 7-8, and the tower was re-reinforced on September 18-23, with data missing.Suspicious data caused by instrument drift, etc., shall be marked in red font.When 10Hz data is missing, the missing data is replaced by the 30-minute data output of the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Li et al. (2013), and for observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
The data set contains the vortex correlativity instrument observation data of zhangye wetland station in the middle reaches of heihe hydrometeorological observation network from January 1, 2017 to December 31, 2017.The site is located in zhangye city, gansu province.The latitude and longitude of the observation point is 100.44640E, 38.97514N, and the altitude is 1460.00m.The height of the vortex correlation instrument is 5.2m, the sampling frequency is 10Hz, the ultrasonic direction is due to the north, and the distance between the ultrasonic wind speed and temperature instrument (Gill) and the CO2/H2O analyzer (Li7500A) is 25cm. The original observation data of vorticity correlativity is 10Hz, and the released data is the data of 30 minutes processed by Eddypro software. The main steps of its processing include: outfield value elimination, delay time correction, Angle correction, coordinate rotation (secondary coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened.(2) data of 1h before and after precipitation were excluded;(3) the missing rate of 10Hz original data is more than 10% every 30min;(4) the observed data of weak turbulence at night were excluded (u* less than 0.1m/s).The average period of observation data was 30 minutes, 48 data a day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked in red. From April 3 to 12, due to instrument calibration, data is missing.When 10Hz data is missing due to a memory card storage problem, the data is replaced by 30min flux data output from the collector. Observations published include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Carbon dioxide flux mass identification QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest are 2).The meaning of data time, such as 0:30 represents the average between 0:00 and 0:30;The data is stored in *.xls format. Please refer to Li et al. (2013) for hydrometeorological network or site information, and Liu et al. (2011) for observation data processing.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This data set contains the eddy correlation-meter observation data from January 1, 2017 to December 31, 2017 at the upper reaches of the heihe hydrometeorological observation network.The station is located in qilian county, qinghai province.The longitude and latitude of the observation point are 98.9406e, 38.8399N and 3739 m above sea level.The frame of the vortex correlator is 4.5m high, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500, replaced with Li7500RS in April 2017) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift, etc., shall be marked in red font.The eddy current system Li7500 was calibrated from April 13 to 15, and the collector's data storage problem occurred from July 8 to 12, resulting in missing data.When 10Hz data is missing due to a problem with the memory card storage data, the data is replaced by the 30min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al. (2018), and for observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This data set contains the eddy correlativity observation data from January 1, 2017 to December 31, 2017 at the super station at the upper reaches of heihe hydrometeorological observation network.The station is located in caoban village, aru township, qilian county, qinghai province.The longitude and latitude of the observation point are 100.4643e, 38.0473n and 3033m above sea level.The rack height of the vortex correlativity meter is 3.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift and other reasons are marked with red font, in which the calibration data of the vortex system Li7500A from April 13 to April 14 is missing;When 10Hz data is missing due to a problem with the storage card (2.17-2.23, 3.3-4.12), the data will be replaced by the 30-min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Liu et al. (2018), and for observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
The data set contains the flux observation data of large aperture scintillator from daman station in the middle reaches of heihe hydrometeorological observation network.Large aperture scintillators of BLS450 and BLS900 models were installed at daman station in the middle reaches of China. The north tower was the receiving end of BLS900 and the transmitting end of BLS450, and the south tower was the transmitting end and the receiving end of BLS900.The observation period is from January 1, 2017 to December 31, 2017.The station is located in dazman irrigation district, zhangye city, gansu province.The latitude and longitude of the north tower is 100.3785 E, 38.8607 N, and the latitude and longitude of the south tower is 100.3685 E, 38.8468 N, with an altitude of about 1556m.The effective height of the large aperture scintillator is 22.45m, the optical diameter length is 1854m, and the sampling frequency is 1min. Large aperture flicker meter raw observation data for 1 min, data released for after processing and quality control of data, including sensible heat flux is mainly combined with the automatic meteorological station observation data, based on similarity theory alonzo mourning - Mr. Hoff is obtained by iterative calculation, the quality control of the main steps include: (1) excluding Cn2 reach saturation data (Cn2 e-13 > 1.43);(2) data with weak demodulation signal strength (Average X Intensity<1000) were eliminated;(3) data at the time of precipitation were excluded;(4) data of weak turbulence under stable conditions were excluded (u* < 0.1m/s).In the iterative calculation process, the stability universal function of Thiermann and Grassl(1992) was selected. Please refer to Liu et al(2011, 2013) for detailed introduction.Due to instrument failure, data of large aperture scintillator was missing from June 6 to July 2, 2017. Some notes on the released data :(1) the middle LAS data is mainly BLS900, the missing time is supplemented by BLS450 observation, and the missing time of both is marked with -6999.(2) data table head: Date/Time: Date/Time (format: yyyy/m/d h:mm), Cn2: structural parameters of air refraction index (unit: m-2/3), H_LAS: sensible heat flux (unit: W/m2).The meaning of data time, such as 0:30 represents the average between 0:00 and 0:30;The data is stored in *.xls format. Please refer to Li et al. (2013) for hydrometeorological network or site information, and Liu et al. (2011) for observation data processing.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
The data set contains the flux observation data of scintillator with large aperture from sidaoqiao station downstream of heihe hydrometeorological observation network.A large aperture scintillator of BLS900 type is installed in the downstream. The north tower is the receiving end and the south tower is the transmitting end.The observation period is from January 1, 2017 to December 31, 2017.The site is located in ejin banner, Inner Mongolia, with tamarix chinensis, populus populus, bare land and cultivated land under it.The latitude and longitude of the north tower is 101.137e, 42.008n, and the latitude and longitude of the south tower is 101.131e, 41.987 N, with an elevation of about 873m.The effective height of the large aperture scintillator is 25.5m, the optical diameter length is 2350m and the sampling frequency is 1min. Large aperture flicker meter raw observation data for 1 min, data released for after processing and quality control of data, including sensible heat flux is mainly combined with the automatic meteorological station observation data, based on similarity theory alonzo mourning - Mr. Hoff is obtained by iterative calculation, the quality control of the main steps include: (1) excluding Cn2 reach saturation data (e-14 Cn2 > 7.58);(2) data with weak demodulation signal strength (Average X Intensity<1000) were eliminated;(3) data at the time of precipitation were excluded;(4) data of weak turbulence under stable conditions were excluded (u* < 0.1m/s).During the iterative calculation, the stability universal function of Thiermann and Grassl(1992) was selected.Please refer to Liu et al(2011, 2013) for detailed introduction.Due to the problem of data storage unit, data of large aperture scintillator was missing from February 21 to March 5, and July 10 to August 18, 2017. A few notes on published data :(1) data missing time is marked by -6999.(2) data table head: Date/Time: Date/Time (format: yyyy/m/d h:mm), Cn2: structural parameters of air refraction index (unit: m-2/3), H_LAS: sensible heat flux (unit: W/m2).The meaning of data time, such as 0:30 represents the average between 0:00 and 0:30;The data is stored in *.xls format. Please refer to Li et al. (2013) for hydrometeorological network or site information, and Liu et al. (2011) for observation data processing.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
The data set contains the flux observation data of large aperture scintillator at areau station upstream of heihe hydrometeorological observation network.Two large aperture scintillation devices of BLS450 and zzlas type were set up in the upstream areau station respectively. The north tower was the receiving end of zzlas and the transmitting end of BLS450, and the south tower was the transmitting end of zzlas and the receiving end of BLS450.The observation time is January 1, 2017, solstice, December 31, 2017.The station is located in the grass daban village, a soft township, qilian county, qinghai province.The latitude and longitude of the north tower is 100.4712e, 38.0568n, and the latitude and longitude of the south tower is 100.4572e, 38.0384 N, with an altitude of about 3033m.The effective height of the large aperture scintillator is 9.5m, the optical diameter length is 2390m, and the sampling frequency is 1min. Large aperture flicker meter raw observation data for 1 min, data released for after processing and quality control of data, including sensible heat flux is mainly combined with the automatic meteorological station observation data, based on similarity theory alonzo mourning - Mr. Hoff is obtained by iterative calculation, the quality control of the main steps include: (1) excluding Cn2 reach saturation data (BLS450: Cn2 > 7.25 e-14, zzlas: Cn2 > 7.84 E - 14).(2) data with weak demodulation signal strength (BLS450: Mininum X Intensity <50) were eliminated;Zzlas: Demod>-20mv);(3) data at the time of precipitation were excluded;(4) data of weak turbulence under stable conditions were excluded (u* < 0.1m/s).In the iterative calculation process, for BLS450, Thiermann and Grassl(1992) stability universal function was selected.For zzlas, select Andreas 1988's stability universal function.Please refer to Liu et al(2011, 2013) for detailed introduction.From April 16 to May 26, 2017, the measurement signal of large aperture scintillator was relatively small, resulting in a large number of missing data. Several notes on the released data :(1) the upstream LAS data is mainly BLS450, the missing time is supplemented by zzlas observation, and the missing time of both is marked by -6999.(2) data table head: Date/Time: Date/Time (format: yyyy/m/d h:mm), Cn2: structural parameters of air refraction index (unit: m-2/3), H_LAS: sensible heat flux (unit: W/m2).The meaning of data time, such as 0:30 represents the average between 0:00 and 0:30;The data is stored in *.xls format, please refer to the references for details. Please refer to Li et al. (2013) for hydrometeorological network or site information, and Liu et al. (2011) for observation data processing.
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This data set contains the eddy correlation-meter observation data of the mixed forest station downstream of heihe hydrometeorological observation network from January 1, 2017 to December 31, 2017.The station is located in Inner Mongolia ejin banner four road bridge, under the surface is populus and tamarix.The longitude and latitude of the observation point are 101.1335e, 41.9903n and 874 m above sea level.The rack height of the vortex correlativity instrument is 22m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500) is 17cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.2m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.April 7 solstice April 8 due to instrument calibration, 3.24-4.08 infrared gas analyzer error, data missing.Suspicious data caused by instrument drift, etc., are identified in red font.When 10Hz data is missing due to a problem with the memory card storage data, the data is replaced by the 30min flux data output by the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Li et al. (2013), and for observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This data set contains the eddy correlativity observation data of huachaizi desert station in the middle reaches of heihe hydrological meteorological observation network from January 1, 2017 to December 31, 2017.The station is located in zhangye city, gansu province.The longitude and latitude of the observation point are 100.3201E, 38.7659N and 1731.00m above sea level.The rack height of the vortex correlator is 4.5m, the sampling frequency is 10Hz, the ultrasonic orientation is due north, and the distance between the ultrasonic wind speed and temperature meter (CSAT3) and CO2/H2O analyzer (Li7500A) is 15cm. The original observation data of the vortex correlativity instrument is 10Hz, and the published data is the 30-minute data processed by Eddypro software. The main processing steps include: outliers, delay time correction, coordinate rotation (quadratic coordinate rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction.Quality assessment for each intercompared to at the same time, mainly is the atmospheric stability (Δ st) and turbulent characteristics of similarity (ITC) test.The 30min pass value output by Eddypro software was also screened :(1) data when instrument error was eliminated;(2) data of 1h before and after precipitation are excluded;(3) remove the data with a missing rate of more than 10% in the original 10Hz data within every 30 minutes;(4) the observation data of weak turbulence at night (u* less than 0.1m/s) were excluded.The average observation period was 30 minutes, 48 data per day, and the missing data was marked as -6999.Suspicious data caused by instrument drift, etc., shall be marked in red font.April 3 solstice on April 4, due to the calibration of vortex correlator Li7500A, data was missing.When the 10Hz data of the vortex correlator is missing (1.1-1.22), the data will be filled by the 30-minute data output of the collector. The published observational data include:Date/Time for the Date/Time, wind Wdir (°), Wnd horizontal wind speed (m/s), standard deviation Std_Uy lateral wind speed (m/s), ultrasonic virtual temperature Tv (℃), the water vapor density H2O (g/m3), carbon dioxide concentration CO2 (mg/m3), friction velocity Ustar) (m/s), Mr. Hoff length L (m), sensible heat flux Hs (W/m2), latent heat flux LE (W/m2), carbon dioxide flux Fc (mg/(m2s)), the quality of the sensible heat flux identifier QA_Hs, the quality of the latent heat flux identifier QA_LE,Quality indicator for co2 flux QA_Fc.The quality of the sensible heat and latent heat, carbon dioxide flux identification is divided into three (quality id 0: (Δ st < 30, the ITC < 30);1: (Δ st < 100, ITC < 100);The rest is 2).The meaning of data time, such as 0:30 represents the average of 0:00-0:30;The data is stored in *.xls format. For information of hydrometeorological network or station, please refer to Li et al. (2013), and for observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
Our project entrust the L band radiosonde sounding encrypt observations to Zhangye National Climate Observatory, and collect regular observation twice a day. The dataset contains three times one day at 8:00, 14:00, 20:00, which can support the remote sensing image atmospheric correction and atmospheric science research. Observation Site: Zhangye National Climate Observatory located in Shajing Town, west of ZhangYe. The coordinates of this site: 39°5′15.68" N, 100°16′39.11" E。 Observation Instrument: China Meteorological Administration Operational L Band radiosonde system. Observation Time: The observation date last from 1 May, 2012 to 31 September, 2012, among which: Three times observations at 7:00-8:00, 13:00-14:00 and 19:00-20:00 during 1 June, 2012 to 31 August, 2012; twice at 7:00-8:00 and 19:00-20:00 during 2012-5-1 to 5-31 and 2012-9-1 to 9-31. Accessory data： Pressure, temperature, relative humidity, wind speed and wind direction profiles data.
MA Mingguo
This dataset contains the flux measurements from the large aperture scintillometer (LAS) at A’rou Superstation in the hydrometeorological observation network of Heihe River Basin between 14 October, 2012, and 31 December, 2013. There were two types of LASs at A’rou Superstation: German BLS450 and China zzlas. The north tower was set up with the zzlas receiver and the BLS450 transmitter, and the south tower was equipped with the zzlas transmitter and the BLS450 receiver. Zzlas has been in use since 14 October, 2012, and the observation period of BLS450 was from 9 August to 10 December, 2013. The site (north: 100.467° E, 38.050° N; south: 100.450° E, 38.033° N) was located in Caodaban village of A’rou town in Qilian county, Qinghai Province. The underlying surface between the two towers was alpine meadow. The elevation is 3033 m. The effective height of the LASs was 9.5 m, and the path length was 2390 m. The data were sampled at 5 Hz and 1 Hz intervals for BLS450 and zzlas, respectively, and then averaged over 1 min. The raw data acquired at 1 min intervals were processed and quality controlled. The data were subsequently averaged over 30 min periods, in which sensible heat flux was iteratively calculated by combining Cn2 with meteorological data according to the Monin-Obukhov similarity theory. The main quality control steps were as follows: (1) The data were rejected when Cn2 exceeded the saturated criterion (BLS450: Cn2>7.25E-14, zzlas: Cn2>7.84E-14). (2) The data were rejected when the demodulation signal was small (BLS450: Average X Intensity<1000; zzlas: Demod>-20 mv). (3) The data were rejected when collected during precipitation. (4) The data were rejected if collected at night when weak turbulence occurred (u* was less than 0.1 m/s). In the iteration process, the universal functions of Thiermann and Grassl, 1992 and Andreas, 1988 were selected for BLS450 and zzlas, respectively. Several instructions were included with the released data. (1) The data were primarily obtained from BLS450 measurements, and missing flux measurements from the BLS450 instrument were substituted with measurements from the zzlas instrument. The missing data were denoted by -6999. Due to the drift of the zzlas signal, data from 10 November to 23 November, 2012, and 14 March to 10 April, 2013, were excluded. Due to the LAS tower’s lean, the data from 10 April to 31 May, 2013, were not collected. (2) The dataset contained the following variables: data/time (yyyy-m-d h:mm), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). In this dataset, a time of 0:30 corresponds to the average data for the period between 0:00 and 0:30, and the data were stored in *.xls format. Moreover, suspicious data were marked in red. For more information, please refer to Li et al. (2013) (for hydrometeorological observation network or sites information), Liu et al. (2011) (for data processing) in the Citation section.
LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
