SCIENTIFIC DATA

<figure><img src="image_1.png" /></figure>

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# Multiyear in-situ L-band microwave radiometry of land surface processes on the Tibetan Plateau

Z. Su<sup>1,9</sup>✉, J. Wen<sup>2</sup>✉, Y. Zeng<sup>1</sup>, H. Zhao<sup>1</sup>, S. Lv<sup>2</sup>, R. van der Velde<sup>1</sup>, D. Zheng<sup>3</sup>, X. Wang<sup>4</sup>, Z. Wang<sup>4</sup>, M. Schwank<sup>5,6</sup>, Y. Kerr<sup>7</sup>, S. Yueh<sup>8</sup>, A. Colliander<sup>8</sup>, H. Qian<sup>9</sup>, M. Drusch<sup>10</sup> & S. Mecklenburg<sup>11</sup>

We report a unique multiyear L-band microwave radiometry dataset collected at the Maqu site on the eastern Tibetan Plateau and demonstrate its utilities in advancing our understandings of microwave observations of land surface processes. The presented dataset contains measurements of L-band brightness temperature by an ELBARA-III microwave radiometer in horizontal and vertical polarization, profile soil moisture and soil temperature, turbulent heat fluxes, and meteorological data from the beginning of 2016 till August 2019, while the experiment is still continuing. Auxiliary vegetation and soil texture information collected in dedicated campaigns are also reported. This dataset can be used to validate the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite based observations and retrievals, verify radiative transfer model assumptions and validate land surface model and reanalysis outputs, retrieve soil properties, as well as to quantify land-atmosphere exchanges of energy, water and carbon and help to reduce discrepancies and uncertainties in current Earth System Models (ESM) parameterizations. Measurement cases in winter, pre-monsoon, monsoon and post-monsoon periods are presented.

## Background & Summary

Microwave remote sensing of land surfaces on a global scale has mainly focused on soil moisture retrieval in the recent past. This is because soil moisture strongly influences hydrological and agricultural processes in controlling runoff generation, drought development, and agroecosystem functioning. As a source of water for evaporation and transpiration over land surfaces, soil moisture is involved in the water, energy and carbon cycles of the Earth system and impacts on the climate system through atmospheric feedbacks. Due to its significance, soil moisture was recognized as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS) in 2010, and several international programs have been established in recent years to produce global soil moisture data. They include the Climate Change Initiative (CCI) of the European Space Agency (ESA), the Soil Moisture and Ocean Salinity (SMOS<fnref n="1" />) satellite mission of ESA, and the Soil Moisture Active Passive (SMAP²,³) satellite mission of the National Aeronautics and Space Administration (NASA). The two dedicated satellite missions for observation of soil moisture SMOS and SMAP provide global soil moisture products at nearly daily temporal resolution and coarse spatial resolution (e.g., 15/25 km for SMOS and 36 km for SMAP).

Despite these advances, the operational retrieval algorithms have relied on zeroth-order radiative transfer theory (the so-called τ-ω model) and empirical assumptions in passive microwave retrievals in the past. Although

<footnote n="1">¹Faculty of Geo-information Science and Earth Observation (ITC), University of Twente, Enschede, The Netherlands.
²College of Atmospheric Sciences, Plateau Atmosphere and Environment Key Laboratory of Sichuan Province, Chengdu University of Information Technology, Chengdu, China.
³Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, China.
⁴Key Laboratory of Land Surface Process and Climate Change in Cold and Arid Regions, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, China.
⁵Swiss Federal Research Institute WSL, Birmensdorf, Switzerland.
⁶Gamma Remote Sensing AG, Gümligen, Switzerland.
⁷CESBIO (CNES/CNRS/UPS/IRD), Toulouse, France.
⁸Jet Propulsion Laboratory, Pasadena, USA.
⁹Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, School of Water and Environment, Chang'an University, Xi'an, 710054, China.
¹⁰European Space Agency, ESTEC, Earth Observation Programmes, Noordwijk, The Netherlands.
¹¹European Space Agency, ESA Climate Office, Harwell Campus, Oxfordshire, UK.
✉e-mail: z.su@utwente.nl; jwen@cuit.edu.cn</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** A schematic overview of the ELBARA-III tower setup (top panel) and the footprints (bottom panel). The footprints vary from 3.31 m² to 43.64 m² for incidence angle from 40° to 70°. The half-axes of the elliptic footprint are indicated as *a* and *b* for a given incidence angle θᵢ (top panel) and the projected ground distances from the radiometer to the closest- and the farthest-side of the elliptic footprints at –3 dB sensitivity of the antenna are indicated as *d*min and *d*max (bottom panel). The locations of the installed *in situ* soil moisture and soil temperature sensors are indicated as SMST_Z and SMST_LC. The fence (25 m × 45 m) is not drawn to scale.</figcaption></figure>

more sophisticated approaches are currently being developed to account for multiple scattering, in particular over forested areas⁴⁻⁷, there is still little consensus about the effects of vegetation interception and litter on observed brightness temperature, which remains difficult to be addressed by studies based on space-borne observations from passive microwave satellite data⁸. As an example, soil moisture retrieval using both SMAP vertical and horizontal polarization does not outperform the use of single vertical polarization⁹. Similarly, the precise nature of the vegetation scattering and emission and its representation in active microwave retrievals have remained unresolved, resulting in the current largely empirical approaches and large uncertainties and inconsistencies among different operational soil moisture products¹⁰,¹¹. As a further complication, when a land surface undergoes freeze-thaw processes, the behaviour of microwave observation abruptly changes in response to changes in the phase of the soil water (i.e., liquid or solid phase) at different soil depths. While such a dynamic process in space and time can be observed and modelled with *in situ* measurements¹²⁻¹⁴, current satellite retrievals can only provide freeze-thaw information (date and range of depths) using passive microwaves at a low resolution¹⁵,¹⁶ or a binary indication of the frozen or thawed surfaces¹⁷. All these results point to a fundamental need to advance knowledge in understanding the precise scattering-emission mechanism of vegetated lands and the need for in-depth investigations of freeze-thaw processes.

To contribute to filling this knowledge gap, an L-band microwave observation system (ELBARA-III; See Fig. 1) was set up at the Maqu site of the Tibetan soil moisture and soil temperature observation network¹⁰ for long-term observation of the land surface processes since the beginning of 2016. ELBARA-III is known as the third generation ESA L-Band Radiometer. The collected dataset contains measurements of L-band brightness temperature in horizontal and vertical polarization, profile soil moisture and soil temperature, turbulent heat fluxes, and meteorological data. Auxiliary vegetation and soil texture information were also collected during intensive field campaigns. This dataset can be used to validate SMOS and SMAP brightness temperature observations and retrievals, verify radiative transfer model assumptions and validate land surface model and reanalysis outputs, retrieve soil properties, as well as to quantify land-atmosphere exchanges of energy, water and carbon that are essential in realistic Earth system modeling. Figure 1 provides a schematic overview of the ELBARA-III radiometry setup.

## Methods

The dataset is collected at the Maqu regional-scale soil moisture and soil temperature (SMST) monitoring network on the north-eastern Tibetan Plateau¹⁰,¹¹ and it mainly contains measurements of L-band brightness temperature in horizontal and vertical polarization by an ELBARA-III microwave radiometer, profile soil moisture and soil temperature by 5TM sensors, turbulent heat fluxes by a CSAT eddy covariance system, and meteorological data by an automatic weather station. Auxiliary vegetation and soil texture information were also collected during intensive field campaigns.

The Maqu SMST monitoring network (33°30′–34°15′ N, 101°38′–102°45′ E) is located in the source region of the Yellow River on the north-eastern part of the Tibetan Plateau at an altitude between 3200 m and 4200 m above mean sea level. The Maqu area has a cold climate with dry winter and warm summer (Dwb) in the updated Köppen-Geiger climate classification¹⁹. Land cover is mainly alpine meadows with grass heights varying from 5 to 15 cm throughout the growing season due to intensive grazing by livestock. The network is equipped with 20 profile measurements of SMST distributed over an area of 40 km by 80 km. The ELBARA-III radiometer is
installed at the center of the SMST monitoring network (hereafter Maqu site) to collect continuous microwave radiometric signatures of the grassland site. The prevailing soil types are sandy loam, silt loam, and organic soil with an average ~30.3% sand, ~9.9% clay, and a maximum of ~39.0% organic matter based on soil sampling²⁰.

**L-band brightness temperature.** ELBARA-III is an L-band (1.4 GHz) Dicke-type radiometer with a dual-polarized conical horn antenna with –3 dB beam width of 12°. It uses a resistive load (RL), and an active cold load (ACL) as internal calibration sources to derive calibrated brightness temperature ($T_b^p$) of the ground footprints in horizontal ($p=H$) and vertical polarization ($p=V$). The 50 Ω RL is kept at a stabilized instrument internal temperature $T_{inst}$ to better than ±0.1 K to ensure the noise temperature $T_{RL} = T_{inst}$. The ACL is a low-noise amplifier with its noise temperature $T_{ACL}$ calibrated by cold sky measurements. To mitigate and detect potential radio frequency interference (RFI), the radiometric signal is split into two sub-bands, one in 1.402–1.413 GHz and the other in 1.414–1.425 GHz, within the protected 1.400–1.427 GHz of the microwave L-band (1–2 GHz). The absolute accuracy of ELBARA-III $T_b^p$ measurement is better than 1 K and the corresponding sensitivity is at least 0.1 K. ELBARA-III has upgraded components and functions compared to ELBARA-II²¹ including a new temperature-controlled Radiometer Microwave Assembly (RMA), a new thermoelectric cooling (TEC) controller for the Peltier elements, a new weather resistant Instrument Computer (IC) enclosure, a new detector assembly, longer antenna feed cables (1 m), addition of a 6 dB attenuator at the input of the sub-band filters to improve isolation between the filters which reduces the leakage from one filter to the other by 12 dB, and a bootloader of the TEC allowing upgrading firmware from the IC.

The ELBARA-III measurement cycle consists of successive measurements of the ambient, hot and cold loads followed by the H and V antenna polarization channels. At each position the total power of the radiometer is measured for a period of approximately 4 seconds. The measurements of the ambient and cold load establish the calibration line for the radiometer. The slope of the line is the radiometer gain and the intercept is related to the noise generated by the radiometer itself. Once the internal noise temperatures are calibrated the brightness temperature $T_b^p$ at polarization $p=H, V$ can be computed from the measured radiometer voltages. The raw data $U_{in}^p$ is measured when the input switch is on the antenna H or V polarization position, which needs to be corrected with the unavoidable losses from the feed cable. The noise temperature at the radiometer input can be computed as $T_{b,in}^p = (T_{RL} - T_{ACL}) (U_{in}^p - U_{ACL}) / (U_{RL} - U_{ACL}) + T_{ACL}$, where $U_{RL}, U_{ACL}$ are the voltages measured when the input switch is set to the resistive load and the active cold load, and $T_{RL}$ and $T_{ACL}$ are noise temperatures related to calibration sources and are internally calculated. After considering the contribution of the feed cable loss, the brightness temperature is calculated from the raw data as $T_b^p = [T_{b,in}^p - (1 - \tau_{fc})T_{fc}] / \tau_{fc}$, where $\tau_{fc}$ is the transmissivity of the feed cable and $T_{fc}$ equals to the ambient temperature of the feed cable.

The ELBARA-III radiometer is installed at Maqu site on a 4.8 m high scaffold tower, making the center of rotation at 6.5 m in height above ground with the antenna pointing towards the south. The continuous measurements include angular scanning of the ground at every 30 min and sky measurements. The angular range of the ground scans is between 40°–70° (relative to nadir) in steps of 5°, and sky measurements are performed at 23:55 local time every day with an observation angle of 155° (relative to nadir). The elliptic footprints at –3 dB sensitivity of the antenna are estimated according to the installation height and the incidence angle²² (Fig. 1).

**Profile soil moisture and soil temperature.** SMST profiles are automatically measured by 5TM ECH₂O probes (METER Group, Inc. USA) next to the radiometer tower at 15-min time-intervals installed at the depths for SMST-Z profile¹³: 5, 10, 20, 40, 80 and 160 cm (one sensor per depth) and for SMST-LC profile²³: 2.5 cm (2 sensors), then one sensor every 2.5 cm in the top 20 cm, every 5 cm between 20–50 cm and every 10 cm between 50–100 cm (19 layers and 20 sensors in total). The 5TM probe is a capacitance sensor operating at 75 MHz measuring the dielectric permittivity of the surrounding soil, and the measured dielectric permittivity values are converted to volumetric liquid water contents using the Topp equation²⁴. Specific calibrations were conducted for the profile soil textures¹⁸,²⁰.

**Turbulent heat fluxes.** The micro-meteorological observations are conducted by an eddy-covariance (EC, CSAT3) system which is installed for measuring the turbulent heat fluxes.

**Meteorological data.** The meteorological observation system consists of an automatic weather station which provides wind speed and direction, air humidity and temperature measurements at 2.5 m above ground, and instrumentations for measuring four radiation components (i.e., up- and downwelling shortwave and long-wave radiations), as well as air pressure and liquid precipitation (Geonor T-200B Series). The albedo is calculated as the ratio of up- and downwelling shortwave radiations, and the surface temperature is derived from the up- and downwelling longwave radiations.

**Soil and vegetation data.** Soil samples are collected around the ELBARA-III field site for laboratory analyses to quantify the soil hydraulic and thermal properties, i.e., soil texture (sand, clay, and silt), organic matter content, bulk density, porosity, soil water potential at air-entry and empirical parameters related to the pore-size distribution of the soil matrix. Field measurements of the saturated hydraulic conductivity are also carried out. The detailed analysis and measured hydraulic properties are given in¹³,²⁰.

The MODIS leaf area index (LAI) product derived from data by the Terra and Aqua satellites (MCD15A2H²⁵) is extracted to represent the vegetation status. The time series of LAI is processed with the harmonic analysis of the time series (HANTS) algorithm²⁶ to remove cloud contaminations. Additional field measurements of fresh and dry above-ground biomass, LAI, and vegetation height were conducted during a field campaign in 2018²⁷.
<table><thead><tr><th rowspan="2">Data Type</th><th colspan="6">2016-2017 period</th></tr><tr><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th></tr></thead><tbody><tr><td>Meteo data*</td><td>25-Mar-2016</td><td></td><td></td><td></td><td></td><td>30-Jan-2017</td></tr><tr><td>Eddy covariance</td><td></td><td></td><td>5-Jun-2016</td><td>4-Sep-2016</td><td>1-Dec-2016</td><td>29-Mar-2017</td></tr><tr><td>SMST LC</td><td></td><td></td><td>7-Aug-2016</td><td></td><td></td><td>29-Mar-2017</td></tr><tr><td>SMST_Z</td><td>1-Jan-2016</td><td>6-Apr-2016</td><td></td><td></td><td></td><td></td></tr><tr><td>ELBARA TB</td><td>1-Jan-2016</td><td>6-Apr-2016</td><td>7-Aug-2016</td><td>30-Nov-2016</td><td>1-Jan-2017</td><td>29-Mar-2017</td></tr><tr><td>In-situ LAI</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MODIS LAI</td><td colspan="6">From Jan-2016 to Dec-2019</td></tr><tr><td>SMAP L1 TB</td><td colspan="6">From 1-Jan-2016 to 31-Dec-2019</td></tr></tbody></table><table><thead><tr><th rowspan="2">Data Type</th><th colspan="6">2017-2018 period</th></tr><tr><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th></tr></thead><tbody><tr><td>Meteo data</td><td>31-Jul-2017</td><td>27-Aug-2017</td><td>22-Oct-2017</td><td>12-Aug-2018</td><td></td><td></td></tr><tr><td>Eddy covariance</td><td>29-Mar-2017</td><td>16-Nov-2017</td><td>9-Dec-2017</td><td>12-Aug-2018</td><td></td><td></td></tr><tr><td>SMST LC</td><td>27-Jul-2017</td><td></td><td></td><td>12-Aug-2018</td><td></td><td></td></tr><tr><td>SMST_Z</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ELBARA TB</td><td>29-Mar-2017</td><td></td><td></td><td>12-Aug-2018</td><td></td><td></td></tr><tr><td>In-situ LAI</td><td>12-Jul-2018, 17-Aug-2018</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MODIS LAI</td><td colspan="6">From Jan-2016 to Dec-2019</td></tr><tr><td>SMAP L1 TB</td><td colspan="6">From 1-Jan-2016 to 31-Dec-2019</td></tr></tbody></table><table><thead><tr><th rowspan="2">Data Type</th><th colspan="6">2018-2019 period</th></tr><tr><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th><th>StartT</th><th>EndT</th></tr></thead><tbody><tr><td>Meteo data</td><td>12-Aug-2018</td><td>30-Oct-2018</td><td>11-Nov-2018</td><td>28-Aug-2019</td><td></td><td></td></tr><tr><td>Eddy covariance</td><td>12-Aug-2018</td><td>10-Nov-2018</td><td>26-Mar-2019</td><td>28-Aug-2019</td><td></td><td></td></tr><tr><td>SMST LC</td><td>15-Aug-2018</td><td></td><td></td><td>31-May-2019</td><td></td><td></td></tr><tr><td>SMST_Z</td><td>—</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ELBARA TB</td><td>12-Aug-2018</td><td>29-Dec-2018</td><td>25-Mar-2019</td><td>28-Aug-2019</td><td></td><td></td></tr><tr><td>In-situ LAI</td><td>12-Jul-2018, 17-Aug-2018</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MODIS LAI</td><td colspan="6">From Jan-2016 to Dec-2019</td></tr><tr><td>SMAP L1 TB</td><td colspan="6">From 1-Jan-2016 to 31-Dec-2019</td></tr></tbody></table>

**Table 1.** Overview of data availability. (StartT: start of a data period, EndT: end of a data period; Meteo data: meteorological data; Eddy covariance: micrometeorological data; SMST LC: Soil moisture and soil temperature data at LC location²³, SMST_Z: Soil moisture and soil temperature at Z location¹³; ELBARA TB: ELBARA-III brightness temperature; In-situ LAI: in-situ measured leaf area index²⁷; MODIS LAI: Leaf area index from the MODIS sensor²⁵; SMAP L1 TB: SMAP Level 1 brightness temperature³). (*The precipitation data from 1-Sep-2016 to 30-Jan-2017 was provided by the Zoige Plateau Wetlands Ecosystem Research Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Science, Lanzhou, China; The precipitation gauge is the same type Geonor T-200B Series installed 200 m to the north-west of Maqu site).

## Data Records

We present here the collected Maqu ELBARA-III radiometry dataset²⁸,²⁹ in detail including L-band brightness temperature, profile soil moisture and soil temperature, turbulent heat fluxes, meteorological data, as well as soil and vegetation data. Online-only Table 1 summarizes the general characteristics of the Maqu site, the included variables in the data, the used instruments and their setups. Table 1 presents the availability of the collected data by different instrumentations.

The dataset as described in Table 1 can be accessed at https://figshare.com. Detailed technical description of the data records can be found in the readme.txt files and codes and procedures for processing and plotting the figures, and for downloading satellite data are also included. The contents of the two figshare data records²⁸,²⁹ are described in Online-only Table 2 and Table 2, respectively.

## Technical Validation

We present here analyses to support the technical quality of the Maqu microwave radiometry dataset. While the used instrumentation and the data collection have been presented in the methods section, this section focuses on the consistency of the different variables that enable further exploration and application of the data.

The quality assurance and quality control (QA/QC) of ELBARA-III measurements are carried out through the 'histogram test' on the voltage samples (raw-data) of the detector output at sampling frequency of 800 Hz. The 'histogram test' is an end-to-end test of the radiometer. Specifically, statistics of the histogram test can show if there is an internal non-Gaussian Radio Frequency Interference (RFI) via excess kurtosis and standard deviation. Furthermore, serious imbalances between the two receiver channels (for each polarization) can be identified via
<table><caption>Table 2. Codes for filtering brightness temperature (TB) outliers and TB data with corrected local time²⁹.</caption><thead><tr><th>Directory</th><th>Subdirectory/File name</th><th>Content</th></tr></thead><tbody><tr><td>Code_update</td><td></td><td>Code for filtering TB outlier</td></tr><tr><td></td><td>TB_FilteringByQuantile.py</td><td>Code for quantile filtering</td></tr><tr><td></td><td>TB_FilteringByHANTS.py</td><td>Code for filtering TB outlier by using HANTS</td></tr><tr><td></td><td>Plot_DailyScale_20170701.py</td><td>Code for plotting Fig. 2</td></tr><tr><td></td><td>AWS_TB_SMSTDisplay_updated.py</td><td>Code for plotting Fig. 4</td></tr><tr><td>Data_update</td><td></td><td>TB data with corrected local time and explanations</td></tr><tr><td></td><td>Data_version1.1</td><td></td></tr><tr><td></td><td>ELBARA-III dataset-2018-2019ELBARA-III TB.csv</td><td>ELBARA TB data, 2018-2019</td></tr><tr><td></td><td>ELBARA-III dataset-2017-2018ELBARA-III TB.csv</td><td>ELBARA TB data, 2017-2018</td></tr><tr><td></td><td>ELBARA-III dataset-2016-2017ELBARA-III TB.csv</td><td>ELBARA TB data, 2016-2017</td></tr><tr><td></td><td>SMAP_ATBD_TimeInfor.pdf</td><td>SMAP Time information</td></tr><tr><td></td><td>README_Version1.1.txt</td><td>README file for version 1.1</td></tr><tr><td></td><td>Note on Filtering brightness temperature caused by solar reflection_v1.2.docx</td><td>Note on filtering TB data</td></tr><tr><td></td><td>Maqu_NoonTime.dat</td><td>Local noon time at Maqu</td></tr></tbody></table>

mean values. Another QA/QC test is via the 'sky-looking' measurement. If the voltages at the antenna ports measured during sky measurements are close to each other, it indicates that the losses in the two antenna cables (i.e., H and V polarization) are almost the same (meaning stable ELBARA-III operations). Other indicators used for QA/QC purposes include the instrument internal temperature, active cold source temperature, ambient temperature, and the angular behaviour of the processed brightness temperatures. These QA/QC procedures are carried out routinely after data downloading and have indicated excellent performance of the ELBARA-III radiometer. There is no potential (non-Gaussian) RFI detected at the Maqu site during the ELBARA-III operation and the ELBARA-III measurements appear reliable and reflect the environmental conditions of the observed footprint areas.

Figure 2 shows an example of angular variations of the Maqu ELBARA-III radiometry dataset for 01/07/2018 (date is given in dd/mm/yyyy). The presented variables in Fig. 2(a) are ground surface temperature (TG), soil temperature at 2.5 cm depth (ST_2.5 cm), and soil moisture at 2.5 cm depth (SM_2.5 cm) (top panel) to assist the interpretation of the observed brightness temperature in horizontal and vertical polarization ($T_b^H, T_b^V$) from 40° to 70° incidence angle in combination with precipitation (Pre) (bottom panel). TG and ST_2.5 cm present a sinusoidal variation with ST_2.5 cm lagging TG by about one hour and having a smaller amplitude. SM_2.5 cm shows a gradual decrease from early morning (6:00) until the heavy precipitation just before 18:00 after which it jumps from ca. 0.26 to 0.29 m³/m³, and then to ca. 0.31 m³/m³ after a second precipitation event around 21:00. $T_b^H$ and $T_b^V$ present typical angular variations (Fig. 2(b)) with $T_b^H$ decreasing from 9:00 till 18:00 and strongly reacting to the two precipitation events, while the changes in $T_b^V$ are in general much smaller throughout the whole period. Despite the fact that the precise contributions to $T_b^H$ and $T_b^V$ from soil and vegetation emission and their interaction need to be quantified by more detailed modeling, we select a few characteristic periods from different years to illustrate the panoply of land-atmosphere conditions affecting the ELBARA-III observations at the Maqu site and provide explanations in detail in the following as examples.

Figures 3–6 present samples of seasonal variations of the Maqu ELBARA-III microwave radiometry dataset for different seasons, including pre-monsoon (late March to late June) (Fig. 3), monsoon (late June to late September) (Figs. 4–6), post-monsoon (early October to late November) (Fig. 7), and winter season (late November to late March) (Fig. 8). For each figure, three panels are plotted to display the ELBARA-III brightness temperature observations in horizontal and vertical polarization at 40° incidence angle, together with precipitation in panel 3 and other most relevant variables that can be used to explain the observed variations in the brightness temperature in panels 1 and 2. Panel 1 displays soil moisture at 2.5 cm depth and albedo, and panel 2 displays ground surface temperature (TG), air temperature (Tair), soil temperature at 2.5 cm depth (ST_2.5 cm), and the nominal freezing point as a reference (273.15 K). Both panels support to appreciate and understand the variations in the ELBARA-III observations in terms of environmental variables.

In the pre-monsoon period in 2018, the gradual increase of ($T_b^H, T_b^V$) (panel 3, the bottom panel in Fig. 3) is due to the increasing of soil temperature (panel 2, resulting in increasing effective temperature³⁰) and increase of vegetation which is indicated by the decreasing difference between $T_b^H$ and $T_b^V$. The soil moisture at 2.5 cm depth (SM_2.5 cm) changed gradually from approximately 0.23 m³/m² in late March to approximately 0.27 m³/m³ (panel 1), due to thawing of the frozen ground in late March and increasing in precipitation after mid-April. Despite the frequent precipitation events from April to June, the strong evaporation prevented further wetting of the soil profile. The strong evaporation is caused by an increase in solar radiation and a slight reduction of albedo from 0.24 to 0.23 due to the greening of the surface by vegetation. $T_b^H, T_b^V$ increased gradually from 140 K to 210 K and from 180 K to 240 K respectively (along the trend line), with some very high values (spikes) in $T_b^H$ and $T_b^V$ the reasons for which will be further explained in the following.

Monsoon period 2016 (Fig. 4): This late monsoon period sees a gradual decrease of ($T_b^H, T_b^V$) from 250 and 260 K to 210 and 225 K respectively (panel 3) mainly due to the increasing of soil moisture from 0.1 on 23/08/2016
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Angular variations of the Maqu ELBARA-III radiometry dataset for 01/07/2018 (date is given in dd/mm/yyyy). (a) plotted are ground surface temperature (TG), soil temperature at 2.5 cm depth (ST_2.5 cm), and soil moisture at 2.5 cm depth (SM_2.5 cm) (top panel), and the brightness temperature in horizontal and vertical polarization ($T_b^H, T_b^V$) from 40° to 70° incidence angle and precipitation (Pre) (bottom panel). (b) angular plot of ($T_b^H, T_b^V$) at 3 hours intervals.</figcaption></figure>

to 0.30 m³/m³ in September and October (panel 1). The role of vegetation is reduced due to senescence as indicated by the increasing difference between $T_b^H$ and $T_b^V$.

Monsoon period 2018 (Fig. 5): This early monsoon period sees a gradual increase of ($T_b^H, T_b^V$) from 225 and 240 K to 240 and 260 K respectively (panel 3) with some higher values at the end of July. Because the level of soil
<figure><img src="image_5.png" /></figure>

**Fig. 3** Seasonal variations of the Maqu ELBARA-III radiometry dataset for pre-monsoon season (late March to late June) in 2018. Plotted are soil moisture at 2.5 cm depth (SM_2.5 cm), albedo, ground surface temperature (TG), air temperature (Tair), soil temperature at 2.5 cm depth (ST_2.5 cm), the nominal freezing point as a reference (273.15 K), and the brightness temperature in horizontal and vertical polarization (T<sub>b</sub><sup>H</sup>, T<sub>b</sub><sup>V</sup>) at 40° incidence angle and precipitation (Pre). Trend lines (dashed lines) are added to SM_2.5 cm and (T<sub>b</sub><sup>H</sup>, T<sub>b</sub><sup>V</sup>) time series to assist interpretation.

<figure><img src="image_6.png" /></figure>

**Fig. 4** Same as Fig. 3 but for monsoon period (August to late September) in 2016.

moisture remained at around 0.30 m³/m³ in June, decreasing gradually in July to below 0.15 m³/m³ and increased to 0.30 after precipitation in August, the increase in (T<sub>b</sub><sup>H</sup>, T<sub>b</sub><sup>V</sup>) may be attributed to the increase of vegetation biomass. This period is characterized by the strong drying down between precipitation events.

Post-monsoon period 2016 (Fig. 6): The precipitation events in the first half of October caused a reduction in (T<sub>b</sub><sup>H</sup>, T<sub>b</sub><sup>V</sup>) correspondingly but both gradually increased afterwards from 160 and 200 K (lower bound) to 200 and 230 K. This period is characterized by a gradual decrease of soil moisture from around 0.30 to 0.24 m³/m³ from 01/10/2016 to 25/11/2016 due to drying and followed by a sharp reduction to 0.11 m³/m³ on 26/11/2016 due to soil freezing. TG, Tair, and ST_2.5 cm all show a gradual decrease in this period, when ST_2.5 cm touched the freezing point on 26/11/2016 indicating the soil layer from surface to the depth of 2.5 cm is frozen.
<figure><img src="image_7.png" /><figcaption>Fig. 5 Same as Fig. 3 but for monsoon period (late June to mid-August) in 2018.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 6 Same as Fig. 3 but for post-monsoon (early October to late November) in 2016.</figcaption></figure>

Winter period 2017–2018 (Fig. 7): $(T_b^H, T_b^V)$ increased till the end of January and decreased again till the end of March 2018 with strong oscillations at the beginning and end of this period. This period is characterized by a gradual decrease in SM_2.5 cm in late November, a very low SM_2.5 cm of around 0.07 m³/m³ in the middle of the period due to frozen soil, and a gradual increase of SM_2.5 cm from late February to March. TG and Tair oscillate around the freezing point in the whole period, but ST_2.5 cm gradually moves from above the freezing point at the beginning of this period to below it and later again above it depicting the freezing processes in the winter period.

Next we present in detail some fascicles of diurnal dynamics of the Maqu ELBARA-III microwave radiometry dataset for different periods. The same variables are plotted as in Figs. 2–7, but for a few typical days in different periods with diurnal characteristics in different seasons highlighted. Similarly characteristic features in each sub-figure are explained in detail.

08/08/2016–04/09/2016, monsoon period (Fig. 8): The gradual increase of $(T_b^H, T_b^V)$ after a precipitation event on 14/8/2016 is due to the drying of the soil profile (panel 1) as a consequence of evaporation when the ground temperature (TG), the air temperature (Tair) and soil temperature (ST_2.5 cm) are all above the freezing level (273.15 K), while the sudden drops in $(T_b^H, T_b^V)$ are caused by a decrease of surface emissivity due to the increase of surface soil moisture by precipitation. The spikes in $T_b^V$ (on 24/08/2016, 30/08/2016, and 03/09/2016 and 04/09/2016) are consequences of surface reflected solar beams into the ELBARA-III antenna horn. The fact that
<figure><img src="image_9.png" /><figcaption>**Fig. 7** Same as Fig. 3 but for winter season (late November to late March) in 2017–2018.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 8** Diurnal dynamics of the Maqu ELBARA-III radiometry dataset for 08/08/2016–04/09/2016 monsoon period. Plotted variables are the same as in Fig. 3, except diurnal characteristics in different seasons are highlighted.</figcaption></figure>

such phenomena occur after precipitation events indicates the presence of water on vegetation and soil surfaces. The geometry for such occurrence requires that the solar elevation is in the range of 44–56° (see Fig. 1 for an illustration of such geometry). More information is provided in next section to guide users wishing to filter out these signals.

11/10/2016–29/10/2016, post-monsoon period (Supplementary Figure S1): The occurrence and amount of precipitation decreased in this period (panel 3). After the two precipitation events on 12/10/2016 and 15/10/2016 which caused a gradual decrease of ($T_b^H, T_b^V$), the soil profile continued the drying (panel 1) as a consequence of evaporation when the ground temperature (TG), the air temperature (Tair) were mostly above the freezing level (273.15 K) during the day but TG and Tair decreased to or below the freezing level at night, while ST_2.5 cm was still above the freezing point both during day and at night. The spikes in $T_b^V$ are still visible but with much reduced amplitude which are mostly below 25 K except on days from 11/10/2016 to 14/10/2016. The much smaller spikes on other days are consequences of melting of surface frozen soil in daytime which is not detected by SM_2.5 cm but can be inferred by the diurnal oscillating TG and Tair around the freezing point in panel 2.
08/08/2017–26/08/2017, monsoon period (Supplementary Figure S2): This monsoon period is similar to that presented in Fig. 8 but showing a delayed first major precipitation event on 20/08/2017 following a prolonged drying down of the soil profile (panel 1).

24/10/2017–14/11/2017, post-monsoon period (Supplementary Figure S3): The period is similar to that in Supplementary Figure S1 but showing a snow event on 29/10/2017–31/10/2017 indicated by the high albedo (0.95 on 29/10/2017 and 0.4 on 31/10/2017). The amplitude of diurnal variations of TG, Tair on these days, and those on 1/11/2017 were much reduced due to the snow event compared to the preceding and the succeeding days. The snow event was a warm snowfall as TG and Tair stayed mostly above the freezing point, and the snow started to melt as seen by the reduced albedo and increased SM_2.5 cm on 29/10/2017 and completely melted on 1/11/2017. The measured mid-day $T_b^H$ jumped by 30 K from 190 K on 28/10/2017 to 220 K on 29/10/2017, while less than 10 K was seen to increase in $T_b^V$ indicating the much bigger sensitivity of $T_b^H$ than $T_b^V$ to the presence of snow on the ground. Other notable features are the freezing events at nights on 13/11/2017 and 14/11/2017 indicated by the reduced SM_2.5 cm (from 0.27 to 0.18 cm³/cm³ on 13/11/2017 and from 0.26 to 0.17 cm³/cm³ on 14/11/2017) and the near freezing point ST_2.5 cm (meaning the soil column above it must be partially frozen as TG and Tair were both below the freezing point). The corresponding $T_b^H$ differ by 30 K and $T_b^V$ by 10 K between mid-day and mid-night. It is important to notice that while both $T_b^H$ and $T_b^V$ increased in the snow and freezing events, the causes of the increases are very different. In the former case, the snow contributed to $T_b^H$ by emission, while in the latter one, it was because of the reduction in soil moisture due to freezing that caused the increase of the emissivity. Similar freezing events can also be seen at the nights on 7/11/2017 and 8/11/2017 but with smaller amplitudes.

14/11/2017–08/12/2017, post-monsoon to winter period (Supplementary Figure S4): The freezing-thawing processes continued from 14/11/2017 till 23/11/2017 when the soil layer above 2.5 cm appeared frozen resulting in reduction of measured SM_2.5 cm from 0.25 at mid-day 21/11/2017 to 0.10 cm³/cm³ at night of 23/11/2017, while it still increased to ca. 0.17 cm³/cm³ because TG (and on some days Tair) increased to above the freezing point during the mid-day due to the strong solar radiation on the Tibetan plateau. $T_b^H$ and $T_b^V$ showed regular patterns modulated by the changes in liquid soil water content governed by freezing-thawing processes (Note TG was always above the freezing point during mid-day).

07/12/2017–31/12/2017, winter period (Supplementary Figure S5): The period is similar to that in Supplementary Figure S4, but for the early winter period. One notable feature is the continuation of the freezing-thawing processes but with reduced amplitude in SM_2.5 cm and the correspondingly reduced amplitude in the variations of $T_b^H$ and $T_b^V$ while their levels keep increasing the whole periods.

04/03/2018–01/04/2018, winter to pre-monsoon period (Supplementary Figure S6): In this period, SM_2.5 cm, TG, Tair, and ST_2.5 cm all started to increase. Except during the few snow events, TG and Tair both reached above freezing point during mid-day. The few brief snow events prevented heat transfer to the ground causing TG and as a result Tair to stay below the freezing point for several days (e.g., on 06-07/03/2018, 19/03/2018, 01/04/2018). The melting of the snow and soil ice after 24/03/2018 by increased solar radiation caused SM_2.5 cm to increase to 0.25 cm³/cm³ during day time. $T_b^H$ and $T_b^V$ briefly increased after the snowfall but restored to regular patterns until about 22/03/2018 when the amplitude of variations increased due to the freezing and thawing processes along with the melting of the snow cover and increase of SM_2.5 cm. This variation was only temporarily interrupted by the six small snow events. The biggest diurnal changes in $T_b^H$ and $T_b^V$ were 100 K and 50 K respectively on 21–22/03/2018.

02/04/2018–30/04/2018, pre-monsoon period (Supplementary Figure S7): In this period, SM_2.5 cm, TG, Tair, and ST_2.5 cm all increased to above the freezing point most of the time. Noteworthy is the snow event on 13/04/2018, which temporarily reduced TG and Tair but ST_2.5 cm to below the freezing point. The precipitation events in this period were likely wet snowfall that melted during the day, as indicated by the sharp reduction of the albedo during day time. $T_b^H$ and $T_b^V$ remained respectively at ca. 150 and 200 K on average in this period except after the snow events when both temporarily increased. The strong variations in both $T_b^H$ and $T_b^V$ in 06–09/04/2018 were due to the freezing of the soil at nights as indicated by the variations in SM_2.5 cm. The similar variations in both $T_b^H$ and $T_b^V$ in 02–04/04/2018 are due to the melting of the snow indicated by reduction of albedo from 1.0 on 01/04/2018 to 0.5 on 02/04/2018 (panel 1) and refreezing at nights and thawing of the soil during daytime as indicated by SM_2.5 cm and TG, Tair and ST_2.5 cm.

06/07/2018–27/07/2018, monsoon period (Supplementary Figure S8): in this period, TG, Tair and ST_2.5 cm were all about 15 K above the freezing point. The more frequent precipitation events were all in the form of rainfall. $T_b^H$ and $T_b^V$ increased respectively to ca. 220 and 240 K on average from 02–12/07/2018. Given that SM_2.5 cm remained at 0.25–0.30 cm³/cm³ and was similar to the last week of April (Supplementary Figure S7), the increase in both $T_b^H$ and $T_b^V$ of 40 and 20 K respectively can only be attributed to vegetation emission because we can assume safely that the surface roughness remained the same from April to July. The increase in $T_b^H$ and $T_b^V$ from 215 and 240 K on 13/07/2018 to 245 and 265 K on 30/07/2018 is caused in part by the drying down in this period, but vegetation emission must have also contributed to the increase because of the favorable condition for vegetation growth in the monsoon period (we do notice the reduction of SM_2.5 cm from 0.30 to 0.15 cm³/cm³ but expect the vegetation have much deeper roots to tap water from deeper soil layers to sustain growth which is also confirmed by the measured leaf area index (LAI) in the fenced area which increased from 3.95 on 12/07/2018 to 7.37 cm²/cm² on 17/08/2018.

28/07/2018–11/08/2018, monsoon period (Supplementary Figure S9): The heavier precipitation events of 5–10 mm per hour in this period caused SM_2.5 cm to increase from 0.15 to 0.30 cm³/cm³ from 28/07/2018 to 03/08/2018. However $T_b^H$ and $T_b^V$ remained respectively at 240 and 265 K indicating the dominant role played by vegetation attenuation and emission. After 03/08/2018, $T_b^H$ and $T_b^V$ decreased slightly by ca. 10 and 5 K respectively.
tively which can be attributed to contributions by the increased SM_2.5 cm to 0.25 cm³/cm³ on average and a further increase in biomass as indicated above by increased LAI.

## Usage Notes

This dataset can be used to validate satellite based observations and retrievals³¹,³², verify radiative transfer model assumptions¹¹ and validate land surface model and reanalysis outputs¹⁴,³⁰, retrieve soil physical properties, as well as to quantify land-atmosphere exchanges of energy, water and carbon and help to reduce discrepancies and uncertainties in current Earth System Models (ESM) parameterizations.

As reported above, there appear some spikes in the observed ELBARA-III brightness temperature (e.g. on 24/08/2016, 30/08/2016, and 03/09/2016 and 04/09/2016, Fig. 8), which are consequences of surface reflected solar beams into the ELBARA-III antenna horn under certain surface conditions. For users who wish to filter out these signals, we suggest two methods which are briefly described as follows: 1) Quantile filtering, and 2) A harmonic analysis of time series (HANTS). These two filtering methods²⁹ are applied to a sample data to demonstrate their effectiveness. The data covers the period from 22/03/2018 to 25/06/2018, with time interval of 30 minutes (see Fig. 3).

**Quantile filtering.** In quantile filtering, we assume that any data point outside of a defined quantile is an outlier. The filtering procedure is as follows:

$$ \text{If } T_b^p(i) > \text{Quantile}[T_b^p(i - K: i + K), q] \\ \text{Then } \text{loc.append(i)} \quad \# \text{ record position} \quad (1) $$

$$ \text{If } PI(i) > \text{Quantile}[PI(i - K: i + K), q] \\ \text{Then } \text{loc.append(i)} \quad \# \text{ record position} \quad (2) $$

where $T_b^p$ is the observed brightness temperature with p (H, V) polarization. q is the quantile (ranging in [0, 1]) to be computed and K is half-time window. In the case of Eq. (1), $q = 0.85$ for $T_b^H$, $q = 0.90$ for $T_b^V$ and K takes the value of 100. For Eq. (2), PI is the polarization index calculated as $((T_b^V - T_b^H)/(T_b^V + T_b^H))$, for which $q = 0.90$ and K takes the value of 100.

**Filtering using the HANTS algorithm.** In applying the HANTS²⁶ algorithm, we define an outlier as

$$ T_b^p(i) > \text{Maxium}[T_b^p\_HANTS(i - K: i + K), q] \quad (3) $$

where $T_b^p\_HANTS$ is the estimated brightness temperature by using the HANTS algorithm. The frequency (nf) as one of inputs to the HANTS algorithm is set at 50 in this case. Descriptions of other inputs such as the valid minimum and maximum values can be found in the script (the uploaded HANTS.py)²⁸. q and K have the same meanings as in Eq. (1), $q = 0.90$ and K takes the value of 150 in this case.

**Filtering results.** Supplementary Figure S10 shows results by quantile filtering of $T_b^p$ and Supplementary Figure S11 shows those by quantile filtering of PI. Supplementary Figure S12 shows seasonal variations of the Maqu ELBARA-III radiometry dataset for pre-monsoon period, in which the displayed $T_b^p$ are based on quantile filtering of $T_b^p$. Supplementary Figure S13 shows results by using the HANTS algorithm and Supplementary Figure S14 shows seasonal variations of the Maqu ELBARA-III radiometry dataset for the same period, in which the displayed $T_b^p$ in the figure are based on filtering using the HANTS algorithm.

The results in Supplementary Figure S10 appear adequate for outliers with big deviations, however the data points represented by black (blue) dots from 26/03/2018 to 09/04/2018 may contain surface freeze/thaw information, and one may want to keep them from being masked out. Supplementary Figure S11 shows extra 'outliers' that have minimum values are also identified while the identified 'outliers' that have maximum values are the same as in Supplementary Figure S10. Therefore, Eq. (1) is sufficient for masking outliers caused by solar reflection. Supplementary Figure S14 shows that $T_b^p$ outliers can be effectively detected by using the HANTS algorithm, and the results are better than those shown in Supplementary Figure S12 by quantile filtering of $T_b^p$ (see also Supplementary Figures S11 and S13 for the identified 'outliers').

## Code availability

The codes for processing the collected data, for plotting the figures, and for downloading SMAP data are included in the dataset²⁸. The codes for filtering the reflected solar signals and for plotting the results are included in the update²⁹ together with the updated brightness temperature that contain corrected local time stamps. Overviews of the codes, data and explanations can be found in Online-only Table 2²⁸ and Table 2²⁹.

Received: 2 April 2020; Accepted: 21 August 2020;
Published online: 30 September 2020

## References

1. Kerr, Y. H. et al. The SMOS soil moisture retrieval algorithm. *IEEE Trans. Geosci. Remote Sens.* **50**, 1384–1403 (2012).

2. Entekhabi, D. et al. The soil moisture active passive (SMAP) mission. *Proc. IEEE.* **98**, 704–716 (2010).

3. Colliander, A. et al. Validation of SMAP surface soil moisture products with core validation sites. *Remote Sens. Environ.* **191**, 215–231 (2017).
4. Kurum, M. et al. A first-order radiative transfer model for microwave radiometry of forest canopies at L-band. *IEEE Trans. Geosci. Remote Sens.* **49**, 3167–3179 (2010).

5. Zhao, Q. & Lang, R. H. Fresnel double scattering from tree branches. *IEEE Trans. Geos. Remote Sens.* **50**, 3640–3647 (2012).

6. Schwank, M., Naderpour, R. & Mätzler, C. "Tau-Omega"-and two-stream emission models used for passive L-band retrievals: Application to close-range measurements over a forest. *Remote Sens.* **10**, 1868 (2018).

7. Kerr, Y. H. et al. Overview of SMOS performance in terms of global soil moisture monitoring after six years in operation. *Remote Sens. Environ.* **180**, 40–63 (2016).

8. Wigneron, J. P. et al. Modelling the passive microwave signature from land surfaces: A review of recent results and application to the L-band SMOS & SMAP soil moisture retrieval algorithms. *Remote Sens. Environ.* **192**, 238–262 (2017).

9. Chaubell, M. J. et al. Improved SMAP dual-channel algorithm for the retrieval of soil moisture. *IEEE Trans. Geosci. Remote Sens.* **58**, 3894–3905 (2020).

10. Su, Z. et al. The Tibetan Plateau observatory of plateau scale soil moisture and soil temperature (Tibet-Obs) for quantifying uncertainties in coarse resolution satellite and model products. *Hydrol. Earth Syst Sci.* **15**, 2303–2316 (2011).

11. Zheng, D. et al. Impact of surface roughness, vegetation opacity and soil permittivity on L-band microwave emission and soil moisture retrieval in the third pole environment. *Remote Sens. Environ.* **209**, 633–647 (2018).

12. Yu, L., Zeng, Y., Wen, J. & Su, Z. Liquid-vapor-air flow in the frozen soil. *J. Geophys. Res. Atmos.* **123**, 7393–7415 (2018).

13. Zheng, D. et al. L-band microwave emission of soil freeze–thaw process in the third pole environment. *IEEE Trans. Geosci. Remote Sens.* **55**, 5324–5338 (2017).

14. Mwangi, S., Zeng, Y., Montzka, C., Yu, L. & Su, Z. Assimilation of cosmic-ray neutron counts for the estimation of soil ice content on the eastern Tibetan Plateau. *J. Geophys. Res. Atmos.* **125**, e2019JD031529 (2020).

15. Rautiainen, K. et al. Detection of soil freezing from L-band passive microwave observations. *Remote Sens. Environ.* **147**, 206–218 (2014).

16. Rautiainen, K. et al. SMOS prototype algorithm for detecting autumn soil freezing. *Remote Sens Environ.* **180**, 346–360 (2016).

17. Entekhabi, D. et al. SMAP Handbook–Soil Moisture Active Passive: Mapping Soil Moisture And Freeze/Thaw From Space. (JPL Publication, Pasadena, CA, 2014).

18. Dente, L., Vekerdy, Z., Wen, J. & Su, Z. Maqu network for validation of satellite-derived soil moisture products. *Int. J. Appl. Earth Observ. Geoinform.* **17**, 55–65 (2012).

19. Peel, M. C., Finlayson, B. L. & McMahon, T. A. Updated world map of the Köppen-Geiger climate classification. *Hydrol. Earth Syst. Sci.* **11**, 1633–1644 (2007).

20. Zhao, H., Zeng, Y., Lv, S. & Su, Z. Analysis of soil hydraulic and thermal properties for land surface modeling over the Tibetan Plateau. *Earth Syst. Sci. Data.* **10**, 1031–1061 (2018).

21. Schwank, M. et al. ELBARA II, an L-band radiometer system for soil moisture research. *Sensors.* **10**, 584–612 (2010).

22. Schwank, M., Matzler, C., Guglielmetti, M. & Fluhler, H. L-band radiometer measurements of soil water under growing clover grass. *IEEE Trans. Geosci. Remote Sens.* **43**, 2225–2237 (2008).

23. Lv, S., Zeng, Y., Wen, J., Zhao, H. & Su, Z. Estimation of penetration depth from soil effective temperature in microwave radiometry. *Remote Sens.* **10**, 519 (2018).

24. Topp, G. C., Davis, J. L. & Annan, A. P. Electromagnetic determination of soil water content: measurements in coaxial transmission lines. *Water Resour. Res.* **16**, 574–582 (1980).

25. Myneni, R., Knyazikhin, Y. & Park, T. MCD15A2H MODIS/Terra + Aqua Leaf Area Index/FPAR 8-day L4 Global 500m SIN Grid V006. NASA EOSDIS Land Processes DAAC (2015).

26. Verhoef, W., Menenti, M. & Azzali, S. Cover a colour composite of NOAA-AVHRR-NDVI based on time series analysis (1981-1992). *Int. J. Remote Sens.* **17**, 231–235 (1996).

27. Hofste, J. G. et al. Year-long, broad-band, microwave backscatter observations of an alpine meadow over the Tibetan Plateau with a ground-based scatterometer. *Earth Syst. Sci. Data Discuss.*, https://doi.org/10.5194/essd-2020-44 (in review, 2020).

28. Su, Z. et al. Multiyear in-situ L-band microwave radiometry of land surface processes on the Tibetan Plateau. figshare https://doi.org/10.6084/m9.fishshare.12058038.v1 (2020).

29. Su, Z. Update: Multiyear in-situ L-band microwave radiometry of land surface processes on the Tibetan Plateau. figshare https://doi.org/10.6084/m9.fishshare.12728444.v1 (2020).

30. Lv, S., Zeng, Y., Su, Z. & Wen, J. A closed-form expression of soil temperature sensing depth at L-band. *IEEE Trans. Geosci. Remote Sens.* **57**, 4889–4897 (2019).

31. Zeng, Y. et al. Blending satellite observed, model simulated, and in situ measured soil moisture over Tibetan Plateau. *Remote Sens.* **8**, 268 (2016).

32. Zhuang, R., Zeng, Y., Manfreda, S. & Su, Z. Quantifying long-term land surface and root zone soil moisture over Tibetan Plateau. *Remote Sens.* **12**, 509 (2020).

## Acknowledgements

This work was supported in part by ESA ELBARA-II/III Loan Agreement EOP-SM/2895/TC-tc, the ESA MOST Dragon IV Program (Monitoring Water and Energy Cycles at Climate Scale in the Third Pole Environment), the Netherlands Organization for Scientific Research under Project ALW-GO/14-29, the National Natural Science Foundation of China (grant no. 41971033) and the Fundamental Research Funds for the Central Universities, CHD (grant no. 300102298307). We thank the Zoige Plateau Wetlands Ecosystem Research Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Science, Lanzhou, China, for providing precipitation data from 1-Sep-2016 to 30-Jan-2017 and we are grateful to colleagues and institutions who have supported our field experiment and campaigns.

## Author contributions

Z.S., J.W., Y.Z. and R.v.d.V. conceptualised the experiment. H.Z., S.L., D.Z., X.W. and Z.W. conducted and maintained the campaigns and experiment. Y.Z., H.Z. and S.L. processed the data. M.S., Y.K., S.Y., A.C., H.Q., M.D. and S.M. provided technical inputs and guidance to the experiment. Z.S. wrote the paper with inputs from all authors.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41597-020-00657-1.
**Correspondence** and requests for materials should be addressed to Z.S. or J.W.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

© The Author(s) 2020