Improving the Seasonal Representation of ASCAT Soil Moisture and Vegetation Dynamics in a Temperate Climate

Previous validation studies have demonstrated the accuracy of the Metop-A ASCAT soil moisture (SM) product, although overand underestimation during different seasons of the year suggest a need for improving the retrieval algorithm. In this study, we analyzed whether adapting the vegetation characterization based on global parameters to regional conditions improves the seasonal representation of SM and vegetation optical depth (τ). SM and τ are retrieved from ASCAT using both a seasonal (mean climatological) and a dynamic vegetation characterization that allows for year-to-year changes. The retrieved SM and τ are compared with in situ and satellite SM, and with vegetation products (SMAP, AMSR2, and SPOT-VGT/PROBA-V). The study region is set in an agricultural area of Lower Austria that is characterized by heterogeneous land cover and topography, and features an experimental catchment equipped with a SM network (HOAL SoilNet). We found that a stronger vegetation correction within the SM retrieval improves the SM product considerably (increase of the Spearman correlation coefficient rs by 0.15 on average, and rs comparable to SMAP and AMSR2). The vegetation product derived with a dynamic vegetation characterization compares well to the reference datasets and reflects vegetation dynamics such as start and peak of season and harvest. Although some vegetation effects cannot be corrected by the adapted vegetation characterization, our results demonstrate the benefits of a parameterization optimized for regional conditions in this temperate climate zone.


Introduction
Soil moisture (SM) plays an important role in the water and carbon cycle and needs to be considered in a number of related applications.SM datasets are used in hydrological model calibration and runoff predictions [1][2][3], irrigation scheduling [4,5], rainfall estimation [6,7], drought monitoring [8][9][10], modeling of groundwater depletion [11], and vegetation and crop growth monitoring [12], amongst many others.These applications require accurate and readily available datasets on different scales.
Over the past 30 years, remote sensing missions have emerged that observe global SM conditions from space [13].The Soil Moisture Ocean Salinity (SMOS) [14] and Soil Moisture Active Passive (SMAP) [15] missions are dedicated SM missions operating in L-band, providing global passive SM observations.Another sensor that can be used for the retrieval of global SM, though not designed for this purpose, is the Advanced Microwave Scanning Radiometer 2 (AMSR2), a passive multi-frequency instrument on board GCOM-W1 [16].An active C-band sensor suitable for SM retrieval is the Advanced Scatterometer (ASCAT) on board the Metop satellites.
At the Vienna University of Technology (TU Wien), a change detection algorithm has been developed to retrieve SM from the scatterometers on board the European remote sensing satellites (ERS) [17], which was later adapted and improved for ASCAT [18].Many studies have evaluated the ASCAT SM product over different regions.High correlations between ASCAT and in situ data have been observed in the Bibeschbach catchment in Luxembourg [19], in southwestern France [20] and for selected networks across Europe [21].However, weaknesses in the representation of the seasonal cycle in ASCAT SM have been described: Wagner et al. [22] observed ASCAT SM values in the summer months that are consistently higher than SM values from SMOS and in situ stations in two watersheds in the United States.Three possible reasons were discussed which can lead to high SM in these catchments: Under very dry conditions, sub-surface scattering can increase the backscatter, and consequently the derived SM [12,23].On the other hand, wet soil surfaces and wetlands can also lead to enhanced backscatter.A third possible reason is a too weak vegetation correction.Barbu et al. [24] found that, over France, temporal correlations between the satellite data and modeled SM increase by 0.07-0.09when applying a seasonal correction to the SM values instead of a bias correction that is static throughout the year.This is due to low ASCAT SM values in May, which are adjusted by the seasonal correction.
The TU Wien SM model uses the multi-angle viewing capacity of the ASCAT sensor to correct for vegetation [25].It has been shown that the vegetation characterization can be converted into vegetation optical depth (τ a ) [26], which is an indicator of the vegetation water content, both the leaf and woody components of the total above-ground biomass [27].Similar to SM, vegetation is an important variable in the water and carbon cycle.It has been shown that τ from microwave remote sensing is suitable for continuously monitoring vegetation dynamics [27,28], and complements datasets from optical satellites such as the normalized difference vegetation index (NDVI) and leaf area index (LAI) [27].
The aim of this study was to improve both the ASCAT SM and τ retrievals by optimizing model parameters linked to the vegetation characterization.The representation of the seasonal cycle in the ASCAT SM and τ datasets is assessed against SM and τ from different satellite and in situ datasets.Since the seasonal bias was observed in temperate climate regions, a region in Lower Austria was chosen as a representative area for this analysis.Since 2013, an in situ SM network is operated in an agricultural catchment in this area.Section 2 gives an overview of the study area, followed by a description of datasets in Section 3. The applied methods are described in Section 4. Section 5 describes the results of the study for SM and τ, which are discussed in Section 6.

Study Site
The study site is situated in Petzenkirchen, Lower Austria (Figure 1).A Hydrological Open Air Laboratory (HOAL) is operated in an agricultural catchment, providing extensive field data used for investigating catchment evapotranspiration, runoff processes, sediment transport, nutrient dynamics, contaminant pathways and spatial patterns in SM [29].The area is characterized by a humid climate with higher precipitation in summer than in winter.From 1990 to 2014, a mean annual temperature of 9.5 °C and a mean annual rainfall of 823 mm year −1 have been observed.Eighty-seven percent of the total area of the HOAL is arable land; the remaining parts are forests (6%), pasture (5%) and paved areas (2%).Two main crop growing seasons are found in the HOAL: winter crops such as wheat, barley and rapeseed are usually planted in autumn and harvested in July.Summer crops (in the HOAL mainly corn) are usually planted in April and harvested in September/October.Between the harvest and seeding of the main crops, green fertilizers are often planted on the fields.Depending on the weather, seeding and harvest dates can vary by a few weeks from year to year.The HOAL catchment is classified as "Cropland, rainfed" by the ESA CCI land cover dataset [30], but it also features evergreen and deciduous forests and grasslands.Figure 2 (left) shows the area that is (approximately) covered by the ASCAT, AMSR2 and SMAP 3-dB footprints.This area is considerably larger than the catchment, but has similar topographic conditions and land cover as the HOAL.The ESA CCI land cover dataset for the ASCAT, SMAP and AMSR2 footprints is shown in Figure 3.The dominant land cover classes found in the study area are rainfed cropland, evergreen needle leaf and deciduous broad leaf forest as well as grasslands.More details about the study area are listed in Table 1.

In Situ Soil Moisture
Since 2013, an in situ SM network has been operated in the Hydrological Open Air Laboratory (HOAL), which measures SM at different depths (5, 10, 20, and 50 cm) every 30 min using the Time Domain Transmission (TDT) method.Twenty permanent stations are installed at selected locations that represent the different hydrologic conditions and land cover types present in the catchment.Additionally, 11 stations are installed temporarily inside agricultural fields, which are removed during farmers' field management practices (planting, harvesting, ploughing, etc.).For this study, the time series of the upper sensors (5 cm) of the in situ stations, averaged over the entire catchment, are used.The network is in the following referred to as HOAL SoilNet. Figure 2 (right) shows the distribution of the permanent and temporary SM stations in the catchment.

ASCAT
Backscatter measurements from the ASCAT sensor on board Metop-A are available since January 2007 and used to retrieve SM and τ a with the TU Wien model [17,18,25,26].A schematic overview of the semi-empirical TU Wien modeling scheme is provided in Figure 4 (top).Model parameters (Wagner et al., 1999;Naeimi et al., 2009;Vreugdenhil et al., 2016;Hahn et al., 2017) Wet conditions

Bare soil
Fully grown vegetation The observed backscattering coefficient σ • is assumed to depend only on the incidence angle (θ), degree of saturation (Θ s ), and vegetation cover (V).Wagner et al. [31] found that both a change in Θ s and in V lead to an increased σ • , and that the effect of V is typically more pronounced at high incidence angles (Figure 4, bottom panel).At two incidence angles, seasonal vegetation changes are assumed to not affect the observed σ • ; those are the so-called dry and wet cross-over angles (θ dry and θ wet ), which were empirically set to 25°and 40°globally [18].The incidence angle dependency of σ • is described by a second-order Taylor polynomial in the TU Wien model.The first and second derivatives of the σ • -θ-relationship are referred to as slope (σ ) and curvature (σ ).Those can be used to estimate the backscatter at any arbitrary incidence angle [25].In the TU Wien model, backscatter observations are normalized to a reference incidence angle of 40°to remove the incidence angle dependency [31]; this backscatter time series is termed σ • 40 .
At θ dry and θ wet , where, as mentioned above, the observed backscatter is assumed to be independent of V, the lowest and highest 10% of backscatter measurements are estimated and averaged to obtain the historically "driest" and "wettest" backscatter values.These two values are then transferred back to the reference incidence angle of 40°using σ and σ .Due to σ and σ varying from day to day, time series instead of single values for the driest and wettest backscatter values are obtained.These time series are termed the dry and wet references (σ • dry and σ • wet ).They describe the effects of different static and dynamic parameters (e.g., soil roughness, soil texture, and V) on the backscatter observations.Since soil roughness and soil texture are assumed to be static over the retrieval period, dynamics in σ • dry and σ • wet are an indicator of the vegetation state [32].The difference between σ • wet and σ • dry , also called the sensitivity, determines the dynamic range of the retrieved SM dataset.A small difference leads to a low dynamic range and consequently higher noise in the retrieval [18].
In a final step of the TU Wien model, σ • 40 is scaled between σ • dry and σ • wet to obtain a SM time series that is corrected for the effect of vegetation (see Equation ( 1)).
ASCAT τ (hereafter referred to as τ a ) is derived from σ • dry and σ • wet using a water cloud model [26].SM and τ a are retrieved at a spatial resolution of 25 km for every 1-3 days and sampled on a discrete global grid with regular 12.5 km point spacing within the TU Wien model [17,18,26].
Traditionally, the vegetation state is modeled using seasonal parameters, i.e., one value for each day of the year, and not accounting for year-to-year variations in the vegetation cycle (referred to as seasonal vegetation characterization VC seas ).In times of limited backscatter data availability, e.g., the scatterometers on board ERS-1 and ERS-2 (ESCAT), this was necessary to obtain robust estimates of the incidence angle dependency of the observed backscatter.Since the operation of the ASCAT sensor, which has a much higher temporal revisit rate than its predecessor ESCAT, this is no longer necessary.Melzer [33] introduced a Kernel smoother in order to obtain dynamic parameters not only for each day of the year, but for all years individually.This is referred to as dynamic characterization (VC dyn ) in the following.Hahn et al. [25] showed the robustness of the dynamic characterization globally and Vreugdenhil et al. [34] showed its sensitivity to inter-annual vegetation dynamics over Australia.

AMSR2
The passive radiometer AMSR2 on-board GCOM-W1 was launched by the Japan Aerospace Exploration Agency (JAXA) in 2012.Similar to its predecessor AMSR-E (Advanced Microwave Scanning Radiometer on the Earth Observing System), it observes brightness temperature at C-and X-band.In this study, SM and τ derived with the Land Parameter Retrieval Model (LPRM) algorithm (version LPRMv05) [16,35] were selected as reference datasets.To be comparable to the ASCAT datasets, only C-band observations are used in this study.The AMSR2 datasets are provided daily at a spatial resolution of 62 km × 35 km (C-band).
Kim et al. [36] compared AMSR2 SM derived with the LPRM algorithm to field measurements from COSMOS stations over the USA and found consistent temporal patterns of the two datasets, although AMSR2 tended to overestimate SM [37].Cho et al. [38] showed AMSR2 to provide a valuable successor to the AMSR-E mission in Australia.More on the validation of the LPRM AMSR2 SM product can be found, e.g., in [39][40][41].To the authors' knowledge, only a brief comparison of the LPRM AMSR2 τ (later referred to as τ a2 ) to other τ products has been carried out by Cui et al. [42]; τ derived from AMSR-E using the LPRM algorithm has however been evaluated against the widely used NDVI [27], applied in agricultural drought monitoring studies [43] and validated in the Sahel [28].Furthermore, it has been applied to the monitoring of global change in the total above ground vegetation water content and biomass over various ecosystems and to the attribution of observed changes to environmental and human drivers [44].

SMAP
The SMAP L-band radiometer measures brightness temperature with a revisit time of 1-3 days globally at a spatial resolution of approximately 40 km [15].The SM product is derived by inverting a τ-ω-model [45,46].Here, the SMAP L3 passive product, version V004, has been selected for a comparison with ASCAT datasets, which is available from 2015 onward.Since SMAP operates in L-band, it is assumed to penetrate deeper into the vegetation than higher frequency bands, and is less affected by radio frequency interference (RFI) [47].
A validation of the SMAP SM products was carried out by Colliander et al. [48], who found that, over 18 core validation sites, the passive product meets the target unbiased root mean square error of 0.04 m 3 /m 3 volumetric SM.A comprehensive assessment of the performance of the passive SM product can be found in [49].The SMAP V004 product includes estimates of vegetation water content (VWC), which is a combined estimate from NDVI (canopy water content) and past field observations and LAI (stem water content) [50].For the comparison with τ from ASCAT and AMSR2, SMAP vegetation opacity (VO; later referred to as τ s ) is used, which is equivalent to the VWC product multiplied by a constant scaling factor.To the authors' best knowledge, the τ s product has so far only been evaluated by Cui et al. [42]; however, an alternate approach for the retrieval of τ from SMAP observations using the multi-temporal dual channel algorithm has been applied and evaluated by Konings et al. [47].

SPOT-VGT and PROBA-V
SPOT-VGT (1999-May 2014) and its successor PROBA-V (June 2014-ongoing) are global vegetation monitoring missions operating in the optical domain [51].In this study, leaf area index (LAI) from the satellites is included as a reference for the evaluation of τ a .The CCI land cover dataset [30] is used to classify the SPOT-VGT/PROBA-V pixels and to obtain LAI time series for croplands, forests and grasslands.

Pre-Processing
In situ measurements from 20 permanent and 11 temporary sensors installed at 5 cm depth are averaged to a mean in situ SM time series.Frozen conditions and surface snow are masked using average Layer 1 soil temperature (0.00-0.10 m), average surface temperature and snow water equivalent data from the global land data assimilation system (GLDAS; version 2.1) [52].The same masking is applied to the satellite datasets.The soil penetration depth of active C-band systems such as ASCAT ranges from 0.5 and 2 cm under normal, not too dry conditions; passive systems observe the brightness temperature emissions which originate from the top 0-2 cm (e.g., C-band sensor AMSR2) to 3-5 cm of the soil (e.g., L-band sensor SMAP) [15,18,53].To make the satellite and in situ datasets comparable, a soil water index (SWI) [54] with a small time scale parameter (T value of two days) is applied to all satellite SM time series.Thus, the satellite SM values better comply with the measurement depths of the HOAL SoilNet stations (5 cm).

Methods
Over-or underestimation of SM during spring or summer are likely to be caused by an insufficient correction of vegetation effects on the backscattered signal.The vegetation correction applied in the TU Wien SM retrieval algorithm depends on two components: the choice of cross-over angles and the use of a seasonal or dynamic vegetation characterization.Similar to the vegetation correction needed for the retrieval of SM, τ a is derived from the dry and wet references.Thus, an improvement of the vegetation characterization should lead to an improvement of both the SM and τ a datasets.Figure 5 shows σ • dry and σ • wet for two selected pairs of incidence angles (10°/30°and 25°/40°) and both VC seas (top) and VC dyn (middle).Metop-A ASCAT backscatter is displayed after its normalization to an incidence angle of 40°(σ • 40 ).To enable a comparison of the time series with meteorological data, air temperature and rainfall and snowfall rates from GLDAS are displayed (bottom).Additionally, the backscatter normalized to an incidence angle of 40°(σ • 40 ) is displayed.(c) GLDAS air temperature, rainfall rate and snowfall rate for the study region.

Type of Vegetation Characterization
Originally, and in the latest released ASCAT SM product (H111), yearly climatologies are used for correction, i.e., every day of the year undergoes the same correction each year (seasonal vegetation correction, VC seas ; Figure 5a).The approach presented by Melzer [33], Vreugdenhil et al. [26] and Hahn et al. [25] included inter-annual variations in the vegetation dynamics (VC dyn ; Figure 5b).In this study, the effect of applying VC dyn instead of VC seas was analyzed.

Selection of Cross-Over Angles
Due to the transferring of the driest and wettest observed values from θ dry and θ wet back to the reference incidence angle of 40°, the amplitude of the vegetation characterization is determined by the choice of the dry and wet cross-over angles; lower cross-over angles lead to a larger amplitude of the vegetation characterization, i.e., the model assumes a stronger effect of vegetation (Figure 5).In the following, this is referred to as a stronger vegetation characterization.Wagner [55] found that the cross-over angle of dry surfaces is lower than for wet surfaces.Globally, the cross-over angles are set to the empirically determined values of 25°(dry soil) and 40°(wet soil; compare Figure 4, bottom) [18].To optimize the retrieval to regional conditions, the effect of modifying the cross-over angles was analyzed in this study.

Evaluation of the Results
The effect of changing the cross-over angles and selecting VC seas or VC dyn was analyzed for both SM (Section 5.1) and τ (Section 5.2).For SM, HOAL in situ data were used as a reference, as well as SM from AMSR2 and SMAP.ASCAT τ (τ a ) was compared to τ from AMSR2 (τ a2 ) and SMAP (τ s ) as well as to LAI retrieved from SPOT-VGT/PROBA-V observations.
For SM, the evaluation was based on the Spearman and Pearson correlation coefficients (r s and r p ) as well as on the unbiased root-mean-square deviation (ubRMSD).In the case of τ, the focus is given to the visual interpretation of the time series, since high r values are mainly associated with smooth τ curves that approximate the mean yearly cycle, but do not necessarily contain interesting features coming from different crop types (see Section 5.2).As a reference for the visual interpretation, photos of the HOAL catchment, information on the crop types, planting and harvesting dates as well as meteorological data (see Figure 5c) were used.
Although the analyzed spatial scales are considerably different (local in situ network vs. sensor footprints of multiple km 2 ), we consider the HOAL catchment appropriate for the comparison with satellite data due to the similar land cover types and topographic conditions described in Section 2.

Soil Moisture
In the following, retrievals of ASCAT SM using different settings for the above described parameters are compared to satellite and in situ SM.

Cross-Over Angle Optimization
Figure 6 shows r s between different versions of ASCAT SM and the reference datasets.r s was calculated for the available period of the respective datasets: 2015-2017 for SMAP, 2012-2017 for AMSR2, and 2013-2017 for the in situ data from the HOAL catchment.For all datasets and for both VC seas and VC dyn , higher correlation coefficients (up to 0.8) are obtained with a stronger vegetation correction, i.e., lower cross-over angles than the original pair of 25°/40°.Only when comparing with SMAP, r s increases with increasing cross-over angles when correlating it with dynamically corrected ASCAT SM (Figure 6b).All correlations are highly significant (p < 0.01).The sensitivity to SM was determined by the selection of cross-over angles.In general, lower cross-over angles (θ dry = 5 • -15 • , θ wet = 30 • ) lead to a lower dynamic range between σ • dry and σ • wet , because it is expected that vegetation has a stronger effect on the overall backscatter intensity.However, this leads to a higher noise in the SM retrieval, because the dynamic range of σ • dry and σ • wet defines how sensitive the (remaining) backscatter signal is to changes in soil moisture.Furthermore, low cross-over angles can create more backscatter outliers (σ • 40 < σ • dry and σ • 40 > σ • wet ), which are typically corrected to 0% (−25% > SM < 0%) and 100% (100% > SM < 125%).However, extreme outliers (i.e., SM < −25% or SM > 125%) are marked as invalid SM estimates and therefore set to NaN.Due to these characteristics of the algorithm, the cross-over angle pair of 10°and 30°is considered the optimal choice for the study region, although other cross-over angle pairs lead to higher correlation coefficients (Figure 6).In the following, all results are shown for the original and optimized cross-over angles pairs of 25°/40°and 10°/30°.

Quantitative Comparison
As stated by Wagner et al. [12], validation results of different types of datasets should also be interpreted in a relative context, for example by comparing different satellite datasets with the same in situ dataset.Therefore, r s and r p as well as the ubRMSD were calculated between each of the satellite datasets and the HOAL SoilNet time series (Figure 7).For each satellite-in situ pair, the time period for which both datasets are available was used; the time periods are limited by the availability of HOAL SoilNet data (from August 2013 onward) and SMAP data (from April 2015 onward).SMAP and AMSR2 show r s and r p of 0.77 and 0.78 (SMAP) and 0.71 and 0.69 (AMSR2) with the in situ data.For the respective periods, ASCAT and in situ SM show correlations ranging from 0.57 to 0.75, with higher correlation coefficients clearly associated with a stronger vegetation correction (lower cross-over angles).The ubRMSD between SMAP and AMSR2 and the in situ data is 0.034 and 0.038, and ranges from 0.036 to 0.048 for the different ASCAT datasets.Using cross-over angles of 10°and 30°leads to ASCAT datasets that have a smaller ubRMSD and that are more similar to SMAP and AMSR2.Figure 7 also shows that the usage of VC seas or VC dyn has a smaller effect on the metrics than the choice of cross-over angle pair.These results show that in the current SM product, the cross-over angles are not optimal and the applied vegetation correction thus not ideal for the study area.In a next step, r s has been calculated between the in situ and satellite datasets for different seasons (spring: March-June; summer: July-September; and autumn and winter: October-February). Figure 8 shows r s and the change obtained from using different cross-over angle pairs.In the spring months, r s between ASCAT and HOAL increases by 0.16 (VC seas ) and 0.12 (VC dyn ) when using lower cross-over angles, i.e., a stronger vegetation correction.In summer, only little (VC seas ) or even negative changes (VC dyn ) are achieved; in autumn, both a stronger VC seas and VC dyn improve ASCAT with respect to the HOAL SoilNet (+0.18 and +0.09, respectively).Over the entire period, the change of r s is positive (+0.15 on average) when applying a stronger vegetation correction, and the HOAL-ASCAT r s becomes comparable to the HOAL-AMSR2 and HOAL-SMAP r s .Whether the correlation improves or deteriorates in different seasons probably depends on several reasons.r s has been calculated over up to six years, all of which experienced different weather conditions.When looking at individual years, larger improvements are made in 2015 and 2017 (rather dry years) than in 2014 and 2016 (rather wet); in 2014 and 2016, r s also increases, but is significantly higher when using VC seas (not shown).Such differences between the years are also observed in the τ a time series (see Section 5.2).A possible reason could be different temperature and rainfall conditions, leading to differences in the vegetation growth that are not sufficiently represented by the vegetation correction.

Qualitative Comparison
Finally, we examined the time series obtained from different model parameters (Figure 9; a moving mean with a 14-day window has been applied on all datasets for better readability).In this figure, it can be seen that, in all years, a stronger vegetation correction for the SM retrieval (blue solid line) leads to a seasonal representation of SM in spring and summer that is closer to the reference datasets, compared to the original SM product (dotted black line).However, not the entire effect is corrected this way; ASCAT SM is still-in some years (e.g., 2014) considerably-lower in spring than the other datasets.This will be further discussed in Section 6.
Apart from the mostly positive effect on the seasonal representation, we observe a negative side effect that can occur in years with high vegetation activity.In the study area, 2015 was a very warm and dry year.When choosing VC dyn and cross-over angles of 10°and 30°, the effect of the vegetation correction may be so strong during summer that no SM can be retrieved anymore, leading to a data gap (see Figure 9b, solid black line in July-August 2015; also visible in Figure 5b).Another issue to be considered are the penetration depths of the C-and L-band sensors (upper few centimeters of the soil), compared to the installation depth of the in situ sensors (5 cm, i.e., deeper than satellite sensors penetrate under normal, not extremely dry conditions).Although this discrepancy is reduced by the application of the SWI with a T value of two days, there can still be effects such as a faster drying out of the soil that is only seen by satellites but not measured by the deeper in situ sensors.Consequently, lower SM values observed by satellite sensors do not necessarily mean erroneous retrievals, but can simply be caused by different reference depths.

Vegetation Optical Depth (τ)
To evaluate the performance of τ a in the study area, it is compared to τ s and τ a2 , and LAI from SPOT-VGT/PROVA-V, resampled to the ASCAT footprint.

Quantitative Comparison
In a first step, linear regressions of LAI and satellite τ datasets have been calculated (Figure 10).The steeper slopes (0.04-0.05) observed with τ a2 and τ s (right column) are in line with the findings of Lawrence et al. [56], who observed a slope of 0.05 between MODIS LAI and τ from SMOS.The slope of τ a and LAI is smaller, whereas larger slopes (0.02) for τ a are obtained when applying a VC seas (left column).Using a VC dyn (middle column) leads to smaller slopes and more outliers, especially when applying a stronger vegetation characterization.A r s of 0.79 and 0.75 is observed for τ time series with small amplitudes and no inter-annual variation (τ s , and τ a derived with a VC seas ).If inter-annual variations are included in a product, r s decreases significantly to values ranging from 0.29 to 0.55 (τ from AMSR2 and τ a derived with a VC dyn ), but corresponds to correlation values found by Jones et al. [57] between MODIS LAI and τ from AMSR-E in temperate climates.R 2 and r s suggest that VC dyn does not add value, but deteriorates the τ a product.However, effects caused by sensitivities to different vegetation characteristics of optical and microwave sensors, which lead to time lags and varying dynamics between the time series [27,58], are neglected or wrongly reflected by r s .Therefore, a qualitative analysis of τ and LAI is given in the following.

Qualitative Comparison
Figure 11 shows the different datasets for 2012-2017.In the left column, τ a is shown as yearly climatologies, i.e., no inter-annual variations [26]; the right column shows τ a calculated as described by Melzer [33] and Vreugdenhil et al. [34] (VC dyn ).The τ a time series retrieved with original (25°/40°) and modified (10°/30°) cross-over angles are displayed as dashed and solid lines, respectively.As can be seen in Figure 11, a change to lower cross-over angles leads to a larger range of values, i.e., a stronger vegetation signal.The value range of all τ and LAI time series is relatively similar over all years, but some inter-annual dynamics are visible, e.g., a variable number of main peaks and different timings of those peaks (Figure 11).This supports the implementation of dynamic parameters for the retrieval of τ a in the study area.The τ a time series shows significant inter-annual differences, which are likely to reflect different weather conditions.In the study area, 2015 was a warm and dry year, thus encouraging early vegetation growth and leading to a high peak in τ a .In contrast, 2016 was rather wet, with light and heavy rainfalls the whole year.This is reflected by a relatively flat τ a time series with some distinct ups and downs throughout the year.
For an interpretation of signal dynamics between years and during each growing season, we examine the reference datasets.Two main peaks occur around June and September in τ s , which does not include inter-annual variations, and in τ a2 in every year.Those two peaks coincide with the peak of season in winter crops (first peak) and summer crops (second peak) and their harvest shortly after.The LAI dataset also shows these two peaks, however not as pronounced in every year.In most of the years, the rise in spring occurs later, but increases more quickly than the microwave datasets.Recurring features can also be identified in the τ a time series: each year, τ a starts to rise in February, shows one peak in April/May, one peak around June/July and one peak in August/September, and then drops again in autumn.The two peaks in June/July and August/September coincide with those visible in the reference datasets, however slightly shifted in time compared to LAI.This is likely due to the sensitivity of the different variables to different vegetation characteristics; since LAI is more sensitive to the green canopy, it observes peaks earlier than τ, which is more sensitive to the water content in the vegetation, including the woody parts of the canopy [34,57,59].To test this for the datasets used in this study, correlations were calculated between the datasets after introducing negative and positive temporal shifts of 0-12 weeks.Figure 12 shows the results of the analysis, which confirm that highest r s values are obtained when applying a positive temporal shift of two weeks (τ s ) and 6-8 weeks (τ a , τ a2 ) to the LAI dataset.The smaller shift between LAI and τ s can be explained by the fact that τ s is based on NDVI, which has been found to have a smaller phase shift to LAI than τ in several land cover regions [57].Between τ a and τ a2 , r s is highest with no introduced time shift.Similar results are obtained for τ s , except when correlating it with τ a retrieved with cross-over angles 10 • and 30 • and using VC seas .However, all r s values obtained between those two time series with positive temporal shifts of 0-6 weeks lie within a range of 0.10, and the maximum r s value is not as distinct as in other dataset combinations.

Land Cover Effect
Land cover in the study region is mostly cropland, but it also features grasslands and forests (mostly evergreen and mixed evergreen/deciduous). Figure 13 shows the LAI time series of up to three selected only-forest (MF, DF, EF), only-cropland (C) and only-grassland (G) areas, as well as the τ a time series of the 25 km ASCAT footprint over the study region.The three selected cropland LAI time series (Figure 13, bottom left) show different peaks due to the planted crop type (winter crops vs. summer crops).As can be seen in the figure, dynamics in LAI and τ a correlate well, and, as expected, features from all selected land cover types can be found in the τ a time series of the study area.Differences in the absolute values of τ a between the years might be caused by different temperature and rainfall conditions.
Every year, a peak is visible in the τ a signal around April/May (Figure 11, right column).This peak is not present in LAI (Figures 11 and 13) and cannot be explained by a particular change in vegetation water content or biomass.To investigate a possible effect of remained snow or wet snow, we looked at photos taken in the HOAL catchment and analyzed meteorological data from GLDAS (air temperature, rainfall and snowfall).However, the GLDAS snowfall rate was zero in most years (see Figure 5c) and no snow was detected on the photos during the periods of interest.Snow can thus most likely be excluded as a reason for the visible peak in spring.A likely interpretation of the peak could be changes in the scattering behavior of the vegetation or bare surface due to structural changes during this period.The size of the peak varies over the years; in 2014, it is very high and wide, whereas, one year later, it is hardly present.The origin of those peaks and their amplitudes are still under investigation.Because we also see inconsistencies in spring in SM, this issue will be addressed in the Discussion Section.

Discussion
This study analyzed the effect of a regionally optimized parameterization of vegetation in the TU Wien algorithm for both the SM and τ a products derived from Metop-A ASCAT backscatter.The analysis of SM shows that, in the study area, a stronger vegetation correction than the globally optimal parameter is needed.Correlations with reference datasets increase, and seasonal differences such as an underestimation of SM in spring and an overestimation in summer are mitigated; however, they are not fully corrected.The application of VC seas or VC dyn does not affect the SM dataset as much as the choice of the cross-over angles.
To transfer these findings to other regions in temperate climates and test if the problems described in Section 1 are mitigated by the application of a stronger vegetation correction, correlations and the ubRMSD between ASCAT and in situ SM have been calculated for the Little River watershed (Georgia, United States; [22]), an area in southwest France [24] and a catchment in western Denmark (hydrological observatory HOBE; [60,61]).At all sites, r p and r s increase from 0.64 to 0.71 on average when applying a stronger vegetation correction; the average ubRMSD decreases from 0.031 to 0.028.Table 2 shows the metrics for every location.The metrics have been calculated for the period 2012-2017 (LR and southwest France) and 2012-2015 (HOBE), and periods with negative temperatures or snowfall have been masked out, as described in Section 3.3.When analyzing the τ a time series, the added value of a dynamic vegetation parameterization becomes obvious.Inter-annual differences due to different weather conditions are visible, and vegetation dynamics such as start and peak of season and harvest of different crop types are reflected in the signal.As for SM, discrepancies between ASCAT and the reference datasets are observed in spring.This discrepancy in spring, where the effect of vegetation is over-and SM underestimated, is currently being investigated.In the course of the SM retrieval, the dependency of σ • on the incidence angle θ is modeled as a second-order Taylor polynomial (see Section 3.2.1).The first component, i.e., the slope of the σ • -θ curve (σ ), decreases (in absolute values) with vegetation growth, since volume scattering from the vegetation increases σ • at larger incidence angles (Figure 4, bottom).Figure 14 shows the time series of σ .In addition to the main seasonal cycle of vegetation growth, we also see a peak in spring (Day of Year 85-130).This means that, during that time, the backscattered component is less dependent on the incidence angle.Behavior such as this is commonly associated with vegetation growth, but can in this case also be due to changes in the structure of the vegetation (e.g., crops growth changes from grass-like to vertical), or due to SM-induced alteration of σ .Structural differences would induce a change in the single scattering albedo (ω), i.e., in the scattering and absorption behavior of the vegetation layer.At the moment, ω is considered to be static over time in the TU Wien retrieval algorithm.The same assumption is made in the retrieval algorithm of τ a2 , where we also see peaks in April.Assuming that ω shows a seasonal cycle and the effect of considering this cycle in the ASCAT SM and τ a retrieval should be explored in further research.

Conclusions
We analyzed the retrieval of surface soil moisture (SM) and vegetation optical depth τ from Metop-A ASCAT over a region in Lower Austria.Both SM and τ have been retrieved from ASCAT backscatter values using a seasonal and a dynamic vegetation characterization, which takes into account year-to-year changes.Furthermore, the effect of adjusting parameters related to the strength of the vegetation characterization has been analyzed.We found that changing the empirically identified cross-over angles used for the global SM retrieval according to regional conditions improves the SM retrieval considerably (increase of r s by 0.15 on average, and r s and ubRMSD comparable to SMAP and AMSR2).The vegetation product derived with a dynamic vegetation characterization compares well to τ observed by AMSR2 and SMAP, and LAI from SPOT-VGT/PROBA-V.Vegetation dynamics such as start and peak of season and harvest are reflected in the signal, however partially shifted in time, which reflects the different sensitivities of LAI and τ.The results highlight the benefits of a parameterization optimized to regional conditions for an area in the temperate climate zone.

Figure 1 .
Figure 1.Location of the Hydrological Open Air Laboratory (HOAL) in Lower Austria.

Figure 2 .
Figure 2. Schematic overview of the 3-dB sensor footprints of ASCAT, SMAP and AMSR2 (left); and distribution of permanent and temporary stations in the HOAL catchment (right).Map data ©2018 Google.

Figure 4 .
Figure 4. Main steps of the TU Wien algorithm for the retrieval of SM and τ (top); and dependency of backscatter (σ • ) on incidence angle θ, soil moisture and vegetation (bottom); after[12].
Seasonal vegetation characterization.
Dynamic vegetation characterization.

Figure 5 .
Figure 5. Dry and wet reference derived from Metop-A ASCAT using two different pairs of cross-over angles (10°/30°and 25°/40°) and seasonal (a) and dynamic vegetation characterization (b).Additionally, the backscatter normalized to an incidence angle of 40°(σ • 40 ) is displayed.(c) GLDAS air temperature, rainfall rate and snowfall rate for the study region.

31 Figure 7 .
Figure 7. Spearman and Pearson correlation coefficients (r s and r p ) and ubRMSD between HOAL in situ SM and different satellite SM datasets.The metrics have been calculated for 20 August 2013 (limited by availability of HOAL SoilNet data) and 1 April 2015 (limited by availability of SMAP data) until the end of 2017.

Figure 8 .
Figure 8. Spearman correlation coefficient (r s ) between HOAL in situ SM and different satellite SM datasets for different seasons.

Figure 9 .
Figure9.SM from ASCAT (SM a ), AMSR2 (SM a2 ), SMAP (SM s ) and the in situ network (HOAL SM) for the years 2012 to 2017; SM a retrieved using: (a) VC seas ; and (b) VC dyn .A moving mean (window: 14 days) has been applied to all datasets for better readability.

Figure 10 .
Figure 10.Satellite τ as a function of SPOT-VGT/PROBA-V LAI; the color indicates the number of observations.The black line shows the linear regression.In each panel, the function of the linear regression, the coefficient of determination (R 2 ) and r s are given.

Figure 14 .
Figure 14.Average annual behavior of the slope of the σ • -θ relationship.The peak in spring is highlighted in red.

Table 1 .
Study area characteristics.

Table 2 .
Metrics between different ASCAT SM versions and three in situ SM datasets.