Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Growth and trend analysis of area, production and yield of rice: A scenario of rice security in Bangladesh

  • Md. Abdullah Al Mamun ,

    Contributed equally to this work with: Md. Abdullah Al Mamun, Sheikh Arafat Islam Nihad, Md. Abdur Rouf Sarkar

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    mamunru4777@gmail.com

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Sheikh Arafat Islam Nihad ,

    Contributed equally to this work with: Md. Abdullah Al Mamun, Sheikh Arafat Islam Nihad, Md. Abdur Rouf Sarkar

    Roles Conceptualization, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Plant Pathology Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Md. Abdur Rouf Sarkar ,

    Contributed equally to this work with: Md. Abdullah Al Mamun, Sheikh Arafat Islam Nihad, Md. Abdur Rouf Sarkar

    Roles Investigation, Resources, Supervision, Writing – review & editing

    Affiliation Agricultural Economics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Md. Abdullah Aziz ,

    Roles Data curation, Resources

    ‡ These authors also contributed equally to this work.

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Md. Abdul Qayum ,

    Roles Data curation, Resources

    ‡ These authors also contributed equally to this work.

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Rokib Ahmed ,

    Roles Data curation, Resources

    ‡ These authors also contributed equally to this work.

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Niaz Md Farhat Rahman ,

    Roles Data curation, Resources

    ‡ These authors also contributed equally to this work.

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Md. Ismail Hossain ,

    Roles Investigation, Methodology, Supervision, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Agricultural Statistics Division, Bangladesh Rice Research Institute, Gazipur, Bangladesh

  • Md. Shahjahan Kabir

    Roles Investigation, Methodology, Supervision, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Director General, Bangladesh Rice Research Institute, Gazipur, Bangladesh

Abstract

Bangladesh positioned as third rice producing country in the world. In Bangladesh, regional growth and trend in rice production determinants, disparities and similarities of rice production environments are highly desirable. In this study, the secondary time series data of area, production, and yield of rice from 1969–70 to 2019–20 were used to investigate the growth and trend by periodic, regional, seasonal and total basis. Quality checking, trend fitting, and classification analysis were performed by the Durbin-Watson test, Exponential growth model, Cochrane-Orcutt iteration method and clustering method. The production contribution to the national rice production of Boro rice is increasing at 0.97% per year, where Aus and Aman season production contribution significantly decreased by 0.48% and 0.49% per year. Among the regions, Mymensingh, Rangpur, Bogura, Jashore, Rajshahi, and Chattogram contributed the most i.e., 13.9%, 9.8%, 8.6%, 8.6%, 8.2%, and 8.0%, respectively. Nationally, the area of Aus and Aman had a decreasing trend with a -3.63% and -0.16% per year, respectively. But, in the recent period (Period III) increasing trend was observed in the most regions. The Boro cultivation area is increasing with a rate of 3.57% per year during 1984–85 to 2019–20. High yielding variety adoption rate has increased over the period and in recent years it has found 72% for Aus, 73.5% for Aman, and 98.4% for Boro season. As a result, the yield of the Aus, Aman, and Boro seasons has been found increasing growth for most of the regions. We have identified different cluster regions in different seasons, indicating high dissimilarities among the rice production regions in Bangladesh. The region-wise actionable plan should be taken to rapidly adopt new varieties, management technologies and extension activities in lower contributor regions to improve productivity. Cluster-wise, policy strategies should be implemented for top and less contributor regions to ensure rice security of Bangladesh.

Introduction

Globally rice (Oryza sativa) is the third major cereal grain [1] and more than half of the world population consumes its as a staple food [2]. Broad adaptive capability in different ecosystems and less cultivation risk, several farmers preferred rice cultivation instead of other crops. World population is increasing and it is assumed that 14,886 million tons (MT) of foods need to be produced in 2050 to meet up the food demand [3]. Worldwide 503.17 MT rice is produced where China produces 29.5% of the total, followed by India (23.8%), Bangladesh (7.0%), Indonesia (6.9%), Vietnam (5.4%), and Thailand (3.7%) [4].

Rice is also the staple food in Bangladesh and accounting for approximately 78 percent of the country’s total net cropped areas cultivation. The country achieves an autarky to meet up the rice demand for its 169.04 million peoples from 11.55 million hectares of cultivated gross area [5, 6]. In Bangladesh, food security is equivalent to rice security [5]. Rice is cultivated in three seasons namely Aus, Aman and Boro throughout the year. Since independence, rice production has been increased three-fold from approximately 11 MT in 1971–72 to about 36.6 MT in 2019–20 [7]. This revolution has transformed the country from so called “Bottomless Basket” to a “Full of Food Basket”. After a long period, rice production in Bangladesh has risen significantly after 1990–1991, especially during two periods: 1996–1997 and 2000–01, as well as from 2009–10 to 2013–14. Improved loan distribution policies (credit deposits directly to farmers’ 10 Taka bank accounts), well-organized fertilizer supplies, availability of high-quality seeds by the public and commercial sectors, and technical interventions (e.g. genetic improvements of varieties for favorable and unfavorable ecosystems) make it possible to make Bangladesh as one of the largest contributors of rice in the world [5, 8]. Bangladesh recently placed the third position worldwide in rice production, behind China and India, with a production volume of 3.6 crore tonnes [1].

In reality, the global food production has increased sharply since the Borlaug and Jennings days, keeping pace with an increasingly higher rate of population growth. To cite a country-specific example, we can easily refer to the Bangladesh scenario. Over the past four decades, Bangladesh succeeded in outpacing the population growth rate (1.3%) with its growth in rice output (2.8%) [7]. To increase or sustain the rice production, it is very important to interpose the rice-based technology in a specific ecosystems or specific locations. Knowledge about region specific rice cultivation scenario will be helpful to disseminate newly release technologies and to take necessary policy for sustainable rice production in Bangladesh. Therefore, this study investigates the growth and trend of area, production and yield of rice in Bangladesh based on periodic, regional, seasonal and total basis.

Materials and methods

Study area

Bangladesh is located in the northeastern part of South Asia. The majestic Himalayas stand some distance to the north, where the Bay of Bengal is in the south. These picturesque geographical boundaries frame a low-lying plain of about 1,47,570 square kilometers, criss-crossed by innumerable rivers and streams. The geographical extent of the country is between latitudes 20°34’ and 26°38’, and longitudes 88°01’ and 92°41’. The country’s borders are as follows: in the west, India (West Bengal); in the east, India (Tripura and Assam); in the south, Myanmar; and in the north, India (West Bengal and Meghalaya) [9]. Fig 1 depicts a conceptual framework of the study.

thumbnail
Fig 1. A conceptual framework of the study.

Shapefile republished from Bangladesh Agricultural Research Council (BARC) database (http://maps.barcapps.gov.bd/index.php) under a CC BY license, with permission from Computer and GIS unit, BARC, original copyright 2014.

https://doi.org/10.1371/journal.pone.0261128.g001

Data used

Seasonal and total rice area, production and yield data were used from 1969–72 to 2019–20 at the national level (aggregate) and 1984–85 to 2019–20 at the regional level (disaggregate).

This study is designed based on secondary sources data published by the different issues of the Bangladesh Bureau of Statistics (BBS). For the analysis, region-wise time-series data were divided into three periods, Period I (1984–85 to 1995–96); Period II (1996–97 to 2007–08), and Period III (2008–09 to 2019–20). The regional variations were analyzed by considering an unchanging regional base of homogeneous environments. Region-wise scenarios would provide a base to explain the effects of specific conditions as well as agricultural development. In the published issues of the Year Book of Agricultural Statistics, district-wise rice data is available after 2006–07 [10]. Earlier area, production, and yield of rice data were published according to 23 crop production regions representing the 30 agro-ecological zones (AEZ) of Bangladesh. In this study, we aggregated all data into 14 agricultural regions because the Department of Agricultural Extension (DAE) conduct their activities according to these regions representing 64 districts of Bangladesh. The distribution of the studied regions is presented in Fig 1.

Geographic Information System (GIS) map was used to describe the regional variations of area, production, yield and total rice in Bangladesh. However, administrative shape file of Bangladesh was downloaded and used from website of the Bangladesh Agricultural Research Council. The specific link of the shape file (“Administrative map”) is: http://maps.barcapps.gov.bd/index.php.

Methods

The autocorrelation detection, regression model for growth estimations, and cluster analysis techniques were used in this study. In time series regression data, there are numerous causes of auto-correlation. Autocorrelation is often caused by the analyst’s inability to incorporate one or more key predictor variables in the model. The existence of autocorrelation in the errors affects the ordinary least-squares regression method in many ways [11]. Although regression coefficients are still unbiased, they are no longer considered minimum-variance estimates. When errors are positively autocorrelated, the residual mean square may underestimate the error variance significantly. Autocorrelation may be detected using various statistical techniques. Durbin and Watson [12] developed a procedure that is frequently utilized worldwide. This test assumes that the errors in the regression model are produced by a first-order autoregressive process observed at evenly spaced time intervals [11]. A significant Durbin-Watson statistic or a suspicious residual plot suggests that autocorrelated model errors may be present. Table 1 represent the value of Durbin-Watson statistics and found that the data have positive autocorrelation. It may be due to a real-time dependency in error or an ’artificial’ time dependence induced by the absence of one or more key predictor variables.

thumbnail
Table 1. Durbin-Watson statistic for detection of autocorrelation.

https://doi.org/10.1371/journal.pone.0261128.t001

The observed autocorrelation in the model errors cannot be eliminated by adding additional predictor variables to the model. The autocorrelative structure in the model must be explicitly considered, and an appropriate parameter estimation technique must be used. The method developed by Cochrane and Orcutt [13] is a very excellent and frequently utilized approach. We used the linear model suggested by Finger [14] for analysis of the production contribution trend and expressed as: (1)

where,

Yt = Expected value at t

t = Time index

β0 = Model intercept

β1 = Regression coefficient

e = Residual

Also, exponential models were employed to assess the growth analysis. The exponential growth trend model and the equation of the model is as follows: (2)

Yt = Area, production, and yield of rice in year t

a = Model intercept

b = Annual rate of change of rice area, production, and yield

e = Residual

The cluster analysis is then used to find geographical groupings of rice-producing areas with comparable features in terms of production growth and adoption of high-yielding varieties. The variables used in the cluster analysis is independent. Thus, the cluster analysis uses the distinct principal components proposed by Huth and Pokorná [15]. Applying the most commonly used non-hierarchical clustering technique, K-means clustering [16, 17] that classify the 14 agricultural regions into K clusters using Euclidean distance as the linkage method. For all of the analyses, we used Microsoft Excel, ArcGIS 10.3, and R programming tools.

Results

Historical (1969–70 to 2019–20) rice area coverage and production contribution by seasons

The area coverage and production contribution of Aus, Aman and Boro rice have been shown in Fig 2. We found that historically the area coverage of Boro rice has increased significantly whereas the area of Aus and Aman season significantly decreased from 1969–70 to 2019–20. On the other hand, the production contribution of the Boro rice has increased from 16.1% to 53.7% and the increasing rate is 0.97% per year. Aus season production contribution has significantly decreased from 25.1% in 1969–70 to 7.5% in 2019–20 and the decreased rate was 0.48% per year. Production contribution of the Aman season has been significantly decreased 58.8% to 38.8% with a decreasing rate of 0.49% per year. Among the rice season, Aus production contribution is drastically reduced compared to the Aman season.

thumbnail
Fig 2.

National area coverage (a) and production contribution (b) of Aus, Aman and Boro rice during 1969–70 to 2019–20.

https://doi.org/10.1371/journal.pone.0261128.g002

Region-wise area and production of rice cultivation in Bangladesh

Thirty-six years (1984–85 to 2019–20) region-wise area and production of rice cultivation in Bangladesh is depicted in Fig 3. Periodic production contributions follow an increasing trend in Dinajpur, Mymensingh, Rajshahi and Rangamati where a periodic decreasing trend was observed in Chattogram and Cumilla regions (Fig 3A). In the recent period, the production contribution of Dhaka, Bogura, Jashore and Faridpur is decreased compared to the earlier period. Based on historical production trend, Mymensingh, Rangpur, Bogura, Jashore, Rajshahi and Chattogram significantly contributed 13.9%, 9.8%, 8.6%, 8.6%, 8.2% and 8.0%, respectively of the total production, and they were positioned as the most rice contributed regions compared to others. On the other hand, Khulna (4.4%) and Faridpur (3.7%) contributed less compared to other regions and Rangamati was the lowest production contributor (0.7%) among all regions. However, the rice cultivation area of Mymensingh has been found as higher compared to other regions where Rangamati had the lowest rice cultivation area in Bangladesh (Fig 3B).

thumbnail
Fig 3.

Periodic production contribution of major rice-growing regions (a), and regional production contribution and area coverage of rice (b) during 1984–85 to 2019–20. Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20.

https://doi.org/10.1371/journal.pone.0261128.g003

Regional trend and growth analysis of rice

Spatial and temporal (periodic) variability has been found in the rice cultivation area of Bangladesh (Fig 4). Over the period, regional variation prevails in terms of rice cultivation area (Fig 5 and Table 2). A continuous decreasing trend was observed in Aus cultivated area and an increasing trend was found in all regions for Boro season. The highest area decreasing rate (-18.43%) was found in Rangpur throughout the period in Aus season but in the recent period, the area is increasing at the rate of 29.62%. However, the Aman rice area fluctuated slightly in all regions, but Dhaka, Faridpur, Khulna and Sylhet regions showed a significantly decreasing trend (-0.43 to -1.54%). Overall Bangladesh, the Aus and Aman area decreased by -3.63% and -0.16% per year, respectively and the Boro area is increasing with an annual rate of 3.57% during 1984–85 to 2019–20. In Period I and Period II, cultivated areas were decreased in most of the regions for the Aus and Aman season, but in the recent period (Period III), i.e., an increasing trend was observed in most of the regions. In total, the highest growth rate of rice area was found in Rangamati (2.03%) and the lowest was observed in Faridpur (-1.25%).

thumbnail
Fig 4.

Spatial (regional) and temporal (periodic) distribution of Aus (a), Aman (b), Boro (c) and total (d) rice cultivation area of Bangladesh. Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20. Shapefile republished from Bangladesh Agricultural Research Council (BARC) database (http://maps.barcapps.gov.bd/index.php) under a CC BY license, with permission from Computer and GIS unit, BARC, original copyright 2014.

https://doi.org/10.1371/journal.pone.0261128.g004

thumbnail
Fig 5.

Trend of Aus (a), Aman (b), Boro (c) and total (d) rice cultivation area in Bangladesh.

https://doi.org/10.1371/journal.pone.0261128.g005

thumbnail
Table 2. Periodic growth rate (%) of Aus, Aman, Boro and total rice cultivation area at the regional level.

https://doi.org/10.1371/journal.pone.0261128.t002

Spatial and temporal (periodic) variability of rice production were observed in Bangladesh (Fig 6). Rice production trends of Aus, Aman, Boro and total are presented in Fig 7 and Table 3. The long-term production trend was found significantly increasing (2.18% to 10.25%) in all regions for the Boro season. In Aman season, production trends of all regions except Cumilla (-0.06%) were significantly increased (0.95% to 3.61%) over the period. Significant decreasing production trend (-15.33% to -4.05%) was found in most of the region for the Aus season but Barishal (2.39%), Rajshahi (3.49%) and Sylhet (1.07%) showed a significant increasing trend. In Aus, Aman and Boro season, the highest production growth rate was found in Rajshahi (3.49%), Rangamati (3.61%) and Dinajpur (10.25%), respectively. In total, the highest production growth rate was found in Dinajpur (4.07%) and the lowest was observed in Chattogram (1.64%).

thumbnail
Fig 6.

Spatial (regional) and temporal (periodic) distribution of Aus (a), Aman (b), Boro (c) and total (d) rice production (metric ton) of Bangladesh. Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20. Shapefile republished from Bangladesh Agricultural Research Council (BARC) database (http://maps.barcapps.gov.bd/index.php) under a CC BY license, with permission from Computer and GIS unit, BARC, original copyright 2014.

https://doi.org/10.1371/journal.pone.0261128.g006

thumbnail
Fig 7.

Trend of Aus (a), Aman (b), Boro (c) and total (d) rice production in Bangladesh.

https://doi.org/10.1371/journal.pone.0261128.g007

thumbnail
Table 3. Periodic growth rate (%) of Aus, Aman, Boro and total rice production at the regional level.

https://doi.org/10.1371/journal.pone.0261128.t003

Spatial and temporal (periodic) variability of rice yield were observed in Bangladesh (Fig 8). The yield trend was found increasing for all the regions and seasons of Bangladesh but region to region yield variations was very high (Fig 9 and Table 4). In most of the region, a significant increasing growth rate was found for Aus, Aman, Boro and total rice yield. In Aus, the highest annual growth rate of yield was found in Rajshahi (3.86%) and the lowest was in Rangamati (0.89%). In Aman, the highest growth rate of yield was observed in Faridpur (2.95%) and the lowest was in Rangamati (1.21%). In Boro, the highest growth rate of yield was found in Sylhet (2.78%) and the lowest was observed in Bogura (1.30%). In total, the highest growth rate of yield was found in Faridpur (4.09%) and the lowest was observed in Rangamati (1.38%).

thumbnail
Fig 8.

Spatial (regional) and temporal (periodic) distribution of Aus (a), Aman (b), Boro (c) and total (d) rice yield (metric ton) of Bangladesh. Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20. Shapefile republished from Bangladesh Agricultural Research Council (BARC) database (http://maps.barcapps.gov.bd/index.php) under a CC BY license, with permission from Computer and GIS unit, BARC, original copyright 2014.

https://doi.org/10.1371/journal.pone.0261128.g008

thumbnail
Fig 9.

Trend of Aus (a), Aman (b), Boro (c) and total (d) rice yield in Bangladesh.

https://doi.org/10.1371/journal.pone.0261128.g009

thumbnail
Table 4. Periodic growth rate (%) of Aus, Aman, Boro and total rice yield at the regional level.

https://doi.org/10.1371/journal.pone.0261128.t004

Modern variety adoption

Periodic modern varieties adoption (%) in different regions are illustrated in Fig 10. In Aus season, high yielding varieties (HYVs) adoption was very low in Periods I and II, but in recent years adoption is much higher than earlier and it is gradually increasing. For the Aman season, the HYVs adoption percentage is gradually increased in all regions, while the lowest adoption was found in Barishal and Faridpur regions. HYVs adoption in Boro season has found always been higher than Aus and Aman. Though adoption is higher in Boro season from the earlier periods, it is still gradually increasing and it reaches almost 100% in most of the areas. However, the adoption rate was comparatively low in Sylhet region during Boro season. Over the period, the modern variety adoption percentage gradually increasing in all over Bangladesh. The countrywide average modern variety adoption (%) in different seasons are shown in Fig 11. In recent period, the adoption percentage of Aus rice reached 72% from 21.4%. In Period I, the HYV adoption percentage for Aman season was 34.15% and it increased by 73.5% in Period III. Similarly, for Boro season adoption is gradually increasing and it reaches 88.6% to 98.4% from Period I to Period III.

thumbnail
Fig 10.

Periodic and regional high yielding varieties adoption (%) of Aus (a), Aman (b), Boro (c) and total (d) rice season in Bangladesh. Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20.

https://doi.org/10.1371/journal.pone.0261128.g010

thumbnail
Fig 11. Seasonal and total high yielding varieties adoption (%) in Bangladesh.

Period I: 1984–85 to 1995–96, Period II: 1996–97 to 2007–08, Period III: 2008–09 to 2019–20.

https://doi.org/10.1371/journal.pone.0261128.g011

Cluster analysis for identifying the dynamics of rice varietal adoption and production growth

The cluster analysis is used to identify groups of identical rice production regions depicting similar characteristics in their modern variety adoption and production growth rate (Fig 12). From the analysis, different classes of cluster have been identified to classify the regions exhibiting significant rising trends, significant decreasing trends and mixed or insignificant trends in the rice production and high yielding adoption growth rates. Fig 12 is prepared by using K-means clustering to reflect the aforesaid distinct characteristics. In Aus season, HYV variety adoption growth rate is positive for all the regions except Rangamati. But in Mymensingh, Bogura, Dinajpur, Rangpur, Dhaka and Faridpur have negative production growth rate and formed the same cluster. For Aman season, we have identified four clusters where Dinajpur, Rajshahi, Sylhet, Barishal formed a similar cluster. Another cluster contains Jashore, Mymensingh, Bogura, Rangpur and Cumilla regions and, Rangamati and Chattogram formed a similar cluster. In the case of Boro season, Chattogram, Rangamati, Bogura, Mymensingh, Faridpur, Dhaka and Mymensingh are in the same cluster environment whereas Rajshahi, Jashore, Rangpur and Barishal are in the same cluster. The other cluster was found in Khulna and Dinajpur regions, whereas Sylhet formed a single separate cluster (Fig 12).

thumbnail
Fig 12. Cluster analysis for grouping the regions based on high yielding variety (HYV) adoption (%) and production growth rate (%).

https://doi.org/10.1371/journal.pone.0261128.g012

Discussion

Growth and trend analysis are necessary to demonstrate the rice scenario by season, area, yield, production and varietal adaptability. Rice production depends on the season, variety, environment and geographical segmentation. Rice is cultivated in three seasons in Bangladesh namely Aus, Aman and Boro. This study revealed that the production contribution of Boro season is much higher compared to other seasons since 1998–99. Whereas production contribution of Aman season was higher from 1969–70 to 1997–98. Boro season is less vulnerable to rice disease, pests and other natural calamities and government incentives on Boro cultivation more specifically on irrigation facilities increase the Boro cultivation area [6, 18]. Moreover, after the release of two Boro mega varieties, BRRI dhan28 and BRRI dhan29 in 1994, the second silent green revolution of rice has occurred in Bangladesh. More than 60% of the area is covered by these two varieties during Boro season [19]. Currently, Boro cultivation area recorded for 61% of total cropped area in the Rabi season, which contributes 55% of total rice production in the country [18]. Our findings are also similar to this results. The maximum growing degree days, long vegetative growth, sunny weather and high amount of fertilizer utilization capacity favors the potential yield of rice during Boro season [18]. Moreover, the price of Boro rice is higher compared to other seasons which is also a key regulator of the intend of farmers to cultivate Boro rice [20]. Boro-Fallow-Fallow is the second dominant cropping pattern (13% of total net cropped area) of Bangladesh [6, 21] and this pattern profoundly found in Haor areas where other season rice cultivation is not possible except Boro due to stagnant water [22, 23]. Haor area is a major contributor (18%) of rice production in Boro season [24] and it might be another cause of yield differentiation from other seasons. However, cold stress at seedling stage in north-west and reproductive stage in north-east (haor) regions, heat wave during flowering stage, biotic stresses mainly blast and brown plant hopper, and flash flood are becoming major challenges for Boro rice cultivation [25]. To overcome these challenges, adoption of biotic and abiotic stress tolerant varieties, precision management and irrigation infrastructure are the key production risk management strategies.

Rainfed Aman rice yield has been found as static compared to other seasons during the study period. Therefore, to ensure the future food demand of overgrowing populations, we have to sustain or increase the Aman rice production in Bangladesh. Though the rice cultivation area of Aman is higher than Aus and Boro, the yield capacity of Aman varieties is lower than the Boro varieties [5]. Climatic conditions i.e., cloudy weather, low uptake capacity of fertilizer, lowest growing degree days and short life duration are the main causes of low rice yield of Aman season. Delay planting due to anomalous rainfall, early flood, drought and other climatic hazards causes low yield of Aman rice [26, 27] which discourage farmers to cultivate rice in Aman season. Late planting of Aman sometimes hampered the rabi crops (potato, wheat, mustard, vegetables, Boro rice etc.) cultivation. Rabi crops cultivation are more profitable and safer than Aman rice and so farmers prefer to keep the land fallow during Aman instead of rice cultivation [28]. To encourage the farmers for cultivating Aman rice, short duration HYVs i.e., BRRI dhan56, BRRI dhan57, BRRI dhan71 and BRRI dhan75 could be a potential technology. Coastal area comprises 20% of the country and it covers 30% of the net cultivation area of Bangladesh [29]. To withstand the challenge of salinity of coastal region during Aman season, farmers tend to cultivate traditional local saline tolerant rice varieties which gives poor yield [30]. Prawn culture is a profitable income source to the farmers of southern part and so they prefer to culture prawn in “Gher (closed flat area)” than the rice cultivation during Aman season [31]. Deficiency of nutrients such as N, P, Cu and Zn in saline soil are also a major drawback for low yield of rice in coastal areas. However, technological advancement and development of high yielding varieties can play a vital role for increasing Aman rice production. Farmers of the southern region tend to cultivate Aman instead of Boro due to availability of soil moisture and low salinity problem. Rainfall of Aman season diminishes the salinity which favors the rice cultivation in this season [32]. Moreover, salinity is a major problem for Boro (dry season) cultivation in the southern part of Bangladesh [33]. Heavy tide is another hinder for rice cultivation in southern-coastal regions. BRRI dhan76, BRRI dhan77 and BRRI dhan78 are modern high yielding rice varieties of Bangladesh Rice Research Institute (BRRI), bred in such a way so that their seedlings will be long to withstand the tidal wetland condition [34]. These varieties give around 5 tonnes yield per hectare which is higher than the local indigenous varieties (2.5 to 3 tonnes/ha) and these varieties shed a light on increasing rice production in the coastal areas. BR23 and BRRI dhan47 are high yielding salt tolerant varieties for Aman and Boro season, respectively which covered a significant area of southern regions. Now-a-days, the newly released HYV i.e., BRRI dhan67 is gaining popularity in saline prone areas of Bangladesh.

Aus is one of the most vulnerable rice growing seasons of Bangladesh [35]. Climatic conditions i.e., hot humid weather favors the outbreak of diseases and insect pest during Aus season. Tungro is one of the severe threats of Aus production and it can cause 100% yield loss of rice under severe outbreak condition [36]. Flash flood, drought, high temperature and low yield are also a major drawback of Aus production. Moreover, late transplanting of Boro rice is one of the reasons for Aus area reduction [37]. High yield and production as well as net return of Boro rice also dampened farmers’ interest to cultivate Aus rice. But intensive irrigated Boro rice cultivation depletes the underground water resulting irrigation water scarcity in the northern part of Bangladesh. So, shifting of Boro area to Aus cultivation is the key concern of the present time. However, adaptation of modern varieties (like BRRI dhan48), disease and pest management technologies and irrigation facilities could be a good option to withstand the challenge of the critical environment of the Aus season.

Mymensingh is the largest rice producing region while the lowest rice producing region was Rangamati in terms of area and production. This is because Mymensingh is favorable for rice cultivation and Rangamati is the hilly disadvantageous areas and less access to modern technologies. We found most of the regions have positive production growth rate for Aman and Boro season. But in Aus season, especially Mymensingh, Bogura, Dinajpur, Rangpur, Dhaka and Faridpur have negative production growth rate. This is due to most of the aforementioned regions have dominated Boro-Fallow-T. Aman cropping pattern. In addition, those regions have intensified with non-rice cash crops and thereby the required growth duration of Aus is insufficient. Noticeably, the released Aus rice varieties are not easily fit into this cropping pattern. Besides, some regions have single Boro dominated cropping patterns due to adverse agro-climatic and geographical conditions. This is why some regions have negative production growth rate. To increase the regional productivity equally, varietal and management interventions is must. Kabir et al., [5] reported that there are five unexplored areas in Bangladesh, where rice cultivation is possible. They mentioned that greater Barishal region, greater Sylhet region, South-west and greater Jashore region, Coastal charland in Barishal and Noakhali, and Chattogram hill tracts have unexplored areas. New HYVs, irrigation facilities, proper drainage system, and re-excavation of canals needs to be applied in these regions to increase the rice production.

Improved varieties and production technologies are the key drivers for enhancing rice production in Bangladesh. So far, we have developed 137 modern rice varieties and more than 300 production technologies by addressing different favorable and fragile ecosystems. Out of 36 stress tolerant varieties, BRRI dhan71 and BINA dhan19 in drought; BRRI dhan51, BRRI dhan52 and BINA dhan12 in submergence; BRRI dhan67, BRRI dhan97 and BINA dhan10 in saline, and BRRI dhan36 and BRRI dhan55 in cold ecosystems are the promising varieties for increasing production of stress prone areas in Bangladesh. Besides, BR11, BRRI dhan28, BRRI dhan29, BRRI dhan48, BRRI dhan49, BRRI dhan50, BRRI dhan58, BRRI dhan63, BRRI dhan81, BRRI dhan87, BRRI dhan89, BRRI dhan92, BRRI dhan96, BRRI dhan98, BRRI hybrid dhan3, BRRI hybrid dhan5 and BRRI hybrid dhan7 are the main players to boost up the rice production in Bangladesh. Use of water saving alternate wet and drying (AWD) techniques for irrigation, use of shallow water tubewell for irrigation from pond or river water, government subsidies for fertilizer and irrigation facilities, development of cropping pattern to relief the abject poverty of the northern part, disease and insect management technologies, and training of farmers to adopt modern rice cultivation techniques are acted as a catalyst to increase the rice production in Bangladesh.

Conclusion

This study examined the area, production and yield trend, and growth rates of rice from 1984–85 to 2019–20; and analytically classified the rice production regions. Trend analyses showed an increasing and decreasing growth rates for the Aus, Aman, Boro, and total rice determinants in different periods. In Aus rice, area and production growth rate had significantly decreased in all the regions, while yield was significantly increased over the period. In the Aman season, the area growth rate was decreased for seven regions, but production and yield growth rates were significantly increased for all regions. Based on area, production, and yield, Boro rice have found a significant increasing trend in all the regions. In the recent period, HYVs adaption rates were found 72% for Aus, 73.5% for Aman, and 98.4% for Boro season. During 1969–70 to 2019–20, the production contribution to the national rice production of Boro rice is significantly increasing at 0.97% per year, where Aus and Aman season production contribution significantly decreased by 0.48% and 0.49% per year, respectively. Aggregated in last 36 years rice production, Mymensingh (13.9%), Rangpur (9.8%), Bogura (8.6%), Jashore (8.6%), and Rajshahi (8.2%) were the top five production contributor regions in national rice production of Bangladesh. We have identified different cluster regions in different seasons, indicating high dissimilarities among the rice production regions. It is recommended that steps need to be taken to increase and sustain the rice production by implementing several specialized approaches. The region-wise actionable plan should be taken to rapidly adopt new technologies and highlight the research and extensions activities for fewer contributor regions to improve productivity. Cluster-wise policy strategies should be implemented for top and less contributor regions to ensure rice security as well as food security in Bangladesh.

Acknowledgments

Authors express their sincere thanks to the Bangladesh Bureau of Statistics (BBS) and Bangladesh Agricultural Research Council (BARC) for making available of relevant rice data and administrative GIS shape file of Bangladesh, respectively. The authors also acknowledge several scientists of the Bangladesh Rice Research Institute for participating discussion at various stages of preparing the manuscript.

References

  1. 1. Rahman M, Islam M, Rahaman M, Sarkar M, Ahmed R, Kabir M. Identifying the threshold level of flooding for rice production in Bangladesh: An Empirical Analysis. J Bangladesh Agric Univ. 2021.
  2. 2. Khush GS. What it will take to feed 5.0 billion rice consumers in 2030. Plant Molecular Biology. 2005. pmid:16217597
  3. 3. Mohammad Fakhrul Islam S, Karim Z. World’s demand for food and water: The consequences of climate change. Desalination -challenges and opportunities. 2020.
  4. 4. USDA (United State of Department of Agriculture). Rice production by country—world agricultural production 2020/2021. World agricultural production. 2020. Available: http://www.worldagriculturalproduction.com/crops/rice.aspx
  5. 5. Kabir MS, Salam M, Islam A, Sarkar MAR, Mamun M, Rahman M, et al. Doubling rice productivity in Bangladesh: A way to achieving SDG 2 and moving forward. Bangladesh Rice J. 2020;24: 1–47.
  6. 6. Nasim M, Khatun A, Kabir M, Mostafizur A, Mamun M, Sarkar M, et al. Intensification of cropping through utilization offallow period and unutilized land resources in Bangladesh. Bangladesh Rice J. 2021;25: 89–100.
  7. 7. BBS (Bangladesh Bureau of Statistics). Yearbook of agricultural statistics of Bangladesh. Government of Bangladesh, Dhaka. 2020.
  8. 8. Rabbi S, Biswas P, Rashid E, Iftekharuddaula K, Rahman N, Rahman M, et al. Increasing rice yield through targeting genetic potentials by rice types. Bangladesh Rice J. 2020;24: 67–82.
  9. 9. Bangladesh—an-introduction—ministry of fForeign affairs-. [cited 22 Sep 2021]. Available: https://mofa.gov.bd/site/page/6dde350b-1ca6-4c69-becd-a3f12cf14ac1/Bangladesh—An-Introduction
  10. 10. BBS (Bangladesh Bureau of Statistics). Yearbook of agricultural statistics of Bangladesh. Government of Bangladesh, Dhaka. 2007.
  11. 11. Douglas C. M, Elizabeth A. P, Geoffrey V G. Introduction to linear regression analysis, Fifth Edition. A JOHN WILEY & SONS, INC., PUBLICATION. 2012.
  12. 12. Durbin J, Watson GS. Testing for serial correlation in least squares regression. III. Biometrika. 1971.
  13. 13. Cochrane D, Orcutt GH. Application of least squares regression to relationships containing auto-correlated error terms. J Am Stat Assoc. 1949.
  14. 14. Finger R. Evidence of slowing yield growth—The example of swiss cereal yields. Food Policy. 2010. pmid:25411518
  15. 15. Huth R, Pokorná L. Simultaneous analysis of climatic trends in multiple variables: An example of application of multivariate statistical methods. Int J Climatol. 2005.
  16. 16. Wilks DS. Statistical methods in the atmospheric sciences, second edition. Meteorological Applications. 2007.
  17. 17. Kannan S, Ghosh S. Prediction of daily rainfall state in a river basin using statistical downscaling from GCM output. Stoch Environ Res Risk Assess. 2011.
  18. 18. Mainuddin M, Alam MM, Maniruzzaman M, Kabir MJ, Mojid MA, Hasan MM, et al. Yield, profitability, and prospects of irrigated Boro rice cultivation in the North-West region of Bangladesh. PLoS One. 2021. pmid:33914816
  19. 19. Salam MU, Mahalder BK, Bhandari H, Kabir MS, Abdur Rouf Sarkar M, Nessa B, et al. Policy directions toward increasing rice productivity—lessons from Bangladesh. Adv Rice Res Abiotic Stress Toler. 2019; 895–913.
  20. 20. Chanda SC, Ali MA, Haque ME, Abdullah MR, Sarwar AKMG. Cost of production and cost benefit analysis of different rice in Sirajganj district. Asian J Crop Soil Sci Plant Nutr. 2019.
  21. 21. Nasim M, Shahidullah SM, Saha A, Muttaleb MA, Aditya TL, Ali MA, et al. Distribution of crops and cropping patterns in Bangladesh. 2017;21: 1–55.
  22. 22. Alam MS, Quayum M, Islam M. Crop production in the haor areas of Bangladesh: Insights from Farm Level Survey. Agric. 1970.
  23. 23. CEGIS. Master plan of haor area, Volume II: Main Report. 2012;II: 73–80.
  24. 24. Kabir MJ, Kabir MS, Salam MA, Islam MA, Omar MI, Sarkar MAR, et al. Harvesting of Boro paddy in haor areas of Bangladesh: Interplay of local and migrant labour, mechanized harvesters and COVID-19 vigilance in 2020. 2020. Available: http://brri.portal.gov.bd/sites/default/files/files/brri.portal.gov.bd/page/6af533b2_278e_4e30_9e27_0626f3448845/2020-06-29-15-48-e0c6a8f2f5d11ff1ee4b79762198e087.pdf
  25. 25. Rahman M, Haque M, Kabir MJ, Islam A, Sarkar M, Mamun M, et al. Enhancing rice Productivity in the unfavourable ecosystems of Bangladesh. Bangladesh Rice J. 2020;24: 83–102.
  26. 26. Mir Kabir MM, Naher UA, Panhwar QA, Shamshuddin J, Khan FH. Effect of transplanting dates on growth and yield of inbred and hybrid rice varieties grown during rainfed season in Bangladesh. Philipp Agric Sci. 2014.
  27. 27. Biswas A, Ahmed MME, Halder T, Akter S, Yasmeen R, Rahman MS. Photosensitive rice (Oryza sativa L.) varieties under delayed planting as an option to minimize rice yield loss in flood affected T. Aman season. Bangladesh Rice J. 2020.
  28. 28. Sarker BC, Mainuddin M, Scientific TC, Bell RW, Bengal W. Effect of sowing dates on yield of wheat grown in excess water and salt affected soils in Southwestern coastal Bangladesh J. of the Indian Societ. of Coast. Agri. Res. 2019;37: 51–59.
  29. 29. Haque SA. Salinity problems and crop production in coastal regions of Bangladesh. Pak. J. Bot. 2006.
  30. 30. Shelley IJ, Takahashi-Nosaka M, Kano-Nakata M, Haque MS, Inukai Y. Rice cultivation in Bangladesh: present scenario, problems, and prospects. J Intl Cooper Agric Dev. 2016.
  31. 31. Fatema K, Miah TH. Rice versus shrimp farming in Khulna district of Bangladesh: interpretations of field-level data. Bangladesh J Agric Econ. 2011.
  32. 32. Mondal MK, Bhuiyan SI, Franco DT. Soil salinity reduction and prediction of salt dynamics in the coastal ricelands of Bangladesh. Agric Water Manag. 2001.
  33. 33. Alam MZ, Carpenter-Boggs L, Mitra S, Haque MM, Halsey J, Rokonuzzaman M, et al. Effect of salinity intrusion on food crops, livestock, and fish species at kalapara coastal belt in Bangladesh. J Food Qual. 2017.
  34. 34. Barua R, Islam MA, Mamin MSI, Karmakar B, Biswash MR, Islam MR, et al. Performance of some tidal submergence tolerant rice genotypes in southern coastal regions of Bangladesh. Sci Agric. 2017;17: 112–117.
  35. 35. Kabir M, Sarkar M, Rahman M, Rahman N, Mamun M, Chowdhury A, et al. Risk of rice cultivation under current and future environment and market. Bangladesh Rice J. 2021.
  36. 36. Nihad SAI, Manidas AC, Hasan K, Hasan MA-I, Honey O, Latif MA. Genetic variability, heritability, genetic advance and phylogenetic relationship between rice tungro virus resistant and susceptible genotypes revealed by morphological traits and SSR markers. Curr Plant Biol. 2021;25: 100194.
  37. 37. Uddin MT, Dhar AR. Government input support on Aus rice production in Bangladesh: Impact on farmers’ food security and poverty situation. Agric Food Secur. 2018. pmid:32983428