Next Article in Journal
Labile Soil Organic Matter Pools Are Influenced by 45 Years of Applied Farmyard Manure and Mineral Nitrogen in the Wheat—Pearl Millet Cropping System in the Sub-Tropical Condition
Next Article in Special Issue
Developing and Testing Remote-Sensing Indices to Represent within-Field Variation of Wheat Yields: Assessment of the Variation Explained by Simple Models
Previous Article in Journal
The Shape Diversity of Olive Stones Resulting from Domestication and Diversification Unveils Traits of the Oldest Known 6500-Years-Old Table Olives from Hishuley Carmel Site (Israel)
Previous Article in Special Issue
Stability and Adaptability of Maize Hybrids for Precision Crop Production in a Long-Term Field Experiment in Hungary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of El Niño on Oil Palm Yield in Malaysia

1
Centre of Disaster Risk Reduction (CDRR), Lee Kong Chian Faculty of Engineering & Science, Civil Engineering Department, Universiti Tunku Abdul Rahman, Jalan Sungai Long, Kajang 43000, Malaysia
2
Centre for Environmental Sustainability and Water Security, Universiti Teknologi Malaysia, Skudai 81310, Malaysia
3
Centre for Mathematical Sciences (CMS), Lee Kong Chian Faculty of Engineering & Science, Department of Mathematical and Actuarial Sciences, Universiti Tunku Abdul Rahman, Jalan Sungai Long, Kajang 43000, Malaysia
4
Department of Liberal Arts and Sciences, American Degree Programme, Taylor’s University, No. 1, Jalan Taylors, Subang Jaya 47500, Malaysia
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(11), 2189; https://doi.org/10.3390/agronomy11112189
Submission received: 6 October 2021 / Revised: 26 October 2021 / Accepted: 26 October 2021 / Published: 29 October 2021
(This article belongs to the Special Issue Crop Yield Prediction in Precision Agriculture)

Abstract

:
Oil palm crop yield is sensitive to heat and drought. Therefore, El Niño events affect oil palm production, resulting in price fluctuations of crude palm oil due to global supply shortage. This study developed a new Fresh Fruit Bunch Index (FFBI) model based on the monthly oil palm fresh fruit bunch (FFB) yield data, which correlates directly with the Oceanic Niño Index (ONI) to model the impact of past El Niño events in Malaysia in terms of production and economic losses. FFBI is derived from Malaysian monthly FFB yields from January 1986 to July 2021 in the same way ONI is derived from monthly sea surface temperatures (SST). With FFBI model, the Malaysian oil palm yields are better correlated with ONI and have higher predictive ability. The descriptive and inferential statistical assessments show that the newly proposed FFBI time series model (adjusted R-squared = 0.9312 and residual median = 0.0051) has a better monthly oil palm yield predictive ability than the FFB model (adjusted R-squared = 0.8274 and residual median = 0.0077). The FFBI model also revealed an oil palm under yield concern of the Malaysian oil palm industry in the next thirty-month forecasted period from July 2021 to December 2023.

1. Introduction

The edible oil and fat sectors consist of 17 varieties of oils and fats that can be categorized into vegetable oils and animal fats. As seen in Figure 1, global edible oil and fat output increased by around 176% in the span of 26 years, from 84.6 million tons (Mt) in 1992/93 to 233.3 Mt in 2018/19. That equates to a mean annual increment of 5.7 Mt [1].
Palm oil made a significant 40% rise in demand, which was an increase of about 2.3 Mt per year in the whole edible oil and fat sectors within the aforementioned period [1]. The said increase was a consequence of the growing demand for palm oil consumption within the food sector, specifically household consumable oil, the cosmetics sector and biodiesel fuel. Typically, vegetable oils are regarded to be an effective replacement for fossil fuels that can lead to a renewable energy source [2].
Palm oil is by far the most manufactured oil relative to many other edible oils and fats, accounting for nearly 31.3% of total edible oils and fats manufactured [1]. As presented in Table 1, oil palm is the most efficient in terms of edible oil yield per hectare with yields of about 4.27 tons per hectare per year (t/ha/yr) (combined palm and palm kernel oil) and at least six times more productive than any other oilseed [3,4]. This claim is supported by the national oil yield of 4.02 t/ha/yr in Malaysia (Appendix A: 3.66 t/ha/yr for palm oil and 0.36 t/ha/yr for palm kernel oil) [5,6,7]. Its input-output ratio is greater than that of other edible oils and fats. As a result, a smaller plantation area is required for oil palms to produce the same quantity of edible oil than other crops. Oil palm cultivation is also extremely alluring given the limited land resources available globally. In addition, the production cost for palm oil is about 700 USD per ton (USD/t), much cheaper compared to other vegetable oils. The next lowest cost is 850 USD/t for rapeseed oil [8,9,10].
Oil palm crop is a crucial income resource in Malaysia, as the revenue from palm oil exports amounted to 5.16% of the country’s Gross Domestic Product (GDP) in 2020 [11,12,13]. Situated between the Indian and Pacific Oceans, Malaysia, being the second largest exporter (after Indonesia) of palm oil in the world [14] is susceptible to El Niño and La Niña occurrences. El Niño (warm phase) and La Niña (cool phase) are phases of a larger phenomenon called the El Niño-Southern Oscillation (ENSO). During non-El Niño conditions (ENSO neutral), trade winds blow westwards across the tropical Pacific Ocean, pushing warm surface water towards the western Pacific (Asia and Australia). During an El Niño event, trade winds which normally blow westwards weaken or reverse direction, now blowing eastwards across the tropical Pacific Ocean, bringing warm water towards the coasts of South America. As warm water accumulates at the coasts of South America, convection above the warm surface water causes warm air to rise and precipitation happens. As a result, drastically increased rainfall would be seen in Ecuador and northern Peru. On the other side of the tropical Pacific Ocean, El Niño brings droughts to Malaysia, Indonesia and Australia [15,16,17].
Scientists use the Oceanic Niño Index (ONI) to measure deviations from normal sea surface temperatures [17]. This study uses ONI to classify El Niño events from January 1986 to June 2021 as shown in Table 2. The ONI is the running three-month mean SST anomaly for the Niño 3.4 region (5° N–5° S, 120° W–170° W). El Niño events are characterized by five consecutive overlapping three-month periods at or above +0.5 °C anomaly [18,19]. In general, El Niño causes lower rainfall, higher temperatures and abnormally dry conditions in Malaysia [20,21,22]. This study analysed the impact of past El Niño events on the oil palm production in Malaysia.
Yield of fresh fruit bundles (FFB) and crude palm oil (CPO) are affected by El Niño as oil palm crop is sensitive to prolonged drought periods like other oil crops [22,23,24]. During the occurrence of El Niño, a high level of water stress is generated for the palm trees due to reduced rainfall and increased temperature [2].
Many researchers have looked into the relationship between rainfall, temperature and oil palm production in Malaysia [2,20,22,25,26,27,28,29]. However, the direct correlation between ENSO index and oil palm production has yet to be analysed. Furthermore, past studies lack the quantification study of actual financial losses due to the impact of El Niño on the Malaysian oil palm industry. Therefore, this study investigated the direct relationship between ONI and the monthly oil palm FFB yield. The objective of this study is to model the impact of El Niño on the palm oil production in Malaysia by proposing a newly derived Fresh Fruit Bunch Index (FFBI). The impacts are then quantified financially in terms of opportunity losses and real Malaysian GDP value. This study also aims to forecast the Malaysian oil palm yields in the near future using the newly proposed FFBI time series model.

2. Materials and Methods

2.1. Study Area

Located in Southeast Asia, Malaysia consists of two geological regions separated by the South China Sea. Peninsula Malaysia (or West Malaysia) shares a land border with Thailand to the north and Singapore to its south. Located on the northern part of the Borneo Island, East Malaysia borders Indonesia and surrounds the Sultanate of Brunei [30]. The total land area is 328,550 km2 [31]. In 2020, the total area planted with oil palm in Malaysia is 58,653 km2 (nearly 18% of the total land area) [7]. Figure 2 shows the location map of oil palm plantation areas in Malaysia.

2.2. The New Fresh Fruit Bunch Index (FFBI) Model

Fresh Fruit Bunch Index (FFBI) is created with the monthly Malaysian oil palm FFB yield data [5,6,7] (Appendix B) through a similar calculation method of ONI (Appendix D) which is based on the monthly sea surface temperatures (SST) tracked by the National Oceanic and Atmospheric Administration (NOAA) [18]. In this study, monthly oil palm FFB yields from January 1986 to July 2021 (total 427 data points) in Malaysia were used [5,6,7]. The national average FFB yield for each month was calculated on the month-to-month basis from 1986 to 2020 in Table 3.
FFBanomaly, the difference between FFB yield and average FFB yield (Table 3) at each corresponding month was calculated using Equation (1):
FFBanomaly = FFB − FFBaverage
The new FFBI (Appendix C) was then calculated by taking the 3-month running mean of FFBanomaly, as shown in Equation (2):
FFBIi = (FFBanomaly, i-1 + FFBanomaly, i + FFBanomaly, i+1)/3,
where FFBanomaly, i-1 and FFBanomaly, i+1 were the FFB anomalies one month before and after the month that was being calculated. Hence, from 427 points of FFB yields, a total of 425 points of FFBI was derived, from February 1986 to June 2021.

2.3. Correlation Test between Malaysian Oil Palm Yields and Oceanic Niño Index (ONI)

The relationships between FFB and FFBI with ONI were tested using correlation tests. Normality tests were conducted beforehand to determine whether parametric or non-parametric correlation tests should be used. If a given dataset has less than 2000 samples, the Shapiro-Wilk test is suggested rather than the Kolmogorov-Smirnov test. This paper used a dataset which has less than 2000 samples, and therefore, the Shapiro-Wilk test was used. If the p-value is greater than 0.05, then the dataset is considered to be normally distributed [33]. Based on the normality tests (Table 4), FFB and ONI data distribution were determined to be non-normally distributed, and FFBI data distribution was normally distributed. As such, non-parametric Spearman’s rho correlation test was used to investigate the relationship between the oil palm yields and ONI in a lag period from 0 to 18 months. A Spearman’s rho correlation is also referred as Spearman correlation or Spearman rank correlation. The strength of association between the variables tested is expressed in a single value between −1 and +1, which is a bivariate correlation analysis. A positive correlation coefficient indicates a positive relationship between the two variables (as values of one variable increase, values of the other variable also increase), while a negative correlation coefficient expresses a negative relationship (as values of one variable increase, values of the other variable decrease). A correlation coefficient with a value of zero indicates that no relationship exists between the variables. The Spearman correlation does not assume that the variables are normally distributed, hence being used in this study as not all the variables were normally distributed [34,35]. All of the tests were conducted using IBM SPSS Statistics version 26.0 [36].

2.4. Time Series Forecasting

FFB and FFBI time series forecasting models with ONI as its predictor were created using the Expert Modeler in IBM SPSS Statistics version 26.0 [36]. The FFBI time series model was also validated repeatedly using a 30-year moving time frame. Six FFBI time series models were created and validated using different block periods (1986–2015, 1987–2016, 1988–2017, 1989–2018, 1990–2019 and 1991–2020). Each model forecasted FFBI one year ahead of the respective block period to be compared and validated with the observed FFBI. Residual analyses for each model were also performed in which the 99% confidence intervals of the residuals’ median, standard deviation and variance were computed using BCa bootstrapping method in this study. Using monthly Malaysian FFB data points from February 1986 to June 2021 (N = 425), the final FFB and FFBI time series forecasting models were created to forecast for another 30 months from July 2021 until December 2023.
R-squared measures the proportion of the variable in the dependent variable explained by all of the independent variables in the model, while adjusted R-squared measures the proportion of variation explained by only those independent variables that really help in explaining the dependent variable. As such, the adjusted R-squared will decrease if independent variables that do not help in predicting the dependent variable are added to the model [37]. The adjusted R-squared (R2adj) was derived from the R-squared produced by the models using Equation (3):
R 2 adj = 1 [ ( 1 R 2 ) ( n 1 ) n k 1 ] ,
where R2 is the R-squared of the models, N is the number of data points and K is the number of independent regressors [38].
To compare the FFB and FFBI time series models, residual analyses were conducted using IBM SPSS Statistics version 26.0. Residual between the predicted and observed data point was calculated [39,40]. Descriptive statistics of the residuals from both models were calculated for model comparison, including residuals’ skewness, range, median, standard deviation and variance. Shapiro-Wilk normality tests were also conducted for the residuals to determine whether the median or mean residuals should be referred for model prediction accuracy comparison. As the residuals from both models were determined to be non-normally distributed, the FFB and FFBI model’s median residuals were used for comparison assessment. The non-parametric inferential statistics of the bias corrected and accelerated (BCa) bootstrapping method was conducted with 2000 random samples (with replacement) to compute the 99% confidence intervals (CI) of the residuals’ median, standard deviation and variance.

2.5. Mann-Kendall Trend Test and Sen’s Slope

The Mann-Kendall trend test was used to evaluate the trends of FFBI during the occurrence of each El Niño events. Mann-Kendall test statistics (S) can be calculated using Equation (4):
S = i = 1 n 1 · j = i + 1 n sign ( x i x k ) ,
where xi and xk are sequential data in the series and
sign ( x i x k ) = { + 1   when   ( x i x k ) > 0 0   when   ( x i x k ) = 0 1   when   ( x i x k ) < 0
The variance of S (Var(S)) is estimated as:
Var ( S ) = s ( n 1 ) ( 2 n + 5 ) p 1 q ( t p 1 ) ( 2 t p + 5 ) 18 ,
where tp defines the ties of the pth value, and q represents the number of the tied value. The standardized test static for the Mann-Kendall test (Z) can be calculated, as shown in Equation (7):
Z = { S 1 Var ( S ) when   S > 0 0   when   S = 0 S + 1 Var ( S ) when   S < 0
The sign of Z indicates the direction of the trend. The negative value of Z indicates a decreasing trend and vice versa. At the 5% significance level, the null hypothesis of no trend is rejected if the absolute value of Z is higher than 1.64 [41]. This study used the excel template developed by Finnish Meteorological Institute [42] to conduct the Mann-Kendall trend tests and Sen’s slope estimation for FFBI during El Niño events.
The Sen’s slope method was used for the estimation of the monthly change in FFBI from February 1986 to June 2021. The Sen’s slope is calculated as the median of all the slopes estimated between all the successive data points of FFBI time series as [41]:
Sen s   Slope = median [ Δ FFBI Δ t ] ,
where ∆FFBI is the change of FFBI in the change of time, ∆t between two subsequent FFBI data.

2.6. Opportunity Loss Modelling of El Niño Events

During or after the occurrence of an El Niño event in the study period, the first clump of continuous monthly negative FFBI was considered to be the oil palm yield losses triggered by El Niño. The oil palm FFB yield losses (FFBL) due to the observed El Niño event were computed by summing up the absolute magnitude of the negative FFB. As such, the opportunity loss caused by the El Niño event was calculated using Equation (9):
OL = FFBL × Area × OER × Price,
where OL is the opportunity loss caused by El Niño in the unit of USD, area is the oil palm matured area (Appendix K), OER is the oil extraction rate (Appendix K) to convert oil palm from fresh fruit bunch to crude palm oil (CPO) and price is the CPO price (Appendix J) [5,6,7].
The opportunity losses for all El Niño events within the study period were projected to December 2021 for fair comparison using a hypothetical discount rate of 6%. This computation was done using FV function in Microsoft Excel software as shown in Figure 3,
where FV is the future value, rate is the discount rate per period, NPER is the total number of payment periods, PMT is the payment made for each period, PV is the present value and type is the representation of the timing of payment (1 for payment at the beginning of the period and 0 for payment at the end of the period).

3. Results

3.1. Predictive Ability of Fresh Fruit Bunch Index (FFBI) Model

3.1.1. Correlation with Oceanic Niño Index (ONI)

The newly proposed FFBI model (Appendix C) has an improved ability to model the impact of El Niño events on the oil palm production in Malaysia compared to FFB model. Based on the non-parametric Spearman’s rho correlation test (Table 5), FFBI has statistically significant correlations with ONI at lag periods from 2 to 16 months, with the highest correlation of −0.399 at 0.01 alpha level. On the other hand, the correlations of FFB and ONI are only significant at lag periods from 6 to 13 months, with the highest correlation of −0.217 at 0.01 alpha level.

3.1.2. Time Series Forecasting

Based on the statistical assessment of the FFB and FFBI time series model (Table 6), the FFBI model (Appendix F) has an adjusted R-squared of 0.9312, which is significantly higher than the FFB model (Appendix E) with adjusted R-squared of 0.8274 only. Compared to the FFB model with residual sum of squares (RSS) of 4.6876, the FFBI model has much lower RSS at 0.5459 only. The FFBI model also has a lower median residual (0.0051) than the FFB model (0.0077). Other than that, residual of the FFBI model has lower skewness, range, standard deviation and variance compared to the FFB model.

3.2. Production, Financial and Economical Loss in the Malaysian Oil Palm Industry due to El Niño

With the application of FFBI, oil palm FFB trends from January 1986 to June 2021 were studied. Based on the Mann Kendall’s (MK) trend test (Table 7), almost every El Niño event in the aforementioned period is accompanied by a significant decreasing trend of FFBI (at least 0.05 alpha level), except for 2002/03′s moderate El Niño and 2014/15′s weak El Niño. Those downward FFBI trends could be observed from 0 to 6 months after an El Niño event (Table 8).
When the FFBI plummets below zero, oil palm production would suffer losses due to El Niño. Those losses could last between 4 to 25 months (Table 8) depending on the severity of the El Niño events, before recovering to normal productions. Both very strong El Niño events caused oil palm yield losses of about 3.0739 ton per hectare (t/ha) in 1997/98 and 3.5197 t/ha in 2015/16. Those losses equate to an average oil palm yield loss of 3.2968 t/ha in Malaysia, causing approximately 3.07 billion USD of opportunity losses when projected to December 2021 using a 6% discount rate (Table 9).
The oil palm production and opportunity losses decrease as the severity of El Niño events decrease (Figure 4). The oil palm yields were estimated to reduce by 1.1250 t/ha during a strong El Niño event, 0.5974 t/ha during a moderate El Niño event and about 0.2953 t/ha during a weak El Niño event on average. In financial terms, a strong El Niño event could cause opportunity losses of approximately 522 million USD, 416 million USD for a moderate El Niño event and 231 million USD for a weak El Niño event in Malaysia (Table 9).
This study then accumulates the opportunity losses at a hypothetical 6% discount rate to show the collective impact of every El Niño event from January 1986 to June 2021 (except for 2018/19 weak El Niño) on the oil palm industry in Malaysia (Figure 5). The total combined loss is estimated to be 9.55 billion USD by the end of 2021, which is approximately 2.84% of Malaysia’s GDP in 2020 [11]. Using SPSS, the cumulative opportunity losses model was modelled with a cubic model with an adjusted R-squared of 0.985, standard error of estimate of 0.331 billion USD and significant at p-value less than 0.0001 (Figure 5).

4. Discussion

4.1. Predictive Ability of the New Fresh Fruit Bunch Index (FFBI)

Based on Spearman’s rho correlation tests (Table 5), both FFB and the proposed FFBI have negative correlations with ONI. This shows that when ONI increased during past El Niño events, water stress was induced in palm trees, causing oil palm production reduction in Malaysia [2,20,21,22,25]. The new FFBI has significantly higher correlation with ONI compared to FFB. By converting monthly FFB data into FFBI, the impact of El Niño on the monthly oil palm production in Malaysia could be modelled with higher precision.
The comparison between the FFB and FFBI time series model (Table 6) shows that FFBI model has better predictive ability than FFB model. Based on the residual analyses, both models have median residuals with 99% confidence intervals which span across zero, indicating their capabilities to produce a monthly oil palm production prediction with low error [40]. However, the FFBI time series model has a higher adjusted R-squared, lower RSS, median residual, residual skewness, residual range, residual standard deviation and residual variance than the FFB model. By having higher correlation with ONI and significantly better monthly oil palm yield forecasting ability, FFBI model is a better predictive modelling tool for the Malaysian oil palm industry stakeholders to better understand the impact of El Niño events. With the application of FFBI model, decision making processes could be enhanced so that appropriate measures could be taken to prepare for the adverse effects of El Niño on oil palm production in the future.
The validation tests using six different FFBI models in different time frames showed that the proposed FFBI model is consistent. The adjusted R-squared of all the validated models are consistent and stable, in the range of 0.916 to 0.9354 (Appendix G). The residuals of all the models have median residuals with 99% confidence intervals which span across zero, further supporting the forecasting ability and stability of the proposed FFBI model. Furthermore, the predicted FFBI produced by the validation models tend to follow the trend of real observed FFBI (derived from historical FFB [5,6,7], see Appendix H). The FFBI model also outperformed the FFB model as the forecasted FFBI (Appendix F) from 2021 to 2023 is closer to the observed FFBI from January 2021 to June 2021 (Appendix I).

4.2. Trends of Oil Palm Yields during El Niño Events

In the Spearman’s rho correlation tests (Table 5), FFBI shows a significant correlation with ONI with lag periods from 2 to 16 months. This is on par with other studies that reported lag periods from 3 to 24 months to witness oil palm production reduction in Malaysia due to El Niño events and drought [2,20,22,25,26,27,28,29]. Since there is no fixed way to explain the nature, the El Niño events could affect the oil palm yields after a period of time and are not constrained to a single fixed period. The FFBI model of this study suggests that the oil palm production could be affected after 2 to 16 months of the occurrence of El Niño events in Malaysia.
Using the FFBI model, trends of oil palm yield in Malaysia during each El Niño event from January 1986 to June 2021 were studied on an event-by-event basis. In general, oil palm yields started to decrease after 0 to 6 months of an El Niño event, according to the Mann Kendall’s trend tests (Table 7). The industry would suffer production losses when the oil palm yields continue to plummet, and this could last for 4 to 25 months depending on the severity of El Niño events (Table 8) before the oil palm yields would recover.

4.3. Impact of El Niño Events on the Malaysian Oil Palm Industry

As shown in Figure 4, very strong El Niño events have a significantly bigger impact on the Malaysian oil palm production compared to other categories of El Niño. Averaged oil palm reduction during very strong El Niño events is about 2.9 times more than the production loss in strong El Niño, 5.5 times more than the moderate El Niño and 11.2 times more than weak El Niño (Table 9). When a very strong El Niño happens, oil palm yields in Malaysia are expected to suffer production losses for about 21 months continuously. For other types of El Niño events, oil palm yields are expected to be negative in terms of FFBI for approximately 6 to 9 months. The total combined opportunity losses from the past 11 El Niño events in Malaysia since 1986 are estimated to be 9.55 billion USD when projected to December 2021. This equates to 2.84% of the Malaysia GDP in 2020 [11] (Figure 5). Undeniably, extreme temperature caused by El Niño is a major factor affecting the oil palm production in Malaysia.
However, the occurrence of El Niño event cannot be controlled. The adverse weather condition caused by El Niño cannot be remedied by concerted intervention [33]. To enhance performance of the oil palm industry in Malaysia, other underlying factors affecting oil palm production should be further investigated. These factors include labour shortage and ageing palm trees [1,43]. Based on the annual Malaysian oil palm FFB yields from 1986 to 2020 [5,6,7] (Appendix K), the FFB yields did not recover to the national average after the 2015/16 very strong El Niño. The Malaysian oil palm FFB yields continued to stay below the mean lower confidence limit (BCa bootstrapped 99% CI) continuously from 2017 to 2020, facing a downward trend, even though there was only a weak occurrence of El Niño in 2018/19 [18,19]. In addition, based on the FFBI validation models, the observed FFBI (derived from historical FFB) from 2016 to 2021 is found to be mostly below the predicted FFBI (Appendix I), showing the possibility that there might be other significant underlying factors affecting oil palm production in Malaysia.
Surely, El Niño would not be the reason to cause the current production downtrend in the Malaysian oil palm industry. One of the reasons might be labour shortage. Oil palm plantations struggle to hire foreign workers for the harvesting of oil palm, while locals shun the idea of working in plantations. The shortage of manpower became worse under the prolonged foreign labour hiring freeze due to the COVID-19 situation and caused more losses to the industry [44]. Another underlying factor is the ageing of palm trees due to the slow replanting process in Malaysia [45]. In 2020, about 35% of the oil palm plantations owned by a private plantation firm (IOI Corporation Berhad) in Malaysia were of prime age, while 37% already passed their prime production phase [46]. Replanting schemes are not attractive to plantation owners due to the lower profit margin of the oil palm sector in the recent years. Backlog of old oil palm trees slow down the improvement of crop yield in Malaysia [47]. According to the United States Department of Agriculture (USDA) [43], aged oil palm trees are to be replaced with new high yield varieties to increase the crop yield. However, many plantation owners are reluctant to replace aged trees in times of high palm oil prices in order to continue gaining profits. Researchers [8] also supported the idea that a proper replanting program is required to improve the oil palm yield in Malaysia. Based on the historical records, oil palm yields dropped during the 1997/98 very strong El Niño and successfully recovered about a year after. However, the oil palm yields maintained at low production since the 2015/16 very strong El Niño. It is possible that the aged palm trees suffered damage from the heat stress of recent very strong El Niño and failed to recover.

4.4. Global Applicability of FFBI

For the previous ten years, more than USD 50 billion has been invested in the Malaysian and Indonesian oil palm sectors [48]. In 2019, almost 88% of the world’s palm oil production comes from the top three palm oil producing countries, Indonesia, Malaysia and Thailand [14]. They are located in the same region (Southeast Asia) with the same climatic pattern [49,50,51]. Hence, the supply of palm oil in the global market will be heavily affected in the event of extreme weather such as El Niño event. Thus, the newly proposed FFBI can be utilized to model the threat of El Niño in the oil palm industry, so that these countries can brace themselves for the adverse effects of intense climate conditions on the oil palm yields in the industry.

5. Conclusions

This study proposes a new oil palm yield index, namely Fresh Fruit Bunch Index (FFBI) that is derived from the monthly oil palm FFB yields. It has shown improvement in the trend and oil palm yield modelling in Malaysia with respect to El Niño events. The FFBI shows a significantly higher correlation with ONI compared to FFB yields based on Spearman’s rho correlation tests at 0.01 alpha level. The FFBI time series model has also shown a significantly higher predictive accuracy when compared to the FFB model. The new FFBI model has an adjusted R-squared of 0.9312, which is significantly higher than the FFB model with adjusted R-squared of 0.8274 only. The FFBI model also has smaller errors based on residual analyses.
The proposed FFBI provides an improved method to model the impact of El Niño events in the Malaysian oil palm industry, as financial losses caused by El Niño events cannot be neglected. A very strong El Niño event is estimated to cause 3.07 billion USD of opportunity losses when projected to December 2021 using a 6% discount rate. The total combined opportunity losses of 13 past El Niño events, starting from 1986 (excluding 2018/19 weak El Niño) is approximately 9.55 billion USD, which is about 2.84% of Malaysia’s 2020 GDP.
The FFBI model reveals that the recent under yield situation of the Malaysian oil palm yield may continue in future. This study further shows that El Niño is not the sole factor affecting the Malaysian oil palm production. Hence, other underlying factors affecting Malaysian oil palm yields such as ageing oil palm trees should be further investigated for the implementation of preventive steps to preserve the sustainable competitive edge of the Malaysian oil palm industry in the international arena. Future work should incorporate Indonesian and Thailand monthly FFB dataset to study the possible global production impact as these countries (Indonesia, Malaysia and Thailand) in the same region produce around 88% of the global palm oil supply.

Author Contributions

Conceptualization, L.L. and Z.Y.; methodology, L.L. and J.F.K.; software, L.L. and J.F.K.; validation, L.L., W.L.T. and Z.Y.; formal analysis, L.L. and J.F.K.; investigation, L.L. E.Z.X.S. and J.F.K.; resources, J.F.K., L.L. and Z.Y.; data curation, L.L. and J.F.K.; writing—original draft preparation, J.F.K.; writing—review and editing, L.L., J.F.K., J.L.L. and Z.Y.; visualization, L.L. and J.F.K.; supervision, L.L., W.L.T. and Z.Y.; project administration, L.L. and Z.Y.; funding acquisition, L.L. and Z.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by Ministry of Higher Education (MoHE) through Fundamental Re-search Grant Scheme (FRGS/1/2021/WAB07/UTAR/02/1), Universiti Tunku Abdul Rahman (IPSR/RMC/UTARRF/2019-C2/L07) and partly supported by the Brunsfield Engineering Sdn. Bhd., Malaysia (Brunsfield 8013/0002 & 8126/0001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Oil palm data can be found in Palm Oil Registration & Licensing Authority (PORLA) reports “PORLA palm oil statistics” from 1986 to 1999 and Malaysian Palm Oil Board (MPOB) reports “Malaysian oil palm statistics” from 2000 to 2020 in MPOB libraries located in Kelana Jaya and Bangi, Malaysia. Recent oil palm data is available at <https://bepi.mpob.gov.my/index.php/en/> (accessed 17 August 2021). ONI data is available at <https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php> (accessed 25 October 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Annual productivity of crude palm oil (CPO) and crude palm kernel oil (CPKO) from 1975 to 2020 in Malaysia [5,6,7].
Table A1. Annual productivity of crude palm oil (CPO) and crude palm kernel oil (CPKO) from 1975 to 2020 in Malaysia [5,6,7].
YearCPO Yield
(t/ha)
CPKO Yield
(t/ha)
YearCPO Yield
(t/ha)
CPKO Yield
(t/ha)
19753.660.1720003.460.41
19763.480.1620013.660.44
19773.540.1820023.590.40
19782.950.1720033.750.43
19793.650.2120043.730.42
19803.780.2220053.800.45
19813.760.2220063.930.47
19823.830.2820073.830.44
19833.430.2920084.080.47
19844.250.3220093.930.45
19854.330.3520103.690.42
19864.410.3620114.010.43
19873.390.3420123.840.43
19883.470.3420133.850.43
19893.880.3820143.840.42
19903.640.4120153.780.40
19913.480.3720163.210.34
19923.430.3720173.530.39
19933.780.4220183.420.39
19943.430.4120193.470.39
19953.500.4120203.330.38
19963.550.41Average:3.660.36
19973.630.40
19983.020.36
19993.580.40

Appendix B

Table A2. Monthly oil palm fresh fruit bunch (FFB) yield in tons per hectare from January 1986 to July 2021 in Malaysia [5,6,7].
Table A2. Monthly oil palm fresh fruit bunch (FFB) yield in tons per hectare from January 1986 to July 2021 in Malaysia [5,6,7].
JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
19861.681.371.381.341.291.291.651.641.951.941.671.33
19871.091.041.081.201.201.491.641.872.051.781.471.17
19880.950.971.191.341.381.671.691.851.781.681.571.41
19891.191.091.201.261.431.611.812.052.142.022.021.68
19901.371.261.351.341.721.551.721.881.931.741.511.22
19910.990.971.251.411.591.431.731.881.981.881.491.23
19921.121.071.201.351.411.421.671.751.921.771.711.41
19931.211.171.201.661.631.731.861.992.312.061.831.53
19941.381.161.231.291.301.411.531.711.921.981.771.69
19951.551.291.391.381.351.501.641.731.871.831.821.58
19961.331.061.431.401.641.641.631.861.901.941.641.43
19971.231.121.421.581.631.681.771.821.951.881.671.30
19981.031.181.151.201.271.381.491.611.651.481.411.10
19990.981.051.291.651.851.761.751.761.881.891.761.61
20001.381.261.261.321.331.351.391.621.882.021.961.56
20011.801.511.461.511.571.481.431.571.781.811.701.52
20021.471.211.341.321.401.431.511.681.871.791.571.38
20031.291.121.401.561.671.731.831.821.831.661.461.62
20041.251.181.251.331.441.531.691.711.971.771.681.74
20051.501.311.521.561.611.501.621.701.821.751.611.40
20061.231.341.511.611.701.641.671.871.971.701.901.43
20071.451.251.301.351.431.391.621.831.891.851.951.68
20081.681.421.511.521.621.651.781.791.761.901.871.68
20091.521.301.421.391.481.521.601.631.742.141.771.69
20101.441.221.441.381.441.491.641.711.661.721.551.34
20111.151.201.511.601.781.791.791.691.961.951.681.59
20121.351.181.211.251.341.431.681.672.001.961.971.85
20131.631.321.301.331.351.371.651.741.901.931.831.67
20141.491.211.391.451.531.431.551.941.801.791.681.37
20151.151.071.331.551.641.611.681.891.811.841.571.34
20161.040.971.101.171.271.451.471.541.571.511.461.36
20171.201.131.321.381.451.371.611.611.621.781.751.67
20181.481.221.381.361.321.171.281.391.621.701.631.61
20191.521.341.421.421.401.341.491.551.571.571.351.22
20201.061.101.191.411.441.651.591.621.651.501.331.19
20211.030.971.201.281.321.351.30

Appendix C

Table A3. Derived monthly Fresh Fruit Bunch Index (FFBI) in tons per hectare from February 1986 to June 2021 in Malaysia.
Table A3. Derived monthly Fresh Fruit Bunch Index (FFBI) in tons per hectare from February 1986 to June 2021 in Malaysia.
JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
1986NULL0.200.06−0.07−0.16−0.13−0.100.000.040.07−0.01−0.13
1987−0.17−0.21−0.20−0.24−0.17−0.100.040.110.10−0.01−0.18−0.29
1988−0.30−0.24−0.14−0.100.000.040.110.03−0.03−0.10−0.10−0.10
1989−0.10−0.12−0.12−0.11−0.030.070.190.260.270.280.250.20
19900.110.050.010.070.070.120.090.100.05−0.05−0.16−0.25
1991−0.27−0.21−0.100.010.010.040.050.120.110.00−0.12−0.21
1992−0.19−0.15−0.10−0.08−0.07−0.04−0.020.040.010.02−0.02−0.05
1993−0.06−0.080.040.090.210.200.230.310.320.290.150.09
19940.03−0.02−0.08−0.13−0.13−0.13−0.08−0.020.070.110.160.18
19950.180.130.05−0.03−0.06−0.05−0.010.000.010.060.090.09
19960.000.00−0.010.090.090.090.080.050.100.050.02−0.06
1997−0.07−0.020.070.140.160.150.130.100.080.05−0.04−0.16
1998−0.16−0.16−0.13−0.20−0.18−0.16−0.14−0.16−0.22−0.27−0.32−0.33
1999−0.28−0.170.020.190.290.240.130.050.040.060.100.09
20000.090.02−0.03−0.10−0.13−0.19−0.18−0.110.040.170.190.28
20010.300.310.190.110.05−0.05−0.14−0.15−0.08−0.020.020.07
20020.070.06−0.02−0.05−0.08−0.10−0.09−0.06−0.02−0.04−0.07−0.08
2003−0.06−0.010.050.140.190.200.160.08−0.03−0.13−0.07−0.05
20040.02−0.05−0.05−0.06−0.030.010.010.050.010.030.080.15
20050.190.170.160.160.090.03−0.02−0.03−0.05−0.05−0.07−0.08
20060.000.080.180.200.180.130.100.090.040.080.020.10
20070.050.06−0.01−0.04−0.08−0.06−0.020.040.050.120.170.28
20080.270.260.180.150.130.140.110.030.010.060.160.20
20090.170.140.060.030.00−0.01−0.05−0.090.030.100.210.14
20100.120.090.040.02−0.03−0.02−0.02−0.07−0.11−0.14−0.12−0.14
2011−0.100.010.130.230.260.240.130.070.060.080.090.05
20120.05−0.03−0.09−0.14−0.13−0.06−0.040.040.070.200.270.33
20130.270.140.01−0.08−0.12−0.09−0.040.020.050.110.160.17
20140.130.090.040.050.00−0.040.010.020.04−0.02−0.04−0.09
2015−0.13−0.090.010.100.130.100.100.050.04−0.04−0.07−0.17
2016−0.21−0.24−0.23−0.22−0.17−0.15−0.14−0.22−0.26−0.27−0.21−0.15
2017−0.10−0.06−0.03−0.02−0.07−0.07−0.10−0.13−0.13−0.060.080.14
20180.130.080.01−0.05−0.18−0.29−0.35−0.31−0.23−0.13−0.010.10
20190.160.150.090.01−0.08−0.13−0.17−0.21−0.24−0.28−0.27−0.28
2020−0.20−0.16−0.07−0.060.030.02−0.01−0.12−0.21−0.29−0.31−0.31
2021−0.26−0.21−0.16−0.14−0.15−0.22NULL

Appendix D

Table A4. Oceanic Niño Index (ONI) from January 1983 to June 2021 [18].
Table A4. Oceanic Niño Index (ONI) from January 1983 to June 2021 [18].
JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
19832.21.91.51.31.10.70.3−0.1−0.5−0.8−1.0−0.9
1984−0.6−0.4−0.3−0.4−0.5−0.4−0.3−0.2−0.2−0.6−0.9−1.1
1985−1.0−0.8−0.8−0.8−0.8−0.6−0.5−0.5−0.4−0.3−0.3−0.4
1986−0.5−0.5−0.3−0.2−0.10.00.20.40.70.91.11.2
19871.21.21.10.91.01.21.51.71.61.51.31.1
19880.80.50.1−0.3−0.9−1.3−1.3−1.1−1.2−1.5−1.8−1.8
1989−1.7−1.4−1.1−0.8−0.6−0.4−0.3−0.3−0.2−0.2−0.2−0.1
19900.10.20.30.30.30.30.30.40.40.30.40.4
19910.40.30.20.30.50.60.70.60.60.81.21.5
19921.71.61.51.31.10.70.40.1−0.1−0.2−0.3−0.1
19930.10.30.50.70.70.60.30.30.20.10.00.1
19940.10.10.20.30.40.40.40.40.60.71.01.1
19951.00.70.50.30.10.0−0.2−0.5−0.8−1.0−1.0−1.0
1996−0.9−0.8−0.6−0.4−0.3−0.3−0.3−0.3−0.4−0.4−0.4−0.5
1997−0.5−0.4−0.10.30.81.21.61.92.12.32.42.4
19982.21.91.41.00.5−0.1−0.8−1.1−1.3−1.4−1.5−1.6
1999−1.5−1.3−1.1−1.0−1.0−1.0−1.1−1.1−1.2−1.3−1.5−1.7
2000−1.7−1.4−1.1−0.8−0.7−0.6−0.6−0.5−0.5−0.6−0.7−0.7
2001−0.7−0.5−0.4−0.3−0.3−0.1−0.1−0.1−0.2−0.3−0.3−0.3
2002−0.10.00.10.20.40.70.80.91.01.21.31.1
20030.90.60.40.0−0.3−0.20.10.20.30.30.40.4
20040.40.30.20.20.20.30.50.60.70.70.70.7
20050.60.60.40.40.30.1−0.1−0.1−0.1−0.3−0.6−0.8
2006−0.9−0.8−0.6−0.4−0.10.00.10.30.50.80.90.9
20070.70.2−0.1−0.3−0.4−0.5−0.6−0.8−1.1−1.3−1.5−1.6
2008−1.6−1.5−1.3−1.0−0.8−0.6−0.4−0.2−0.2−0.4−0.6−0.7
2009−0.8−0.8−0.6−0.30.00.30.50.60.71.01.41.6
20101.51.20.80.4−0.2−0.7−1.0−1.3−1.6−1.6−1.6−1.6
2011−1.4−1.2−0.9−0.7−0.6−0.4−0.5−0.6−0.8−1.0−1.1−1.0
2012−0.9−0.7−0.6−0.5−0.30.00.20.40.40.30.1−0.2
2013−0.4−0.4−0.3−0.3−0.4−0.4−0.4−0.3−0.3−0.2−0.2−0.3
2014−0.4−0.5−0.30.00.20.20.00.10.20.50.60.7
20150.50.50.50.70.91.21.51.92.22.42.62.6
20162.52.11.60.90.4−0.1−0.4−0.5−0.6−0.7−0.7−0.6
2017−0.3−0.20.10.20.30.30.1−0.1−0.4−0.7−0.8−1.0
2018−0.9−0.9−0.7−0.5−0.20.00.10.20.50.80.90.8
20190.70.70.70.70.50.50.30.10.20.30.50.5
20200.50.50.40.2−0.1−0.3−0.4−0.6−0.9−1.2−1.3−1.2
2021−1.0−0.9−0.8−0.7−0.5−0.4

Appendix E

Figure A1. FFB time series, ARIMA (2,0,2) (1,1,1) model with ONI as predictor (adjusted R-squared = 0.8274). Modelled period: February 1986 to June 2021 (N = 425). 30 months forecasted period: July 2021 to December 2023 (N = 30, highlighted area of Figure A1). Note: Dash lines represent BCa 99% confidence interval range of monthly oil palm yield in Malaysia. FFB model predicts monthly oil palm yield to fluctuate around the BCa 99% confidence interval range in next 30 months.
Figure A1. FFB time series, ARIMA (2,0,2) (1,1,1) model with ONI as predictor (adjusted R-squared = 0.8274). Modelled period: February 1986 to June 2021 (N = 425). 30 months forecasted period: July 2021 to December 2023 (N = 30, highlighted area of Figure A1). Note: Dash lines represent BCa 99% confidence interval range of monthly oil palm yield in Malaysia. FFB model predicts monthly oil palm yield to fluctuate around the BCa 99% confidence interval range in next 30 months.
Agronomy 11 02189 g0a1

Appendix F

Figure A2. The new FFBI time series, ARIMA (3,1,3) (0,0,1) model with ONI as predictor (adjusted R-squared = 0.9312). Modelled period: February 1986 to June 2021 (N = 425). 30 months forecasted period: July 2021 to December 2023 (N = 30, highlighted area of Figure A2). Note: Dash lines represent BCa 99% confidence interval range of monthly oil palm yield in Malaysia. Unlike FFB model (Figure A1), FFBI model revealed a future monthly under yield concern of Malaysian oil palm industry.
Figure A2. The new FFBI time series, ARIMA (3,1,3) (0,0,1) model with ONI as predictor (adjusted R-squared = 0.9312). Modelled period: February 1986 to June 2021 (N = 425). 30 months forecasted period: July 2021 to December 2023 (N = 30, highlighted area of Figure A2). Note: Dash lines represent BCa 99% confidence interval range of monthly oil palm yield in Malaysia. Unlike FFB model (Figure A1), FFBI model revealed a future monthly under yield concern of Malaysian oil palm industry.
Agronomy 11 02189 g0a2

Appendix G

Table A5. FFBI time series forecasting validation models’ summary with residual analyses using descriptive and inferential statistics at α = 0.01.
Table A5. FFBI time series forecasting validation models’ summary with residual analyses using descriptive and inferential statistics at α = 0.01.
Modelled Period1986–20151987–20161988–20171989–20181990–20191991–2020
ModelARIMA (1,1,5) (0,0,1)ARIMA (0,0,4) (0,0,0)ARIMA (3,1,3) (0,0,1)ARIMA (0,0,4) (0,0,0)ARIMA (0,0,4) (0,0,1)ARIMA (3,1,3) (0,0,1)
R−Squared0.92850.91640.92810.91980.92380.9354
Adjusted R−Squared0.92830.91620.92790.91960.92360.9352
Residual Skewness−0.5314−0.1992−0.7146−0.3753−0.5032−0.7131
Residual Range0.23100.33720.22880.27410.24050.2265
Residual Median0.0001−0.00020.00010.00200.00150.0036
(BCa 99% CI)[−0.0044, 0.0063][−0.0058, 0.0083][−0.0052, 0.0063][−0.0039, 0.0064][−0.0048, 0.0081][−0.0028, 0.0079]
Residual Standard Deviation0.03420.03750.03450.03740.03670.0343
(BCa 99% CI)[0.0304, 0.0382][0.0331, 0.043][0.0303, 0.0387][0.0333, 0.0415][0.0326, 0.0408][0.0305, 0.0382]
Residual Variance0.00120.00140.00120.00140.00130.0012
(BCa 99% CI)[0.0009, 0.0015][0.0011, 0.0018][0.0009, 0.0015][0.0011, 0.0017][0.0011, 0.0017][0.0009, 0.0015]

Appendix H

Figure A3. FFBI time series forecasting validation models (total of 6 models combined). Modelled period: moving 30 years starting from 1986 (N = 360) with 12 months forecasted period. Note: Predicted FFBI (Red dash line) from January 2016 to December 2016 are values predicted from 1986–2015 model and so on with a moving 30 years validation.
Figure A3. FFBI time series forecasting validation models (total of 6 models combined). Modelled period: moving 30 years starting from 1986 (N = 360) with 12 months forecasted period. Note: Predicted FFBI (Red dash line) from January 2016 to December 2016 are values predicted from 1986–2015 model and so on with a moving 30 years validation.
Agronomy 11 02189 g0a3

Appendix I

Figure A4. Predicted values from 2021 to 2023 (N = 36) using FFB model (Figure A1) and proposed FFBI model (Figure A2) were compared to observed FFBI from January 2021 to June 2021 (N = 6).
Figure A4. Predicted values from 2021 to 2023 (N = 36) using FFB model (Figure A1) and proposed FFBI model (Figure A2) were compared to observed FFBI from January 2021 to June 2021 (N = 6).
Agronomy 11 02189 g0a4

Appendix J

Table A6. Monthly crude palm oil (CPO) price in Malaysian Ringgit (MYR) per ton from January 1986 to June 2021 [52].
Table A6. Monthly crude palm oil (CPO) price in Malaysian Ringgit (MYR) per ton from January 1986 to June 2021 [52].
JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
1986637476653553589573504448504723718712
19879007507486907257516787218007608831093
198810688959069161017122010119819791073985967
1989955942903922904810668699779787633659
1990696683695705700710701708707752863896
1991848850770775755768867791807883884914
1992888891905848848859826841868917938953
1993960968933875851820847852780787857992
1994970943974108911071086118512781404147816471580
1995144915511489131112961345141913211355141713921354
1996118511781214125011881077109912381142118112071255
1997129612641258128112301157116313051573170917842075
1998239221172145229522922323242824992234233323022156
199918541563152315861355123298813101239123512451176
200011331040119611491035103010331003901852822804
20017807518807667648651236106994099010971146
2002115511501152120414261355150214701361151015881644
2003161915941433135514161408128013391460176217371774
2004177819471937188816101550141915151413143814251387
2005127414031461142914121408137713701460144513911415
2006144315071437148114441491164115781559166119401995
2007190019602070221425812427260824202644288029303050
2008323240053395339534983598305026202090151516321695
2009177918952000259525602230218923702105220824722663
2010244525952556255824362373251725702730306134123788
2011380934723326327033933072309630092905293830183175
2012307832703433347131013020298030192546249623702438
2013255723972378228623972344223624042320259326542659
2014255928002634262324232426225719292217230621722266
2015214623052165210222162229212019912375236323442485
2016244325482725259326202353231625262636276330733109
2017302927702646250824992459267527062695281526032503
2018249225592425236224292326219422482174215020402121
2019229921212106209520691951207022342135248527443052
2020260423192402208822922297267727382714301133053600
2021349037423612386839193599

Appendix K

Table A7. Annual oil palm matured area, oil extraction rate (OER) and fresh fruit bunch (FFB) yield from 1986 to 2020 in Malaysia [5,6,7].
Table A7. Annual oil palm matured area, oil extraction rate (OER) and fresh fruit bunch (FFB) yield from 1986 to 2020 in Malaysia [5,6,7].
YearMatured Area (ha)OER (%)FFB Yield (t/ha)
19861,360,57919.6222.15
19871,373,14719.8717.10
19881,530,90619.8717.52
19891,672,09619.7819.57
19901,746,05419.6418.53
19911,826,26719.4717.85
19921,890,26819.2117.83
19932,020,51618.6720.26
19942,144,08018.6318.42
19952,243,06518.5118.93
19962,353,14718.7118.95
19972,513,18319.0319.10
19982,638,02018.9115.98
19992,856,70118.6019.26
20002,941,79118.8618.33
20013,005,26719.2219.14
20023,188,30719.9117.97
20033,303,13319.7518.99
20043,450,96020.0318.60
20053,631,44020.1518.88
20063,703,25420.0419.60
20073,764,38920.1319.03
20083,915,92420.2120.18
20094,075,70220.4919.20
20104,202,21320.4518.03
20114,281,83720.3519.69
20124,352,87220.3518.89
20134,526,08920.2519.02
20144,689,32120.6218.63
20154,859,39720.4618.48
20165,001,43820.1815.91
20175,110,71319.7217.89
20185,189,34419.9517.16
20195,216,82220.2117.19
20205,231,74319.9216.73
Annual FFB yields highlighted in red are the ones below the mean 99% confidence interval [18.02, 19.08]. All FFB yields below the mean confidence limit are associated with El Niño years, except for 2017 and after. Although there was a weak El Niño in 2018/19, the annual FFB yields never fall below the mean confidence limit due to weak El Niño based on the historical 35 years records.

References

  1. Mielke, T. Global supply, demand and price outlook of oils and fats in 2018/19. In Proceedings of the GLOBOIL, Mumbai, India, 28 September 2018. [Google Scholar]
  2. Oettli, P.; Behera, S.K.; Yamagata, T. Climate based predictability of oil palm tree yield in Malaysia. Sci. Rep. 2018, 8, 2271. [Google Scholar] [CrossRef] [Green Version]
  3. Malaysian Palm Oil Industry. Available online: http://www.palmoilworld.org/about_malaysian-industry.html (accessed on 17 August 2021).
  4. World Energy Council. Biofuels: Policies, Standards and Technologies; World Energy Council: London, UK, 2010. [Google Scholar]
  5. Palm Oil Registration & Licensing Authority (PORLA). PORLA Palm Oil Statistics; PORLA: Kelana Jaya, Malaysia, 1986–1999. [Google Scholar]
  6. Malaysian Palm Oil Board (MPOB). Malaysian Oil Palm Statistics; MPOB: Bangi, Malaysia, 2000–2020. [Google Scholar]
  7. Economics and Industry Development Division. Available online: https://bepi.mpob.gov.my/index.php/en/ (accessed on 17 August 2021).
  8. Wicke, B.; Sikkema, R.; Dornburg, V.; Junginger, M.; Faaij, A. Drivers of Land Use Changes and the Role of Palm Oil Production in Indonesia and Malaysia; Universiteit Utrecht Copernicus Institute: Utrecht, The Netherlands, 2008. [Google Scholar]
  9. Ling, A.H. Global palm oil trade—Prospects and outlook. In Proceedings of the Malaysian-China Business Forum, Kuala Lumpur, Malaysia, 4 March 2019. [Google Scholar]
  10. Barcelos, E.; Rios, S.D.A.; Cunha, R.N.V.; Lopes, R.; Motoike, S.Y.; Babiychuk, E.; Skirycz, A.; Kushnir, S. Oil palm natural diversity and the potential for yield improvement. Front. Plant. Sci. 2015, 6, 190. [Google Scholar] [CrossRef] [PubMed]
  11. World Development Indicators. Available online: https://databank.worldbank.org/reports.aspx?source=2&country=MYS# (accessed on 25 August 2021).
  12. Covert 1 US Dollar to Malaysian Ringgit. Available online: https://www.xe.com/currencyconverter/convert/?Amount=1&From=USD&To=MYR (accessed on 25 August 2021).
  13. Overview of the Malaysian Oil Palm Industry 2020. Available online: https://bepi.mpob.gov.my/images/overview/Overview_of_Industry_2020.pdf (accessed on 28 August 2021).
  14. Which Countries Produce the Most Palm Oil? Available online: https://www.statista.com/chart/23097/amount-of-palm-oil-produced-in-selected-countries/ (accessed on 25 October 2021).
  15. What is El Niño-Southern Oscillation (ENSO)? Available online: https://www.weather.gov/mhx/ensowhat (accessed on 17 August 2021).
  16. What are El Niño and La Niña? Available online: https://oceanservice.noaa.gov/facts/ninonina.html (accessed on 25 October 2021).
  17. El Niño. Available online: https://www.nationalgeographic.org/encyclopedia/el-nino/ (accessed on 25 October 2021).
  18. Cold & Warm Episodes by Season. Available online: https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php (accessed on 25 October 2021).
  19. El Niño and La Niña Years and Intensities. Available online: https://ggweather.com/enso/oni.htm (accessed on 17 August 2021).
  20. Shanmuganathan, S.; Narayanan, A. Modelling the climate change effects on Malaysia’s oil palm yield. In Proceedings of the 2012 IEEE Symposium on e-Learning, e-Management and e-Services (IS3e 2012), Kuala Lumpur, Malaysia, 21–24 October 2012. [Google Scholar]
  21. Tawang, A.; Tengku Ahmad, T.A.; Abdullah, M.Y. Stabilization of Upland Agriculture under El Nino Induced Climatic Risk: Impact Assessment and Mitigation Measures in Malaysia; United Nation: New York, NY, USA, 2001. [Google Scholar]
  22. Kamil, N.N.; Omar, S.F. Climate variability and its impact on the palm oil industry. Oil Palm Ind. Econ. J. 2016, 16, 18–30. [Google Scholar]
  23. Wójtowicz, M.; Wójtowicz, A. The effect of climate change on linolenic fatty acid in oilseed rape. Agronomy 2020, 10, 2003. [Google Scholar] [CrossRef]
  24. Sobko, O.; Stahl, A.; Hahn, V.; Zikeli, S.; Claupein, W.; Gruber, S. Environmental effects on soybean (Glyciune Max (L.) Merr) production in Central and South Germany. Agronomy 2020, 10, 1847. [Google Scholar] [CrossRef]
  25. Kamil, N.N.; Omar, S.F. The impact of El Niño and La Niña on Malaysian palm oil industry. Oil Palm Bull. 2017, 74, 1–6. [Google Scholar]
  26. Harun, M.H.; Mohammad, A.T.; Noor, M.R.M.; Din, A.K.; Latiff, J.; Sani, A.R.A.; Abdullah, R. Impact of El Niño occurrence on oil palm yield in Malaysia. Planter 2010, 86, 837–852. [Google Scholar]
  27. Shanmuganathan, S.; Narayanan, A.; Mohamed, M.; Ibrahim, R.; Khalid, H. A hybrid approach to modelling the climate change effects on Malaysia’s oil palm yield at the regional scale. In Recent Advances on Soft Computing and Data Mining; Herawan, T., Ghazali, R., Deris, M., Eds.; Springer: Basel, Switzerland, 2014; pp. 335–345. [Google Scholar]
  28. Puah, P.W.; Sidik, M.J. Impacts of rainfall, temperature and recent El Niños on fisheries and agricultural products in the west coast of Sabah [2000–2010]. Borneo Sci. 2011, 28, 73–85. [Google Scholar]
  29. Rahman, A.K.A.; Abdullah, R.; Balu, N.; Shariff, F.M. The impact of La Niña and El Niño events on crude palm oil prices: An econometric analysis. Oil Palm Ind. Econ. J. 2013, 13, 38–51. [Google Scholar]
  30. Malaysia. Available online: https://www.nationsonline.org/oneworld/malaysia.htm (accessed on 25 October 2021).
  31. Land Area (sq. km)—Malaysia. Available online: https://data.worldbank.org/indicator/AG.LND.TOTL.K2?locations=MY (accessed on 25 October 2021).
  32. GeoRSPO—RSPO Mapbuilder App. Available online: https://rspo.org/members/georspo (accessed on 25 October 2021).
  33. Ling, L.; Yusop, Z.; Yap, W.S.; Tan, W.L.; Chow, M.F.; Ling, J.L. A calibrated, watershed-specific SCS-CN method: Application to Wangjiaqiao watershed in the three gorges area. China. Water 2019, 12, 60. [Google Scholar] [CrossRef] [Green Version]
  34. Conduct and Interpret a Spearman Rank Correlation. Available online: https://www.statisticssolutions.com/free-resources/directory-of-statistical-analyses/spearman-rank-correlation/ (accessed on 27 August 2021).
  35. Artusi, R.; Verderio, P.; Marubini, E. Bravais-Pearson and Spearman correlation coefficients: Meaning, test of hypothesis and confidence interval. Int. J. Biol. Mark. 2002, 17, 148–151. [Google Scholar] [CrossRef]
  36. Downloading IBM SPSS Statistics 26. Available online: https://www.ibm.com/support/pages/downloading-ibm-spss-statistics-26 (accessed on 27 August 2021).
  37. Difference between Adjusted R-squared and R-squared. Available online: https://www.listendata.com/2014/08/adjusted-r-squared.html (accessed on 28 August 2021).
  38. Adjusted R2/Adjusted R-Squared: What is it Used for? Available online: https://www.statisticshowto.com/probability-and-statistics/statistics-definitions/adjusted-r2/ (accessed on 27 August 2021).
  39. Ling, L.; Yusop, Z.; Chow, M.F. Urban flood depth estimate with a new calibrated curve number runoff prediction model. IEEE Access 2020, 8, 10915–10923. [Google Scholar] [CrossRef]
  40. Ling, L.; Yusop, Z.; Ling, J.L. Statistical and Type II Error Assessment of a Runoff Predictive Model in Peninsula Malaysia. Mathematics 2021, 9, 812. [Google Scholar] [CrossRef]
  41. Ali, R.; Kuriqi, A.; Abubaker, S.; Kisi, O. Long-term trends and seasonality detection of the observed flow in Yangtze River using Mann-Kendall and Sen’s innovative trend method. Water 2019, 11, 1855. [Google Scholar] [CrossRef] [Green Version]
  42. MAKESENS-Application for Trend Calculation. Available online: https://en.ilmatieteenlaitos.fi/makesens (accessed on 20 June 2021).
  43. Foreign Agricultural Service (FAS). Commodity Intelligence Report—Malaysia: Stagnating Palm Oil Tields Impede Growth; United States Department of Agriculture (USDA): Washington, DC, USA, 2012. [Google Scholar]
  44. Malaysia’s Palm Oil Yield to Continue Declining on Labour Shortage. Available online: https://www.theedgemarkets.com/article/malaysias-palm-oil-yield-continue-declining-labour-shortage (accessed on 28 August 2021).
  45. Palm Oil Analyst Raises 2019 Malaysia Output Forecast to 20m Tonnes. Available online: https://www.thestar.com.my/business/business-news/2019/04/30/palm-oil-analyst-raises-2019-malaysia-output-forecast-to-20mil-tonnes (accessed on 28 August 2021).
  46. IOI Group Oil Palm Plantations Malaysia 2020 by Age. Available online: https://www.statista.com/statistics/1097997/ioi-group-oil-palm-plantations-in-malaysia-by-age/ (accessed on 28 August 2021).
  47. Organisation for Economic Co-operation and Development (OECD)/Food and Agriculture Organization (FAO). OECD-FAO Agricultural Outlook 2019–2028; OECD/FAO: Paris, France, 2019. [Google Scholar]
  48. Banks and Palm Oil. Available online: https://www.banktrack.org/campaign/banks_and_palm_oil (accessed on 25 October 2021).
  49. Indonesia. Available online: https://climateknowledgeportal.worldbank.org/country/indonesia/climate-data-historical (accessed on 25 October 2021).
  50. Climate—Thailand. Available online: https://www.climatestotravel.com/climate/thailand (accessed on 25 October 2021).
  51. Malaysia Information. Available online: https://www.malaysia.gov.my/portal/content/144 (accessed on 25 October 2021).
  52. Palm Oil. Available online: https://tradingeconomics.com/commodity/palm-oil (accessed on 25 October 2021).
Figure 1. Market shares of 17 edible oils and fats produced internationally in: (a) 2018/19 and; (b) 1992/93. Created with data from [1].
Figure 1. Market shares of 17 edible oils and fats produced internationally in: (a) 2018/19 and; (b) 1992/93. Created with data from [1].
Agronomy 11 02189 g001
Figure 2. Location map of oil palm plantation areas in Malaysia. Created with data from [32].
Figure 2. Location map of oil palm plantation areas in Malaysia. Created with data from [32].
Agronomy 11 02189 g002
Figure 3. FV function in Microsoft Excel.
Figure 3. FV function in Microsoft Excel.
Agronomy 11 02189 g003
Figure 4. Oil palm yield loss and opportunity loss due to El Niño according to category.
Figure 4. Oil palm yield loss and opportunity loss due to El Niño according to category.
Agronomy 11 02189 g004
Figure 5. Monthly cumulative opportunity losses in the Malaysian oil palm industry due to El Niño projected to December 2021.
Figure 5. Monthly cumulative opportunity losses in the Malaysian oil palm industry due to El Niño projected to December 2021.
Agronomy 11 02189 g005
Table 1. Average productivity of various major oil crops. Created with data from [3,4].
Table 1. Average productivity of various major oil crops. Created with data from [3,4].
CropOil Yield (t/ha)
Palm Oil3.82
Palm kernel0.45
Palm + Palm kernel oil4.27
Rapeseed0.69
Sunflower0.52
Groundnut0.45
Coconut0.34
Cottonseed0.19
Soybean0.38
Corn0.15
Table 2. Categories of El Niño events occurred during the period of 1986 to 2021. Created with data from [19].
Table 2. Categories of El Niño events occurred during the period of 1986 to 2021. Created with data from [19].
Very StrongStrongModerateWeak
1997/981987/881986/872004/05
2015/161991/921994/952006/07
2002/032014/15
2009/102018/19
Table 3. Average FFB yield for each month from 1986 to 2020 in Malaysia. Created with data from [5,6,7].
Table 3. Average FFB yield for each month from 1986 to 2020 in Malaysia. Created with data from [5,6,7].
MonthAverage FFB Yield (t/ha)
January1.3194
February1.1903
March1.3234
April1.4049
May1.4829
June1.5109
July1.6329
August1.7420
September1.8543
October1.8146
November1.6746
December1.4743
FFB yields from January 2021 to July 2021 were not included in calculating these averages as data points in the full year of 2021 are still not available.
Table 4. Shapiro-Wilk normality test results for FFB, FFBI and ONI.
Table 4. Shapiro-Wilk normality test results for FFB, FFBI and ONI.
Variablep-ValueNormality
FFB0.0103 (<0.05)Non-normal
FFBI0.1922 (>0.05)Normal
ONI0.0000 (<0.05)Non-normal
Table 5. Spearman’s rho correlation between FFB and FFBI with ONI at lag periods from 0 to 18 months.
Table 5. Spearman’s rho correlation between FFB and FFBI with ONI at lag periods from 0 to 18 months.
Lag Period
(Months)
FFB & ONIFFBI & ONI
Spearman’s rho Coefficientp-ValueSpearman’s rho Coefficientp-Value
0−0.0270.578−0.0640.188
1−0.0290.557−0.0870.074
2−0.0300.540−0.1060.029 *
3−0.0340.479−0.1250.010 **
4−0.0490.317−0.1540.001 **
5−0.0790.104−0.1970.000 **
6−0.1200.013 *−0.2540.000 **
7−0.1650.001 **−0.3130.000 **
8−0.2000.000 **−0.3640.000 **
9−0.2170.000 **−0.3940.000 **
10−0.2110.000 **−0.3990.000 **
11−0.1880.000 **−0.3810.000 **
12−0.1590.001 **−0.3480.000 **
13−0.1280.008 **−0.3030.000 **
14−0.0950.050−0.2480.000 **
15−0.0600.214−0.1830.000 **
16−0.0300.542−0.1120.021 *
17−0.0070.878−0.0390.424
180.0090.8490.0230.642
* Correlation is significant at α = 0.05. ** Correlation is significant at α = 0.01. Lag period refers to the number of months FFB and FFBI are delayed corresponding to the ONI at that month.
Table 6. FFB and FFBI time series forecasting models’ summary with residual analyses using descriptive and inferential statistics at α = 0.01.
Table 6. FFB and FFBI time series forecasting models’ summary with residual analyses using descriptive and inferential statistics at α = 0.01.
Dependent VariableFFB (tons/Hectare)New FFBI (tons/Hectare)
PredictorONIONI
Modelled PeriodFeb 1986–Jun 2021 (N = 425)Feb 1986–Jun 2021 (N = 425)
Forecasted PeriodJul 2021–Dec 2023 (N = 30)Jul 2021–Dec 2023 (N = 30)
ModelARIMA (2,0,2) (1,1,1)ARIMA (3,1,3) (0,0,1)
Adjusted R-Squared0.82740.9312
Residual Sum of Squares (RSS)4.68760.5459
Residual Skewness−0.6903−0.4162
Residual Range0.76440.2790
Residual Median0.00770.0051
Residual Median (BCa 99% CI)[−0.0076, 0.0258][−0.0001, 0.0096]
Residual Standard Deviation0.10670.0359
Residual Standard Deviation (BCa 99% CI)[0.0941, 0.1192][0.0320, 0.0396]
Residual Variance0.01140.0013
Residual Variance (BCa 99% CI)[0.0088, 0.0142][0.0010, 0.0016]
Table 7. Mann Kendall’s trend test and Sen’s Slope Test during each El Niño event from January 1986 to June 2021.
Table 7. Mann Kendall’s trend test and Sen’s Slope Test during each El Niño event from January 1986 to June 2021.
CategoryEl Niño EventsMK Trend of FFBISen’s Slope
Very Strong1997/98Decreasing ***−0.0348
2015/16Decreasing ***−0.0438
Strong1987/88Decreasing **−0.0986
1991/92Decreasing *−0.0970
Moderate1986/87Decreasing **−0.0479
1994/95Decreasing *−0.0527
2002/03No Trend
2009/10Decreasing ***−0.0379
Weak2004/05Decreasing ***−0.0267
2006/07Decreasing *−0.0360
2014/15No Trend
2018/19Decreasing ***−0.0439
* Trend is significant at α = 0.05. ** Trend is significant at α = 0.01. *** Trend is significant at α = 0.001.
Table 8. Oil palm yield loss and opportunity loss due to each El Niño event from January 1986 to June 2021.
Table 8. Oil palm yield loss and opportunity loss due to each El Niño event from January 1986 to June 2021.
CategoryEl Niño EventsLagNegative FFBI PeriodOil Palm Yield LossOpportunity LossProjected to December 2021
(Months)(Months)(t/ha)(USD)(USD)
Very Strong1997/980163.0739799,512,628.043,118,809,718.44
2015/162253.51972,223,629,847.243,026,040,134.75
Strong1987/88481.270286,205,737.87621,682,504.89
1991/92390.979874,237,749.64422,833,751.88
Moderate1986/87181.227959,553,760.49451,416,000.07
1994/95340.137317,843,504.2983,902,729.92
2002/03360.284468,135,456.56206,588,187.64
2009/10490.7400484,337,375.58925,974,111.78
Weak2004/05670.301674,053,854.49189,770,606.98
2006/07350.203589,253,638.73208,552,751.05
2014/15050.3808197,547,898.49296,149,942.22
2018/194122.15761,308,607,860.031,481,072,202.50
USD = 4.2189178 as of 24 August 2021 [12]. Opportunity losses were projected using 6% discount rate. Opportunity loss for 2018/19 weak El Niño is excluded in the following analyses because the event has unexplained high production loss.
Table 9. Averaged oil palm yield loss and opportunity loss due El Niño according to category.
Table 9. Averaged oil palm yield loss and opportunity loss due El Niño according to category.
CategoryLagNegative FFBI PeriodOil Palm Yield LossOpportunity LossProjected to December 2021
(Months)(Months)(t/ha)(USD)(USD)
Very Strong1.020.53.29681,511,571,237.643,072,424,926.59
Strong3.58.51.125080,221,743.75522,258,128.39
Moderate2.86.80.5974157,467,524.23416,970,257.35
Weak3.05.70.2953120,285,130.57231,491,100.09
USD = 4.2189178 as of 24 August 2021 [12]. Opportunity losses were projected using 6% discount rate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khor, J.F.; Ling, L.; Yusop, Z.; Tan, W.L.; Ling, J.L.; Soo, E.Z.X. Impact of El Niño on Oil Palm Yield in Malaysia. Agronomy 2021, 11, 2189. https://doi.org/10.3390/agronomy11112189

AMA Style

Khor JF, Ling L, Yusop Z, Tan WL, Ling JL, Soo EZX. Impact of El Niño on Oil Palm Yield in Malaysia. Agronomy. 2021; 11(11):2189. https://doi.org/10.3390/agronomy11112189

Chicago/Turabian Style

Khor, Jen Feng, Lloyd Ling, Zulkifli Yusop, Wei Lun Tan, Joan Lucille Ling, and Eugene Zhen Xiang Soo. 2021. "Impact of El Niño on Oil Palm Yield in Malaysia" Agronomy 11, no. 11: 2189. https://doi.org/10.3390/agronomy11112189

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop