Next Article in Journal
Human-Based New Approach Methodologies in Developmental Toxicity Testing: A Step Ahead from the State of the Art with a Feto–Placental Organ-on-Chip Platform
Previous Article in Journal
Effectiveness of an E-Book App on the Knowledge, Attitudes and Confidence of Nurses to Prevent and Care for Pressure Injury
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Climate Change Impacts Assessment Using Crop Simulation Model Intercomparison Approach in Northern Indo-Gangetic Basin of Bangladesh

by
Md Rafique Ahasan Chawdhery
1,*,
Murtuza Al-Mueed
2,
Md Abdul Wazed
3,
Shah-Al Emran
4,
Md Abeed Hossain Chowdhury
5 and
Sk Ghulam Hussain
6
1
Department of Agroecology and Crop Production, Faculty of Agrobiology, Food and Natural Resources, Czech University of Life Sciences Prague, Kamýcká 129, Suchdol, 16500 Prague, Czech Republic
2
Ministry of Public Administration, Dhaka 1000, Bangladesh
3
Department of Chemical and Materials Engineering, University of Auckland, Private Bag 92019, Auckland 1142, New Zealand
4
Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
5
Bangladesh Agricultural Research Council, Crop Zoning Project, Farmgate, Dhaka 1000, Bangladesh
6
International Maize and Wheat Improvement Center (CIIMMYT), Gulshan-1, Dhaka 1212, Bangladesh
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(23), 15829; https://doi.org/10.3390/ijerph192315829
Submission received: 7 October 2022 / Revised: 18 November 2022 / Accepted: 23 November 2022 / Published: 28 November 2022

Abstract

:
The climate change impacts of South Asia (SA) are inextricably linked with increased monsoon variability and a clearly deteriorating trend with more frequent deficit monsoons. One of the most climate-vulnerable nations in the eastern and central Indo-Gangetic Basin is Bangladesh. There have been numerous studies on the effects of climate change in Bangladesh; however, most of them tended to just look at a small fraction of the impact elements or were climatic projections without accounting for the effects on agriculture. Additionally, simulation studies using the CERES-Rice and CERES-Wheat models were conducted for rice and wheat to evaluate the effects of climate change on Bangladeshi agriculture. However, up to now, Bangladesh has not implemented farming system ideas by integrating cropping systems with other income-generating activities. This study was conducted as part of the Indo-Gangetic Basin (IGB) regional evaluations using the protocols and integrated assessment processes of the Agricultural Model Intercomparison and Improvement Project (AgMIP). It was also done to calibrate crop models (APSIM and DSSAT) using rice and wheat. To assist policymakers in creating national and regional plans for anticipated future agricultural systems, our work on the integrated evaluation of climate change impacts on agricultural systems produced realistic predictions. The outcome of this research prescribes a holistic assessment of climate change on future production systems by including all the relevant enterprises in the agriculture sector. The findings of the study suggested two major strategies to minimize the yield and increase the profitability in a rice–wheat cropping system. Using a short-term HYV (High Yielding Variety) of rice can shift the sowing time of wheat by 7 days in advance compared to the traditional sowing days of mid-November. In addition, increasing the irrigation amount by 50 mm for wheat showed a better yield by 1.5–32.2% in different scenarios. These climate change adaptation measures could increase the per capita income by as high as 3.6% on the farm level.

1. Introduction

Climate change impacts are increasingly visible in South Asia (SA), with greater variability of the monsoon, noticeably, a declining trend with more frequent deficit monsoons [1]. Bangladesh is one of the most climate-vulnerable countries in the central and eastern Indo-Gangetic Basin. The rapidly growing population of the country puts tremendous pressure on its scarce natural resources. The country is vulnerable to many climatic hazards, including frequent floods, droughts, cyclones, and storm surges that damage life, property, and agricultural production [2]. Climate change could make it more difficult than it already is to accelerate agricultural production to meet the growing demands in countries in the region and more specifically for Bangladesh. The above-mentioned types of disasters make the problems more complicated. In the foreseeable future, Bangladesh is likely to be one of the most vulnerable countries in the world in the event of climate change [3]. Global warming due to the increase in greenhouse gas concentrations in the earth’s atmosphere and the consequent sea level rise is adding fuel to the fire. Almost every sector of socio-economic life in Bangladesh is likely to be affected by climate change. The 6th Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) has recently reported that “Warming of the climate system is unequivocal, and since the 1950s, many of the observed changes are unprecedented over decades to millennia. The atmosphere and ocean have warmed, the amounts of snow and ice have diminished, sea level has risen, and the concentrations of greenhouse gases have increased” [4]. Climatic factors (e.g., temperature, precipitation, wind, etc.) regulate agricultural activities and security in production [5,6]. In developing countries, such as Bangladesh, poor households mostly live in rural areas and are largely dependent on the agricultural production system, and therefore, climate change has a major impact on their livelihood and food security [7,8].
There have been several studies of climate change impacts in Bangladesh, but they have provided climate projections without quantifying agricultural impacts or tended to examine just a subset of the impact factors. For example, some studies of climate change impacts on Bangladesh rice have focused on temperature and carbon dioxide effects [9,10]. Additionally, simulation studies were carried out for rice and wheat using the APSIM-Oryza and APSIM-Wheat models and CERES-Rice and CERES-Wheat for DSSAT to assess the impact of climate change on Bangladesh agriculture [11,12]. The Agricultural Production Systems Simulator (APSIM) is a software package that facilitates sub-modules combined to simulate agricultural systems [13]. The APSIM simulator has several modules categorized as Plant-Soil-Atmosphere and crop management operations. The model can simulate crop growth stage developments, soil processes, and crop management options. The APSIM model requires an input dataset of weather, soil, and crop management development on a daily basis to run the simulation outputs. The DSSAT cropping model is a computer-based crop growth application that can simulate growth, development and yield of crop growing under defined management over time [14,15]. The DSSAT model requires databases of daily basis weather data, soil, crop management conditions and cultivar-specific genotypes to simulate and compare results with observed estimates. Moreover, a number of investigations have been made in the Indo-Gangetic area to achieve sustainable land use. For instance, a conceptual modeling framework could be adopted to discover dynamic land consequences and understand their responses for sustaining ecosystem services in the dominant agricultural region [16] and land transition information for sustainable management and development of land resources [17]. It has also been reported that the input data uncertainty has a higher influence on land use output and hence high-resolution geospatial inputs were recommended for crop and land management in this dominant agricultural region [18,19]. In most cases, the detrimental effect of temperature rise was observed even with elevated CO2 levels. Considerable spatial and temporal variations were also noted. Resilience development through climate change adaptation strategies is also a strong part of the challenges for crop production in Bangladesh [20]. Wheat was more susceptible to high temperature than rice and temperature increase generally reduces production across all scenarios [21]. Precipitation changes can have either a positive or a negative impact, with a high degree of uncertainty across Global Circulation Models (GCMs) [22]. Carbon dioxide impacts on crop production are positive and depend on the emissions pathway. Precipitation uncertainties from different GCMs and emissions scenarios are reduced when integrated across the large Ganges-Brahmaputra-Meghna Basins’ hydrology [23]. Irrigation requirements are likely to increase with regional variations in agricultural production [24].
Agriculture in southern Bangladesh is severely affected by sea level rise even when cyclonic surges are not fully considered, with impacts increasing under the higher emissions scenario. The impact of climate change on cereals and food security has been studied by several authors; however, no effort in using farming system approaches integrating cropping systems with other income-generating activities has been made for Bangladesh until now. This study was undertaken as a component of the IGB (Indo-Gangetic Basin) regional assessments following AgMIP protocols and integrated assessment procedures. It is expected that this integrated assessment of climate change impacts on agricultural systems would generate reasonable estimates to help policymakers to develop national and regional plans for projected future agricultural systems.

2. Materials and Methods

2.1. Description of Study Area

The study area is located in the Dinajpur district (25°22’–26°06′ N, and 89°31′–88°38′ E), which is part of the Eastern Gangetic Plain region in Bangladesh (Figure 1). The total area of the Dinajpur district is 3438 km2 of which 78.87 km2 is forested [25]. Dinajpur is one of the most intensive agricultural areas among the crop production hubs in Bangladesh. The total cultivable cropped land was around 2707 km2 in 2014–2015 [26]. The general weather type of this region is humid and wet with subtropical climate variations identified as summer, monsoon, and winter seasons. The yearly mean highest temperature of the study location was around 33.5 °C and the lowest one was nearly 10.5 °C. The annual total rainfall was commonly 2540 mm. The dominant soil types of this region are divided into three categories, e.g., Piedmont plan, Tista Floodplain, and Barind Tract [27]. Based on flooding during the monsoon, the lands of Dinajpur can be divided into three distinct land types: High Land-HL (e.g., land higher than the normal flooding point), Medium High Land-MHL (e.g., land flooded up to 90 cm for a minimum of two weeks), and Medium Low Land-MLL (e.g., the land flooded between 90 to 180 cm for over two weeks) [27]. The texture of the topsoil (usually 0–15 cm) is primarily loam and silt loam, although it can also be sandy loam, silty clay loam, or clay loam [28].

2.2. Farming System Diagram

Rice (Oryza sativa) is the most important commodity in terms of livelihood and food in the study area, followed by wheat (Triticum aestivum L.). In Bangladesh, rice is grown in three distinct seasons: Aus (Kharif-I), Aman (Kharif-II), and Boro (Rabi), while wheat is grown only during the winter (Rabi). The most dominant farming system of the study area was rice–wheat covering about 67% of the net sown area. The rice and wheat acreage in Dinajpur were 0.42 million hectare (Mha) and 0.019 Mha, respectively, with almost the same level of average productivity between rice and wheat (3.2 and 3.1 tha−1, respectively). Though livestock is an integral part of the farming system, this study focused on crops enterprise only. The cultivar used for rice was BR11 (T. Aman) cultivar which is grown during the wet season (Kharif-II). BR11 is the most popular variety for T. Aman season which is grown from 15th to 30th July. Usually, 25–35-day-old seedlings are transplanted with 20 cm row spacing. The crop’s lifecycle is 145 days. For wheat, BARI-Gom21 (Shatabdi) cultivar was used, which is grown during the winter season. The optimum sowing date for wheat is 15th November to 7th December with a lifecycle of 109–112 days. The yield potentials of BR11 and Shatabdi is 5.5 tha−1 and 3.5–5.0 tha−1, respectively [29]. For integrated assessment, a field survey was carried out during May–September 2019 to collect data from 50 randomly selected farms from four Upazilas (sub-districts) of the Dinajpur district to cover the major cropping pattern, e.g., the rice (Transplanted Aman)–wheat cropping pattern. All data were collected through a questionnaire survey directly from fifty sampled farms. The collected data were then cleaned, edited, standardized, and organized as per the prescribed format in spreadsheets to fulfil the requirements for running the crop and socio-economic models (Section 2.5). The farming system diagram of the study location is presented in Figure 2.

2.3. Representative Agricultural Pathways (RAPs) and Adaptation Package

The Representative Agricultural Pathways (RAPs) is an overall narrative description of a plausible future development pathway, and also contains key variables with qualitative storylines and quantitative trends, consistent with higher-level pathways (e.g., SSPs, global RAPs developed by the AgMIP Global Modeling Group). RAPs are translated into one or more scenarios (parameterizations) for The Tradeoff Analysis Model for Multi-dimensional Impact Assessment (TOA-MD) and crop models. These scenarios represent a set of technology and management adaptations to climate change. These scenarios, developed for specific RAPs, typically include changes in the types of crops or livestock produced and the way they are managed (e.g., use of fertilizers and improved crop cultivars). Procedures for RAP development are based on the stepwise process as described earlier [30].
The future trends in the agricultural system, and the RAP narratives were developed by the project team and were shared with the stakeholders. Several interactions were made with different stakeholders that included academicians from agricultural universities and scientists from the national agricultural research system to prepare the RAPs. Moreover, during the farm survey, interactions were held with 50 farmers from four sub-districts of the Dinajpur district, and their opinions about current climatic variability and how they perceive climate change and their anticipated impact on agriculture were recorded.
A combination of population growth, government subsidy on fertilizers, and improved economic performance is expected to cause a shift from agriculture to the service industry. From the historic data, it is observed that farm and family size is decreasing over time. The stakeholders agreed that farm size would decrease by 20 percent and family size would decrease by 10 percent during the next 30-year period [31]. The RAPS parameters and adaptation strategy used in the TOA-MD analysis are given below (Table 1).
To minimize the adverse impact of climate change on wheat production, several adaptation options were proposed by the participants in the stakeholders’ workshop. However, the following adaptation strategies were finalized: Modification of irrigation dates and amount of water for wheat. Advancement of sowing date by one week to avoid terminal heat stress for wheat. Improved fertilizer management and soil test-based application in both rice and wheat. Use of short-duration rice varieties so that advanced wheat planting can be accommodated. Out of the suggested adaptation measures, increasing the irrigation amount by 50 mm per irrigation for wheat was tested in this integrated assessment.

2.4. Data Collection and Methods of Study

2.4.1. Soil Data Preparation

The soil profile-wise data were gathered from Reconnaissance Soil Survey Reports (Soil Resource Development Institute, Bangladesh). A total number of 45 soil series profiles were created as input for DSSAT models (e.g., AGMIP13001 to AGMIP13045). Examples from the Amnura soil profile are presented in Supplementary Tables S1 and S2. Four soil profiles, e.g., Amnura (AGMIP13028), Jamun (AGMIP13035), Ekdala (AGMIP13032), and Gangachara (AGMIP13033), were used for the simulation runs for Dinajpur.

2.4.2. Observed Trends in Temperature and Precipitation

Historical time series (1990–2020) daily agro-meteorological data (daily maximum and minimum temperature, daily precipitation, and daily sunshine hour) were collected from the Bangladesh Meteorological Department (BMD, Supplementary Table S3). The sunshine hour data were converted to solar radiation (MJ m−2 d−1) as per the input format of the crop models as described by Allen et al. [32]. The input weather dataset was structured as per the model requirement. These were used as a historical climate dataset. The trends in temperature and precipitation observed data are presented in Figure 3.

2.4.3. Climate Projections

Using the AgMIP Climate Scenario Generation Tools with R, total 20 Global Circulation Model (GCM) data were generated for RCP 8.5 (2040–2069, Mid-Century) based on the observed (OXXX) climate (30 years weather data (1990–2020-baseline)) data. Based on GCM climate change scenarios, five more scenarios, namely IEXA, IIXA, IKXA, IOXA, IRXA, for the median future scenarios were also generated for model runs. Here, the first letter I stands for RCP 8.5 (2040–2069, Mid-Century); E = CCSM4; I = GFDL-ESM2M; K = HadGEM2-ES; O = MIROC5; R = MPI-ESM-MR; type of scenario- X = Observations (no scenario); A = the mean change from GCM. In the case of RCP8.5 Mid-Century 2040–2069 scenarios, CO2 concentration was considered at 571 ppm [33]. The winter period (December, January, and February) showed a similar pattern regardless of GCMs. In this study, all the GCMs forecast higher summer temperature. Moreover, a reduced amount of rainfall was predicted during the winter season and higher during the monsoon period (mid-June to mid-October). Predicted maximum and minimum temperatures and rainfall in Dinajpur are presented in Figure 4.
From Figure 5, it is clear that all 20 Mid-Century scenarios projected higher maximum and minimum temperatures in February and March compared to the baseline time period (blue X). Again, March is warmer than February. The highest February maximum temperature (29.69 °C) is predicted by IPSL-CM5A-LR (M), and the highest March maximum temperature (34.58 °C) is predicted by CCSM4 (E); the rest of the GCMs are predicted above 30 °C.

2.5. Socio-Economic Analysis

Socio-economic data of 50 representative farm households were collected from Dinajpur in 2019. The average household size was 5.74 persons per house against the national average of 5.31 persons per house [34]. The average farm size was 0.90 ha. The average annual non-agricultural income was Bangladeshi taka (BDT) 59,600 farm−1, and the contribution of crop component was BDT 40,690 farm−1 year−1. The average plot size of rice was 0.65 ha. The majority of the farmers cultivated rice in medium high land (58%) followed by high land (26%) and medium low land (16%). Farmers applied 150.15 kg nitrogen, 69.14 kg TSP, 51.37 kg MoP (Muriate of potash) per hectare, which is higher than the recommended doses of fertilizers for T. Aman rice (90–120 kg N, 80–100 kg TSP and 80–120 kg MoP, and 50–72 kgha−1 gypsum (FRG, 2012)). The crop was grown in rainfed conditions. The farmers obtained 3826 kgha−1 grain yield which was a little higher than the five-years mean yield of rice (3801 kgha−1) in Dinajpur. The average plot size of wheat was 0.40 hectares. Most of the farmers cultivated wheat in medium high lands (58%) followed by high land (28%) and medium low land (16%). The input use per hectare by the farmers was 2139 kg farmyard manure (FYM), 86.5 kg N, 84.3 kg TSP, 78.2 kg MoP, 20.2 kg gypsum, 1.83 kg zinc sulphate, and 9.09 kg borax. The recommended fertilizer doses for wheat were 81–99 kg N, 140–180 kg TSP, 40–45 kg MoP, and 110–120 kg gypsum per hectare (BARI, 2019). The farmers obtained 3346 kgha−1 grain yield, which was much higher than the five-year mean yield of wheat (2497 kgha−1) in Dinajpur.

2.6. Crop-Model Calibration (APSIM and DSSAT)

In this study, two crops (rice and wheat) were modeled. The crop models used are CERES-Rice and CERES-Wheat for DSSAT Ver. 4.5.1.023 –Stub [35] and APSIM 7.5-Oryza and APSIM-Wheat [36]. The Agricultural Production Systems Simulator (APSIM) is a software package that facilitates sub-modules combined to simulate agricultural systems [13]. The APSIM simulator has several modules categorized as Plant-Soil-Atmosphere and crop management operations. The model can simulate the crop growth stage developments, soil processes, and crop management options. The APSIM model requires the input dataset of weather, soil, and crop management development on a daily basis to run the simulation outputs. The DSSAT cropping model is a computer-based crop growth application that can simulate growth, development, and yield of crop growing under defined managements over time [14,15]. The DSSAT model requires databases of daily basis weather data, soil, crop management conditions, and cultivar specific genotypes to simulate and compare results with observed estimates. Estimated genetic coefficients for the selected cultivars are presented in Table 2.
The calibration outputs for rice (APSIM and DSSAT) and wheat (APSIM and DSSAT) are presented in Table 3 and Table 4, respectively. Supplementary Figures S1 and S2 depict the calibration and validation outputs (observed and simulated) using the APSIM model for the rice (BR11-T. Aman) and wheat (Shatabdi) cultivars, respectively.

2.7. Validation of Crop Simulation Models

The crop model calibration is the process of comparison between the observed and simulated values that confirmed the acceptable estimates of the model outputs. Crop models require cultivar specific genetic coefficients to simulate performance of diverse genotypes under different soil, weather and management conditions (crop growth, development, and grain production). The calibration outputs and the validation of different crop simulation models need to operate in different ways. Two statistical criteria were used to validate the model: the normalized root mean square error (NRMSE) and the root mean square error (RMSE) [37].
RMSE = i = 1 n ( S i O i ) 2 n
NRMSE = 100 1 n i = 1 n ( S i O i ) 2 Omax Omin
where s are the simulated values, o are the observed values, and n are the number of data points. Omax and Omin are the maximum and minimum values of the observations, respectively.

3. Results

3.1. Evaluation of Crop Models (APSIM and DSSAT)

Intercomparison was made between two rice models (APSIM-Rice and DSSAT-CERES-Rice). Compared to the observed yields, higher yields were predicted by both models for the Dinajpur location. The ambiguity in yields related to different farms was greater with APSIM model predictions compared to the DSSAT model. However, both the models overvalued the simulated crop yields. Higher crop yields were predicted by the DSSAT-Wheat model in comparison with the APSIM-Wheat model. The APSIM-wheat model predicted less than 3500 kgha−1 at 55% cumulative probability, which was lower than the observed yield. The uncertainty of cumulative probability was higher with APSIM in simulated crop yields among different farms compared to DSSAT (Figure 6). Higher simulated yields compared to farm survey yields are expected since neither models considered pest and disease effects.
To evaluate the impact of future climate scenarios with APSIM and DSSAT rice and wheat models, the yield performance of the 30-year baseline weather dataset was used to simulate 30 years of future climate scenarios. There was a total of six sets of climate data, namely (OXXX) for the current climate, and IEXA, IIXA, IKXA, IOXA, and IRXA future scenarios were predicted for CO2 at 571 ppm level.
The input files for the simulation were generated based on the following parameters.
Rice (rainfed): The simulation runs were made for Dinajpur location using DSSAT CERES-Rice and APSIM-Rice models. The rice cultivar was BR11 (T. Aman), which was sown from 31 May to 10 August of 2018. The cultivation depth was 5 cm. Final crop density was 25 plants m−2 and three plantshill−1. Nitrogen was applied as basal at sowing @27–40.5 kgha−1 with a placement depth of 10 cm. Most of the farmers applied nitrogen in two splits @39–46 kgha−1 and 39–46 kgha−1 as top-dressing at 15 and 45 days after sowing, respectively.
Wheat (irrigated): The simulations were carried out for the Dinajpur location with DSSAT CERES-Wheat and APSIM-Wheat models. The wheat cultivar was BARI Gom21 (Shatabdi). The sowing period of the wheat cultivar was from 11 November to 07 December of 2018 with a cultivation depth of 5 cm and crop density of 240 plantsm−2. Nitrogen was applied as basal at sowing @30–40 kgha−1 with a placement depth of 10 cm. Most of the farmers applied nitrogen in two splits @25–35 kgha−1 and 25–40 kgha−1 as top-dress at 25 and 55 days after sowing, respectively; and three irrigations (flood method) at 20, 45 and 75 days after sowing to maintain the soil moisture at field capacity level.
The simulated yields are presented with box and whisker plots. The box represents data between the 25th and 75th percentiles. The center black horizontal line across the box is the median at 50% of yields, and the whiskers (error bars) above and below the box indicate the 95th and 5th percentiles.
For the single treatment, the DSSAT CERES-Rice and APSIM-Rice models were simulated with 30 seasons of baseline weather (base) and five different climate scenarios (IEXA, IIXA, IKXA, IOXA, and IRXA). APSIM and DSSAT rice models predicted lower yields for all the GCM scenarios compared to the baseline, but the magnitude of yield decreases was small in the case of DSSAT, and a higher yield variability was noted in the case of APSIM (Figure 7). Only higher yields were predicted by the IOXA scenario. The lowest yields were predicted by the IKXA scenario for both the models (Figure 7). On the other hand, APSIM predicted lower yields for all the GCM scenarios compared to the baseline, but the magnitude of yield decreases was quite high. The variability observed in the yields in different farms may be due to variable fertilizer and water management caused by variable availability of rainwater and temperature increase.
The maximum rice yield values predicted by APSIM for OXXX (Historical), IEXA, IIXA, IKXA, IOXA, and IRXA scenarios were 8557, 7408, 7267, 5880, 7538, and 5492 kgha−1, while the minimum values were 6607, 5071, 4920, 3697, 5586, and 3718 kgha−1 with average values of 5904, 7472, 6003, 4749, 6533, and 4829 kgha−1, respectively (Figure 7). In the case of DSSAT, the predicted maximum values were 7832, 7053, 7005, 6584, 7095, and 6483 kgha−1 for OXXX (Historical), IEXA, IIXA, IKXA, IOXA, and IRXA scenarios, respectively. However, the minimum yields were 3756, 3703, 3640, 3652, 3744, and 4112 kgha−1 with average yield values of 7185, 6630, 6503, 6107, 6710, and 6152 kgha−1 (Figure 7).
Similarly, the APSIM-Wheat and the DSSAT CERES-Wheat models were simulated with 30 seasons of baseline weather (base) and five different climate scenarios (IEXA, IIXA, IKXA, IOXA, IRXA). For both DSSAT and APSIM models, predicted yields were higher for all climate scenarios (Figure 8). In the case of the DSSAT model, future climate scenarios IEXA, IKXA, and IRXA have more adverse climatic effects on wheat yield compared to the baseline scenario. In contrast, the future climate scenarios showed positive effects on wheat yield for the APSIM model (Figure 8).
The APSIM model projected maximum wheat yield for OXXX, IEXA, IIXA, IKXA, IOXA, and IRXA scenarios were 5692, 5624, 6223, 5879, 6308, and 6054 kgha−1, and the minimum were 2878, 2825, 3547, 3175, 3621, and 3824 kgha−1 with average yields of 3464, 3398, 4286, 3780, 4369, and 4471 kgha−1, separately (Figure 8).
At the same time, the DSSAT model anticipated that the maximum wheat yields for OXXX, IEXA, IIXA, IKXA, IOXA, and IRXA were 5507, 5299, 5369, 5204, 5548, and a minimum of 5009 kgha−1, and 4036, 3470, 3452, 3380, 3649, and 3338 kgha−1 with average projected yields were 4637, 4481, 4691, 4485, 4885, and 4426 kgha−1 (Figure 8).
To offset the impact of climate change and reduce the variability in wheat yields, an adaptation strategy was incorporated in this study. Although there were several adaptation strategies suggested in the RAP development workshop, only one, e.g., increasing the irrigation amount by 50 mm per irrigation for wheat, was used. Accordingly, both APSIM and DSSAT wheat models were used to generate outputs under five different CGMs.
The APSIM model showed reduced variability in yields with adaptation strategy irrespective of GCMs. In the case of the baseline (OXXX) scenario, the yields were 2878 to 5692 kgha−1 without adaptation and 2998 to 6434 kgha−1 with adaptation. Among all scenarios, the IRXA scenario showed the maximum average yield in both without and with the adaptation strategy (Figure 9).
Correspondingly, the DSSAT model indicated that projected median yields were raised with the adaptation in all GCMs and reduced variability in yields. In the case of DSSAT, the predicted minimum yields for the historic baseline (OXXX) scenario without adaptation was 4036 kgha−1 and with adaptation was 4112 kgha−1 (Figure 9). The maximum yield without adaptation was 5507 kgha−1 and with adaptation was 5514 kgha−1, respectively. The IRXA scenario showed the lowest median yields without and with adaptation of DSAAT compared to APSIM. The simulated minimum and maximum yield predictions for the IRXA scenario differed between 3338 and 5009 kgha−1 without adaptation, and between 3338 and 5034 kgha−1 with adaptation (Figure 10).

3.2. Impact of Climate Change on Crop Production

The sensitivity of the current production system under five different climatic scenarios is presented in Table 5. The current production system, under all five climate change scenarios, shows that the mean yield of rice would decline by 13.6–38.2 percent (with APSIM) and 7.4–14.7 percent (with DSSAT), respectively (Table 5). On the other hand, mean wheat yield is likely to increase (0.8–21.8%) under all five climate scenarios (with APSIM). However, DSSAT estimates show mixed results—a decline in wheat yield (2.9–8.8%) under four climate scenarios (IEXA, IIXA, IKXA and IRXA), and a marginal increase (1%) under climate scenario IOXA.
The current production system is vulnerable to climate change as it shows decline (4.6–17.8%) in mean net farm returns (APSIM model) under all but one climate scenario (IOXA) which shows marginal increase (0.74%) in net farm returns. The extent of losses in net farm returns varies from 4–11.6% (DSSAT model). Consequently, the per capita income under climate change declines (0.8–3.9%) in all the scenarios except one scenario (IOXA for APSIM model) which shows a marginal increase (0.16%). Though changes in the poverty levels are minimal, there is a large proportion of the population (48–80%) which is adversely affected by climate change (Table 5).

3.3. Impact of Climate Change on Future Agricultural Systems

Using RAP parameters and country-specific estimates of productivity and price trends for rice and wheat from global models, TOA-MD analysis was done to study the impact of climate change on future agricultural production systems. In the case of rice, larger yield decline (14–38%) was predicted by the APSIM model than that by the DSSAT model (7–15%). Interestingly, wheat yield was likely to increase (0.8–22%) under all climate scenarios (APSIM model); however, DSSAT estimates show decline in wheat yields (3–9%), except under one climate scenario (IOXA) which shows a marginal increase (1%) in wheat yields (Table 6). Interestingly, the mean returns per farm were projected to increase under all but one climate scenario (IKXA, APSIM model) in future in the range of 2–29% and 0.5–26% for APSIM and DSSAT, respectively (Table 6). It appears that the decline in rice yields is more compensated by the increase in wheat yields (APSIM) and higher growth trend in yield and prices of rice and wheat along with a doubling of off-farm income in future agricultural production systems. The cumulative effect is strong enough to offset the decline in productivity (rice) and acreage under the rice–wheat system. Though increase in per capita income is quite low (0.07–4.8%), except IKXA- APSIM model (−0.78%), it substantially reduces the proportion of the population which is adversely affected by climate change. Total gains, losses, and net impact on farm population are shown in Figure 11A,B. Though it appears a small change in per capita income, it brings about 37–39% of the population out of poverty (Table 6).

3.4. Benefits of Climate Change Adaptations

Though parameters derived from the RAPs scenario and the growth trend in price and productivity of the rice–wheat system resulted in overall increase in net farm returns, yet more than half of the population remains below the poverty line. So, it would be imperative to assess the impact of adaptation measures on net farm income and livelihoods. In this study, the adaptation strategy was tested for wheat crop only, wherein irrigation amount was increased by 50 mm per irrigation. The adaptation package for wheat results in higher increased (10–32%) in wheat yields with APSIM model in comparison to DSSAT model (1–3% only) under five climate scenarios (Table 7). Consequently, the increase in mean net farm returns was larger (5–19%) with APSIM than that for DSSAT (7–10%). Overall, the per capita income increased up to 3.6% only, and the poverty rate declined marginally (Table 7). It is obvious because wheat acreage is less than that of rice, which is the most important crop and affects the livelihoods of most of the population.

4. Discussion

Crop simulation model evaluation is a fundamental activity to compare the observed and simulated growth and development of crops performed by different models [24]. APSIM and DSSAT crop models are widely used for agronomic crop simulation performance for development and yields predictions [30]. Based on the input dataset and the climatic stresses and uncertainty, both the APSIM and DSSAT model can perform different simulation results for rice and wheat from the same location [32]. In the case of wheat, APSIM showed higher sensitivity and lower yield prediction than the observed farm yields and DSSAT estimated higher yield in our study which also studied in different parts of China [33]. Moreover, the ambiguity in yield prediction with different sites was more in the case of APSIM compared to DSSAT in Dinajpur. The variations in yield predictions between the APSIM and DSSAT models could be due to the exclusion of pest and disease infestation and damage [31]. In the current study, the sensitivity of both the APSIM and DSSAT crop models estimated the lower yield predictions by 13.6–38.2% (APSIM) and 7.4–14.7% (DSSAT) for rice under different climate scenarios (Table 7). The current crop production system in most parts of the Indo-Gangetic belt showed a decreasing rice yield trend for both APSIM and DSSAT models [36]. Moreover, increasing CO2 concentration (571 ppm) will influence the total global average temperature, which will affect rice and wheat yields [34]. Climate change has a significant impact on agricultural land use diversity, cropping area, and production based on the shifting pattern of rainfall and its fluctuation, temperature, flooding risk, agroecological features of different regions, humidity, and sunshine hours [38,39]. In response to climate change, farmers may utilize more inputs, which could have an impact on future agricultural productivity and production efficiency as well as on human health and the environment [40]. Other studies noted that due to the use of pesticides and chemical fertilizers, water extraction for irrigation, crop intensification, and CO2 emissions can enhance the uncertainty in sustainable agricultural production (rice–wheat cropping system) in Bangladesh [41]. This declined rice production will affect the net farm returns under climate change scenarios [35]. In contrast, mean wheat yield is likely to increase by 0.8–21.8 per cent (APSIM) and mixed results–decline in wheat yield by 2.9–8.8 per cent (DSSAT). A recent study also suggested that both APSIM and DSSAT yield variabilities for wheat can be managed by strategic management of fertilizer applications in some cases [36]. Socio-economic analysis for the current study estimated the impact of climate change on future agricultural production systems on farm level. In the case of rice, yield decline was predicted by APSIM and DSSAT model in all future climate scenarios. Remarkably, wheat yield is estimated to rise under all climate scenarios by the APSIM model, while DSSAT predicted decline in wheat yields. The cumulative effect is enough to balance the decline in productivity (rice) and land use under the rice–wheat system. The impact of future climate change events on the rice–wheat cropping system has mixed magnitude on rural farm level production [42]. Using crop models can help farmers make decisions about how to manage their crops to reduce climatic risk and make the most use of limited resources [43]. Crop yield predictions, using future climate scenarios and various crop models from the same site where all other necessary parameters for cropping are known, can bring the chance to compare the selection of crop model for predictable crop yield based on the impacts of climate change [14,44,45,46]. In this study, the adaptation strategy was analyzed for wheat crop with irrigation amount increased by 50 mm per irrigation which resulted 10–32 percent higher yield with APSIM compared to the same amount applied to project with DSSAT. The planting time, irrigation schedule, and amount of irrigation along with climate change also influenced the growth and development of the crops. Improved fertilizer method, sowing density, and climate resilient cultivar can also be considered as major adaptation strategies for rice and wheat cropping systems [24,47]. To combat the diverse impact of climate change, maintaining the irrigation water use efficiency and introducing drought resistant cultivars plays a significant role [48]. The current study revealed that the benefits of the climate change adaptation strategy can increase the per capita income by up to 3.6% on the farm level and thus could potentially reduce the poverty rate. It is to be noted that adaptation is crucial and responses to climate change need to be well-integrated with poverty reduction strategies to achieve sustainable development. Moreover, climate change adaptation strategies exert direct and immediate impacts on the poor by affecting factors that condition poverty reduction, such as economic growth.

5. Conclusions

Since rice is the most important crop for the population, there is a need to test the adaptation strategy for rice for a more realistic integrated assessment. This study confirmed the development and potential application of an integrated approach of climate change impact assessment. From the climatic perspective, the findings of the study suggested that the rice–wheat cropping system can be adjusted by some adaptation strategies. The duration of rice cropping can be minimized by using the short-term high yielding rice varieties to adjust the wheat sowing timing 7 days earlier than the traditional mid-November to 1st week of November. Shifting the wheat sowing date will minimize the monsoon threat on the wheat ripening stage, which will ensure the quality of grain and the storage life. Additionally, increasing the irrigation amount by 50 mm in wheat cultivation increased the yield. However, this study could not include other enterprises (livestock, fishery, etc.) in the assessment. Therefore, a holistic assessment of climate change on future production systems would only be possible by including all the relevant enterprises of the agriculture sector and testing the most effective adaptation strategy. Lastly, a methodological framework should be developed for upscaling the integrated assessment for the entire Indo-Gangetic Basin.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph192315829/s1.

Author Contributions

Conceptualization, M.R.A.C., M.A.H.C., and S.G.H.; formal analysis, M.R.A.C. and M.A.H.C.; funding acquisition, M.R.A.C., M.A.H.C. and S.G.H.; investigation, M.R.A.C., M.A.-M. and M.A.W.; project administration, M.R.A.C., M.A.H.C. and S.G.H.; resources, M.R.A.C., M.A.H.C. and S.G.H.; supervision, S.-A.E., M.A.H.C. and S.G.H.; validation, M.R.A.C., M.A.H.C. and S.G.H.; visualization, M.R.A.C., M.A.H.C. and S.G.H.; writing—original draft, M.R.A.C., M.A.-M. and M.A.W.; writing—review and editing, M.R.A.C., S.-A.E., M.A.H.C. and S.G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding, however internally supported by research project of Technology Agency of the Czech Republic: ‘Water systems and water management in the Czech Republic in conditions of climate change’ TAČR “Prostředí pro život” (SS02030027). Moreover, the APC was kindly funded by MDPI.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

Authors would like to acknowledge AGMIP-Bangladesh and Global for technical support and data. The useful comments of Roberto Valdivia and contributions from Harbir Singh and N Subash for reviewing this work are gratefully acknowledged. We also thank Vera Potopová for her enormous support while preparing the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pal, L.; Ojha, C.S.P.; Dimri, A.P. Characterizing Rainfall Occurrence in India: Natural Variability and Recent Trends. J. Hydrol. 2021, 603, 126979. [Google Scholar] [CrossRef]
  2. Islam, F.; Alam, G.M.; Begum, R.; Sarker, M.N.I.; Bhandari, H. Vulnerability, Food Security and Adaptation to Climate Change of Coastal Rice Farmers in Bangladesh. In Climate Vulnerability and Resilience in the Global South; Springer: Cham, Switzerland, 2021; pp. 187–197. [Google Scholar]
  3. Filho, W.L.; Modesto, F.; Nagy, G.J.; Saroar, M.; YannickToamukum, N.; Ha’apio, M. Fostering Coastal Resilience to Climate Change Vulnerability in Bangladesh, Brazil, Cameroon and Uruguay: A Cross-Country Comparison. Mitig. Adapt. Strateg. Glob. Chang. 2018, 23, 579–602. [Google Scholar] [CrossRef]
  4. Intergovernmental Panel on Climate Change (IPCC). Climate Change 2022: Impacts, Adaptation and Vulnerability; Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climaste Change; Pörtner, H.O., Roberts, D.C., Tignor, M.M.B., Poloczanska, E.S., Mintenbeck, K., Alegría, A., Craig, M., Langsdorf, S., Löschke, S., Möller, V., et al., Eds.; Cambridge University Press: Cambridge, UK, 2022. [Google Scholar]
  5. Cattiaux, J.; Yiou, P.; Vautard, R. Dynamics of Future Seasonal Temperature Trends and Extremes in Europe: A Multi-Model Analysis from CMIP3. Clim. Dyn. 2012, 38, 1949–1964. [Google Scholar] [CrossRef]
  6. Berry, P.M.; Rounsevell, M.D.A.; Harrison, P.A.; Audsley, E. Assessing the Vulnerability of Agricultural Land Use and Species to Climate Change and the Role of Policy in Facilitating Adaptation. Environ. Sci. Policy 2006, 9, 189–204. [Google Scholar] [CrossRef]
  7. Arunrat, N.; Wang, C.; Pumijumnong, N.; Sereenonchai, S.; Cai, W. Farmers’ Intention and Decision to Adapt to Climate Change: A Case Study in the Yom and Nan Basins, Phichit Province of Thailand. J. Clean. Prod. 2017, 143, 672–685. [Google Scholar] [CrossRef]
  8. Masud, M.M.; Azam, M.N.; Mohiuddin, M.; Banna, H.; Akhtar, R.; Alam, A.F.; Begum, H. Adaptation Barriers and Strategies towards Climate Change: Challenges in the Agricultural Sector. J. Clean. Prod. 2017, 156, 698–706. [Google Scholar] [CrossRef]
  9. Basak, J.; Ali, M.; Islam, M.; Alam, J. Assessment of the Effect of Climate Change on Boro Rice Production in Bangladesh Using DSSAT Model. INT Cent. 2009, 38, 9. [Google Scholar]
  10. Timsina, J.E. Humphreys Applications of CERES-Rice and CERES-Wheat in Research, Policy and Climate Change Studies in Asia: A Review. Int. J. Agric. Res. 2006, 1, 202–225. [Google Scholar] [CrossRef] [Green Version]
  11. Hussain, S. Assessing Impacts of Climate Change on Cereal Production and Food Security in Bangladesh. In Climate Change and Food Security in South Asia; Springer: Dordrecht, The Netherlands, 2011; pp. 459–476. ISBN 978-90-481-9515-2. [Google Scholar]
  12. Karim, Z.; Hussain, S.G.; Ahmed, M. Assessing Impacts of Climatic Variations on Foodgrain Production in Bangladesh. In Climate Change Vulnerability and Adaptation in Asia and the Pacific; Springer: Dordrecht, The Netherlands, 1996. [Google Scholar]
  13. McCown, R.; Hammer, G.; Hargreaves, J.; Holzworth, D.; Freebairn, D. APSIM: A novel software system for model development, model testing and simulation in agricultural systems research. Agric. Syst. 1996, 50, 255–271. [Google Scholar] [CrossRef]
  14. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  15. Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Shelia, V.; Wilkens, P.W.; Singh, U.; White, J.W.; Asseng, S.; Lizaso, J.I.; Moreno, L.P.; et al. The DSSAT crop modeling ecosystem. In Advances in Crop Modeling for a Sustainable Agriculture; Burleigh Dodds Science Publishing: Cambridge, UK, 2019; pp. 173–216. [Google Scholar] [CrossRef]
  16. Palmate, S.S.; Pandey, A.; Pandey, R.P.; Mishra, S.K. Assessing the land degradation and greening response to changes in hydro-climatic variables using a conceptual framework: A case-study in central India. Land Degrad. Dev. 2021, 32, 4132–4148. [Google Scholar] [CrossRef]
  17. Pandey, A.; Palmate, S.S. Assessments of spatial land cover dynamic hotspots employing MODIS time-series datasets in the Ken River Basin of Central India. Arab. J. Geosci. 2018, 11, 1–8. [Google Scholar] [CrossRef]
  18. Palmate, S.S.; Wagner, P.D.; Fohrer, N.; Pandey, A. Assessment of Uncertainties in Modelling Land Use Change with an Integrated Cellular Automata–Markov Chain Model. Environ. Model. Assess. 2022, 27, 275–293. [Google Scholar] [CrossRef]
  19. Swain, S.; Mishra, S.K.; Pandey, A.; Pandey, A.C.; Jain, A.; Chauhan, S.K.; Badoni, A.K. Hydrological modelling through SWAT over a Himalayan catchment using high-resolution geospatial inputs. Environ. Chall. 2022, 8, 100579. [Google Scholar] [CrossRef]
  20. Emran, S.A.; Krupnik, T.J.; Aravindakshan, S.; Kumar, V.; Pittelkow, C.M. Factors Contributing to Farm-Level Productivity and Household Income Generation in Coastal Bangladesh’s Rice-Based Farming Systems. PLoS ONE 2021, 16, e0256694. [Google Scholar] [CrossRef] [PubMed]
  21. Ruane, A.C.; Cecil, L.D.; Horton, R.M.; Gordón, R.; McCollum, R.; Brown, D.; Killough, B.; Goldberg, R.; Greeley, A.P.; Rosenzweig, C. Climate Change Impact Uncertainties for Maize in Panama: Farm Information, Climate Projections, and Yield Sensitivities. Agric. For. Meteorol. 2013, 170, 132–145. [Google Scholar] [CrossRef]
  22. Tarek, M.; Brissette, F.; Arsenault, R. Uncertainty of Gridded Precipitation and Temperature Reference Datasets in Climate Change Impact Studies. Hydrol. Earth Syst. Sci. 2021, 25, 3331–3350. [Google Scholar] [CrossRef]
  23. Janes, T.; McGrath, F.; Macadam, I.; Jones, R. High-Resolution Climate Projections for South Asia to Inform Climate Impacts and Adaptation Studies in the Ganges-Brahmaputra-Meghna and Mahanadi Deltas. Sci. Total Environ. 2019, 650, 1499–1520. [Google Scholar] [CrossRef] [PubMed]
  24. Potopová, V.; Trnka, M.; Vizina, A.; Semerádová, D.; Balek, J.; Chawdhery, M.R.A.; Musiolková, M.; Pavlík, P.; Možný, M.; Štěpánek, P.; et al. Projection of 21st Century Irrigation Water Requirements for Sensitive Agricultural Crop Commodities across the Czech Republic. Agric. Water Manag. 2022, 262, 107337. [Google Scholar] [CrossRef]
  25. BBS (Bangladesh Bureau of Statistics). Zila Statistics 2011; Statistics and Informatics Division, Ministry of Planning, Government of the People’s Republic of Bangladesh: Dhaka, Bangladesh, 2012.
  26. BBS (Bangladesh Bureau of Statistics). Yearbook of Agricultural Statistics of Bangladesh 2015, 35th ed.; Statistics and Informatics Division, Ministry of Planning, Government of the People’s Republic of Bangladesh: Dhaka, Bangladesh, 2015; ISBN 978-984-475-047-0.
  27. Huq, S.I.; Shoaib, J.M. The Soils of Bangladesh; Springer: Dordrecht, The Netherlands, 2013; Volume 1, ISBN 978-94-007-1128-0. [Google Scholar]
  28. Hassan, A.A.; Jahiruddin, M.; Shamsun, N.; Jalal, U.S.; Abdul, L.S.; Monowar, K.K.; Bokhtiar, S.M.; Quddus, A.; Nazmul, H.; Sultana, R.; et al. Fertilizer Recommendation Guide-2012; Bangladesh Agricultural Research Council: Dhaka, Bangladesh, 2012; ISBN 978-984-500-000-0. [Google Scholar]
  29. Chowdhury, M.A.H.; Hassan, M.S. Hand Book of Agricultural Technology; Bangladesh Agricultural Research Council: Farmgate, Dhaka, 2013; p. 230. ISBN 978-984-500-018-5. [Google Scholar]
  30. Antle, J.M.; Stoorvogel, J.J.; Valdivia, R.O. New parsimonious simulation methods and tools to assess future food and environmental security of farm populations. Philos. Trans. R. Soc. B Biol. Sci. 2014, 369, 20120280. [Google Scholar] [CrossRef] [PubMed]
  31. BBS (Bangladesh Bureau of Statistics). Yearbook of Agricultural Statistics of Bangladesh 2020, 40th ed.; Statistics and Informatics Division, Ministry of Planning, Government of the People’s Republic of Bangladesh: Dhaka, Bangladesh, 2020; ISBN 978-984-475-047-0.
  32. Allen, R.; Pereira, L.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper; Food and Agriculture Organization of the United Nations: Rome, Italy, 1998; Volume 56. [Google Scholar]
  33. Ruane, A.C.; Major, D.C.; Yu, W.H.; Alam, M.; Hussain, S.G.; Khan, A.S.; Hassan, A.; Al Hossain, B.M.T.; Goldberg, R.; Horton, R.M.; et al. Multi-Factor Impact Analysis of Agricultural Production in Bangladesh with Climate Change. Glob. Environ. Chang. 2013, 23, 338–350. [Google Scholar] [CrossRef] [Green Version]
  34. Holzworth, D.P.; Huth, N.I.; Devoil, P.G.; Zurcher, E.J.; Herrmann, N.I.; Mclean, G.; Chenu, K.; Van Oosterom, E.J.; Snow, V.; Murphy, C.; et al. APSIM—Evolution towards a New Generation of Agricultural Systems Simulation. Environ. Model. Softw. 2014, 62, 327–350. [Google Scholar] [CrossRef]
  35. Anser, M.K.; Hina, T.; Hameed, S.; Nasir, M.H.; Ahmad, I.; Naseer, M.A.U.R. Modeling Adaptation Strategies against Climate Change Impacts in Integrated Rice-Wheat Agricultural Production System of Pakistan. Int. J. Environ. Res. Public Health 2020, 17, 2522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wajid, A.; Hussain, K.; Ilyas, A.; Habib-ur-Rahman, M.; Shakil, Q.; Hoogenboom, G. Crop Models: Important Tools in Decision Support System to Manage Wheat Production under Vulnerable Environments. Agriculture 2021, 11, 1166. [Google Scholar] [CrossRef]
  37. Muntean, N.; Chawdhery, R.A.; Potopová, V.; Tűrkott, L. The ability of CROPGRO-Tomato model to simulate the growth characteristics of Thomas F1 tomato cultivar grown under open field conditions. J. Agric. Sci. 2021, 159, 473–487. [Google Scholar] [CrossRef]
  38. Amin, M.; Zhang, J.; Yang, M. Effects of climate change on the yield and cropping area of major food crops: A case of Bangladesh. Sustainability 2015, 7, 898–915. [Google Scholar] [CrossRef] [Green Version]
  39. Rahman, S. Impacts of climate change, agroecology and socio-economic factors on agricultural land use diversity in Bangladesh (1948–2008). Land Use Policy 2016, 50, 169–178. [Google Scholar] [CrossRef] [Green Version]
  40. Boxall, A.B.; Hardy, A.; Beulke, S.; Boucard, T.; Burgin, L.; Falloon, P.D.; Haygarth, P.M.; Hutchinson, T.; Kovats, R.S.; Leonardi, G.; et al. Impacts of climate change on indirect human exposure to pathogens and chemicals from agriculture. Environ. Health Perspect. 2008, 117, 508–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Sabiha, N.; Rahman, S.; Salim, R. Bangladesh agricultural sustainability: Economic, environmental and social issues. In Bangladesh Economic, Political and Social Issues, Chapter# 1; Alam, K., Ed.; Nova Science Publishers, INC.: New York, NY, USA, 2018; pp. 1–24. ISBN 9781536142105. [Google Scholar]
  42. Hussain, J.; Khaliq, T.; Ahmad, A.; Akhtar, J. Performance of Four Crop Model for Simulations of Wheat Phenology, Leaf Growth, Biomass and Yield across Planting Dates. PLoS ONE 2018, 13, e0197546. [Google Scholar] [CrossRef]
  43. Saddique, Q.; Ji, J.; Ajaz, A.; Jiatun, X.; Yufeng, Z.; He, J.; Cai, H. Performance Comparison of the APSIM and CERES-Wheat Models in Guanzhong Plain, China; ASABE Annual International Meeting: Boston, MA, USA; American Society of Agricultural and Biological Engineers: St. Joseph, MI, USA, 2019; p. 1. [Google Scholar]
  44. Bregaglio, S.; Willocquet, L.; Kersebaum, K.C.; Ferrise, R.; Stella, T.; Ferreira, T.B.; Pavan, W.; Asseng, S.; Savary, S. Comparing Process-Based Wheat Growth Models in Their Simulation of Yield Losses Caused by Plant Diseases. Field Crops Res. 2021, 265, 108108. [Google Scholar] [CrossRef]
  45. Wang, J.; Hasegawa, T.; Li, L.; Lam, S.K.; Zhang, X.; Liu, X.; Pan, G. Changes in Grain Protein and Amino Acids Composition of Wheat and Rice under Short-Term Increased [CO2] and Temperature of Canopy Air in a Paddy from East China. New Phytol. 2019, 222, 726–734. [Google Scholar] [CrossRef] [PubMed]
  46. Khan, M.A.; Tahir, A.; Khurshid, N.; Husnain, M.I.U.; Ahmed, M.; Boughanmi, H. Economic Effects of Climate Change-Induced Loss of Agricultural Production by 2050: A Case Study of Pakistan. Sustainability 2020, 12, 1216. [Google Scholar] [CrossRef] [Green Version]
  47. Munir, Z.; Shrestha, S.; Zaman, M.; Khan, M.I.; Akram, M.M.; Tahir, M.N. Climate Change Impacts on Wheat Yield: A Multi-Modeling Case Study of Central Punjab, Pakistan. Clim. Res. 2022, 87, 13–37. [Google Scholar] [CrossRef]
  48. Kruseman, G.; Mottaleb, K.A.; Tesfaye, K.; Bairagi, S.; Robertson, R.; Mandiaye, D.; Frija, A.; Gbegbelegbe, S.; Alene, A.; Prager, S. Rural Transformation and the Future of Cereal-Based Agri-Food Systems. Glob. Food Secur. 2020, 26, 100441. [Google Scholar] [CrossRef]
Figure 1. Location of the study.
Figure 1. Location of the study.
Ijerph 19 15829 g001
Figure 2. System diagram for study location.
Figure 2. System diagram for study location.
Ijerph 19 15829 g002
Figure 3. Monthly observed average (1990–2020) maximum and minimum temperature (A) and rainfall (B) of Dinajpur.
Figure 3. Monthly observed average (1990–2020) maximum and minimum temperature (A) and rainfall (B) of Dinajpur.
Ijerph 19 15829 g003
Figure 4. Predicted maximum and minimum temperatures (A) and rainfall (B) in Dinajpur.
Figure 4. Predicted maximum and minimum temperatures (A) and rainfall (B) in Dinajpur.
Ijerph 19 15829 g004
Figure 5. Maximum and minimum temperatures from 20 Mid-Century (RCP8.5) GCMs for Dinajpur (the letters in red are selected for the simulation runs, blue is the baseline, and the rest are different GCMs for RCP8.5).
Figure 5. Maximum and minimum temperatures from 20 Mid-Century (RCP8.5) GCMs for Dinajpur (the letters in red are selected for the simulation runs, blue is the baseline, and the rest are different GCMs for RCP8.5).
Ijerph 19 15829 g005
Figure 6. Cumulative distribution functions for observed rice and wheat yields, and APSIM and DSSAT simulated rice and wheat yields.
Figure 6. Cumulative distribution functions for observed rice and wheat yields, and APSIM and DSSAT simulated rice and wheat yields.
Ijerph 19 15829 g006
Figure 7. Simulated rice grain yields (APSIM and DSSAT) compared with baseline and future scenarios.
Figure 7. Simulated rice grain yields (APSIM and DSSAT) compared with baseline and future scenarios.
Ijerph 19 15829 g007
Figure 8. Simulated wheat yields (APSIM and DSSAT) compared with baseline and future scenarios.
Figure 8. Simulated wheat yields (APSIM and DSSAT) compared with baseline and future scenarios.
Ijerph 19 15829 g008
Figure 9. Simulated wheat yields (APSIM) without and with adaptation.
Figure 9. Simulated wheat yields (APSIM) without and with adaptation.
Ijerph 19 15829 g009
Figure 10. Simulated wheat yields (DSSAT) without and with adaptation.
Figure 10. Simulated wheat yields (DSSAT) without and with adaptation.
Ijerph 19 15829 g010
Figure 11. Comparison between gains, losses, and net impacts under current agricultural production system and future agricultural production system by different climate scenarios and crop models (C1–C5 = Climate models; C1 = CCSM4, C2 = GFDL–ESM2M, C3 = HadGEM2–ES, C4 = MIROC5, C5 = MPI–ESM–MR; M1 = crop model (APSIM), M2 = crop model (DSSAT)). (A) Gains, losses, and net impact under current climate. (B) Gains, losses, and net impact under future climate.
Figure 11. Comparison between gains, losses, and net impacts under current agricultural production system and future agricultural production system by different climate scenarios and crop models (C1–C5 = Climate models; C1 = CCSM4, C2 = GFDL–ESM2M, C3 = HadGEM2–ES, C4 = MIROC5, C5 = MPI–ESM–MR; M1 = crop model (APSIM), M2 = crop model (DSSAT)). (A) Gains, losses, and net impact under current climate. (B) Gains, losses, and net impact under future climate.
Ijerph 19 15829 g011
Table 1. RAPs parameters for Bangladesh used in the model and adaptation strategy for wheat.
Table 1. RAPs parameters for Bangladesh used in the model and adaptation strategy for wheat.
ParametersDirectionMagnitude (%)
Farm sizeDecrease20
Family sizeDecrease10
Variable costs of productionIncrease140
Off farm incomeIncrease100
Rice yield trend without CCIncrease112.5
Rice price trend without CCIncrease127.0
Rice price trend with CCIncrease150.3
Wheat yield trend without CCIncrease111.2
Wheat price trend without CCIncrease126.8
Wheat price trend with CCIncrease156.8
Adaptation strategy for wheat
Irrigation in wheatIncrease50 mm per farm
Table 2. Experimental dataset for genetic coefficient calculation for rice and wheat.
Table 2. Experimental dataset for genetic coefficient calculation for rice and wheat.
RiceWheat
Sowing date: 31 May 2018Sowing date: 24 November 2018
Transplanting: 5 July 2018Sowing density (plants m−2): 240
Sowing depth (cm): 5 cmSowing depth (cm): 5
Row spacing (cm): 20 cmRow spacing (cm): 20
Fertilizer and irrigation management:Fertilizer and irrigation management:
N fertilizer: urea NN fertilizer: urea N
Applied N (kgha−1): 90Applied N (kgha−1): 120 (2 split)
1st N application amount (kgha−1): 301st N application amount (kgha−1): 80 (as basal)
2nd N application amount (kgha−1): 302nd N application amount (kgha−1): 40 (14 December 2018)
3rd N application amount (kgha−1): 301st irrigation: 14. December 2018
2nd irrigation: 31 December 2018
3rd irrigation: 24 January 2019
4th irrigation: 7 February 2019
Table 3. Calibration outputs (Rice-APSIM and DSSAT) for cultivar BR11.
Table 3. Calibration outputs (Rice-APSIM and DSSAT) for cultivar BR11.
Stage Code Date (Comparison)
StageObservedSimulated (APSIM)Simulated (DSSAT)
1Sowing31 May 201831 May 20182 June 2018
2Transplanting5 July 20185 July 201815 July 2018
3End of juvenile10 August 201810 August 201821 August 2018
4Floral initiation7 October 20183 October 201820 September 2018
5Flowering25 October 201825 October 201825 October 2018
6Maturity24 November 201824 November 201823 November 2018
7Grain yield (kg ha−1)461047436811
8Biomass yield (kg ha−1)10,31013,21214,171
Table 4. Calibration outputs (Wheat-APSIM and DSSAT) for cultivar Shatabdi.
Table 4. Calibration outputs (Wheat-APSIM and DSSAT) for cultivar Shatabdi.
Stage Code Date (Comparison)
StageObservedSimulated (APSIM)Simulated (DSSAT)
1Sowing24 November 201824 November 201824 November 2018
2Germination25 November 201825 November 201825 November 2018
3Emergence3 December 201829 November 201827 November 2018
4End of juvenile15 December 201815 December 201816 December 2018
5Floral initiation17 January 201917 January 201920 January 2019
6Flowering31 January 201931 January 201931 January 2019
7Start grain fill6 February 20196 February 20199 February 22019
8End of grain fill13 March 201911 March 201911 March 2019
9Maturity14 March 201913 March 201913 March 2019
Grain yield (kg ha−1)520048204809
Biomass yield (kg ha−1)10,000820310,378
Table 5. Climate sensitivity for rice–wheat system in Dinajpur, Bangladesh.
Table 5. Climate sensitivity for rice–wheat system in Dinajpur, Bangladesh.
ParametersIKXAIOXAIRXAIEXAIIXA
APSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSAT
Observed average crop yield (Rice) (kg ha−1)3681372535983637352435893614368135533634
Average crop yield change (Rice) (%)
[defined as: (average relative crop yield − 1) × 100]
−37.40−14.70−13.60−7.40−38.17−13.40−20.30−8.80−21.70−10.10
Observed average crop yield (Wheat) (kg ha−1)2789282728522889279828412837289428172910
Average crop yield change (Wheat) (%)
[defined as: (average relative crop yield − 1) × 100]
9.00−6.0020.001.0021.80−8.800.80−6.0018.70−2.90
Losers (%)80.1071.3448.5158.0073.3972.0371.5465.3159.3164.25
Gains (% average net farm returns)11.7812.4816.2314.3312.2412.3512.5613.1914.1013.35
Losses (% average net farm returns)−25.23−20.79−15.69−17.22−21.54−20.96−21.04−18.87−17.47−18.60
Observed net earnings without climate change (BDT farm−1)16,90316,90316,90316,90316,90316,90316,90316,90316,90316,903
Observed net earnings with climate change (BDT farm−1)13,88715,00117,02916,23314,78314,93614,96415,59516,12215,690
Observed per-capita earnings without climate change (BDT)14,79514,79514,79514,79514,79514,79514,79514,79514,79514,795
Observed per-capita earnings with climate change (BDT)14,21114,42714,81914,66514,38414,41414,41914,54114,64314,560
Observed poverty rate without climate change (%)96.0796.0796.0796.0796.0796.0796.0796.0796.0796.07
Observed poverty rate with climate change (%)96.5896.4096.0596.1996.4396.4196.4196.3096.2196.28
Note: Exchange rate: USD 1 = BDT 84.86. Poverty line = USD 1.25.
Table 6. Impact of climate change in rice–wheat system in Dinajpur, Bangladesh.
Table 6. Impact of climate change in rice–wheat system in Dinajpur, Bangladesh.
ParametersIKXAIOXAIRXAIEXAIIXA
APSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSAT
Projected average yield (Rice) (kgha−1)4537458345614625468147114516458445924648
Average yield change (Rice) (%)
[defined as: (average relative yield − 1) × 100]
−37.40−14.70−13.60−7.40−38.17−13.40−20.30−8.80−21.70−10.10
Projected average yield (Wheat) (kgha−1)3248335632543315331433573311340634513487
Average yield change (Wheat) (%)
[defined as: (average relative yield − 1) × 100]
9.00−6.0020.001.0021.80−8.800.80−6.0018.70−2.90
Losers (%)58.2748.1524.3333.5946.8449.2619.0315.6012.8915.03
Gains (% average net farm returns)16.6818.0525.0321.3118.2017.8228.1332.5036.2333.11
Losses (%average net farm returns)−20.17−17.31−13.34−14.50−16.94−17.52−12.78−13.13−13.19−13.10
Projected net earnings without climate change (BDTfarm−1)23,00023,00023,00023,00023,00023,00023,00123,00623,01723,008
Projected net earnings with climate change (BDTfarm−1)21,89723,23526,60925,13423,40123,09427,68028,83629,86429,017
Projected per-capita earnings without climate change (BDT)30,53130,53130,53130,53130,53130,53130,53130,53230,53530,533
Projected per-capita earnings with climate change (BDT)30,29430,58231,30630,98930,61730,55131,53631,78432,00531,823
Projected poverty rate without climate change (%)58.7458.7458.7458.7458.7458.7458.7458.7458.7458.74
Projected poverty rate with climate change (%)59.1158.6857.6058.0758.6258.7257.2556.8856.5556.82
Note: Exchange rate: USD 1 = BDT 84.86. Poverty line = USD 1.25.
Table 7. Impacts of adaptation in rice–wheat system in Dinajpur, Bangladesh.
Table 7. Impacts of adaptation in rice–wheat system in Dinajpur, Bangladesh.
ParametersIKXAIOXAIRXAIEXAIIXA
APSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSATAPSIMDSSAT
Projected average yield without adaptation (Rice) (kg farm−1)2261304130933306222530853605413835574076
Average yield change (Rice) (%)
[defined as: (average relative yield − 1) × 100]
0000000000
Projected average yield without adaptation (Wheat) (kgfarm−1)3558295332002930319529104284332536333265
Average yield change (Wheat) (%)
[defined as: (average relative yield − 1) × 100]
22.302.1010.301.3010.100.7032.203.3012.501.50
Adoption rate (%)63.4658.1741.0450.1237.6959.3574.1458.3444.5153.91
Projected net earnings without adaptation (BDTfarm−1)21,89723,23526,60925,13423,40123,09427,68228,83629,86429,017
Projected net earnings with adaptation (BDT farm−1)24,72325,61128,14427,10824,60625,56232,98631,84631,82331,620
Projected per-capita earnings without adaptation (BDT)30,29430,58231,30630,98930,61730,55131,53631,78432,00531,823
Projected per-capita earnings with adaptation (BDT)30,90131,09231,63531,41330,87631,08132,67532,43032,42532,382
Projected poverty rate without adaptation (%)59.1158.6857.6058.0758.6258.7257.2556.8856.5556.82
Projected poverty rate with adaptation (%)58.2157.9257.1157.4458.2557.9455.5855.9355.9456.00
Note: Exchange rate: USD 1 = BDT 84.86. Poverty line = USD 1.25.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chawdhery, M.R.A.; Al-Mueed, M.; Wazed, M.A.; Emran, S.-A.; Chowdhury, M.A.H.; Hussain, S.G. Climate Change Impacts Assessment Using Crop Simulation Model Intercomparison Approach in Northern Indo-Gangetic Basin of Bangladesh. Int. J. Environ. Res. Public Health 2022, 19, 15829. https://doi.org/10.3390/ijerph192315829

AMA Style

Chawdhery MRA, Al-Mueed M, Wazed MA, Emran S-A, Chowdhury MAH, Hussain SG. Climate Change Impacts Assessment Using Crop Simulation Model Intercomparison Approach in Northern Indo-Gangetic Basin of Bangladesh. International Journal of Environmental Research and Public Health. 2022; 19(23):15829. https://doi.org/10.3390/ijerph192315829

Chicago/Turabian Style

Chawdhery, Md Rafique Ahasan, Murtuza Al-Mueed, Md Abdul Wazed, Shah-Al Emran, Md Abeed Hossain Chowdhury, and Sk Ghulam Hussain. 2022. "Climate Change Impacts Assessment Using Crop Simulation Model Intercomparison Approach in Northern Indo-Gangetic Basin of Bangladesh" International Journal of Environmental Research and Public Health 19, no. 23: 15829. https://doi.org/10.3390/ijerph192315829

APA Style

Chawdhery, M. R. A., Al-Mueed, M., Wazed, M. A., Emran, S. -A., Chowdhury, M. A. H., & Hussain, S. G. (2022). Climate Change Impacts Assessment Using Crop Simulation Model Intercomparison Approach in Northern Indo-Gangetic Basin of Bangladesh. International Journal of Environmental Research and Public Health, 19(23), 15829. https://doi.org/10.3390/ijerph192315829

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