Derivation of Vegetation Optical Depth and Water Content in the Source Region of the Yellow River using the FY-3B Microwave Data

This study uses the brightness temperature at the given microwave frequency (18.7 GHz) from the Microwave Radiation Imager (MWRI) on-board the Fengyun-3B (FY-3B) satellite to improve the τ-ω model by considering the radiative contribution from waterbody in the pixels over the wetland of the Yellow River source region, China. In order to retrieve vegetation optical depth (VOD), a dual-polarization slope parameter is defined to express the surface emissivity in the τ-ωmodel as the sum of soil emissivity and waterbody emissivity. In the regions with no waterbody, the original τ-ω model without considering waterbody impact is used to derive VOD. With use of the field observed vegetation water content (VWC) in the source region of the Yellow River during the summer of 2012, a regression relationship between VOD and VWC is established and then the vegetation parameter b is estimated. The relationship is employed to derive the spatial VWC during the entire vegetation growing period. The VOD retrieved is invalid and failed in some part of the study area by using the previous τ-ωmodel, while the results from the improved τ-ωmodel indicate that the VOD is in the range of 0.20 to 1.20 and the VWC is in the range of 0.20kg/m2 to 1.40kg/m2 in the entire source region of the Yellow River in 2012. Both VOD and VWC exhibit a pattern of low values in the west part and high values in the east part. The largest regional variations appear along the Yellow River. The comparison between the remote-sensing-estimated VWC and the ground-measured VWC gives the root mean square error of 0.12kg/m2. These assessments reveal that with considering the fractional seasonal wetlands in the source region of the Yellow River, the microwave remote sensing measurements from the FY-3B MWRI can be successfully used to retrieve the VWC in the source region of the Yellow River.


Introduction
Vegetation water content, vegetation canopy biomass, and soil moisture strongly affect fluxes of the surface water vapor and energy, control the proportion of surface sensible heat, and latent heat allocated to the usable surface energy, as well as the proportion of rainfall among the surface runoff, infiltration, and evapotranspiration, and thus can have significant impacts on numerical weather, climate and hydrologic predictions [1,2].Optical-infrared satellite remote sensing, such as the Normalized Difference Vegetation Index (NDVI) and Leaf-Area Index (LAI), is commonly used to determine vegetation water content.However, clouds, aerosols, and solar illumination effects degrade the ability to monitor vegetation information.Satellite microwave observations are insensitive to clouds at lower frequencies and independent of solar illumination [3].Currently, soil moisture retrievals using microwave radiometers have a high accuracy over bare soil or sparse vegetation coverage area, whereas ongoing challenges remain in high vegetated regions, due to lack of knowledge about the properties of vegetation structure, vegetation optical depth (VOD), vegetation water content (VWC), and biomass of canopy [3][4][5].To accurately retrieve soil moisture, passive soil moisture retrievals must properly account for the effect of VOD, so developing a feasible VOD algorithm is necessary for soil moisture retrievals.
The VOD is a variable that characterizes the contribution of vegetation canopy decay to soil radiation and increases with vegetation density.Most VOD retrievals are conducted simultaneously with soil moisture retrievals [6,7].As a useful indicator for understanding vegetation state and variability, VOD has also been widely used in a proxy of biomass [8], phenology metrics like growing season length and start of season [9,10], forest and agriculture monitoring [11], land evaporation estimation [12], and vegetation water stress [13].
Many studies revealed that VOD is less prone to saturation effect than optical vegetation index which is sensitive to the green components [14].Meesters et al. [15] proposed a nonlinear iterative algorithm, which can simultaneously derive both soil moisture and VOD.In order to minimize the atmosphere impacts and soil surface emission signals, Shi et al. [16,17] proposed the passive microwave vegetation indices (MVIs) index that depends only on vegetation properties such as vegetation fractional coverage, biomass, water content, temperature, and geometry of vegetation canopy.In the microwave bands (>5.0 GHz), the brightness temperatures are obviously affected by absorption and scattering of vegetation, whereas the L-band can almost penetrate most of the dense natural vegetation.However, VOD retrievals at L-band can increase potential sensitivity at higher biomass levels, and have severe radio frequency interference (RFI) contamination and are more sensitive to soil moisture at lower VOD levels [18].High microwave frequencies (>18.7 GHz), such as K-band, X-band, or C-band, have also been used in previous studies [8,19].For higher frequencies >18.7 GHz, however, water vapor and cloud liquid water in various atmospheric layers absorb and emit scattered radiation [14].Therefore, the 18.7 GHz frequency has low RFI and can reduce atmosphere and soil moisture sensitivities to other channels [6].Jones et al. [14,19] retrieved VOD based on the τ-ω model from both the 18.7 GHz and 23 GHz channels (K-band and X-band frequencies), by determining the influence of open water fraction of each pixel and with given albedo.The VOD parameter contributes to a more comprehensive view of land surface phenology.
Vegetation water content (VWC) is another important parameter that affects microwave signals.Since the 1990s, several visible spectral indexes have been proposed to describe VWC, including normalized difference vegetation index(NDVI), water band index (WBI), normalized difference water index (NDWI), and moisture stress index (MSI) [20].There are several factors that may contribute to limitation of the global VWC study.Products with spectral reflectance and index usually contain noise from the background (e.g., soil and exposed waterbodies).Measurements at the visible/near-infrared band mostly reflect greenness of vegetation because of the limitation of penetrating capacity, which directly affects accuracy of leaf water content estimation.Therefore, studies of the global VWC are limited [21][22][23][24][25].Furthermore, the visible/near-infrared vegetation index reflects only vegetation information on the top of the canopy [26], so it cannot be applied for accurately monitoring VWC.Measurements using passive microwave radiometers are more sensitive to variations in the characteristics of vegetation canopy, and can be directly used to retrieve canopy VWC without affected by weather conditions.These observations have the potential to be used for all-weather monitoring of regional vegetation canopy.VWC is proportional to microwave VOD of the above-ground biomass canopy and essentially a linear function of VOD (VOD = b*VWC) [26].The parameter b depends on the vegetation structure and sensors.However, the parameter b was usually assumed constant because of the limitation of vegetation information or satellite remote sensing data.Consequently, the results showed the deficiencies in both temporal and spatial scales [26].
The source region of the Yellow River is located over the Yellow River Basin in the northeast part of the Tibetan Plateau.The vegetation coverage is dominated by alpine meadow with a broad distribution of seasonal marsh and wetland [27].Study of VWC in the source region of the Yellow River is necessary and can provide accurate microwave retrievals of soil moisture and first-hand data of soil and water loss and desertification in this region.However, due to the nonuniform land surface and seasonal wetland in this region, remote sensing signal may be subject to strong interference from the waterbodies in addition to the vegetation and soil.
The primary objectives of this paper are to (a) obtain more accurate VOD information on seasonal vegetation growth in the source region of the Yellow River based on theτ-ω model and open water fraction (fw) by using the brightness temperatures (Tb) at 18.7 GHz (K band) of the Fengyun-3B (FY-3B) Microwave Radiation Imager (MWRI) and (b) use canopy VWC from the field measurements and the correspondence between VOD record to obtain spatiotemporal distribution of VWC in the whole source region of the Yellow River.The paper is organized as follows.In Section 2 the algorithm is proposed.The study area and data source are described in Section 3. Section 4 presents results from the retrievals with the proposed algorithm, and key implications of the results are also discussed.Finally, a summary of the results and conclusions of the study is reported in Section 5.

Classical Retrieval Approach for VOD and VWC
Passive microwave remote sensing retrieval for both VOD and soil moisture is based on the physical mechanism of the classical microwave transfer equation [28].The simplified microwave radiation transfer equation is expressed as follows Equation (1) indicates that the brightness temperature of surface microwave radiation, T p b , is composed of three parts: The first is emissivity of vegetation-covered soil which is sum of microwave emission of the underlying soil with emissivity e p s attenuated by overlying vegetation with transmissivity Г.The second is direct canopy emission with single scattering albedo, ω p veg .The third is downwelling canopy emission reflected by the soil and attenuated by the canopy.Superscript p represents the polarization type; soil moisture, soil temperature, and soil conditions (e.g., clay, sand, and silt content) can be measured from the Maqu regional scale soil moisture and soil temperature monitoring network.The soil emissivity depends on soil dielectric properties according to the Fresnel coefficients modified for surface roughness [29].The soil bulk dielectric properties are determined from a dielectric mixing model [30].ω p veg is the vegetation single scattering albedo of 0.05, Γ p is the transmittance of vegetation, and T soil and T veg represent the thermodynamic temperatures of surface soil and canopy vegetation, respectively, and are considered equal to surface temperature in this study [14].In this study, the method of Holmes [31] is used to estimate the surface temperature with the 37GHz vertically-polarized brightness temperature.The vegetation transmittance, Γ p , can be given as where τ p veg is the microwave optical thickness of vegetation and µ is the observation angle of the sensor.Equations ( 1) and ( 2) are usually called the τ-ω model.Given that the other parameters in Equations ( 1) and ( 2) are known, one can derive VOD.
Jackson and O'Neill [32] found that for a particular type of canopy, the relationship between VOD (τ veg ) and VWC can be described by the following equation.
where b is the empirical coefficient of the linear regression and depends on the vegetation structure, canopy type, and wavelength [33].
In some researches, the parameter b between VOD and canopy VWC for different underlying vegetation surfaces is considered as a constant, and the regression equation between NDVI and VWC based on field observation experiments is also established [34,35].In addition to the microwave brightness temperature, these approaches requires an estimate of the physical temperature of the surface, the VOD, and a vegetation parameter b.Physical temperature and VOD can be determined using the existing algorithms and sensor systems.However, the dependence of the vegetation parameter b on various physical and sensor variables is not determined exactly.Moreover, the value of this parameter is only suitable for large-scale canopy vegetation in the previous researches, and thus is not suitable for seasonal wetland.In this study, with the assumption that VOD is closely related with vegetation canopy density, both the retrieved VOD from the microwave remote sensing data and VWC from the field observations are used to estimate the parameter b, and consequently, the spatial VWC can be obtained.

Determination of Influence of Open Water Fraction
Under conditions of coverage of seasonal wetland in the source region of the Yellow River, the brightness temperatures from the microwave radiometer measurements contain the signal of waterbodies, except for those of exposed soil and vegetation.The emissivity of waterbodies can vary with wavelength and salinity, although it is generally believed that waterbodies on land surface are stable at the low-frequency band (≤18.7 GHz).Jones et al. [14] presumed that mixed pixels are composed of waterbodies and vegetation, with the vegetation component containing both vegetation canopy and soil covered by vegetation.In this case, the emissivity of each pixel can be expressed as the sum of the composition accounted for by the emissivity of waterbodies and that of vegetation: where e p is the surface emissivity of mixed pixels.e l,p is the effective emissivity of land surface with vegetation coverage and its value is 0.90 for high biomass vegetation cover in this study, which is also consistent with previous studies [14,36].e w,p is the emissivity of waterbodies and f w is the fraction of open water.The parameter e p can be simplified and linearized from its usual quadratic form by ignoring atmospheric emissions reflected by the surface [14]: where T bp is the brightness temperature observed by satellite, T s is the ancillary surface temperature, t a is the atmospheric transmittance (a function of atmospheric water vapor content V), and ∆T is the ground-air temperature difference.In this study a similar algorithm from Jones et al. [19] is used to estimate the ancillary T s and V separately.f w can be estimated with Equations ( 4) and (5).
A dual-polarization slope parameter is defined, as it represents the ratio (α) between the V polarization of all pixels (e v ) subtracted by the V polarization of the waterbodies (e w,v ) and the H polarization of all pixels (e h ) subtracted by the H polarization of waterbodies (e w,h ): where the parameter α is not a function of f w .By combining the τ-ω model and the dual-polarized slope parameter (Equations ( 1), (2), and ( 6)), the effective optical depth of land fraction is determined by considering the influence of waterbodies as follows [14]: Remote Sens. 2019, 11, 1536 5 of 18 where A, B, and C are given by The vegetation single scattering albedo (ω) fixs to 0.05 [14].r s,h and r s,v are the Fresnel reflectivities of the soil and determined by soil moisture that define the bounds of minimum and maximum VOD.In addition, in order to improve the accuracy of VOD, the microwave emissivity of surface soil is obtained by the measurements from the Maqu regional scale soil moisture and soil temperature monitoring network by using the Fresnel reflection equation (e s,p = 1 − r s,p ).

Study Area
The Yellow River source region is located in the northeastern Tibetan Plateau of 95 • 30 -103 • 30 E and 31 • 30 -36 • 30 N (Figure 1) and has an area of approximately 122,000 km 2 .The six stations represent different catchments in the Yellow River source region.MaduoX and Maduo represent the headwaters of the Yellow River source region; Dari, Jiuzhi, and Maqu are considered the primary runoff stations; and Tangnaihai is a typical station located at the outlet of the Yellow River source region.The terrain in the region is generally high in the west and low in the east, and the most areas are above the elevation of 3000-5000 m.The annual averaged temperatures in the source region are ~1.2 • C, and generally low in the west and high in the east, representing the typical continental climate of the plateau.The annual rainfall is 250-750 mm and mainly concentrated in June-September [27].The land use in the region is generally a weathered and broken layer of sandy bedrock covered by a thick layer of black soil and a large amount of pebbles.The soil type includes alpine meadow soil and brown calcareous soil.The two blue regions in the upper left corner of Figure 1 denote Gyaring Lake and Ngoring Lake (97 • 02 -97 • 55 E and 34  The vegetation type (Figure 2) in the entire source region of the Yellow River is dominated by grass with sporadic shrubs on the plateau, which are characterized by a uniform land cover of short grassland used for grazing by sheep and yaks.The terrestrial pattern can be identified according to the digital elevation model.Sporadic shrubs dominate the areas in the top and middle elevations above 4200 m.

Field Observation Data
Ground measurements are critical, and a large number of point measurements can be able to provide adequate and accurate vegetation information at relevant scales.However, there are only a few such datasets available from specific field campaigns for a limited range of land cover types.Due to the coarse resolution passive microwave measurements at tens of kilometers, and twenty sampling sites in the Maqu regional scale soil moisture and soil temperature (SMST) monitoring network are selected to collect vegetation samples.Figure 1 also gives the locations and names of the monitoring sites, and the 20 sites chosen for vegetation samples which are spread over the large valley of the Yellow River.There are 11 monitoring stations installed around the valleys of the Yellow and Black Rivers, three stations in the valleys between the hills in the eastern part of the area, four stations on hill slopes, and two stations on wetlands.The vegetation samples covers an area of approximately 80 km×80 km and the site locations are selected in order to monitor the area extensively for various altitudes and slopes.Usually, the measurements obtained at points or at very high resolutions are not comparable to the spatial scale of the passive microwave satellite measurements.In this study, the total number of stations installed in the study area are considered to be a good match between the need for sample vegetation characteristics over a wide area (covering at least one resolution cell of the satellite-derived products), and practical limitations, such as the constraints of cost and time for installation and site revisits [37].To exclude the interference of snow and frozen water from the microwave signal, the period from early-June to mid-October is used as the experimental period.Each site is collected five times so it can basically reflect the vegetation growth in the source region from summer to early-fall.Unfortunately, five samples were scorched by misoperation.So, a total of 95 samples are obtained during the entire experiment

Field Observation Data
Ground measurements are critical, and a large number of point measurements can be able to provide adequate and accurate vegetation information at relevant scales.However, there are only a few such datasets available from specific field campaigns for a limited range of land cover types.Due to the coarse resolution passive microwave measurements at tens of kilometers, and twenty sampling sites in the Maqu regional scale soil moisture and soil temperature (SMST) monitoring network are selected to collect vegetation samples.Figure 1 also gives the locations and names of the monitoring sites, and the 20 sites chosen for vegetation samples which are spread over the large valley of the Yellow River.There are 11 monitoring stations installed around the valleys of the Yellow and Black Rivers, three stations in the valleys between the hills in the eastern part of the area, four stations on hill slopes, and two stations on wetlands.The vegetation samples covers an area of approximately 80 km × 80 km and the site locations are selected in order to monitor the area extensively for various altitudes and slopes.Usually, the measurements obtained at points or at very high resolutions are not comparable to the spatial scale of the passive microwave satellite measurements.In this study, the total number of stations installed in the study area are considered to be a good match between the need for sample vegetation characteristics over a wide area (covering at least one resolution cell of the satellite-derived products), and practical limitations, such as the constraints of cost and time for installation and site revisits [37].To exclude the interference of snow and frozen water from the microwave signal, the period from early-June to mid-October is used as the experimental period.Each site is collected five times so it can basically reflect the vegetation growth in the source region from summer to early-fall.Unfortunately, five samples were scorched by misoperation.So, a total of 95 samples are obtained during the entire experiment period, including a total of 87 samples from the monitoring network, and eight samples from the source region in Maduo and MaduoX (Figure 1) required for validation.In order to improve measurement accuracy, convenience of transportation to the sites must be considered to ensure that the collected vegetation can be immediately sent back to the base station to be weighed and dried to avoid any loss of VWC.
The transit time of FY-3B over the study region was approximately 13:40 to 14:00 at local time, and the time of sampling was synchronous with the transit time of the satellite.For all the observation sites, we measured the vegetation wet weight, coverage, and height within one area unit.However, the rainy season for the study area is summer, and when it rains, the stomata of leaves open and raindrops would adhere to vegetation, causing inaccurate measurements.Therefore, the surface observations in this kind of experiments were processed under clear skies and good weather conditions when the satellite orbited over the sites.As a result, the field observations have strong representativeness.
Destructive above-ground biomass samples were collected in 30 cm × 30 cm.The location coordinates, name, number of strains, height, coverage, the soil dry-wet and weather conditions, and growth of the plant strain were also recorded.The surface dew was processed with a silica-gel desiccant and then rapidly loaded into a sealed bag with negligible weight.The bag labeled with the vegetation type, sampling time, and sampling site location was immediately returned to the laboratory for measuring fresh weight.The temperature was first fixed at 105 • C for 15 minutes and then adjusted to 80 • C for 1-2 days, and then the above-ground biomass were oven dried to a constant weight after which the dry weight was determined.
At all sampling sites, the meadow grew rapidly during the experimental period.The temporal variation of VWC within a square meter region is shown in Figure 3.The overall change of water content indicates increase first and then decrease.From late-May to mid-July, the VWC changes from 0.2kg/m 2 to approximately 0.8 kg/m 2 , and it reaches the maximum of approximately 1.6 kg/m 2 in late-July and then continuously decreases from early-August to early-October.Consequently, VWC becomes approximately 0.5 kg/m 2 in early-October, although the value of certain individual measurements continues to decrease.For example, the value on June 30 is obviously lower than that on several days before.The log indicates that the weather in the study region was continuously dry, so the surface was dry.Dry weather over three continuous days produced dry soil.Soil moisture has a significant impact on the absorption systems of vegetation and the stomata closure of leaves.The overall VWC first increases and then decreases, and the vegetation growth in later months were slow because the vegetation became matured.
The transit time of FY-3B over the study region was approximately 13:40 to 14:00 at local time, and the time of sampling was synchronous with the transit time of the satellite.For all the observation sites, we measured the vegetation wet weight, coverage, and height within one area unit.However, the rainy season for the study area is summer, and when it rains, the stomata of leaves open and raindrops would adhere to vegetation, causing inaccurate measurements.Therefore, the surface observations in this kind of experiments were processed under clear skies and good weather conditions when the satellite orbited over the sites.As a result, the field observations have strong representativeness.
Destructive above-ground biomass samples were collected in 30 cm x 30 cm.The location coordinates, name, number of strains, height, coverage, the soil dry-wet and weather conditions, and growth of the plant strain were also recorded.The surface dew was processed with a silica-gel desiccant and then rapidly loaded into a sealed bag with negligible weight.The bag labeled with the vegetation type, sampling time, and sampling site location was immediately returned to the laboratory for measuring fresh weight.The temperature was first fixed at 105°C for 15 minutes and then adjusted to 80°C for 1-2 days, and then the above-ground biomass were oven dried to a constant weight after which the dry weight was determined.
At all sampling sites, the meadow grew rapidly during the experimental period.The temporal variation of VWC within a square meter region is shown in Figure 3.The overall change of water content indicates increase first and then decrease.From late-May to mid-July, the VWC changes from 0.2kg/m 2 to approximately 0.8 kg/m 2 , and it reaches the maximum of approximately 1.6 kg/m 2 in late-July and then continuously decreases from early-August to early-October.Consequently, VWC becomes approximately 0.5 kg/m 2 in early-October, although the value of certain individual measurements continues to decrease.For example, the value on June 30 is obviously lower than that on several days before.The log indicates that the weather in the study region was continuously dry, so the surface was dry.Dry weather over three continuous days produced dry soil.Soil moisture has a significant impact on the absorption systems of vegetation and the stomata closure of leaves.The overall VWC first increases and then decreases, and the vegetation growth in later months were slow because the vegetation became matured.
Few studies have been conducted about the vegetation detection by using the FY-3B microwave data sets [26,29,37].To examine if the FY-3B MWRI has a capacity to catch the vegetation signals and obtain more accurate information on seasonal vegetation growth in the source region of the Yellow River, we use FY-3B MWRI 18.7 GHz (K band) brightness temperatures provided by the National Satellite Meteorological Center and other ancillary data inputs to derive the VOD and VWC.The daily VOD retrievals are temporally composited into 16-day median values to coincide with MODIS NDVI and MVI (see Section 4.3.2),respectively.VOD retrievals at lower (e.g., 6.9 and 10.7 GHz) microwave frequencies are reported to increase potential VOD sensitivity at higher biomass levels and sensitive to soil moisture at lower VOD levels.However, the brightness temperatures at lower frequency (≤18.7 GHz) are susceptible to radio frequency interference (RFI) mainly from surface and space-based telecommunications, so they were masked by employing frequency and polarization thresholds [18].The channel at 18.7 GHz used in this study has low RFI and reduced atmosphere and soil moisture sensitivity relative to other channels [6].

Spatial Distribution of 16-Day Mean NDVI and MVI
The 16-day composited NDVI is selected from Moderate Resolution Imaging Spectroradiometer (MODIS) data.The FY-3B MWRI is derived from MVI by using the frequency pair in 10.65 GHz and 18.7GHz [40].Mean values of the NDVI are approximately 0.30 and 0.65 during the spring and summer seasons, respectively.For both seasons, the NDVI exhibits high values in the southern parts and low values in the northwest of the study region.The patches in the east and south parts of the study region with high NDVI values are identified as the meadow.The minimum value is located in the northwest of the study region where the Gyaring Lake and Ngoring Lake located.MVI is affected only by the vegetation properties and is negatively correlated to NDVI.Compared with the NDVI (Figure 4a,b), the corresponding MVI (Figure 4c,d) shows a great negative which exhibits lower values in the southern part and higher values in the northwest of the study region.However, the MVI exhibits obviously spatial difference especially during the summer season (Figure 4d).This demonstrates that the influence of vegetation on microwave signals also increases as the vegetation density increases.MVI can provide information on vegetation water content, structure and biomass which is complementary to that provided by NDVI.
Few studies have been conducted about the vegetation detection by using the FY-3B microwave data sets [26,29,37].To examine if the FY-3B MWRI has a capacity to catch the vegetation signals and obtain more accurate information on seasonal vegetation growth in the source region of the Yellow River, we use FY-3B MWRI 18.7 GHz (K band) brightness temperatures provided by the National Satellite Meteorological Center and other ancillary data inputs to derive the VOD and VWC.The daily VOD retrievals are temporally composited into 16-day median values to coincide with MODIS NDVI and MVI (see Section 4.2.3),respectively.VOD retrievals at lower (e.g., 6.9 and 10.7 GHz) microwave frequencies are reported to increase potential VOD sensitivity at higher biomass levels and sensitive to soil moisture at lower VOD levels.However, the brightness temperatures at lower frequency (≤18.7 GHz) are susceptible to radio frequency interference (RFI) mainly from surface and space-based telecommunications, so they were masked by employing frequency and polarization thresholds [18].The channel at 18.7 GHz used in this study has low RFI and reduced atmosphere and soil moisture sensitivity relative to other channels [6].

Fractional Coverage of Open Waterbodies
The underlying surface in the source region of the Yellow River includes seasonal wetland coverage.Therefore, to investigate spatial variations in VOD and VWC with considering seasonal wetland coverage, Equations ( 3) and ( 4) are combined to derive spatial distribution of fractional coverage of open water on June 5, July 30, August 30, September 9, and October 15, separately (Figure 5). Figure 5 shows that the waterbodies during the experimental period were mainly distributed in the area from the source region of the Yellow River to Maqing.In particular, the variation in fractional coverage of open water in the surrounding area of Gyaring Lake and Ngoring Lake was large (0.20-0.40), whereas the variation in the area along the Yellow River was in the range of 0.15 to 0.20.Over the whole of study region, the fractional coverage of open water in early-June and mid-October was lower than that in the other months, and the value was the highest and most obvious in late-August.The algorithm assumes that a given grid cell is composed of a mixture of open water and vegetated land, but does not specifically account for inundated pixels with overlying vegetation.Thus, the accuracy of the effective surface emissivity is improved by considering the fractional coverage of open waterbodies.Moreover, this method can be used for monitoring the composition of surface-exposed waterbodies on the scale of low-frequency microwave radiometer, and providing a basis for studying variation in the seasonal wetland in the source region of the Yellow River.The 16-day composited NDVI is selected from Moderate Resolution Imaging Spectroradiometer (MODIS) data.The FY-3B MWRI is derived from MVI by using the frequency pair in 10.65 GHz and 18.7GHz [40].Mean values of the NDVI are approximately 0.30 and 0.65 during the spring and summer seasons, respectively.For both seasons, the NDVI exhibits high values in the southern parts and low values in the northwest of the study region.The patches in the east and south parts of the study region with high NDVI values are identified as the meadow.The minimum value is located in the northwest of the study region where the Gyaring Lake and Ngoring Lake located.MVI is affected only by the vegetation properties and is negatively correlated to NDVI.Compared with the NDVI (Figure 4a,b), the corresponding MVI (Figure 4c,d) shows a great negative which exhibits lower values in the southern part and higher values in the northwest of the study region.However, the MVI exhibits obviously spatial difference especially during the summer season (Figure 4d).This demonstrates that the influence of vegetation on microwave signals also increases as the vegetation density increases.MVI can provide information on vegetation water content, structure and biomass which is complementary to that provided by NDVI.

Spatial Variation of Retrieved VOD
To determine whether the VOD model is beneficial for analyzing waterbodies, Figures 6 and 7 present the spatial distribution of VOD in the source region of the Yellow River in summer, including the results with and without the influence of waterbodies.VODr represents the VOD without the influence of waterbodies, namely, the reference VOD, and VODe represents the VOD considering the influence of waterbodies.Unlike VWC or soil moisture, VOD can not be observed or verified directly by field work.As shown in the Figures 6 and 7, the results from both methods exhibit obvious seasonal changes.High-value distributions occurred along the Yellow River, where they varied in the range of 0.60 to 1.00 and generally exhibited a pattern of low in the west and high in the east.Over some part of the area, compared with these two figures, VOD without considering the contribution of waterbodies cannot obtain any results in Figure 7b-d.This demonstrates that during the period from July to September, the results from the previous method are poor in some part of the region, while the algorithm with considering waterbodies gives better results.This also means the VOD method without considering the influence of waterbodies is invalid and failed in some area at least in this study.

Temporal Variation of Retrieved VOD
To analyze the seasonal variation trend of VOD, two representative regions are selected in the study area: One is the Maqu area representing the southern part and the other is the Maduo area representing the northwest part.The fresh vegetation bag labeled with the vegetation type, sampling time must be immediately returned to the laboratory and both the sites have the laboratories for measuring and drying fresh vegetation.In each representative area, two representative sampling sites separated by 30 km are selected.Figure 8 shows the temporal distribution of the VOD at four selected sampling sites.The VODr (gray trend line) and VODe (dot-dashed red line) represent the variation trends of the algorithm with and without the influence of waterbodies, respectively.The blue line represents the temporal variation trend of water-body coverage.In the Maqu area or the Maduo area of Sister Lake, the VOD variation with including the influence of waterbodies is small and stable, indicating less interference from the waterbodies.In the Maqu area, during a period from early-June to late-September, the VOD change with considering waterbodies is smaller than that without considering waterbodies, so the former is stable.The maximum VOD appears in early-August, and is approximately 0.60.Because of high elevation of the Maduo area, the temperature and the vegetation density of the meadow are low, VOD is small, and the maximum VOD does not exceed 0.40.Moreover, the variation in the Maduo area is less obvious than that in the Maqu area.The seasonal trend of water-body coverage is consistent in the same region.In the source region, because the water coverage in Sister Lake is abundant, the water-body coverage shows higher than that in the Maqu area.In addition, the algorithm fails obtaining VOD at one site of Maqu without considering the influence of waterbodies.For the spatial distribution of VOD with considering the influence of waterbodies, the high-value region mainly occurred in the area along the Yellow River in the eastern source region and in the Tibetan Autonomous Prefecture of Aba in the southeastern source region, whereas the low-value region could be seen in the western Yellow River source region, especially in the west of Maqing.In particular, in the region from the west of Maqing to the east of Ngoring Lake, the VOD variation was obvious, with strong monthly variation.Another region with dramatic variation appeared in the east of the source region where the water content and temperature conditions were the decisive factors for vegetation growth, which indicates that the vegetation in this area was sensitive to surface temperature and moisture.The distribution of the VOD without considering the influence of waterbodies indicates that in addition to the derivation error of optical depth caused by errors in the soil moisture and temperature data, waterbodies could also interfere with the results if there was surface water in the mixed pixels, especially near Gyaring Lake and Ngoring Lake in the source region of the Yellow River.If the τ-ω model is used to retrieve VOD, the waterbodies in the mixed pixels would produce a small VOD because the application of a simplified surface radiation transfer equation cannot exclude the influence of waterbodies.
Since the surface cover in the source region of the Yellow River is dominated by alpine meadow and the coverage type is simple, the spatiotemporal distribution type of the VOD that includes the influence of waterbodies may be characterizes the seasonal timing of vegetation growing seasons, canopy growth and senescence.The timing, rate and duration of these events are vital on vegetation photosynthesis, carbon sequestration and land-atmosphere water and energy exchange.Therefore, in addition to acting as a criterion for evaluating the optical depth in the derivation of soil moisture, VOD may be used as a characteristic variable of the vegetation to characterize comprehensive features of vegetation canopy.

Temporal Variation of Retrieved VOD
To analyze the seasonal variation trend of VOD, two representative regions are selected in the study area: One is the Maqu area representing the southern part and the other is the Maduo area representing the northwest part.The fresh vegetation bag labeled with the vegetation type, sampling time must be immediately returned to the laboratory and both the sites have the laboratories for measuring and drying fresh vegetation.In each representative area, two representative sampling sites separated by 30 km are selected.Figure 8 shows the temporal distribution of the VOD at four selected sampling sites.The VODr (gray trend line) and VODe (dot-dashed red line) represent the variation trends of the algorithm with and without the influence of waterbodies, respectively.The blue line represents the temporal variation trend of water-body coverage.In the Maqu area or the Maduo area of Sister Lake, the VOD variation with including the influence of waterbodies is small and stable, indicating less interference from the waterbodies.In the Maqu area, during a period from early-June to late-September, the VOD change with considering waterbodies is smaller than that without considering waterbodies, so the former is stable.The maximum VOD appears in early-August, and is approximately 0.60.Because of high elevation of the Maduo area, the temperature and the vegetation density of the meadow are low, VOD is small, and the maximum VOD does not exceed 0.40.Moreover, the variation in the Maduo area is less obvious than that in the Maqu area.The seasonal trend of water-body coverage is consistent in the same region.In the source region, because the water coverage in Sister Lake is abundant, the water-body coverage shows higher than that in the Maqu area.In addition, the algorithm fails obtaining VOD at one site of Maqu without considering the influence of waterbodies.

Correlations between Retrieved VOD and Vegetation Indexes
The VOD parameter is an indicator of canopy density and VWC, while NDVI is directly related to the photosynthetic canopy cover and mainly depends on green leaf material of the vegetation cover.MVI includes information of both leafy and woody parts of vegetation because of greater penetration and sensitivity.The ranked Spearman correlation coefficients [20] between retrieved VOD against NDVI or MVI during the experimental period are given in Figure 9.These

Correlations between Retrieved VOD and Vegetation Indexes
The VOD parameter is an indicator of canopy density and VWC, while NDVI is directly related to the photosynthetic canopy cover and mainly depends on green leaf material of the vegetation cover.MVI includes information of both leafy and woody parts of vegetation because of greater penetration and sensitivity.The ranked Spearman correlation coefficients [20] between retrieved VOD against NDVI or MVI during the experimental period are given in Figure 9.These metrics are potentially complementary and together provide a more comprehensive VOD assessment.While we might expect a favorable, positive correlation between VOD and NDVI time series, negative correlation between VOD and MVI time series in Figure 9. Here, the interferences from Gyaring and Ngoring Lake on these statistics are not considered.As shown in Figure 2, the distribution patterns of the vegetation index derived from the optical and microwave sensors exhibits a good correlation with the land cover type.Land cover categories with characteristically high biomass levels show the highest correlations, with decreasing correspondence for land cover classes with characteristically lower canopy density and biomass.It shows that there is a significant positive correlation between NDVI and VOD in the eastern part of the source region of the Yellow River, but the positive correlation becomes gradually weak in the northwest of the study area.The correlation becomes zero near the boundary areas.In the east part of the study area, the values of correlation between NDVI and VOD can reach 0.8 to 1.0, but the correlation is going to descend in the northwest especially around the two lakes.The strong negative correlations between MVI and VOD are occurred in the east and north part, and the values can reach −0.5 to −0.9.There are many factors that may contribute to correlations between VOD and vegetation index in the study.The fractions of uneven terrain, wet and mixed snow, ice cover and frost cover are found to reduce the correlation.The MVI shows a significant large negative correlation with the VOD in the northwest and southeast of the source region area.

VWC Retrievals
Although VOD is usually assumed to be linearly proportional to VWC [32], a good logarithmic correspondence between these two parameters is obtained in this study.Figure 10 (black line represents logarithm trend line) presents logarithmic function between VOD and VWC, there is a good corresponding relationship between VWC and VOD in the source region during the summer season.According to the Figure 10, the parameter b is determined by the logarithm trend line y = 1.5524 ln(x) + 1.5566.With this relationship the spatial VWC is obtained during the entire growth period (Figure 11).
Since the iterative solution of the VOD directly retrieved based on the τ-ω model cannot be obtained in certain regions, we retrieve spatial distribution of the VWC in the source region with including the impact of waterbodies, as shown from Figure 11.The VWC in the entire source region indicates a variation in the range of 0.20 to 1.20kg/m 2 , and the high VWC appears along the Yellow River and in the southeastern part where the Aba-Tibetan Autonomous Prefecture located.The VWC exhibits less seasonal variation in the western source region, especially in the area from the west of Maqing to the east of Ngoring Lake.This result indicates that the vegetation is sparse in this region, which is dominated by exposed soil.Because the water content and temperature are the most important factors for vegetation growth in the source region of the Yellow River, a large variation of 0.30 to 0.80kg/m 2 appears in the Maqu area to the east, indicating that the vegetation in the Maqu area is more sensitive to surface temperature and moisture.In the northern region near the Qaidam Basin, the vegetation coverage is low and the VWC value is small as well.Such of those results were found for an AMSR-E daily brightness temperature by Chen et al. [41].That study did not specify analysis the quantitative relationship between the VOD and VWC, although a field observation and an inversion model of vegetation index is mentioned for VWC.In the current study, no statistically significant phase lags are found for any of the vegetation classes.Although VOD is usually assumed to be linearly proportional to VWC [32], a good logarithmic correspondence between these two parameters is obtained in this study.Although VOD is usually assumed to be linearly proportional to VWC [32], a good logarithmic correspondence between these two parameters is obtained in this study.Since the iterative solution of the VOD directly retrieved based on the τ-ω model cannot be obtained in certain regions, we retrieve spatial distribution of the VWC in the source region with including the impact of waterbodies, as shown from Figure 11.The VWC in the entire source region indicates a variation in the range of 0.20 to 1.20kg/m 2 , and the high VWC appears along the Yellow River and in the southeastern part where the Aba-Tibetan Autonomous Prefecture located.The VWC exhibits less seasonal variation in the western source region, especially in the area from the west of Maqing to the east of Ngoring Lake.This result indicates that the vegetation is sparse in this region, which is dominated by exposed soil.Because the water content and temperature are the most important factors for vegetation growth in the source region of the Yellow River, a large variation of 0.30 to 0.80kg/m 2 appears in the Maqu area to the east, indicating that the vegetation in the Maqu area is more sensitive to surface temperature and moisture.In the northern region near the Qaidam Basin, the vegetation coverage is low and the VWC value is small as well.Such of those results were found for an AMSR-E daily brightness temperature by Chen et al. [41].That study did not specify analysis the quantitative relationship between the VOD and VWC, although a field observation and an inversion model of vegetation index is mentioned for VWC.In the current study, no statistically significant phase lags are found for any of the vegetation classes.

Validation of VWC
Maqu regional scale soil moisture and soil temperature monitoring network is set up in the water source region of the Yellow River to validate the accuracy of the satellite-derived soil moisture products over the Tibetan Plateau.This network has the capability to monitor the spatial and temporal soil moisture variability of the area with a high degree of accuracy [37,42].The VWC sampling is collected according to the network site.Dente et al [42] proposed that the weighted spatial average of measured soil moisture was successfully used as the ground reference for validation of the AMSR-E soil moisture products and ASCAT soil wetness index products.The results from Dente's study are very encouraging for our applications.So, the match problem can be ignored in our opinion.
Figure 12 presents a comparison between the measured VWC at the sampling sites and the derived VWC from the remote sensing (red line represents reference line of zero bias).During the

Validation of VWC
Maqu regional scale soil moisture and soil temperature monitoring network is set up in the water source region of the Yellow River to validate the accuracy of the satellite-derived soil moisture products over the Tibetan Plateau.This network has the capability to monitor the spatial and temporal soil moisture variability of the area with a high degree of accuracy [37,42].The VWC sampling is collected according to the network site.Dente et al. [42] proposed that the weighted spatial average of measured soil moisture was successfully used as the ground reference for validation of the AMSR-E soil moisture products and ASCAT soil wetness index products.The results from Dente's study are very encouraging for our applications.So, the match problem can be ignored in our opinion.
Figure 12 presents a comparison between the measured VWC at the sampling sites and the derived VWC from the remote sensing (red line represents reference line of zero bias).During the experimental period from early-June to mid-October, 95 samples were collected in the southern area.The validation data, however, were collected from the northwest of the region.The altitude of the northwest of the study area is more than 4200 m.Compared with the southern part, the northwest area shows a higher latitude and less vegetation cover (Figures 1 and 2).Because the validation region was less affected by human activity and the road was rough, we only collected samples from four sites.The root mean square error is 0.12 kg/m 2 .At one sampling, the value for the remote sensing retrieval was ~0.15 kg/m 2 , but the value for the surface observation was ~0.55 kg/m 2 .This discrepancy might have been caused by the sampling site location, which was close to the lake.Thus, the sensor likely received more signals from the waterbody than the vegetation.Figure 12 shows the retrieved VWC from the remote sensing is smaller than the measured one.The difference is resulted from the model systematic error that the non-exposed waterbodies are excluded.However, the difference may also be contributed by the fixed values of the parameters used during the iterative processes, such as the effective emissivity of vegetation.Moreover, the observation angle and transient field of the view of the sensor can also affect the model accuracy.
the remote sensing retrieval was ~0.15 kg/m , but the value for the surface observation was ~0.55 kg/m 2 .This discrepancy might have been caused by the sampling site location, which was close to the lake.Thus, the sensor likely received more signals from the waterbody than the vegetation.Figure 12 shows the retrieved VWC from the remote sensing is smaller than the measured one.The difference is resulted from the model systematic error that the non-exposed waterbodies are excluded.However, the difference may also be contributed by the fixed values of the parameters used during the iterative processes, such as the effective emissivity of vegetation.Moreover, the observation angle and transient field of the view of the sensor can also affect the model accuracy.

Summary and Conclusions
In this study, the τ-ω model was improved by using the dual-polarized brightness temperature observations at the K-band obtained from the MWRI microwave radiometer on board the FY-3B meteorological satellite and removing the influence of waterbodies.By constructing a microwave VOD model that is suitable for the source region of the Yellow River, one can be able to quantitatively retrieve the canopy VWC.
The VOD scheme considering the influence of waterbodies demonstrated superior to the one directly using the τ-ω model for the inversion.The VOD with considering the influence of waterbodies showed a large-value region occurring primarily in the area along the Yellow River and in the southeastern source region, where seasonal variation was dramatic.In particular, in the region from the west of Maqing to the east of Ngoring Lake, the VOD exhibited strong monthly variation.In the Maqu area, the maximum VOD of approximately 0.60 occurred in early-August.Due to low vegetation density in the source region of the Yellow River, the maximum VOD was less than 0.40.
The VWC in the entire source region was in the range of 0.20 to 1.40kg/m 2 , and the large-value area was mainly located along the Yellow River in the eastern source region.Large variation in the source region in the range of 0.20 to 0.50 kg/m 2 occurred in the Maqu area to the east, whereas in the

Summary and Conclusions
In this study, the τ-ω model was improved by using the dual-polarized brightness temperature observations at the K-band obtained from the MWRI microwave radiometer on board the FY-3B meteorological satellite and removing the influence of waterbodies.By constructing a microwave VOD model that is suitable for the source region of the Yellow River, one can be able to quantitatively retrieve the canopy VWC.
The VOD scheme considering the influence of waterbodies demonstrated superior to the one directly using the τ-ω model for the inversion.The VOD with considering the influence of waterbodies showed a large-value region occurring primarily in the area along the Yellow River and in the southeastern source region, where seasonal variation was dramatic.In particular, in the region from the west of Maqing to the east of Ngoring Lake, the VOD exhibited strong monthly variation.In the Maqu area, the maximum VOD of approximately 0.60 occurred in early-August.Due to low vegetation density in the source region of the Yellow River, the maximum VOD was less than 0.40.
The VWC in the entire source region was in the range of 0.20 to 1.40kg/m 2 , and the large-value area was mainly located along the Yellow River in the eastern source region.Large variation in the source region in the range of 0.20 to 0.50 kg/m 2 occurred in the Maqu area to the east, whereas in the western source region it did not show seasonal variation obviously, especially in the east of Ngoring Lake.
Although the VWC was derived by using the field observations, the number of samples is still insufficient.An analytic experiment about vegetation canopy samples with an even broader coverage range is necessary if these data are employed to derive the parameter b that is suitable for variation in vegetation coverage in the source region of the Yellow River to retrieve canopy VWC.
The error analysis of VWC indicated that the τ-ω model considering the influence of waterbody coverage was reasonable for derivation of vegetation canopy characteristics in the source region of the Yellow River.This study has shown that the brightness temperature data from the MWRI sensor on board the FY-3B could be used for land surface parameters retrieval in the source region of the Yellow River.

Figure 1 .Figure 1 .
Figure 1.Location of the source region of the Yellow River catchment (left) and Maqu soil moisture and soil temperature (SMST) monitoring network (right) [37].(Left: red five-stars represent the sampling stations: Maqu 87 times and Maduo 8 times.) Figure 1.Location of the source region of the Yellow River catchment (left) and Maqu soil moisture and soil temperature (SMST) monitoring network (right) [37].(Left: red five-stars represent the sampling stations: Maqu 87 times and Maduo 8 times.).

Figure 1 .Figure 2 .
Figure 1.Location of the source region of the Yellow River catchment (left) and Maqu soil moisture and soil temperature (SMST) monitoring network (right) [37].(Left: red five-stars represent the sampling stations: Maqu 87 times and Maduo 8 times.)

Figure 2 .
Figure 2. Land cover map in the source region of the Yellow River.

Figure 3 .
Figure 3. Variation of observed vegetation water content (VWC) during the experimental period.

Figure 3 .
Figure 3. Variation of observed vegetation water content (VWC) during the experimental period.

Figure 5 .
Figure 5. Spatial distribution of the fractional coverage of open waterbodies (fw) during summer of 2012.

Figure 5 .
Figure 5. Spatial distribution of the fractional coverage of open waterbodies (fw) during summer of 2012.

Figure 6 .
Figure 6.Vegetation optical depth (VOD) considering waterbody influence during the summer of 2012 in the source region of the Yellow River, China.

Figure 6 .
Figure 6.Vegetation optical depth (VOD) considering waterbody influence during the summer of 2012 in the source region of the Yellow River, China.

Figure 6 .Figure 7 .
Figure 6.Vegetation optical depth (VOD) considering waterbody influence during the summer of 2012 in the source region of the Yellow River, China.

Figure 7 .
Figure 7. VOD without waterbody influence during the summer of 2012 in the source region of the Yellow River, China.

Figure 8 .
Figure 8. Temporal distribution of VOD in the Yellow River source region during summer 2012.

Figure 8 .
Figure 8. Temporal distribution of VOD in the Yellow River source region during summer 2012.
Figure 9 indicates that both NDVI and MVI have a strong correlation with VOD, so they can be considered as characteristic variables for monitoring VOD variation during the vegetation growing season.Nevertheless, in the middle of source region NDVI reflects the VOD more reasonable variation than MVI does because of the impact of land cover.Remote Sens. 2019, 11, x FOR PEER REVIEW 14

Figure 9 .
Figure 9. Correlations between retrieved VOD against (a) NDVI or (b) MVI during the experimental period.

Figure 9 .
Figure 9. Correlations between retrieved VOD against (a) NDVI or (b) MVI during the experimental period.

Figure 9 .
Figure 9. Correlations between retrieved VOD against (a) NDVI or (b) MVI during the experimental period.

Figure 10 .
Figure 10.Relationship between measured VWC and microwave-retrieved VOD in the source region of the Yellow River.
Figure 10 (black line represents logarithm trend line) presents logarithmic function between VOD and VWC, there is a good corresponding relationship between VWC and VOD in the source region during the summer season.According to the Figure 10, the parameter b is determined by the logarithm trend line 1.5524 ln( ) 1.5566 y x = +.With this relationship the spatial VWC is obtained during the entire growth period (Figure11).

Figure 10 .
Figure 10.Relationship between measured VWC and microwave-retrieved VOD in the source region of the Yellow River.

Figure 9 .
Figure 9. Correlations between retrieved VOD against (a) NDVI or (b) MVI during the experimental period.

Figure 10 .
Figure 10.Relationship between measured VWC and microwave-retrieved VOD in the source region of the Yellow River.
Figure 10 (black line represents logarithm trend line) presents logarithmic function between VOD and VWC, there is a good corresponding relationship between VWC and VOD in the source region during the summer season.According to the Figure 10, the parameter b is determined by the logarithm trend line 1.5524 ln( ) 1.5566 y x = +.With this relationship the spatial VWC is obtained during the entire growth period (Figure11).

Figure 11 .
Figure 11.Spatial distribution of VWC during summer of 2012 in the Yellow River source region, China.

Figure 11 .
Figure 11.Spatial distribution of VWC during summer of 2012 in the Yellow River source region, China.

Figure 12 .
Figure 12.Comparison between the remote-sensing derived VWC and the measured VWC.

Figure 12 .
Figure 12.Comparison between the remote-sensing derived VWC and the measured VWC.
Spatial distribution of 16-day mean normalized difference vegetation index (NDVI) and microwave vegetation index (MVI).(a,b) NDVI from 23 April to 8 May and 13 to 28 July, 2012, respectively.(d,e) MVI from 23 April to 8 May and 12 to 28 July, 2012, respectively.