Vegetation Effects on Soil Moisture Retrieval from Water Cloud Model Using PALSAR-2 for Oil Palm Trees

In oil palm crop, soil fertility is less important than the physical soil characteristics. It is important to have a balance and sufficient soil moisture to sustain high yields in oil palm plantations. However, conventional methods of soil moisture determination are laborious and time-consuming with limited coverage and accuracy. In this research, we evaluated synthetic aperture radar (SAR) and in-situ observations at an oil palm plantation to determine SAR signal sensitivity to oil palm crop by means of water cloud model (WCM) inversion for retrieving soil moisture from L-band HH and HV polarized data. The effects of vegetation on backscattering coefficients were evaluated by comparing Leaf Area Index (LAI), Leaf Water Area Index (LWAI) and Normalized Plant Water Content (NPWC). The results showed that HV polarization effectively simulated backscatter coefficient as compared to HH polarization where the best fit was obtained by taking the LAI as a vegetation descriptor. The HV polarization with the LAI indicator was able to retrieve soil moisture content with an accuracy of


Introduction
Soil moisture content is a critical input variable in a wide variety of scientific studies in the field of agriculture. Soil moisture is an environmental element that connects the Earth's surface and the atmosphere. When soil moisture levels are balanced, agricultural yields improve, yield losses due to drought are reduced, and groundwater levels are recharged, ensuring the continuity of rivers and stream flows [1]. Oil palm trees require a reasonably steady high temperature, and continuous precipitation all year. Furthermore, prolonged dry periods of more than 2-3 months do not directly harm vegetative growth but have a significant impact on the yield and quality of fruit bunches [2]. The yield of oil palm is highly dependent on the availability of water during the sex differentiation of its inflorescences, which occurs approximately 28 months before bunch harvest. Soil moisture is therefore critical for optimal production in the oil palm crop [3]. Plants may be stressed if there is a water shortage. When water is not a constraint, potential evapotranspiration is the quantity of water that might be evaporated or transpired at a given temperature and humidity [4]. Water is a key component of plant tissue and is used to transport metabolites and minerals inside the plant. Water is also required for cell expansion, whereby it increases from VH polarization data [30]. Similarly, the biomass estimation of oil palm plantations using a regression analysis of HV polarized PALSAR was carried out [31].
A few SAR backscatter modelling techniques have been carried out to explore the influence of surface-related variables on the backscatter coefficient, such as topography and vegetation. These models are generally divided into theoretical, semi-empirical and empirical methods. Theoretical techniques such as the Integral Equation Model (IEM) and advanced IEM are complex, where the main challenge in assembling model parameters that accurately characterize the canopy is tough [32] as it requires a large number of parameters [33,34]. In addition to that, the empirical methods, for example the Dubois model [35], are simpler than the IEM models. However, the Dubois model has been reported to be dependent on experimental site and data conditions [36], where the saturation points were reported at lower NDVI values [37].
When the semi-empirical method was used, effects of vegetation cover using the WCM [38] were initially developed and later on were modified [39][40][41]. WCM is preferred for its simplicity [42] where it defines the overall backscatter coefficient obtained by the sensor over vegetated surfaces as the incoherent sum of the effects of vegetation and soil [43]. In this model, the canopy is often represented as a set of exact variables such as plant density per area, leaf size and orientation, which complicates and makes the model hard to interpret [44]. Many factors influence the backscatter coefficients from vegetation canopies, e.g., size of disperses in a canopy; the shape of scatters in a canopy; the orientation in a canopy; and the geometry of the canopy cover on ground [32]. To begin with, the WCM was modelled using plant height and water content of the vegetation layer to allow for the retrieval of soil moisture [38], and Leaf Area Index (LAI) was subsequently evaluated [45]. Other vegetation descriptors in the WCM were used from various combinations such as leaf water area index (LWAI), Normalized Plant Water Content (NPWC), vegetation water mass, and biomass [43,[46][47][48].
The goal of this research is to optimize soil moisture retrieval using the WCM, thereby reducing the effect of vegetation on the crop-covered soil moisture backscatter coefficient. The WCM attempts to express vegetation cover scattering and attenuation terms in the concept through plain vegetation descriptors. Furthermore, a comparison analysis is performed on the use of the three primary vegetation descriptors, LAI, NPWC, and LWAI, singly or in combination, in the retrieval of soil moisture from PALSAR-2 data from Malaysian oil palm trees.

Materials and Methods
To recover soil moisture, a basic WCM was used in this investigation. Assuming that the effect of soil surface roughness on observed backscatter is consistent over a short time period for a given site, the temporal change in SAR backscattering only reflects changes in vegetation and soil moisture [49]. As a result, for this investigation, multi-temporal SAR data were used in the WCM. The vegetation descriptors and actual soil moisture were determined using field collected data. The Root Mean Square Error (RMSE) and mean absolute error coefficient were utilized to analyze the soil moisture retrieval accuracy. The next sections describe WCM, vegetation descriptors used in the model and evaluation of data processing methods.

Water Cloud Model (WCM)
WCM was developed assuming that the canopy "cloud" contains identical water droplets randomly distributed within the canopy [38]. In a water cloud model, the expression developed incorporated the soil moisture and the vegetation parameters in the equation. WCM for a given polarization (pp) is given as where σ • tot,pp is the total backscatter coefficient, σ • veg,pp is the backscatter contribution of the vegetation cover, σ • soil+veg, pp is the multiple scattering involving vegetation elements and the soil surface, σ • soil,pp is the backscatter contribution of the soil surface and τ 2 pp is the two-way vegetation attenuation. The second component in Equation (1) reflects the interaction of incident radiation with the underlying soil. Because the interaction is not a dominant factor in co-polarized radiation, it may be ignored [43,50]. Subsequent to the model development, there were studies that modified the mentioned model [45]. WCM, Equation (2) [43], with four empirical coefficients, A pp , B pp , C pp , and D pp where A pp and B pp are vegetative characteristics and C pp and D pp are soil parameters, was presented for the given polarization where pp is either HH or HV polarization. The parameter A pp corresponds to the albedo of the vegetation, and B pp is an attenuation factor seen in (3) and (4). The parameter D pp indicates the sensitivity of the radar signal to soil moisture, and C pp can be considered as a calibration constant in (5). Hence, the equation written for the given polarization where pp is HH and HV is modified to where V 1 and V 2 describe the effect of canopy water content and its geometry on backscatter coefficients. Vegetation descriptors have different effects on the WCM model. Therefore, several experimental studies on different combination of vegetation variables to quantify V 1 and V 2 in WCM were conducted by using plant height, LAI, LWAI and NPWC [51]. In Equations (3) and (4), θ is the incident angle of the image used, A pp and B pp are the vegetation parameter. M v is the volumetric soil moisture and C pp and D pp are the soil parameter in Equation (5).
In this study, the LAI of palm fronds is defined by the amount of leaflet surface area per unit ground area [3]. LWAI is a product of LAI of palm fronds multiplied by the amount of water (W), expressed as the ratio of the difference between wet and dry mass to wet mass as in Equations (6) and (7).
where m w and m d are field records for freshly plucked and oven-dried mass of vegetation samples, respectively. Normalised plant water content, NPWC, also plays a dominant part in attenuating backscatter [38]. NPWC is like Equation (8) but divided with dry mass instead of wet mass.
In this study, 5 combinations were used to understand the vegetation effects of WCM by using Case 1 where V 1 = 1 and V 2 = LAI Case 2 where V 1 = LAI and V 2 = 1 Case 3 where V 1 = LAI and V 2 = LAI Case 4 where V 1 = LWAI and V 2 = LWAI Case 5 where V 1 = NPWC and V 2 = NPWC

Estimating Parameters of A, B, C and D in the WCM
The estimation of parameters C pp and D pp are solved using a linear model fitting, following which, the values of C pp and D pp are substituted into Equation (2), which allows for parameter A pp and B pp to be solved using the Nonlinear Least Squares Method (NLSM) [51][52][53][54]. It was noted that using Levenberg-Marquardt (LM) optimization in NLSM, an estimation of A pp and B pp can be made [55]. According to the optimization, LM is a common approach for addressing nonlinear least square issues that emerge from fitting a parameterized function to a collection of observed data points by minimizing the sum of the squares of the errors between the observed data points and the function output [56]. The LM algorithm is an iterative process which combines the Gauss-Newton method and the gradient descent method [57].

Evaluating WCM
To validate the WCM prior to inversing the model for accuracy and estimation performance metrics which include the coefficient of determination (R 2 ), the RMSE and the mean absolute error (MAE) were examined. The RMSE derived using Equation (9) was tested in most studies [49,55,57,58]. RMSE is a frequently used measure of the difference between values predicted by a model and the values observed from the environment that is being modelled. RMSE evaluation is given in Equation (9); X sim is the simulated σ • tot,pp , X obs is the observed σ • tot,pp and n is the sample count. These individual differences are also called residuals, and the RMSE serves to aggregate them into a single measure of predictive power [59].
Together with R 2 and RMSE, the MAE was explored to further evaluate the model [60]. MAE is the average magnitude of the errors in a set of predictions, without considering their direction [61]. It is the average over the test sample of the absolute differences between prediction and actual observation where all individual differences have equal weight, as seen in (10)

Study Area
The study area of this research was in Chuping district, Perlis state, which comprises a region of flat-terrain oil palm cultivation. The area mentioned is about 28 ha. It covers 4 years of old palm stands which have just started fruiting. The center point coordinates of the study area are 6 • 31 07.2 N 100 • 19 07.7 E. In Chuping district, there are seven sub-districts namely Panggas, Sungai Buloh, Kubang Perun, Guar Nangka, Felda Chuping, Sungai Buloh, and Kilang Gula Chuping. This research took place at the Kilang Gula Chuping subdistrict as shown in Figure 1. The field data collected in this research were soil moisture, leaf moisture content, and leaf area index from frond 17 of the palms. Field data collection was planned to match the acquisition date of PALSAR-2 as shown in Section 3.2.
In the research site, the soil type was identified as Chuping and Dampar soil series, where the Chuping soil type is categorized as sandy clay loam with a hue of 7.5-10 YR (Yellow to Red) following the Munsell color chart standards, and the latter soil type was identified as clay loam with a hue of 7.5 YR [62]. The site was considered flat topography with an elevation of 21.6 m while the slope class was identified as 4-12% [63]. As per the study site, the months of February to March of the calendar year are regarded as the driest season of the year; with a maximum of 28 days of dry spell every year and an average of 1362.38 mm of precipitation per year in the years of 2015-2017 [64]. In the PALSAR-2 image acquisition months, Table 1 shows the monthly meteorological data obtained from Malaysian Meteorological Department. The observation dates in Table 1 represent the annual weather of the study area where the month of January was observed to be drier than April and July; however, the Mean Evaporation and Mean Radiation were similarly recorded. In the research site, the soil type was identified as Chuping and Dampar soil series, where the Chuping soil type is categorized as sandy clay loam with a hue of 7.5-10 YR (Yellow to Red) following the Munsell color chart standards, and the latter soil type was identified as clay loam with a hue of 7.5 YR [62]. The site was considered flat topography with an elevation of 21.6 m while the slope class was identified as 4-12% [63]. As per the study site, the months of February to March of the calendar year are regarded as the driest season of the year; with a maximum of 28 days of dry spell every year and an average of 1362.38 mm of precipitation per year in the years of 2015-2017 [64]. In the PALSAR-2 image acquisition months, Table 1 shows the monthly meteorological data obtained from Malaysian Meteorological Department. The observation dates in Table 1 represent the annual weather of the study area where the month of January was observed to be drier than April and July; however, the Mean Evaporation and Mean Radiation were similarly recorded.

Data and Processing
PALSAR-2 signals were used for extracting the backscatter coefficient from the oil palm trees. High resolution scenes of PALSAR-2 were obtained under our collaboration with the Earth Observation Research Announcement 2 (EO-RA-2) of the Japanese Aerospace Exploration Agency (JAXA). We used a total of 3 PALSAR-2 scenes in this research, details of which are given in Table 2.
The acquired PALSAR-2 data were first converted from their Digital Number (DN) to backscatter coefficients in decibels, dB for both polarizations. Once the σ • HH and σ • HV for each of the field points was available, georeferencing of both images using SRTM (3 Sec) was done. The Lee filter was applied with a 5 × 5 window size to remove speckle and noise from the image. All the rectification process was completed in the SNAP opensource software. Additionally, Unmanned Aerial Vehicle (UAV) imagery was acquired on 17 January 2019 with the spatial resolution of 8 by 8 cm consisting of the red, blue, green, red edge and near infrared bands. The imagery from the UAV platform allowed us to compute the NDVI to identify bare soil areas in the field for the calibration process, with an NDVI value less than 0.2 taken to be bare soil. Ground truthing was done to validate the bare soil areas.

In-Situ Data Collection
The determination of soil moisture content at soil depth 0-5 cm was done using the soil gravimetric method. Fresh weight of soil moisture collected in field was recorded whereas the dry weight was determined after oven drying of the soil samples in the laboratory. A total of 32 points of soil sample were taken at each observation dated 17 January, 19 April and 9 July 2019, totalling to 96 soil samples altogether.
For the LAI for oil palm trees, a conventional method of destructive sampling was carried out using Equation (11). The frond to be chosen must not only be a good indicator of the nutrient status of the palm, but must also be one which is easy to identify and is consistent. By convention, frond 17 is used commonly in the oil palm fields. It is reported that frond 17 provides satisfactory results from its use as an input in the LAI formula [65]. According to the localised LAI formula, LAI was determined by the following equation: where A f is the leaf area per frond (m 2 ), F n is the number of fronds per palm; P DEN is the planting density where the number of palms per hectare is identified. Leaflets from the one-sided leaflet area were multiplied by two to obtain the total leaflet area of the frond in this research [66]. All leaflets were brought to the laboratory to determine the leaflet area using the LI-3100C, LI-COR Inc., which is USA-made equipment. Samples from the leaflets were taken and oven dried at 70 • C for 72 h until constant weight was achieved. An electronic balance was used to weigh the oven dried leaflets upon the completion of drying step. Similarly, the LWAI and NPWC were used from the same samples and methods as seen in Equations (6)-(8).

In-Situ Results
Field data collection for the three observation periods was based on the availability of PALSAR-2 orbit path data acquisition. In situ information was simultaneously collected for the soil moisture. LAI, LWAI and NPWC were calculated to understand the range values and average of each component in field. In Table 3, it can be noted that during the January and April data collection, the soil moisture average was 0.240 m 3 /m 3 whereas the July data acquisition showed an increase in soil moisture recording at 0.273 m 3 /m 3 . Overall, for the palm trees cultivated areas, the soil moisture at the depths of 0-5 cm was averaged at 0.251 m 3 /m 3 . LAI for the palm fronds was averaged at 1.845 m 2 /m 2 where it was seen to increase at every observation date, being 1.748, 1.784 and 2.005 m 2 /m 2 , respectively, during the January, April, and July observations dates. For the water content based on LAI in the palm fronds, an increase on each observation date was observed as well; the January observation date was at 0.205, 0.353 for the April observation date and 0.396 in the July observation date. This gave an average of 0.368 for all three observation dates. It can be noted that the study area faced drought stress in the earlier months of the calendar year where January to March recorded a lower rainfall [64]. Normalized plant water content was similar in all the observations at an average of 0.199%.

Water Cloud Model Parameterization
The WCM model parameterization values were calibrated using bare soils and vegetation points from the NDVI values, derived from the RGB-NIR data obtained from UAV data acquisition, to fit the WCM model with the ground range values. Table 4 shows the vegetation parameters, A pp and B pp and the soil parameters C pp and D pp for the given polarization together with the combination indicator that represents the simulated scenarios of Cases 1-5. Using these vegetation and soil parameters together with LAI, LWAI, and NPWC as the input values, the WCM model was able to replicate the backscatter coefficients.

Backscatter Simulation Based on the Proposed Vegetation Descriptors
Backscatter coefficients at frequency range of C and X bands are dominated by scattering activities in the crown layer of branches and foliage in the canopies, whereas scattering processes involving substantial trunks and branches would be dominated at lower frequencies like P and L bands [67]. In this study, the L band from PALSAR-2 was considered to give a better penetration in oil palm tree structure where backscatter coefficient characterizes the nature of oil palm structure. It was reported that oil palm trees using L band were capable to penetrate to the basal trunk in the oil palm plantation using both HH and HV polarizations [68]. Using the WCM model, the backscatter coefficients were simulated using a combination of vegetation descriptors. These simulated backscatter coefficients were then compared with the image-derived backscatter where the observed backscatter values were extracted from PALSAR-2. Table 5 shows the comparison of the backscatter coefficients observed and simulated for the study site. Both HH and HV polarizations showed a positive indication of a fitted WCM model using the respective vegetation descriptors where R 2 values ranged from 0.823 to 0.998, indicating the regression model fits the observed data well. In most HV polarization, the regression was found to be higher than HH polarization using the WCM model. On the evaluation of the RMSE and MAE, it was found that the errors were minimized in the HV polarization. In the Case 1 where V 1 = 1 and V 2 = LAI, the HV polarization was more efficient to the simulate the backscatter coefficient with RMSE 0.222 dB and R 2 = 0.997. However, the HH polarization, even though it shows good R 2 of 0.962, had higher RMSE with 2.259 dB. Using the LAI indicator as V 1 and V 2 as 1, both HH and HV polarizations showed a good R 2 with 0.956 and 0.965, respectively, but the RMSE was lower by 1.584 dB in the HV polarization. The MAE gave lower values for the HV polarization with 0.212. In Case 4 using LWAI, the model accurately predicted the simulated backscatter with R 2 = 0.995 in the HV with 0.387 dB RMSE and MAE of 0.366. When the NPWC vegetation descriptors was evaluated for HH polarization, its R 2 value of 0.969 was higher than the HV polarization R 2 of 0.951, but the MAE in HV polarization was computed to be lower than HH polarization by 1.531 and RMSE by 1.905 dB, which means the latter polarization is more accurately simulated. It is important to note that in the LAI combinations of the vegetation descriptors, dual use of LAI was found to be the most accurate parameter to simulate backscatter coefficient values. Furthermore, it is critical to produce an accurate simulation of the backscatter coefficient since it can convey the actual field range values obtained during the in-situ collection utilizing the model. This understanding is important to allow future work to be carried out with minimal calibration values and adopted into large areas of plantations where it is remotely challenging to access the collected field data and senses.

Vegetation Effects on Soil Moisture Retrieval Based on Polarization
Vegetation descriptors play an equally important role in the estimation of the backscatter coefficient as well as providing reliable information to showcase the ground parameters in the WCM model. The backscattering coefficient mechanism interaction with soil moisture is complex [69]. Recently, in an attempt to reduce complexity in the modified IEM, the soil parameters were being simplified from three to two soil parameters [70]. Earlier studies have demonstrated the necessity to eliminate the effects of vegetation on soil moisture retrieval [71]. It was discovered that soil moisture retrieval is much more accurate when vegetation cover is considered [72]. To further explore the backscatter coefficients derived in Section 4.3 with different vegetation descriptors, the soil moisture was retrieved from the backscatter to allow comparison of observed soil moisture in the oil palm study area. It was seen that the simulated backscatter coefficient was predicted from the ground soil moisture values accurately given the observation period variability of lesser rainfall in the January observation when compared to the other observations. By making comparisons with the soil moisture retrieved from PALSAR-2 images, the model can be cross-evaluated to see the model fitting using the WCM. The vegetation descriptors evaluated as mentioned by Cases 1-5 can be seen in Table 6. In this comparison between the retrieved and observed soil moisture values, we have the best case with the highest R 2 value of 0.805 and the lowest RMSE of 0.046 m 3 /m 3 . When looking at Case 1 for HV polarization, the soil moisture retrieved was plotted against simulated backscatter coefficient where the R 2 was demonstrated at 0.916, as seen in Figure 2. Similar accuracy was achieved using Radarsat-2 utilizing the LAI descriptor, whereby RMSE of 0.069 m 3 /m 3 was reported [73]. The initial combination introduced in the WCM [38], with V 1 as 1 and V 2 as LAI providing the most promising soil moisture retrieval in the oil palm trees in this study, as seen in Figure 2. LAI parameter was evaluated as one of the best parameters when VV polarization was compared to LWAI and NPWC, with 4.19-4.43% of the RMSE [51]. Signal backscattering is influenced by canopy structure where it is very sensitive to plant water content, a variable highly correlated with LAI during the vegetative phase [74]. It can be observed that the Case 4 using LWAI as vegetation descriptors, using HH polarization gave an R 2 value of 0.727 where the RMSE was found to be 0.085 m 3 /m 3 , the HV polarization recorded a lower error with an RMSE value of 0.066 m 3 /m 3 with R 2 value of 0.727. In contrast, the lowest R 2 value of 0.301 was from the NPWC case where RMSE was observed at 0.075 m 3 /m 3 . Figure 3 shows the observed soil moisture based on polarization of PALSAR images. It can be noted that the HV polarization gives better lower RMSE values as compared to HH polarization, considering all the vegetation descriptors [75]. It was reported that HH polarization could provide a low RMSE value of 0.049 m 3 /m 3 [71] using remote sensing-based vegetation descriptors, in this case Normalized Difference Infrared Index (NDII). Similarly, when canopy water content based on LAI was applied as a vegetation descriptor, an RMSE of 0.039 was recorded [53]. Overall, when RMSE and MAE are carefully evaluated, the cross-polarized HV backscatter coefficient is revealed to be more vulnerable than the co-polarized backscatter HH in terms of polarization response in all Cases 1-5 observed.  It can be observed that the Case 4 using LWAI as vegetation descriptors, using HH polarization gave an R 2 value of 0.727 where the RMSE was found to be 0.085 m 3 /m 3 , the HV polarization recorded a lower error with an RMSE value of 0.066 m 3 /m 3 with R 2 value of 0.727. In contrast, the lowest R 2 value of 0.301 was from the NPWC case where RMSE was observed at 0.075 m 3 /m 3 . Figure 3 shows the observed soil moisture based on polarization of PALSAR images. It can be noted that the HV polarization gives better lower RMSE values as compared to HH polarization, considering all the vegetation descriptors [75]. It was reported that HH polarization could provide a low RMSE value of 0.049 m 3 /m 3 [71] using remote sensing-based vegetation descriptors, in this case Normalized Difference Infrared Index (NDII). Similarly, when canopy water content based on LAI was applied as a vegetation descriptor, an RMSE of 0.039 was recorded [53]. Overall, when RMSE and MAE are carefully evaluated, the cross-polarized HV backscatter coefficient is revealed to be more vulnerable than the co-polarized backscatter HH in terms of polarization response in all Cases 1-5 observed.  Our research results need to be interpreted with respect to the age of the oil palm stands which were 4 years old and our terrain type, whereby the terrain in our study was considered flat and the backscatter intensity represented in the acquired images was a composition of interactions with the crown, the trunk, and the ground surface of the current palm stands. However, it is known that oil palm crop growth or biomass increases over years, hence it is recommended that larger data sets need to be collected from the field to represent greater diversity of palm ages, possibly using fully polarimetric SAR data. Specifically, an understanding of the effects of backscattering on undulating grounds and peatlands in oil palm growing areas will be helpful in future. This could help in understanding the influence of SAR on other variables of importance to the oil palm industry, such as estimation of vegetation water content of oil palm canopy and yield. Benefits of this research can be useful for the plantation industry especially with Our research results need to be interpreted with respect to the age of the oil palm stands which were 4 years old and our terrain type, whereby the terrain in our study was considered flat and the backscatter intensity represented in the acquired images was a composition of interactions with the crown, the trunk, and the ground surface of the current palm stands. However, it is known that oil palm crop growth or biomass increases over years, hence it is recommended that larger data sets need to be collected from the field to represent greater diversity of palm ages, possibly using fully polarimetric SAR data. Specifically, an understanding of the effects of backscattering on undulating grounds and peatlands in oil palm growing areas will be helpful in future. This could help in understanding the influence of SAR on other variables of importance to the oil palm industry, such as estimation of vegetation water content of oil palm canopy and yield. Benefits of this research can be useful for the plantation industry especially with the increase in available high-resolution multi-polarization satellite SAR sensors. This allows for the possibilities in exploiting the use of oil palm crop-related parameters supported with satellite data input to a decision-making platform for oil palm plantations.

Conclusions
In this study, L band from PALSAR-2 observations and in-situ data collected were used to evaluate a WCM model for soil moisture retrieval in oil palm cultivated regions. The WCM was calibrated using NDVI values, then fitted to a WCM model by different instances of vegetation descriptors to evaluate the best fit for the model. Using the WCM model and LAI specifically as V 1 and V 2 , it was found to provide the best simulated backscatter coefficient in the HV polarization, secondly followed by V 1 which is set as constant with a value of 1 and V 2 as LAI. For HV polarization, LAI specifically as V 1 and V 2 showed an R 2 of 0.998 in simulating the backscatter coefficient accurately with RMSE of 0.158 dB, whereas HH polarization showed a lower R 2 of 0.823 with a higher RMSE of 2.384 dB. However, when LWAI and NPWC were considered, the model fitted better with co-polarized backscatter, where HH backscatter coefficient records R 2 of 0.727 and 0.558, respectively. It can be said that the co-polarized backscatters have shown lower RMSE values in the model fit for all cases where it can be predicted to be more useful in retrieving soil moisture in oil palm cultivated areas. In this respect, this research makes a useful and novel contribution on soil moisture for the benefit of the oil palm cultivation, and we found the PALSAR-2 sensor data beneficial for this purpose. The findings of this research will eventually help the oil palm growers to have systems in place to address the abrupt droughts brought about by climate change. Future work of this research can be explored with SAR-based vegetation descriptors, e.g., Radar Vegetation Index and optical derived indices such as NDVI, Normalised Difference Water Index and NDII. Funding: This research was funded by the Japanese Aerospace Agency (JAXA), using the Earth Observation-Research Announcement (EO-RA2) collaboration with PI number of ER2A2N180 to obtained PALSAR-2 data. Data Availability Statement: PALSAR-2 data was requested via proposal where restrictions apply for high resolution data.