Remote Sensing of Complex Permittivity and Penetration Depth of Soils Using P-Band SAR Polarimetry

A P-band SAR moisture estimation method is introduced for complex soil permittivity and penetration depth estimation using fully polarimetric P-band SAR signals. This method combines eigenand model-based decomposition techniques for separation of the total backscattering signal into three scattering components (soil, dihedral, and volume). The incorporation of a soil scattering model allows for the first time the estimation of complex soil permittivity and permittivity-based penetration depth. The proposed method needs no prior assumptions on land cover characteristics and is applicable to a variety of vegetation types. The technique is demonstrated for airborne P-band SAR measurements acquired during the AirMOSS campaign (2012–2015). The estimated complex permittivity agrees well with climate and soil conditions at different monitoring sites. Based on frequency and permittivity, P-band penetration depths vary from 5 cm to 35 cm. This value range is in accordance with previous studies in the literature. Comparison of the results is challenging due to the sparsity of vertical soil in situ sampling. It was found that the disagreement between in situ measurements and SAR-based estimates originates from the discrepancy between the in situ measuring depth of the top-soil layer (0–5 cm) and the median penetration depth of the P-band waves


Introduction
Soil moisture is one of the essential climate variables (ECVs) [1] "that critically contribute [ . . . ] to the characterization of [the] Earth's climate" [2] and its changes. Hence, it plays a crucial role within the hydrological and biogeochemical cycles [1,[3][4][5][6][7]. Many approaches exist to retrieve soil moisture from air-or space-borne remote sensing observations, taking advantage of the sensitivity of active (radar) and/or passive (radiometer) microwave signals to soil moisture [8][9][10][11][12]. These studies allow an estimation of soil moisture predominantly from L-band measurements: on the one hand, because of the larger penetration capabilities compared to higher frequency bands, such as C-or X-band, and on the other hand, due to the availability of global L-band satellite data, e.g., from Soil of [40,41,47], concluding that at P-band, influences of soil surface roughness and vegetation canopy are diminished. Nevertheless, numerical results from [48] revealed limited influence of small scale roughness on P-band backscatter simulations, but an increasing complexity between P-band signals and soil roughness at spatially larger scales.
In the literature, P-band is assumed to be able to estimate soil moisture across the rootzone of the soil (up to approx. 2 m depth) [5,49,50]. Consequently, a wide range of P-band penetration and sensing depths can be found in the literature, from some centimeters up to one meter or more [5,39,46,51]. Following the assumption that soil depths are accessible "in the order of one half to one tenth of the [employed] wavelength [52]" [45,47,53], P-band penetration depths of 6.97 to 34.86 cm at 430 MHz (λ = 69.72 cm) are potentially realistic. This value range was estimated, e.g., by [46] in a ten-year study based on P-band data near the Ameriflux site at Vaira Ranch, CA, USA. Additionally, within the AirMOSS P-band rootzone, soil moisture (RZSM) product sensing depths of typically 45 cm are considered [51]. Here, RZSM is estimated based on AirMOSS measurements together with simulations from a vegetation, surface, and subsurface scattering model by assuming a second-order polynomial function of the profile [54].
In this study, an adapted non-iterative hybrid P-band SAR decomposition method for moisture estimation is proposed to separate scattering mechanisms in fully polarimetric SAR data (Section 3.2). The proposed method utilizes a soil scattering model, which is suitable for P-band wavelengths (Section 3.1). The decomposition results are used in a consecutive step for the estimation of the complex permittivity, which in turn is applied for penetration depth calculations (Section 3.3). Thus, this study seeks to estimate complex permittivity from decomposed SAR signals and to determine permittivity-based P-band penetration capabilities to better understand P-band microwave behavior in soil along depth. For that, the proposed P-band SAR moisture estimation method continues the efforts of [3,24], and is the first of its kind, to the best of our knowledge, for complex soil permittivity estimation from (P-band) SAR data.

Materials
The polarimetric P-band SAR observations employed in this study were acquired during the AirMOSS campaign (a NASA Earth Venture-1 project) from 2012 to 2015 over nine different biomes across ten sites in Northern and Central America. This campaign was the first P-band airborne mission designed to estimate RZSM. The monitoring sites comprise varying land cover classes from bare soils to tropical rain forests. Each site covers an area of approximately 25 × 100 km at~100 m spatial resolution and was revisited at least two to four times every year during the campaign. The P-band instrument operated at a center frequency of 430 MHz (λ ≈ 69.72 cm). With a high radiometric calibration accuracy (0.5 dB) and a noise equivalent σ o of −40 dB, the AirMOSS dataset provides P-band SAR measurements, which are very well suited for estimating soil and vegetation parameters [3]. Detailed information on the campaign can be found in, e.g., [3,16].
The AirMOSS dataset provides measurements for a total of 168 dates at ten sites, each of them covering an area of 25 × 100 km. Hence, in this study, only a subset of the dataset was processed and analyzed, with focus on the three AirMOSS monitoring sites "MOISST" in Oklahoma, "Walnut Gulch" in Arizona, and "Harvard Forest" in Massachusetts, USA. By choosing these three monitoring sites, varying vegetation types and climatic conditions are analyzed. Furthermore, every study area is filtered for pixels classified as water, snow, wetlands, developed ground, bare soil, and pasture/hay, as well as for incidence angles smaller than 30 • and greater than 50 • , to narrow down the amount of data to be processed and to focus on vegetated soils.
Lastly, only four dates per year and monitoring site are evaluated, in the case data for more than four dates are available ( Figure 1). Walnut Gulch in Arizona is characterized by an arid to semi-arid climate of a hot desert (BWh), MOISST in Oklahoma by a temperate climate with hot summers (Cfa), and Harvard Forest in Massachusetts by a cold humid continental climate with warm summers (Dfb) [55]. as for incidence angles smaller than 30° and greater than 50°, to narrow down the amount of data to be processed and to focus on vegetated soils.
Lastly, only four dates per year and monitoring site are evaluated, in the case data for more than four dates are available ( Figure 1). Walnut Gulch in Arizona is characterized by an arid to semi-arid climate of a hot desert (BWh), MOISST in Oklahoma by a temperate climate with hot summers (Cfa), and Harvard Forest in Massachusetts by a cold humid continental climate with warm summers (Dfb) [55]. In Figure 2, the land cover types according to the National Land Cover Database (NLCD) [56], as provided within the AirMOSS datasets, are shown for every monitoring site. The percentages in Table 1 represent the respective number of pixels per land cover class in comparison to all processed pixels (after filtering). Walnut Gulch (Figure 2A) is mostly covered by shrub/scrub (92.5%) while other classes represent less than 5% of all pixels ( Table 1). The monitoring site MOISST ( Figure 2B) is mainly covered by grassland/herbaceous (51.5%) and cultivated crops (38.3%), while Harvard Forest ( Figure  2C) is covered by deciduous (49%), mixed (30.4%), and evergreen (25.1%) forests (Table  1). In Figure 2, the land cover types according to the National Land Cover Database (NLCD) [56], as provided within the AirMOSS datasets, are shown for every monitoring site. The percentages in Table 1 represent the respective number of pixels per land cover class in comparison to all processed pixels (after filtering). Walnut Gulch (Figure 2A) is mostly covered by shrub/scrub (92.5%) while other classes represent less than 5% of all pixels ( Table 1). The monitoring site MOISST ( Figure 2B) is mainly covered by grassland/herbaceous (51.5%) and cultivated crops (38.3%), while Harvard Forest ( Figure 2C) is covered by deciduous (49%), mixed (30.4%), and evergreen (25.1%) forests (Table 1)     For a later comparison of SAR-based estimates, the daily averages of in situ soil moisture measurements from five different networks at 17 stations are used as shown in Figure 2. At Walnut Gulch, five stations from the FLUXNET2015 dataset [57] and four stations from the Cosmic-ray Soil Moisture Observing System (COSMOS) [58,59] are employed. At MOISST, one station from the Soil moisture Sensing Controller And oPtimal Estimator (SoilSCAPE) project [60], and COSMOS, two stations from the US Climate Reference Network (US-CRN) [61], and three stations from the Plate Boundary Observatory to study the water cycle (PBO H 2 O) project [62] are used. Lastly, at Harvard Forest, one station from COSMOS is employed ( Figure 2).

Methods
In this study, the proposed P-band SAR moisture estimation method consists of a hybrid decomposition method and a subsequent moisture estimation procedure. The hybrid decomposition method combines eigen-and model-based decomposition techniques to separate P-band SAR signals into individual scattering mechanisms. For the determination of the scattering components, the method is supported by a soil scattering model, the first-order multi-layer small perturbation method (SPM) [63]. Subsequently, based on the decomposed scattering components, the complex soil permittivity is estimated. The individual processing steps of the method are shown in Figure 3 and explained in detail in the following.
Remote Sens. 2022, 14, x FOR PEER REVIEW 7 of 26 model, the first-order multi-layer small perturbation method (SPM) [63]. Subsequently, based on the decomposed scattering components, the complex soil permittivity is estimated. The individual processing steps of the method are shown in Figure 3 and explained in detail in the following.

Soil Scattering Based on Multi-Layer SPM
The first-order multi-layer SPM [63] is employed to model the scattering angle , needed within the proposed hybrid decomposition method (cf. Section 3.2). As shown on the right side of Figure 3, modeled complex backscatter channels,

Soil Scattering Based on Multi-Layer SPM
The first-order multi-layer SPM [63] is employed to model the scattering angle α Model s , needed within the proposed hybrid decomposition method (cf. Section 3.2).
As shown on the right side of Figure 3, modeled complex backscatter channels, denoted by σ pp [-] for horizontal or vertical polarization, are required to estimate the model-based α Model s . In this study, σ pp are modeled with the first-order solution of the multi-layer SPM [63], which computes backscatter coefficients based on soil characteristics for multiple soil layers with depth. The ability of this model to consider backscattering from multiple subsurface layers qualifies it to be more suitable for analyzing P-band soil interactions than a soil surface-only scattering model such as the SPM, also called Bragg model [27,64], as shown in [63]. The multi-layer SPM computes first-order scattering from layered surfaces by considering "multiple scattering processes between the boundaries" [63]. Since the model was originally designed to compute backscattering coefficients σ o pp , it was adapted in order to simulate complex backscatter channels σ pp , where σ o pp = σ pp 2 holds. The complex backscatter channels σ pp are estimated by: pp k s ⊥ as coefficients related to the incident field, W f 1 (k ⊥ ) and W f 2 (k ⊥ ) as the power spectral densities of the rough boundaries, and with W f 1 f 2 (k ⊥ ) as their joint spectral density [63]. In detail, if the incident electric field is p-polarized, the first-order p-polarized scattered electric field in the direction k s is a linear function of the Fourier transforms of the rough boundary functions f 1 and f 2 , with the corresponding coefficients denoted by α f 1 pp k s ⊥ and α f 2 pp k s ⊥ . Lastly, the square root of a complex number is calculated according to Moivre's theorem [65]. The reader is referred to [63] for detailed information on backscatter coefficient calculations.
In Table 2, the required input parameters of the multi-layer SPM for simulating σ pp are listed together with their respective values used in this study. Please note that the chosen values are one possible approximation for the acquisition scenario in this study and may need to be adapted for different scenarios and other studies.  40], ε s ∈ [0, 10] ε s 2 = ε s 1 + (10 + j0.5) Table 3. Surface roughness parameter sets of each layer (vertical RMS height s, horizontal correlation length l) as input for the multi-layer SPM based on TanDEM-X derived roughness indicator R TDX .

Roughness Indicator from TanDEM-X [m] Input Roughness Parameters [cm]
R TDX < 5 s 1 = 0.5; l 1 = 30; s 2 = 0.25; l 2 = 60 5 ≤ R TDX < 10 s 1 = 1.5; l 1 = 25; s 2 = 0.75; l 2 = 50 10 ≤ R TDX < 15 Following findings by [48] (cf. Section 1), the surface roughness parameters, required for modeling P-band backscatter channels, are not fixed to static values in this study. However, since no information on surface roughness is available for AirMOSS monitoring sites with comparable spatial resolution, the DLR TanDEM-X DEM at 90 m resolution [66], downloadable at https://download.geoservice.dlr.de/TDM90/ (accessed on 5 March 2021), is used to derive a first-order roughness indicator for each site. For this, TanDEM-X elevations were resampled to the AirMOSS resolution, converted to roughness values with the GDAL DEM utility algorithm [67] in QGIS© [68], and scaled according to the employed wavelength (λ = 69.72 cm). The resulting roughness values, giving the degree of irregularity of the surface, serve as roughness indicators in this study. Hence, depending on this roughness indicator R TDX (Table 3, left column), typical surface roughness parameter sets for each layer (Table 3, right column) were fixed as model input to account for varying roughness during simulations (from smooth to rather rough). It can be seen that the chosen values for vertical RMS height s vary from 0.25 to 3 cm, and for horizontal correlation length l from 15 to 60 cm. These values were fixed based on several sensitivity analyses and after reviewing the literature.

Polarimetric Hybrid Decomposition Method
The hybrid decomposition method was originally introduced by [24] as an iterative approach (cf. Section 1) specifically designed for L-band data. However, it is adapted here to be applied in a non-iterative way to P-band observations, and for estimating complex permittivity. As illustrated in Figure 3, the hybrid decomposition method separates polarimetric SAR observations into individual scattering components (surface, dihedral, and volume).
First, the polarimetric coherency [T] matrix is defined as: The rank-3, model-based vegetation volume component [T v ] is described by: with f v as volume scattering intensity, and A p [-] as particle anisotropy, which describes the shape of vegetation volume scatterers from vertical dipoles (A p = 0) and spheres (A p = 1) to horizontal dipoles (A p = ∞). f v is given by: and the polarimetric scattering angle α s as model-based α Model s [24,31] (cf. Section 3.1). Further, the parameters for estimating the vegetation volume component are given by: with ∆ψ [ • ] as width of the orientation angle distribution, which describes the degree of orientation of the vegetation volume from oriented (∆ψ = 0 • ) to random (∆ψ = 90 • ). Hence, (4) can be rearranged together with (5) in order to determine the soil parameters from the individual soil scattering components ([T s ], [T d ]), by subtracting the volume component: As already mentioned, [24] originally proposed an iterative estimation procedure to find the most suited volume descriptions, namely f v , A p , and ∆ψ. However, since this approach is computationally expensive and needs certain assumptions (bare soil regions in the observation area) to be initialized, a non-iterative way for estimating f v is proposed here, independent of bare areas, by employing the model-based α Model s and a realistic parameter space for vegetation volume parameters, with Ap ∈ [0, 1] and ∆ψ ∈ [0 • , 90 • ].
Thus, together with α Model s , the SAR measurements are decomposed into the two eigen-based scattering angles α s and α d [ • ], as well as the scattering intensities of surface f s , dihedral f d , and volume f v [-] for multiple Ap − ∆ψ combinations. The ability of the method to estimate five variables (α s , α d , f s , f d , and f v ) out of four SAR observations (|S HH | 2 , |S VV | 2 , |S HV | 2 , |S V H | 2 ) is possible, since the ambiguity for α s and α d is solved by assuming an orthogonality condition with α s = π/2 − α d [27], meaning, if the decomposed polarimetric scattering angle α is smaller than or equal to 45 • , surface scattering is assumed, denoted by α s , whereas in the case α is greater than 45 • , dihedral scattering is assumed, denoted by α d [24]. This way, in contrast to model-based-only approaches, the proposed method needs no prior assumptions on vegetation characteristics, and provides the advantage of incorporating actual SAR measurements.
The surface [T s ] and dihedral [T d ] scattering components are defined in (12) and (13), respectively, assuming that both are orthogonal rank-1 components [3,24]: Finally, (11) can be reformulated to: with the surface scattering intensities given by [24]: and the eigen-based scattering angle α to separate contributions from surface and dihedral scattering, given by [24]: In summary, from polarimetric SAR observations and simulated α Model s , the decomposition results for α s , α d , f s , f d , and f v can be estimated for varying ε s and individual Ap − ∆ψ combinations. Hence, the proposed method overcomes the computationally expensive iterative procedure as proposed in [24], and provides the possibility of including multiple vegetation characteristics (Ap, ∆ψ) as suggested by [3] (cf. Section 1). In this study, not all Ap − ∆ψ combinations are applied. Firstly, to ensure non-negative decomposed powers, scattering intensities for any Ap − ∆ψ combination that are smaller than zero are eliminated for subsequent analyses in order to avoid possible overestimation of the volume scattering component (cf. Section 1.). Secondly, the radar vegetation index (RV I) is used to exclude unrealistic Ap − ∆ψ combinations. Here, the RV I model is modeled according to [31,69] for all Ap − ∆ψ combinations used within the approach. Afterwards, the data-based RV I data is calculated for every pixel based on the respective SAR backscatter observations. Finally, all Ap − ∆ψ combinations, where the RV I model deviates from RV I data around the average value (RV I model ) of all RV I model (RV I model / ∈ RV I data ± RV I model ), are excluded from further analyses.

Complex Permittivity Estimation
The P-band SAR moisture estimation method proposed in this study provides the advantage of complex permittivity estimation. The real part of the permittivity ε s is associated with energy or heat storage, and the imaginary part of the permittivity ε s Remote Sens. 2022, 14, 2755 11 of 25 is associated with energy or heat loss, often called dielectric loss factor [35,70]. Since model-based α Model s is used to estimate f v during the decomposition of polarimetric SAR observations, decomposition results can be calculated for a realistic range of soil permittivity values. Subsequently, by determining the closest fit between decomposed data-based α Data s and model-based α Model s for individual Ap − ∆ψ combinations (cf. Section 3.2.), the corresponding ε s value can be estimated.
For that, the absolute value of the modulus r [-] and phase angle φ [-] of the real (α s ) and imaginary (α s ) parts of α Data s and α Model s , respectively, are calculated: Afterwards, the smallest sum of the absolute differences between data-and modelbased r α s and φ α s is used to determine ε s : Since complex ε s are estimated for multiple Ap − ∆ψ combinations to cover realistic vegetation volume characteristics (cf. Section 3.2), the final complex ε s is the average value of ε s (Ap, ∆ψ) for all remaining realistic Ap − ∆ψ combinations (non-negative surface scattering powers, modeled RV I model around data-based RV I data ) (cf. Section 3.2).
For comparison, permittivity can be converted to soil moisture and vice versa, by employing dielectric mixing models, such as the one by Topp et al., [71], Dobson et al., [72], or Mironov et al., [73]. An extended review of dielectric mixing models for soils can be found in [74]. In this study, ε s results are converted to soil moisture θ [vol.%] according to the dielectric mixing model of [71].
Based on estimated complex ε s from decomposed polarimetric SAR observations, one possible application is the estimation of the penetration depth δ p [cm]. In this study, the well-known formulation of [75] is adapted to calculate δ p [44]: This formulation was chosen because almost all available equations to estimate δ p are comparable and mainly dependent on the attenuation factor α, a measure for the penetration depth of an electromagnetic wave in a medium. However, this is one of the first to be published and well known.
Since δ p is only dependent on complex permittivity and wavelength (frequency) (cf. Equation (20)), it is considered in this study as permittivity-based penetration depth to analyze potential depths for P-band signals. For a better understanding of how δ p behaves across varying permittivity levels at P-band (430 MHz), modeling results are displayed in Figure 4. It can be seen that δ p decreases with increasing permittivity, as expected. Further, δ p estimates significantly decrease with increasing ε s , proving the importance of ε s for δ p calculations.

Since
is only dependent on complex permittivity and wavelength (frequency) (cf. Equation (20)), it is considered in this study as permittivity-based penetration depth to analyze potential depths for P-band signals. For a better understanding of how behaves across varying permittivity levels at P-band (430 MHz), modeling results are displayed in Figure 4. It can be seen that decreases with increasing permittivity, as expected. Further, estimates significantly decrease with increasing ′′ , proving the importance of ′′ for calculations. Equation (20)).

Results
In this section, the results of the P-band SAR moisture estimation method are presented, comprising filtered and decomposed SAR signals, estimated complex permittivity, and soil moisture, as well as permittivity-based penetration depths. Furthermore, the results are related to in situ measurements from various networks (cf. Section 2) to compare the complex permittivity estimates.  (20)).

Results
In this section, the results of the P-band SAR moisture estimation method are presented, comprising filtered and decomposed SAR signals, estimated complex permittivity, and soil moisture, as well as permittivity-based penetration depths. Furthermore, the results are related to in situ measurements from various networks (cf. Section 2) to compare the complex permittivity estimates.

Decomposition Results
Based on the proposed hybrid decomposition method, the relative indices P s /P t , P d /P t , and P v /P t for soil, dihedral, and volume are estimated to analyze the normalized contribution of each scattering mechanism. The indices are calculated based on the decomposition results for scattering angles and intensities according to [3].
In Figure 5, only the dominant scattering mechanisms, characterized by their relative indices, are shown for each AirMOSS monitoring site. Apparent data gaps (white areas in Figure 5) in the results originate from filtering for specific land cover classes and incidence angles prior to the processing (cf. Section 2). Here, brown regions represent dominant soil scattering, meaning the normalized contribution of soil scattering is the highest of all three contributions. Accordingly, beige regions represent dominant dihedral scattering, while dark green regions represent dominant volume scattering.
It can be seen that Walnut Gulch is mainly characterized by volume scattering, with dihedral scattering in the eastern part of the site, where cultivated crops are the main land cover class, as well as soil and dihedral scattering in the western part of the site near the Green Valley south of Tucson (Figure 2, Section 2). In this particular region, the influence of the incidence angle is evident since dominant volume scattering is mixed with dominant soil scattering, despite the same land cover class (shrublands). Here, the incidence angles vary across the entire possible range, from 30 • to 50 • , and the terrain heights are lowest, which leads to dominant soil scattering components instead of dominant volume scattering contributions ( Figure 5A). At the monitoring site MOISST, dihedral scattering is dominant in the eastern part around the city Stillwater ( Figure 5B), where the land is predominantly covered by grassland/herbaceous and smaller evergreen forests (Figure 2, Section 2). In contrast, regions mainly characterized by cultivated crops show dominant volume scattering, e.g., in the center of the site. Again, the results are influenced by the varying incidence angles, for instance, in the center of the site, where higher incidence angles lead to rather dominant soil scattering ( Figure 5B). Lastly, Harvard Forest shows primarily dihedral scattering with only some regions in the northern part of the monitoring site, where soil or volume scattering are dominant ( Figure 5C).
/ , and / for soil, dihedral, and volume are estimated to analyze the normalized contribution of each scattering mechanism. The indices are calculated based on the decomposition results for scattering angles and intensities according to [3].
In Figure 5, only the dominant scattering mechanisms, characterized by their relative indices, are shown for each AirMOSS monitoring site. Apparent data gaps (white areas in Figure 5) in the results originate from filtering for specific land cover classes and incidence angles prior to the processing (cf. Section 2). Here, brown regions represent dominant soil scattering, meaning the normalized contribution of soil scattering is the highest of all three contributions. Accordingly, beige regions represent dominant dihedral scattering, while dark green regions represent dominant volume scattering. It can be seen that Walnut Gulch is mainly characterized by volume scattering, with dihedral scattering in the eastern part of the site, where cultivated crops are the main land cover class, as well as soil and dihedral scattering in the western part of the site near the Green Valley south of Tucson (Figure 2, Section 2). In this particular region, the influence of the incidence angle is evident since dominant volume scattering is mixed with dominant soil scattering, despite the same land cover class (shrublands). Here, the incidence angles vary across the entire possible range, from 30° to 50°, and the terrain heights are lowest, which leads to dominant soil scattering components instead of dominant volume scattering contributions ( Figure 5A). At the monitoring site MOISST, dihedral scattering is dominant in the eastern part around the city Stillwater ( Figure 5B), where the land is predominantly covered by grassland/herbaceous and smaller evergreen forests (Figure 2, Section 2). In contrast, regions mainly characterized by cultivated crops show dominant volume scattering, e.g., in the center of the site. Again, the results are influenced by the varying incidence angles, for instance, in the center of the site, where higher incidence angles lead to rather dominant soil scattering ( Figure 5B). Lastly, Harvard Forest shows primarily dihedral scattering with only some regions in the northern part of the monitoring site, where soil or volume scattering are dominant ( Figure 5C). For a comparison of all decomposition results for varying years, Figure 6 displays the average probability distribution function (PDF) of the normalized scattering contributions at each monitoring site. It can be seen that the PDF for soil scattering peaks at the lowest value but with the highest probability at Harvard Forest, while at Walnut Gulch and MOISST the normalized soil scattering is located between dihedral and volume scattering as second dominant scattering contribution, respectively ( Figure 6, first row).
The dihedral scattering contribution is lowest at Walnut Gulch and MOISST but highest at Harvard Forest ( Figure 6, second row), while the volume scattering contribution is highest at Walnut Gulch but lowest at Harvard Forest, with MOISST in between ( Figure 6, third row). Moreover, there are hardly any differences in the estimated PDFs of all three relative indices between the years 2013 and 2015. Only the PDF for soil scattering contribution at Harvard Forest slightly decreases from 2013 to 2015, whereas the PDF for dihedral and volume scattering slightly increases ( Figure 6, right column).
In summary, the monitoring site mainly covered by shrub/scrub (Walnut Gulch) displays the lowest dihedral and the highest volume scattering, while the monitoring site covered by forests (Harvard Forest) displays almost exclusively dihedral and the lowest vegetation scattering. The mixed monitoring site covered by grassland/herbaceous or cultivated crops (MOISST) shows almost balanced contributions of soil or volume scattering with dihedral scattering being predominant in regions of grassland/herbaceous and volume scattering in regions with extensive agriculture.
For a comparison of all decomposition results for varying years, Figure 6 displays the average probability distribution function (PDF) of the normalized scattering contributions at each monitoring site. It can be seen that the PDF for soil scattering peaks at the lowest value but with the highest probability at Harvard Forest, while at Walnut Gulch and MOISST the normalized soil scattering is located between dihedral and volume scattering as second dominant scattering contribution, respectively ( Figure 6, first row). The dihedral scattering contribution is lowest at Walnut Gulch and MOISST but highest at Harvard Forest ( Figure 6, second row), while the volume scattering contribution is highest at Walnut Gulch but lowest at Harvard Forest, with MOISST in between ( Figure  6, third row). Moreover, there are hardly any differences in the estimated PDFs of all three relative indices between the years 2013 and 2015. Only the PDF for soil scattering contribution at Harvard Forest slightly decreases from 2013 to 2015, whereas the PDF for dihedral and volume scattering slightly increases ( Figure 6, right column).
In summary, the monitoring site mainly covered by shrub/scrub (Walnut Gulch) displays the lowest dihedral and the highest volume scattering, while the monitoring site covered by forests (Harvard Forest) displays almost exclusively dihedral and the lowest vegetation scattering. The mixed monitoring site covered by grassland/herbaceous or cultivated crops (MOISST) shows almost balanced contributions of soil or volume scattering with dihedral scattering being predominant in regions of grassland/herbaceous and volume scattering in regions with extensive agriculture.

Results for Complex Soil Permittivity and P-Band Penetration Depth
In Figure 7, the real ′ and imaginary ′′ part of the estimated complex soil permittivity , determined with the proposed P-band SAR moisture estimation method

Results for Complex Soil Permittivity and P-Band Penetration Depth
In Figure 7, the real ε s and imaginary ε s part of the estimated complex soil permittivity ε s , determined with the proposed P-band SAR moisture estimation method (cf. Section 3.3), are displayed for the main land cover classes for each monitoring site and date (Table 1). Walnut Gulch, covered by shrub/scrub (cf. Section 2), shows the overall lowest permittivity results with a median value of ε S = 11.2 − j1.5. At the monitoring site MOISST, the main land cover class grassland/herbaceous shows slightly higher permittivity results with a median value of ε S = 14.9 − j1.7. The latter is in contrast to the second dominant land cover class cultivated crops (cf. Section 2) with a median value of ε S = 13.8 − j1.5. However, the results for cultivated crops span the largest value ranges of all classes and monitoring sites with the real part of permittivity varying from 6 to 31 and the imaginary part of permittivity varying from 1 to 3.2. The three main land cover classes, deciduous, mixed, and evergreen forest (cf. Section 2), at Harvard Forest display the overall smallest value ranges and the highest permittivity results with comparable median values of~ε S = 16 − j1.8. Hence, the forest site Harvard Forest, characterized by a cold humid continental climate (cf. Section 2), shows the highest permittivity with the smallest deviations, while the driest monitoring site, mainly covered by shrub/scrub due to the arid climate (Walnut Gulch), shows the lowest permittivity results. The monitoring site characterized by a temperate climate and grassland/herbaceous or cultivated crops (MOISST) displays permittivity results in between.  Walnut Gulch, covered by shrub/scrub (cf. Section 2), shows the overall lowest permittivity results with a median value of = 11.2 − 1.5. At the monitoring site MOISST, the main land cover class grassland/herbaceous shows slightly higher permittivity results with a median value of = 14.9 − 1.7. The latter is in contrast to the second dominant land cover class cultivated crops (cf. Section 2) with a median value of = 13.8 − 1.5. However, the results for cultivated crops span the largest value ranges of all classes and monitoring sites with the real part of permittivity varying from 6 to 31 and the imaginary part of permittivity varying from 1 to 3.2. The three main land cover classes, deciduous, mixed, and evergreen forest (cf. Section 2), at Harvard Forest display the overall smallest value ranges and the highest permittivity results with comparable median values of ~ = 16 − 1.8. Hence, the forest site Harvard Forest, characterized by a cold humid continental climate (cf. Section 2), shows the highest permittivity with the smallest deviations, while the driest monitoring site, mainly covered by shrub/scrub due to the arid climate (Walnut Gulch), shows the lowest permittivity results. The monitoring site characterized by a temperate climate and grassland/herbaceous or cultivated crops (MOISST) displays permittivity results in between.
Converted to soil moisture, Walnut Gulch revealed values in a range from 12.  Based on the complex permittivity results, the penetration depth δ P for each monitoring site is estimated according to Equation (20) (cf. Section 3.3). In Figure 8, areal results for δ P are displayed for every AirMOSS monitoring site. At all sites, δ P varies in comparable ranges from 5.7 cm, 5.4 cm, or 7.7 cm to 35 cm, with the majority of values from 20 cm to 35 cm. Walnut Gulch and MOISST show higher variations in results according to varying land cover classes, for instance, lower depths in regions where grassland/herbaceous or forests are dominant, and greater depths where shrublands or cultivated crops are dominant ( Figure 2, Section 2). In contrast, at Harvard Forest only small deviations in results are observed with the majority of δ P values being at around 25 cm. The two slightly apparent vertical stripes in the results with overall higher δ P values originate from varying incidence angles of the sensing system across the monitoring site ( Figure 8C).
In detail, variations in δ P are largest, where land cover is rather heterogeneous. For instance, the monitoring site Harvard Forest, which is almost fully covered by forests (Figure 2, Section 2), displays the smallest deviations in results and a homogeneous map of δ P across the entire site ( Figure 8C). In contrast, in the eastern part of the monitoring site Walnut Gulch around the city Elfrida, highest variations in δ P are estimated because in that region the land cover is a rather heterogeneous mixture of cultivated crops, grassland/herbaceous, and shrub/scrub. Further, in the center of the site Walnut Gulch at the Whetstone Mountains, which are covered by forests, the estimated δ P are lower than in the surrounding areas, where shrub/scrub is dominant (Figure 2, Section 2). In the western part of Walnut Gulch, at the dip east of the Green Valley (Figure 2, Section 2), δ P are lowest ( Figure 8A). small deviations in results are observed with the majority of values being at around 25 cm. The two slightly apparent vertical stripes in the results with overall higher values originate from varying incidence angles of the sensing system across the monitoring site ( Figure 8C). In detail, variations in are largest, where land cover is rather heterogeneous. For instance, the monitoring site Harvard Forest, which is almost fully covered by forests ( Figure 2, Section 2), displays the smallest deviations in results and a homogeneous map of across the entire site ( Figure 8C). In contrast, in the eastern part of the monitoring site Walnut Gulch around the city Elfrida, highest variations in are estimated because in that region the land cover is a rather heterogeneous mixture of cultivated crops, grassland/herbaceous, and shrub/scrub. Further, in the center of the site Walnut Gulch at the Whetstone Mountains, which are covered by forests, the estimated are lower than in the surrounding areas, where shrub/scrub is dominant (Figure 2, Section 2). In the western part of Walnut Gulch, at the dip east of the Green Valley ( Figure 2, Section 2), are lowest ( Figure 8A).
For detailed analyses, results are displayed in Figure 9 for the main land cover classes of the three monitoring sites. As can be seen in the upper plot of Figure 9, the land cover classes for forests (deciduous, evergreen, and mixed) at Harvard Forest have the lowest values of all landcover classes at ~24.6 cm, with the smallest variations in results, varying from ~22.8 cm to ~26.4 cm only. In contrast, at Walnut Gulch the highest median value is estimated at 27 cm, varying in the range from 16.8 to 35 cm. The results at the monitoring site MOISST display a median of 25.1 cm for grasslands/herbaceous and 25.7 cm for cultivated crops with the largest deviations in the latter land cover class (varying from 13.6 to 35 cm). Accordingly, PDFs in the lower plot of Figure 9 display the highest densities with the narrowest distributions for forest land cover classes, peaking from 24 to 25 cm, while other land cover classes reveal more distributed results, covering a broader range of estimates at lower densities ( Figure 9). For detailed analyses, δ P results are displayed in Figure 9 for the main land cover classes of the three monitoring sites. As can be seen in the upper plot of Figure 9, the land cover classes for forests (deciduous, evergreen, and mixed) at Harvard Forest have the lowest δ P values of all landcover classes at~24.6 cm, with the smallest variations in results, varying from~22.8 cm to~26.4 cm only. In contrast, at Walnut Gulch the highest median δ P value is estimated at 27 cm, varying in the range from 16.8 to 35 cm. The results at the monitoring site MOISST display a median δ P of 25.1 cm for grasslands/herbaceous and 25.7 cm for cultivated crops with the largest deviations in the latter land cover class (varying from 13.6 to 35 cm). Accordingly, PDFs in the lower plot of Figure 9 display the highest densities with the narrowest distributions for forest land cover classes, peaking from 24 to 25 cm, while other land cover classes reveal more distributed results, covering a broader range of δ P estimates at lower densities ( Figure 9).

Comparison of Permittivity Estimates with In Situ Measurements
For a comparison of the estimated soil permittivity with in situ soil moisture measurements at single locations (cf. Section 2), the results for ε s are converted to soil moisture θ [vol.%] according to the dielectric mixing model of [71].
In Figure 10, retrieved θ, based on the proposed P-band SAR moisture estimation method and P-band SAR measurements, are compared with in situ measured θ at the driest AirMOSS site, Walnut Gulch in Arizona. While retrieved values are in the range from 12.6 vol.% to 34 vol.% with a highest PDF from 12.5 vol.% to 19 vol.%, in situ measurements vary from 2.6 vol.% to 15.4 vol.% with the PDF peaking at 6.7 vol.%. Conducted statistical tests between retrieved and in situ measured θ values with a correlation coefficient (R 2 ) [-] and medium root-mean square error (RMSE) [vol.%] of R 2 = 0.2 and RMSE = 14.1 vol.% confirm a clear overestimation of retrieved θ compared to in situ measurements. The main reason for the low correlation could be the discrepancy between measuring depths since available in situ values are measured near the soil surface at~0-5 cm, while for retrieved θ, based on decomposed P-band SAR observations, penetration depths proved to vary from 16.8 cm to 35 cm with median δ P at 27 cm (Figure 9).

Comparison of Permittivity Estimates with In Situ Measurements
For a comparison of the estimated soil permittivity with in situ soil moisture measurements at single locations (cf. Section 2), the results for ′ are converted to soil moisture [vol.%] according to the dielectric mixing model of [71].
In Figure 10, retrieved , based on the proposed P-band SAR moisture estimation method and P-band SAR measurements, are compared with in situ measured at the driest AirMOSS site, Walnut Gulch in Arizona. While retrieved values are in the range from 12.6 vol.% to 34 vol.% with a highest PDF from 12.5 vol.% to 19 vol.%, in situ measurements vary from 2.6 vol.% to 15.4 vol.% with the PDF peaking at 6.7 vol.%. Conducted statistical tests between retrieved and in situ measured values with a correlation coefficient ( 2 ) [-] and medium root-mean square error (RMSE) [vol.%] of 2 = 0.2 and = 14.1 vol. % confirm a clear overestimation of retrieved compared to in situ measurements. The main reason for the low correlation could be the discrepancy between measuring depths since available in situ values are measured near the soil surface at ~0-5 cm, while for retrieved , based on decomposed P-band SAR observations, penetration depths proved to vary from 16.8 cm to 35 cm with median at 27 cm (Figure 9).  Figure 10. Comparison of retrieved soil moisture values, estimated with the proposed P-band SAR method and converted from permittivity according to the dielectric mixing model from [71], with in situ measured soil moisture values from various networks (cf. Section 2) at the AirMOSS monitoring site Walnut Gulch, AZ. � represents the respective median value, n is the total amount of values per boxplot.
In Figure 11, the comparison is shown for the AirMOSS test site Harvard Forest. While retrieved vary in the range from 25.6 vol.% to 32.6 vol.%, with the PDF peaking at 29.3 vol.%, in situ measurements vary from 25.6 vol.% and 38.8 vol.%, with a median value of 32.4 vol.%. Figure 10. Comparison of retrieved soil moisture values, estimated with the proposed P-band SAR method and converted from permittivity according to the dielectric mixing model from [71], with in situ measured soil moisture values from various networks (cf. Section 2) at the AirMOSS monitoring site Walnut Gulch, AZ. x represents the respective median value, n is the total amount of values per boxplot.
In Figure 11, the comparison is shown for the AirMOSS test site Harvard Forest. While retrieved θ vary in the range from 25.6 vol.% to 32.6 vol.%, with the PDF peaking at 29.3 vol.%, in situ measurements vary from 25.6 vol.% and 38.8 vol.%, with a median value of 32.4 vol.%. monitoring site Walnut Gulch, AZ. � represents the respective median value, n is the total amount of values per boxplot.
In Figure 11, the comparison is shown for the AirMOSS test site Harvard Forest. While retrieved vary in the range from 25.6 vol.% to 32.6 vol.%, with the PDF peaking at 29.3 vol.%, in situ measurements vary from 25.6 vol.% and 38.8 vol.%, with a median value of 32.4 vol.%. Figure 11. Comparison of retrieved soil moisture values, estimated with the proposed P-band SAR method and converted from permittivity according to the dielectric mixing model from [71], with in situ measured soil moisture values from various networks (cf. Section 2) at the AirMOSS monitoring site Harvard Forest, MA. � represents the respective median value, n is the total amount of values per boxplot.
Hence, with estimated 2 = 0.5 and = 4.6 vol. % , results show that the humid forest site displays a higher correlation and lower error between retrievals and in situ measurements compared to the arid site Walnut Gulch (Figure 10), despite the same discrepancy between measuring depths. In situ measurements are again only available from the soil surface at ~0-5 cm, while retrieval results showed average penetration depths around 24.6 cm, overall varying from 22.8 to 26.4 cm (Figure 9). Hence, with estimated R 2 = 0.5 and RMSE = 4.6 vol.%, results show that the humid forest site displays a higher correlation and lower error between retrievals and in situ measurements compared to the arid site Walnut Gulch (Figure 10), despite the same discrepancy between measuring depths. In situ measurements are again only available from the soil surface at~0-5 cm, while retrieval results showed average penetration depths around 24.6 cm, overall varying from 22.8 to 26.4 cm (Figure 9).
Lastly, the comparison of θ values is conducted for MOISST ( Figure 12). Only for this site, in situ measurements for comparisons are available from the soil surface (~0-5 cm) and also for greater soil depths at~5-10 cm and~10-20 cm. The median values of in situ measured θ are increasing with greater measuring depths, from 16.4 vol.% near the soil surface to 24 vol.% at a soil depth of~10-20 cm. Concurrently, the in situ value ranges are decreasing with greater soil depths, displaying the smallest deviations in measurements at a depth of~10-20 cm, varying from 19.3 vol.% to 32.3 vol.%. In comparison, retrieved θ vary between 12.6 vol.% and 47.5 vol.% with a median value of 26.4 vol.%. Conducted statistics show R 2 = 0.09 and RMSE = 13.5 vol.% for 0-5 cm, R 2 = 0.13 and RMSE = 12 vol.% for 5-10 cm, and R 2 = 0.09 and RMSE = 12.2 vol.% for 10-20 cm. Hence, retrieval results correlate slightly worse with in situ measurements from the top-soil at~0-5 cm compared to other measuring depths, and display the smallest difference in median values of 2.4 vol.% with in situ measurements at the greatest measuring depth (~10-20 cm). However, estimated θ at MOISST showed a median penetration depth of 25.7 cm, with most values varying from 13.6 to 35 cm for landcover class cultivated crops ( Figure 9). Hence, there is still a discrepancy between measuring depths, which explains a remaining difference between in situ (10-20 cm) and SAR-based soil moisture estimates. measurements from the top-soil at ~0-5 cm compared to other measuring depths, and display the smallest difference in median values of ~2.4 vol.% with in situ measurements at the greatest measuring depth (~10-20 cm). However, estimated at MOISST showed a median penetration depth of 25.7 cm, with most values varying from 13.6 to 35 cm for landcover class cultivated crops ( Figure 9). Hence, there is still a discrepancy between measuring depths, which explains a remaining difference between in situ (10-20 cm) and SAR-based soil moisture estimates.

Figure 12.
Comparison of retrieved soil moisture values, estimated with the proposed P-band SAR method and converted from permittivity according to the dielectric mixing model from [71], with in situ measured soil moisture values from various networks (cf. Section 2) at the AirMOSS monitoring site MOISST, OK. � represents the respective median value, n is the total amount of values per boxplot.

Discussion
In comparison to model-or eigen-based decomposition methods, the hybrid decomposition technique proposed in this study combines both methods and further employs a soil scattering model suitable for P-band. Hence, the proposed P-band SAR moisture estimation method continues the efforts of [24], an iterative approach for soil moisture estimation applied to L-band observations, and of [3], the first hybrid decomposition method applied to P-band SAR data by assuming volume scattering as a cloud of randomly oriented dipoles (cf. Sections 1 and 3.2). Furthermore, the method proposed here is non-iterative with reduced algorithm complexity and, hence, computationally less expensive. It needs no prior assumptions on initial vegetation conditions since it allows the application of a wide range of vegetation cover characteristics in comparison to [3]. Lastly, it is the first of its kind, to the best of our knowledge, for complex soil permittivity estimation from (P-band) SAR data. The well- Figure 12. Comparison of retrieved soil moisture values, estimated with the proposed P-band SAR method and converted from permittivity according to the dielectric mixing model from [71], with in situ measured soil moisture values from various networks (cf. Section 2) at the AirMOSS monitoring site MOISST, OK. x represents the respective median value, n is the total amount of values per boxplot.

Discussion
In comparison to model-or eigen-based decomposition methods, the hybrid decomposition technique proposed in this study combines both methods and further employs a soil scattering model suitable for P-band. Hence, the proposed P-band SAR moisture estimation method continues the efforts of [24], an iterative approach for soil moisture estimation applied to L-band observations, and of [3], the first hybrid decomposition method applied to P-band SAR data by assuming volume scattering as a cloud of randomly oriented dipoles (cf. Sections 1 and 3.2). Furthermore, the method proposed here is non-iterative with reduced algorithm complexity and, hence, computationally less expensive. It needs no prior assumptions on initial vegetation conditions since it allows the application of a wide range of vegetation cover characteristics in comparison to [3]. Lastly, it is the first of its kind, to the best of our knowledge, for complex soil permittivity estimation from (P-band) SAR data. The well-known problem of decomposition methods to overestimate the volume scattering components (cf. Section 1), e.g., [9,27,33], is solved by excluding any vegetation characterization from further analyses, where scattering intensities become negative (cf. Section 3.2). However, this is only feasible since multiple vegetation characterizations can be considered.
By applying the proposed P-band SAR moisture estimation method to SAR observations at three AirMOSS monitoring sites with varying land cover classes, the scattering contributions of surface, dihedral, or volume mechanisms could be estimated. Results show dominant dihedral scattering over forests at the site Harvard Forest in Massachusetts, and dominant volume scattering over shrub/scrub at the site Walnut Gulch in Arizona as well as over cultivated crops at the site MOISST in Oklahoma. Further, it is interesting to notice that cultivated crops at Walnut Gulch in Arizona show dominant dihedral scattering while cultivated crops at MOISST in Oklahoma display dominant volume scattering ( Figure 5). This may be due to the different types and phenological stages of cultivated vegetation. Hence, results in this study suggest dominant dihedral scattering for increasing vegetation volume (height, density) at P-band. These results concur with previous studies. For example, [77] found strong dihedral scattering for P-band measurements across a boreal forest in Canada. They concluded that at P-band the ground and trunks contribute with more relevance to the SAR signal than the branches and leaves, since with "smaller scatterers compared to the wavelength" [77], these layers become more transparent to the incident wave compared to higher frequencies such as L-or C-band. This is supported by findings of [19] across the AirMOSS monitoring site MOISST, where results also demonstrated that vegetation is more transparent at P-band than at L-band, and that woody vegetation is needed for dihedral scattering to appear in P-band. Similarly, Ref. [78] estimated dominant dihedral scattering contributions due to mainly trunk-ground interactions in P-band SAR signals across eucalyptus species in woodlands on the east coast of Australia. Analogue results at L-or C-band, however, revealed predominant volume scattering from branches and leaves (C-band) or a mix of dihedral and volume scattering, depending on the polarization (L-band). They concluded that with multiple frequencies and polarizations, individual subsets of forest structures could be analyzed [78]. In contrast, findings from [3], which are based on the same AirMOSS dataset employed in this study, showed increasing volume scattering with increasing vegetation. In the case of the monitoring site Harvard Forest, they found mainly volume scattering and less dihedral scattering, contrary to the results presented. One reason may be that the applied method in [3] only considers one type of vegetation (randomly oriented dipoles) across the entire study area, and did not employ any soil scattering model for the determination of the volume intensity component prior to the decomposition of SAR signals. However, the results in this study are supported by the fact that overall lower penetration depths into the soil could be estimated across forests compared to shrub/scrub (Figure 9, Section 4.2.). Over forests, such as at Harvard Forest, P-band microwaves are able to penetrate through vegetation but less into soils, and interact more with trunks. In smaller and less densely vegetated areas such as shrub/scrub at Walnut Gulch, P-band microwaves are able to penetrate deeper into soils but interact more with present branches and leaves due to the absence of trunks (2nd scattering center for dihedral scattering), and thus display a higher volume scattering contribution in contrast to forests. Hence, the characterization of vegetation types within the decomposition method should be carefully addressed. Choosing only one type of representative will likely lead to a biased estimation of volume scattering in P-band SAR signals.
Based on the best fit between simulations and decomposition results, the complex soil permittivity was estimated. Analyses concur with the investigated land cover classes and the climatic circumstances at each monitoring site. Lowest soil permittivity was estimated across shrub/scrub in the hot deserts of Arizona, while the highest permittivity with the smallest deviations in results was determined over forests in Massachusetts, characterized by cold humid continental climates. Accordingly, results at the monitoring site in Oklahoma with a temperate climate displayed soil permittivity in between the former ranges. Lastly, the smallest deviations in results are found in the forest land cover classes, while the largest deviations in results are found for pixels classified as cultivated crops. This fits varying phenological stages of different vegetation types and irrigation situations on agricultural managed lands (Figure 7, Section 4.2).
Converted permittivity to soil moisture was then compared at single locations with in situ measurements. Unfortunately, in situ soil moisture values mainly of the top-soil layer (0-5 cm) are available for comparison, while P-band SAR penetration depth estimates indicate deeper penetration into the soil (mainly around 20-30 cm). Thus, the highest correlation with the lowest deviations was found at Harvard Forest, and the lowest correlations at MOISST, but with slightly increasing correlations for measurements at greater soil depths. Overall, correlations are low and RMSEs are high. The main reason for this is the discrepancy between sensing depths. While in situ soil moisture values are mainly measured near the soil surface at~0-5 cm, the estimated soil moisture is based on Pband SAR observations, which proved to have penetration depths mainly from 20 cm to 35 cm. Thus, retrievals correlate less with near-surface soil moisture measurements but will more likely correlate with measurements at deeper soil depths, as indicated at MOISST (Figure 12, Section 4.3). At Walnut Gulch, for instance, retrievals are clearly overestimating in situ measurements, showing the highest error (cf. Section 4.3). Several studies revealed increasing soil moisture with greater soil depth until a certain point, e.g., [51,[79][80][81]. For example, [80] showed increasing soil moisture values with increasing soil depth from 5 cm to 30 cm at Kendall, Arizona (Figure 2, Section 2). Further, in hot desert regions such as Walnut Gulch, the soil surface layer as a land-atmosphere boundary is highly influenced by climate and weather conditions and, hence, the surface tends to dry out faster than the soil layers below [79]. Another source for error is the spatial discrepancy between field measurements and SAR-based observations. While the in situ measurements are point values from single locations (cf. Figure 2, Section 2), the estimated results based on the AirMOSS dataset have a spatial resolution of~100 m.
Based on the estimated complex soil permittivity, permittivity-based penetration depths were calculated, leading to overall P-band depths from 5 to 35 cm. Highest variations occur in regions, where the land cover is rather heterogeneous or where influences of different phenological stages of the vegetation or irrigation situations lead to varying permittivity estimates. In detail, at the humid forest site the smallest penetration depths at a median value of 24.5 cm with the smallest deviations in results were estimated. Correspondingly, the driest monitoring site shows the greatest variations in penetration depths over cultivated crops, varying mainly from 16.8 to 35 cm around a median value of 27 cm. This is in accordance with previous studies. The authors of [40] showed that P-band can provide soil moisture retrievals for a depth greater than 10 cm when using a frequency lower than 0.5 GHz, while [46] estimated similar P-band penetration depths, varying from 15 to 36 cm, showing significant annual seasonality depending on the soil moisture content.
Lastly, following the assumption of [52], P-band penetration depths between 6.97 and 34.86 cm are realistic (cf. Section 1), which is almost exactly the value range estimated in this study (5 to 35 cm). Hence, depending on soil permittivity and vegetation cover characteristics, P-band is able to penetrate several decimeters into the soil, enabling the estimation of RZSM.

Conclusions and Outlook
In this study, a method for complex soil permittivity and permittivity-based penetration depth estimation based on P-band SAR data is proposed and implemented. For that, a hybrid decomposition technique is set up for separating fully polarimetric P-band SAR observations into individual scattering mechanisms (soil, dihedral, volume).
The proposed method extracts the soil scattering component by applying a hybrid (eigen-and model-based) decomposition method and incorporating a soil scattering model suitable for P-band frequency, the multi-layer SPM. To the best of our knowledge, this method is the first of its kind for complex soil permittivity estimation from P-band SAR data. The approach is only limited to fully polarimetric SAR observations as co-and cross-polarized coherency matrix elements are needed within the decomposition. Further, although the proposed approach provides the ability to consider many different types of vegetation structures, the characterization of varying vegetation types within the hybrid decomposition method should be carefully addressed.
Comparisons of retrieved soil moisture with in situ measurements show lower to medium correlations, proving the overall larger penetration of P-band microwaves into soils, since results match less with top-soil layer measurements but more with deeper (~20 cm) measurements. The lack of possibility to analyze retrieval results, presented in this study, with comparable soil moisture measurements, originating at least from similar soil depths, proves the need for measurements at greater soil depths. Overall, P-band penetration depths from 5 to 35 cm were estimated, depending on regional climate and land cover characteristics and heterogeneity.
In summary, the proposed method was further adapted and revised in order to overcome many drawbacks of previous decomposition methods (e.g., overestimation of volume scattering component, consideration of one single vegetation type, missing multi-layer structures for soil scattering at P-band), and to be applicable for complex soil permittivity estimation. Presented decomposition and complex permittivity results agree well with climate and soil conditions at the three different monitoring sites, proving the feasibility of the method. Further, estimated P-band penetration depths are in accordance with previous studies and confirm the overall higher penetration ability of P-band compared to L-band microwaves. In addition, strong dihedral scattering contributions in the SAR signal at P-band over forests, as shown in Figure 8 (cf. Section 4.2), are encouraging the development of an approach for soil moisture estimation under forests based on the dihedral scattering component instead of the soil scattering component, as proposed by [82].
Finally, since the proposed method only needs fully polarimetric P-band SAR measurements, it can be used to estimate complex soil permittivity and potential P-band penetration depths from space based on the SAR data of the upcoming BIOMASS mission (cf. Section 1) [13].