Abstract

To promote effective oil palm plantation and environmental sustainability, this research assessed influence of climatic parameters on physicochemical properties of Thai acid sulfate soils (ASSs). ASSs under oil palm planting areas (Topsoil: Ap, Subsoil-1: Ap-60 cm, Subsoil-2: 60–100 cm, and Rootzone: sum of the three depth levels) and historical climate data in tropical savanna and tropical monsoon were investigated. Stepwise approach of multiple regression analysis from component defining variable of principal component analysis revealed that, in tropical savanna, daily solar radiation influenced clay content () in Topsoil, cation exchange capacity () in Subsoil-1, soil pH by water () in Subsoil-2, and silt content () in Rootzone. In tropical monsoon, daily solar radiation influenced exchangeable magnesium () in Subsoil-2, and exchangeable sodium () and silt content () in Rootzone. We concluded that daily solar radiation is the most influential climatic parameter on soil properties due to the transport of heat in soils, whereas particle size distribution is mostly influenced by climatic parameters due to their clay and silt fractions. OPP in ASSs under these climates should include management of water usage by using raised beds with irrigation canals, considering the rooting depth and depth of acidic horizon before applying fertilizer or amendment and liming along with integrated organic material management to raise soil pH.

1. Introduction

Oil palm (Elaeis guineensis Jacq.) is a native of Central Africa, but it can thrive successfully as a commercial crop in tropical lowland areas [1], and it has been considered as the most efficient oilseed crop in the world [2]. In addition, the global market of biodiesel and aviation biofuel has been increasing continually that creates opportunities in Southeast Asia, including Thailand to be major producers, consumers, and exporters [3]. Unfortunately, a plan in Thailand to extend oil palm planting areas in the Southern and Eastern regions for alternative energy was restricted because most area has already been used under productive agriculture [4]. Regarding this, the Land Development Department of Thailand had proposed the main target planting area which was not under productive agriculture [5].

Acid sulfate soils (ASSs) are mostly formed in coastal plains in the tropics and generally less in the temperate zones [6]. The soils are found in South Asia, Southeast Asia, and East Asia of the Tropics [7] with the larger part in Indonesia, Thailand, and Vietnam [8]. Thai ASSs are widespread in some parts of country, especially in the Central Plain (CP) with other small areas in the Southeast Coast (SC) and peninsular regions [9].

Acid sulfate soils have been posing problems for crop production management because of their severe acidity. These soils affect oil palm yields by depth of acidic horizon where it occurs within 60 centimeters of soil surface [6, 10]. Nevertheless, oil palm is moderately tolerant to soil acidity and can be grown on ASS under proper management practices [11] which consider parameters affecting growth performance [12].

The Department of Alternative Energy Development and Efficiency of Thailand estimated the potential ability of palm oil for biodiesel production at 14 million liters per day in 2036 [13]. The planting areas have been expanding into the northern and northeastern regions [14] even though these regions were never planted in the past due to their low climate suitability that requires an intensive treatment and irrigation supply [15]. These lead to high greenhouse gas emissions and water scarcity footprint per ton of fresh fruit bunch (FFB) products [15]. Oil palm plantation areas are somewhat limited by climatic conditions [15, 16], especially oil palm in ASS can be affected by droughts in an El Nino year. The yields are reduced due to severe acidity from low water level [17]. Therefore, improving oil palm productivity should take into consideration the control and management of climatic parameters [18]. However, it is difficult to determine climatic limitations in different regions or countries by comparing the FFB yields because many other parameters may influence this comparison [17]. The main climatic elements affecting growth and yield of oil palms are temperature, solar radiation intensity and duration, rainfall, atmospheric vapor pressure deficit, evaporation rate, and wind speed [19, 20].

The suitability of soil and climate combinations for oil palm development can be calculated by using Surre [21] method [22] which has been continuously used in oil palm industry to calculate water deficit [2325]. Seasonal water deficit is considered as the most important climatic parameter affecting oil palm yields [17]. The yields will be lost about 10% for every 100 mm increase in water deficit [26] and 20% if it exists for 3 years continually [27]. The FFB yield will be reduced between 10% and 20% by water deficit depending on each soil’s quality [23, 28]. Therefore, the promotion for expanding oil palm planting areas must consider environmental sustainability based on suitability of land and climate along with good practices for productivity improvement [15]. In addition, the rooting depth should be at least 100 cm as it is important for nutrition, water uptake, and anchoring in the ground [17]. To support the increasing demand and consider environmental sustainability, there is a need to analyze the suitability of climate and ASS and their influence in CP and SC so as to devise an efficient plan for sustainable oil palm cultivation in tropical setting including Thailand.

2. Materials and Methods

2.1. Soil Sampling

Twenty study sites of ASS under oil palm planting areas in CP (10 sites) and SC (10 sites) of Thailand were examined (Table 1 and Figure 1). These soils have developed under tropical savanna and tropical monsoon climates, respectively. The soils in CP were classified by soil taxonomy [29] as Typic Sulfaquepts and Sulfic Endoaquepts, whereas the soils in SC were classified as Typic Sulfaquepts. These soils in both CP and SC were classified by World Reference Base [30] as Thionic Gleysols. The soils were collected during September and October in 2013 by hand auger to a maximum depth of 100 cm consisting of Topsoil (Ap), Subsoil-1 (Ap–60 cm), and Subsoil-2 (60–100 cm), whereas Rootzone (0–100 cm) was the sum of the three depth levels. The Topsoil depth varies in the range of 18–30 cm in CP and 14–26 cm in SC. In order to minimize oxidation of sulfidic materials, soil samples were kept in sealed plastic bags and cooled to 4°C in an ice container during transportation from field to laboratory.

2.2. Climatic Data

Historical climate data downloaded from the “WorldClim version 2” database (http://www.worldclim.org/version2) [31] for the period 1970–2000 consist of monthly minimum temperature, maximum temperature, mean temperature, precipitation, solar radiation, wind speed, and water vapor pressure. The spatial resolutions with 30 seconds (∼1 km2) containing 12 GeoTiff files, one for each month of the year, were used in this study.

2.3. Soil Analysis

Physicochemical properties of soil samples were determined for field-moist conditions but are reported on an oven-dried (105°C) basis. The methods of physicochemical analysis are summarized in Table 2.

2.4. Data Analysis

Monthly potential evapotranspiration (ETP) was calculated based on Hargreaves and Semani [43] equation following Läderach et al. [44]:where ETP = evapotranspiration (mm·day−1), Ra = extraterrestrial solar radiation expressed in water equivalent (mm·day−1), T − t = different between monthly maximum and minimum mean temperature (°C), and tm = mean air temperature (°C). Ra was converted into water equivalent (mm·day−1) with 1 mm·day−1 = 2.45 MJ·m−2·day−1 [45].

ETP was converted to monthly water balance (B) based on the Surre (1968) equation by following Rhebergen et al. [25]:where B = monthly water balance (mm), Sres = residual soil water from previous month (mm), PPTm = monthly precipitation (mm), and ETP = monthly potential evapotranspiration (mm).

Available water capacity (AWC) for each soil texture class (Topsoil, Subsoil-1, and Subsoil-2) was calculated based on the Epebinu and Nwadialo equation [33]:where S = silt (%) and C = clay (%). The sum of AWC was the maximum amount of available soil water in the top 100 cm. The excess of water was assumed to be lost as drainage or runoff; water deficit was indicated by the negative values of BSres of the following month was set to zero; Sres was set to the maximum (total possible AWC) if it was greater than the sum of AWC, the sum of all negative monthly water balances was annual water deficit [25].

2.5. Climate and Soil Suitability Evaluation for Oil Palm

We used climatic suitability classification of Goh [19], water deficit suitability classification of Rhebergen et al. [25], and ideal climate requirement of Hartley [46] to evaluate climatic suitability (Table 3). Also, we used soil limitation of Shamshuddin et al. [47] and soil nutrient status of Goh [48] to evaluate soil suitability (Table 4).

2.6. Statistical Analysis

The climate and ASS data were analyzed by using principal components analysis (PCA) and multiple regression analysis (MRA) from SPSS (version 24) for Windows. The PCA was performed to reduce the set of climate and soil variables; component defining variable (CDV) selected the variables with the highest principal component loading from the rotated principal components, and stepwise approach of MRA determined the relative effects of climatic parameters on physicochemical properties of ASSs.

3. Results

3.1. Climatic Parameters

The climatic parameters of CP and SC are shown in Table 5. Climate suitability for sustainable oil palm was considered by climatic suitability classification of Goh [19]. Figure 2 shows the climatic conditions of CP and SC classified as currently unsuitable and suitable annual precipitations, highly suitable annual mean temperatures, suitable daily solar radiations, and highly suitable annual wind speeds, respectively. Also, by water deficit suitability classification of Rhebergen et al. [25], Figure 3(a) shows CP is favorable, while SC is optimal. In addition, Figure 3(b) shows CP and SC are in the ranges of ideal climate requirement of Hartley [46] with their annual maxima temperatures and minima temperatures, respectively.

3.2. Physicochemical Properties of ASS

Physicochemical properties of ASS under oil palm planting areas in CP and SC are shown in Tables 6 and 7. The median pH1:1-H2O (soil pH by water) of CP are low at all depth levels, whereas the median pH1:1-H2O of SC are low but slightly increased with depth from Topsoil to Subsoil-2 (Table 7). The median TN (total nitrogen) of CP are high but slightly decreased with depth from Topsoil to Subsoil2, whereas the median TN of SC are high and slightly increased with depth from Topsoil to Subsoil2 (Table 7). The median Avail-P (available phosphorus) of CP are high in Topsoil but low in Subsoil1 and Subsoil2, whereas the median Avail-P of SC are low at all depth levels (Table 7). The median Exch-K (exchangeable potassium) of CP are high at all depth levels, whereas the median Exch-K of SC are low but slightly increased with depth from Topsoil to Subsoil2 (Table 7). The median Exch-Mg (exchangeable magnesium) of CP are high at all depth levels, whereas the median Exch-Mg of SC are high and slightly increased with depth from Topsoil to Subsoil2 (Table 7).

3.3. Result of PCA on Climatic Parameters

The result of PCA based on the correlation matrix analysis with Varimax rotation indicated 2 principal components (PCs) with the extraction of eigenvalues greater than 1 of CP and SC (Table S1). The extracted components had overall cumulative variances about 88.3% and 90.1%, respectively, based on the significant variables which factor loading ≥0.7 and the concept of CDV which selected the highest component loading of each principal component. The result indicated the main climate parameters that influence the physicochemical properties of ASSs (Table S1). The main climatic parameters of CP were daily solar radiation and annual water deficit, while those for SC were annual water vapor pressure, annual minimum temperature, and daily solar radiation (the values of annual water vapor pressure and annual minimum temperature were equal).

3.4. Results of PCA on Physicochemical Properties of ASSs

Likewise, the PCA results of CP and SC had overall cumulative variances of the extracted components of Topsoil about 90.4% and 95.4%, Subsoil-1 about 93.0% and 95.6%, Subsoil-2 about 91.3% and 92.6%, and Rootzone about 90.7% and 90.4%, respectively. The CDV results indicated the main soil properties of each area. The Topsoil of CP were Exch-Ca (exchangeable calcium), clay content, CEC, and Acid-BC (acid buffering capacity), while the Topsoil of SC were Extr-SO42− (extractable sulfate), Exch-K, silt content, clay content, and Exch-Ca (Table S2). The Subsoil-1 of CP were Exch-Ca, CEC, Exch-Na (exchangeable sodium), pH1:1-H2O2 (soil pH by hydrogen peroxide), and EC, while the Subsoil-1 of SC were Acid-BC, Exch-Na, Avail-P, silt content, and clay content (Table S3). The Subsoil-2 of CP were pH1:1-H2O, TN, Exch-Mg, Extr-SO42−, and sand content, while the Subsoil-2 of SC were Exch-Mg, AWC, Acid-BC, Avail-P, and EA (Table S4). The Rootzone of CP were Exch-Ca, TN, silt content, EC, sand content, and Acid-BC, while the Rootzone of SC were Acid-BC, Exch-Na, silt content, clay content, Avail-P, and TN (Table S5).

3.5. Influence of Climate on Physicochemical Properties of ASSs

Stepwise MRA was used to select variables of climatic parameters and physicochemical properties from CDV results in order to determine the effects of climatic parameters on physicochemical properties of ASS under oil palm planting areas in CP and SC (Table 8).

In CP, the results showed (i) clay content (F = 5.46, ) in Topsoil, CEC (F = 9.58, ) in Subsoil-1, pH1:1-H2O (F = 12.79, ) in Subsoil-2, and silt content (F = 31.59, ) in Rootzone, were significantly influenced by daily solar radiation; (ii) Acid-BC in Rootzone was significantly influenced by annual water deficit. In other words, (i) daily solar radiation accounted for 33% of clay content in Topsoil, 49% of CEC in Subsoil-1, 57% of pH1:1-H2O in Subsoil-2, and 51% of silt content in Rootzone; (ii) annual water deficit accounted for 11% of Acid-BC in Rootzone.

In SC, the results showed (i) clay content (F = 8.55, ) in Topsoil, AWC (F = 7.75, ) in Subsoil-2 and clay content (F = 28.23, ) in Rootzone were significantly influenced by annual minimum temperature; (ii) clay content (F = 7.29, ) in Subsoil-1 and Acid-BC (F = 6.53, ) in Rootzone were significantly influenced by annual water vapor pressure; (iii) Exch-Mg (F = 20.30, ) in Subsoil-2 and Exch-Na (F = 10.05, ) and silt content (F = 7.85, ) in Rootzone were significantly influenced by daily solar radiation. These indicated that (i) annual minimum temperature accounted for 46% of clay content in Topsoil, 43% of AWC in Subsoil-2, and 48% of clay content in Rootzone; (ii) annual water vapor pressure accounted for 41% of clay content in Subsoil-1 and 16% of Acid-BC in Rootzone; (iii) daily solar radiation accounted for 68% of Exch-Mg in Subsoil-2, and 24% of Exch-Na and 19% of silt content in Rootzone.

4. Discussion

4.1. Suitability of Climate for Oil Palm

Due to the climatic suitability classification [19], oil palm yields in SC may decrease because the median annual precipitation (Table 5) is greater than 2,500 mm (Table 3). Too high continuous rains reduce the population of pollinating weevils and subsequently decrease the yields of oil palm because of more precipitation and associated cloudy conditions affecting less radiation can reach the surface [49]. Due to the climatic suitability classification [19] and water deficit suitability classification [25], oil palm plantation in CP should have to compensate precipitation rate and decrease water deficit because the median annual precipitation (Table 5) is lower than 2,000 mm (Table 3) and the median annual water deficit (Table 5) is greater than 150 mm (Table 3). In Ghana, a previous study by Rhebergen et al. [25] demonstrates that optimal areas (Table 3) for oil palm are in the south of the Western Region and a smaller area west of Koforidua in the Eastern Region which estimated at 580,000 ha or 2% of the total land areas.

It would be better for the oil palm growers in CP and SC to manage water by using raised beds with irrigation canals to compensate unsuitable precipitation and decrease water deficit in CP. This also include draining water excess before flooding, and controlling water table level carefully to prevent aluminum toxicity [12].

4.2. Suitability of ASSs for Oil Palm

In both CP and SC, their pH1:1-H2O2 of Topsoil, Subsoil-1, and Subsoil-2 (Table 7) are extremely lower than their pH1:1-H2O because the lower pH1:1-H2O2 values are caused by oxidation of iron sulfides [29, 50] and soil organic matter [51]. Hydrogen peroxide causes not only the oxidation of OM but also the degradation of clay minerals, especially 2 : 1 clay minerals [51]. Therefore, both ASS areas require soil pH adjustment so as to reduce aluminum toxicity affecting oil palm roots. In addition, considering the median of pH1:1-H2O in each depth level of CP and SC (Table 7) by soil limitation [48], their soil pH is not in the desirable range except Subsoil-2 and Rootzone of SC (Figure 4). As a result, the pH level should be increased to at least pH 4.3 [52] and not to exceed about pH 6.0, with optimal in the range of pH 5.0–5.5 [53] by liming at an appropriate rate and applying organic materials, especially green manure [54]. Liming improve soil conditions and increase oil palm yields in ASS [12], whereas pH neutralization is faster by sulfate reducing bacteria when adding high-organic carbon in anoxic condition [55, 56]. Besides, high acid buffering capacity (Acid-BC) of both areas (Table 7) was good to maintain their soil pH when the pH was raised to the desired level.

The very high Extr-Al (extractable aluminum) contents of both CP and SC are due to their very low soil pH levels (Table 7). At the low pH, toxic aluminum species can form that inhibit the growth of plant roots and subsequently limit ecosystem productivity [57]. Aluminum toxicity is the most important limiting factor for growing plant in ASS [58] and has different effects on the length of primary roots of oil palm in different varieties [59]. Therefore, before applying fertilizer or amendment in ASSs, it is essential to consider (i) the rooting depth of oil palm which can be impeded by dense layers of lateritic or other gravel, solid rock, or unfavorable chemical conditions [17], and (ii) the depth of acidic horizon which oil palm can grow without acidity limitation when the presence of acidic layers occur below 100 cm from soil surface [60, 61]. Considering the rooting depth covering the total 100 cm from soil surface with the interpretation of soil nutrient status for fertilizer recommendations [48] (Figure 5). The growers in CP should maintain N and P fertilizer rate, and not input K and Mg fertilizer except to correct nutrient imbalance because the Rootzone had moderate TN, low Avail-P, and very high Exch-K and Exch-Mg. The input K fertilizer rate is different from the recent study by Palykaew et al. [62] that the rate should be sufficient due to the low to very high of Avail-K (available potassium) which tends to increase with depth levels, whereas the growers in SC should decrease N and K fertilizer rates, maintain P fertilizer rate, and not input Mg fertilizer except to correct nutrient imbalance because the Rootzone had high TN and Exch-K, very low Avail-P, and very high Exch-Mg. The P fertilizer rate of both ASSs areas should not be increased because the reactions that fix P in relatively unavailable forms are related to soil pH [63].

4.3. Influence of Climate on Physicochemical Properties of ASSs

Physical properties of soils are grossly static in their response to change in climatic conditions, whereas chemical properties are more dynamic to changes in climatic conditions [64, 65]. The findings of this study indicated that (Table 8) (i) in CP, daily solar radiation was the most significant climate parameter that influenced ASS properties in all depth levels; it influenced clay content in topsoil, CEC in the range of topsoil-60 cm, soil pH by water in the range of 60–100 cm, and silt content covering the overall 100 cm from soil surface; (ii) in SC, annual minimum temperature was one of the most significant climate parameters that influenced ASS properties in almost all of the depth levels. It influenced clay content in topsoil, AWC in the range of 60–100 cm, and clay content covering the overall 100 cm from soil surface, whereas daily solar radiation influenced Exch-Mg in the range of 60–100 cm, Exch-Na and silt content covering the total 100 cm from soil surface. Therefore, it is important to emphasize the daily solar radiation in both areas because it was the major influential one of climatic parameters on soil properties due to the transport of heat in soils which is mostly governed by solar radiation during day time and the soil radiation to atmosphere during night time [66]. In addition, annual minimum temperature was the major influential of climatic parameters on soil properties in SC. These results indicated influence of soil temperature from daily solar radiation and annual minimum temperature on ASS properties due to their effect on the rates of physical, chemical, and biological reactions and process in the soil [67, 68].

On the other hand, particle size distribution was the most significant ASS property in both areas which was influenced by climatic parameters due to their clay and silt fractions (Table 8). Clay fraction is the most active and important fraction of soils because it retains moisture and nutrients, including its minerals that determine the potential fertility of soils via their CEC [69], whereas the silt fraction has not only a larger surface area but also a faster weathering rate and release of soluble nutrients for plant, compared with the sand fraction [70].

Climate can have an important influence on soil properties and their dynamics, especially precipitation and temperature [71]. Many significant correlations reported about the increasing precipitation and/or decreasing temperature with increment in soil organic carbon [7275]. In addition, contents of Mg, Fe, Ca, and Al decreased with increasing precipitation while contents of Na and K increased with increasing temperature [76]. Our study indicated the importance of soil temperature on changing of ASS properties affecting directly or indirectly effective oil palm plantation. For example, magnesium (Table 8), which is influenced by daily solar radiation, is one of the key nutrients affecting oil palm growth and yield [22]. Soil temperature also affects the production of CO2 in soils which influence the rates of microbial activity [77] and roots respiration [78].

5. Conclusion

The overall results of Thai ASSs under oil palm planting areas in CP and SC indicated that daily solar radiation not only was the most significant climate parameter but it also influenced silt content (mainly quartz) covering the overall 100 cm from soil surface of both ASS areas. Particle size distribution was the most significant ASS property in both areas that was influenced by climatic parameters. These results indicated the importance of soil temperature on changing ASS properties affecting directly or indirectly on effective oil palm plantation and suggested that sustainable oil palm plantation in ASS under tropical savanna and tropical monsoon climates should include management of water use by devising raised beds with irrigation canals to compensate precipitation, decrease water deficit, draining water excess before flooding, and controlling water table level carefully to prevent aluminum toxicity. Besides, the rooting depth of oil palm and depth of acidic horizon should be considered before applying fertilizer or amendment. There is a need to improve soil conditions and increase yields by liming along with integrated organic material management to raise the soil pH to at least 4.3 but not exceed 6.

Data Availability

All data generated or analyzed during this study are included within the article as the supplementary information files.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors would like to thank Chaipat Khongkeaw for his assistance on soil analyses and Yutika Chanlaksana for her advice on chemical analyses. Funding for this research project was provided by the Thailand Research Fund (TRF).

Supplementary Materials

Table S1: results of principal component analysis with concept of component defining variable on climate parameters (n = 10 and 30). Table S2: results of principal component analysis with concept of component defining variable on soil properties of topsoil (n = 10). Table S3: results of principal component analysis with concept of component defining variable on soil properties of subsoil-1 (n = 10). Table S4: results of principal component analysis with concept of component defining variable on soil properties of subsoil-2 (n = 10). Table S5: results of principal component analysis with concept of component defining variable on soil properties of rootzone (n = 30). (Supplementary Materials)