Next Article in Journal
A Case Analysis of Dust Weather and Prediction of PM10 Concentration Based on Machine Learning at the Tibetan Plateau
Next Article in Special Issue
Impact of Climate and Land-Use Change on Groundwater Resources, Study of Faisalabad District, Pakistan
Previous Article in Journal
Trends of Tropical Cyclone Translation Speed over the Western North Pacific during 1980−2018
Previous Article in Special Issue
Possible Relationships between the Interdecadal Anomalies of Heavy Rainfall under Northeastern China Cold Vortexes and the Sea Surface Temperature (SST)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Future Climate Projections Using SDSM and LARS-WG Downscaling Methods for CMIP5 GCMs over the Transboundary Jhelum River Basin of the Himalayas Region

1
Department of Geography, University of Gujrat, Gujrat 50700, Pakistan
2
Department of Civil Engineering, College of Ocean Science, Jeju National University, 102 Jejudaehakro, Jeju 63243, Korea
3
Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan 430079, China
4
Department of Civil, Environmental and Natural Resources Engineering, Lulea University of Technology, 97187 Lulea, Sweden
5
Institute of Applied Technology, Thu Dau Mot University, Thu Dau Mot City 7500, Binh Duong Province, Vietnam
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(6), 898; https://doi.org/10.3390/atmos13060898
Submission received: 2 April 2022 / Revised: 27 May 2022 / Accepted: 27 May 2022 / Published: 1 June 2022

Abstract

:
Climate change is one of the leading issues affecting river basins due to its direct impacts on the cryosphere and hydrosphere. General circulation models (GCMs) are widely applied tools to assess climate change but the coarse spatial resolution of GCMs limit their direct application for local studies. This study selected five CMIP5 GCMs (CCSM4, HadCM3, GFDL-CM3, MRI-CGCM3 and CanESM2) for performance evaluation ranked by Nash–Sutcliffe coefficient (NSE) and Kling–Gupta Efficiency (KGE). CCSM4 and HadCM3 large-scale predictors were favored based on ranks (0.71 and 0.68, respectively) for statistical downscaling techniques to downscale the climatic indicators Tmax, Tmin and precipitation. The performance of two downscaling techniques, Statistical Downscaling Methods (SDSM) and Long Ashton Research Station Weather Generator (LARS-WG), were examined using the Mean Absolute Error (MAE), Root Mean Square Error (RMSE), bias, NSE and KGE with weights (Wi) for the validation period. The results of statistical measures proved SDSM more efficient (0.67) in comparison to the LARS-WG (0.51) for the validation time for the Jhelum River basin. The findings revealed that the SDSM simulation for Tmax and Tmin are more comparable to the reference data for the validation period except simulation of extreme events by precipitation. The 21st century climatic projections exhibited a significant rise in Tmax (2.37–4.66 °C), Tmin (2.47–4.52 °C) and precipitation (7.4–11.54%) for RCP-4.5 and RCP-8.5, respectively. Overall, the results depicted that winter and pre-monsoon seasons were potentially most affected in terms of warming and precipitation, which has the potential to alter the cryosphere and runoff of the Jhelum River basin.

1. Introduction

Climate change is becoming a leading issue for the 21st century due to its devastating environmental and socioeconomic impacts. In the last few decades, the frequency and magnitude of extreme climatic events increased subsequently in response to the anthropogenic activities [1]. Anthropogenic activities, primarily socioeconomic (fossil fuels burning and land use/land cover changes), have influenced the amount of greenhouse gases which trigger climate change and extreme climate events. The occurrences of extreme events are not uniform across the globe; some regions are more susceptible to climate change. Particularly, Pakistan, has faced frequent heatwaves and floods in the last few years [2,3,4]. To cope with these extreme climatic events, the timely and effective monitoring of climate change is required to make policies for adaptation and mitigation. The impacts of climate change on water resources and the hydrological cycle are of extreme importance because all socioeconomic and natural systems ultimately depend on water resources. The climatic changes can directly impact the availability and changing patterns of water resources such as flooding and droughts [3,5], and some indirect impacts on food, agriculture and energy production [6]. These climate change impacts may be worse for the transboundary river’s basins such as Jhelum where management of the basin depends on different economic, political, and social interest of the countries. Jhelum River basin is an integral part of the Himalayas region where an increasing trend of temperature has been observed that increases glacial melting and precipitation, and affects the availability of water resources [7,8]. Climatic change studies in the Jhelum River basin is still at its infancy due to a lack of significant weather station data. To understand the impacts of climate change on the transboundary Jhelum River basin of the Himalayas region, GCMs have been used to assess the present and future climatic changes.
The GCMs provide projections of climate at a global scale for policymakers to adapt better strategies to cope with climate change. GCMs represent significant outputs at the global, hemispherical and continental scales by incorporating the complexity of the global system, however, these global dynamics cannot be represented at the local sub-grid level [9]. The efficiency of GCMs to project the future climate has been debated due to their uncertainties during the validation processes at the regional scale. Despite improvements in the GCMs to represent climate processes in better ways, these uncertainties cater to produce better climate projection but still remain a subject of ample concern at the regional/local level [10]. GCMs are widely used tools to assess the climate change impacts but their coarse spatial resolution restricts direct use for the sustainable management at the regional or local scale [11]. The downscaling techniques are essential to transform GCMs’ spatial resolution from coarse to fine to allow their direct use at the local/regional scale [12]. The two widely used downscaling methods, statistical downscaling and dynamical downscaling applied to relate the GCMs’ coarse resolution and local climatic variables [13]. Dynamic downscaling of GCMs are employed as a Regional Climate Model (RCM) at finer spatial resolution (10–50 km) to simulate regional climate by incorporating local features such as topography. Dynamic downscaling is an emerging and advanced method, but the advanced computational requirement and heavy data storage limit their use at regional scale [2]. Statistical downscaling is both a flexible and computationally efficient approach to downscale GCMs and to use fine resolution data for a climate impact assessment at the local/regional level [13]. Statistical downscaling has directly built a relationship between local observation, climatic variables and GCMs’ output without requiring the physical knowledge of the local region [14]. Therefore, statistical downscaling methods have been extensively applied by the researchers to simulate climate projections at the local/regional scale for the climate impact studies [15].
The main theme of statistical downscaling is to develop the relationship among predictors (GCMs variables) and predictands (local scale variables) through statistical and mathematical techniques such as linear and non-linear regression models, and weather generators [13]. Among linear regression models, the Statistical Downscaling Method (SDSM) is a renowned statistical model developed by [16] that is frequently used by research to downscale GCMs [17]. SDSM is a hybrid model that employed the weather generators and regression models to downscale climatic variables. It facilitates the downscaling of long-term, low-cost, and rapid development of multiple daily weather parameters. The weather generator’s technique, Long Ashton Research Station Weather Generator (LARS-WG), is a well-known stochastic weather generator technique used to simulate the weather data for a single weather station in the form of time series data for both present and future climatic conditions. The long-term time series data of a climate variables group e.g., Tmax, Tmin and precipitation are simulated for the single weather station using the LARS-WS method [6].
SDSM and LARS-WG have been widely used techniques by the researchers to downscale the GCMs’ data for local/regional basins [6,9]. These techniques have been used for three GCMs (BCC-CSM1-1, CanESM2 and MICROC5) and future projections depicted that mean annual temperature and precipitation for the future are on the rise [18]. These studies demonstrated statistical downscaling methods as a vigorous tool to analyze the futuristic climatic changes for the regional/local basin level. SDSM and LARS-WG have been used for the better assessment of climate changes in the Jhelum River basin using different GCMs. The study area is of key importance as it is part of the Indus basin and greater Himalayas that have permafrost mountain tops. The climate changes ultimately trigger the melting of snow/glacier at mountain tops. The recent study focused on examining the efficiency of these statistical downscaling techniques SDSM as a regression model and LARS-WG as weather generators for downscaling the Tmax, Tmin and precipitation data for the Jhelum River basin. The basin is the transboundary and conflicted region located at the greater Himalayas, therefore, future projections of Tmax, Tmin and precipitation will help to study the dynamics of hydrometeorological changes in the basin. The study designed a methodology to incorporate multiple GCMs for the basin based on local conditions. The selection of GCMs helped to downscale the long-term time series climatic data for the 21st century by using two different statistical downscaling techniques (SDSM and LARS-WG). The meteorological station’s data of Jhelum River basin was applied to evaluate the accuracy of SDSM and LARS-WG. After evaluation of the statistical downscaling techniques, climate change projections were simulated for RCP 4.5 and RCP 8.5 using six GCMs for the 21st century.

2. Materials and Methods

2.1. Study Area

The Jhelum River is the major tributary of the Indus basin located in the north of Pakistan. The basin is the transboundary conflicted region divided by the line of control (LoC) between India and Pakistan. The geographical extent of the region exists between 73–75.63° E and 33–35.1° N, covering a total area of about 34,475 km2 (Figure 1). The river is fed by the glacier/snow melting of the drains from the top of the Himalayas mountains. The basin has a diverse altitudinal variation from 233 to 6178 m that covers the permanent snow-covered area in the north. The river basin has unique geomorphology and heterogeneous lithology with varied hydrological conditions which add up to the basin being more susceptible to climate change [19]. The climate of the basin has diversity in terms of spatio-temporal variability. The precipitation in the basin is dominant by the two seasonal rainfall regimes, monsoons during summer and western disturbances during winter. The details of the eleven meteorological stations of the basin are described in Table 1.

2.2. Data Description

The predictands of daily time series data for Tmax and Tmin and precipitation were acquired from the Pakistan Meteorological Department (PMD) from 1976–2020. Few stations have some missing values that were replaced by using interpolation technique [15]. The datasets for 1976–2005 of eleven meteorological stations were applied for the evaluation of the GCMs and selection of the GCMs for the downscaling techniques. The thirty years of data (1976–2005) were used as the baseline data for statistical downscaling and bias correction of GCMs of the daily time series for the period 1976–2099.
The projected climate changes are highly uncertain and rely on technological innovations and socioeconomic development to adapt clean technology against greenhouse gas emissions in the atmosphere. The Coupled Model Intercomparison Project Phase 5 (CMIP5) in accordance of the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC) used RCPs that provide a wider picture of future climate change including a mitigation scenario (RCP-2.6), two medium stabilization scenarios (RCP-4.5 and RCP-6.0) and an extreme scenario (RCP 8.5). The five GCMs of CMIP5 for one medium stabilization scenario (RCP-4.5) and an extreme scenario (RCP-8.5) were downloaded (Table 2) from the Royal Netherlands Meteorological Institute’s (KNMI) Climate Explorer [20].
Twenty different NCEP predictors applied for the screening purposes (Table 3) and based on the relationship among the NCEP predictors and local basin predictands (Tmax and Tmin and precipitation), following predictors selected for the SDSM.

2.3. Mann–Kendall Trend Model

The Mann–Kendall trend test (MK) is a non-parametric model applied to assess the climatic trends in the long-term time series data [21]. MK trend model was applied on the Tmax, Tmin and precipitation data for the baseline period using the following equations.
S = i = 1 n 1 j = i + 1 n s g n x j x i ,  
where 𝑥𝑖 and 𝑥j are the orderly data records in the 𝑖 and 𝑗 years; 𝑛 is the span of the time series data.
w h e r e   s g n x j x i , = + 1 , x j x i > 0 0 , x j x i = 0 1 , x j x i < 0
V a r   S = 1 18 n n 1 2 n + 5 p = 1 q t p t p 1 2 t p + 5
In Equation (3), q is the number of tied groups whereas tp is the number of observations in the pth group.

2.4. Screening of the GCMs

The five GCMs were further evaluated based on the accuracy assessment to select two GCMs for the downscaling purposes because the study focused on the evaluation of the downscaling processes SDSM and LARS-WG that already made it intensive. The accuracy assessment of GCMs were carried out by using gridded Tmax, Tmin and precipitation datasets cross-validated against the observed station data sets [22]. The gridded Tmax, Tmin and precipitation are interpolated using the kriging geospatial technique. The correspondence values of GCMs for climatic indicators (Tmax, Tmin and precipitation) at the meteorological stations were used for the performance evaluation of GCMs. The indices were based on the simulated GCMs’ values and the observed monthly time scale climatic indicators (Tmax, Tmin and precipitation). The performance indicators between reference datasets and GCMs modeled data of climatic indicators by using Pearson’s Correlation Coefficient (r), Kling–Gupta Efficiency (KGE) [19] and Nash–Sutcliffe coefficient (NSE) [20]. The detailed formulas of r, KGE and NSE are as follows:
r = x x i ¯ ( x i x i ) ¯ x i x i ¯ 2 ( x i x i ) ¯ 2
KGE = 1 ( r 1 ) 2 + β 1 2 + ( γ 1 ) 2
NSE = 1 i = 1 n ( x o   x s ) 2 i = 1 n ( x s   x o ) 2
The description of ‘x’ are the climatic indicators (Tmax, Tmin and precipitation) while subscripts s and o denoted the observed referenced values and simulated modeled values, respectively. Equation (4)represents the Pearson’s correlation coefficient to measure the linear relationship between observed and simulated values. Equation (5) describes the KGE equation including values for г, β and ϒ while Equation (3) describes the NSE. The resultant values were normalized and rescaled between 0 and 1 and then finally summed up to generate rankings for each GCM.

2.5. Bias Correction

The bias correction method mean-based biased correction method (MB-BC) applied to eliminate errors from the modeled simulated data. The MB-BC method utilized the mean observed data and GCM simulation for the baseline period [21]. Following, the different Equations (7) and (8) were used for the bias correction of temperature and precipitation to avoid negative values for precipitation.
T de biased =   T sim 2020 2099 × T sim ( 1976 2005 ) T obs 1976 2005
P de biased =   P sim 2020 2099 × P obs 1976 2005 P sim 1976 2005

2.6. Statistical Downscaling Method

Two downscaling methods (SDSM and LARS-WG) were applied for the eleven meteorological stations of Jhelum River basin. The study aimed to compare the efficiency of the downscaling schemes for the selected GCMs. The three predictands (Tmax, Tmin and precipitation) were considered against the GCM predictors to analyze the climate changes in the river basin for the 21st century.

2.6.1. SDSM

SDSM is the hybrid scheme, a combination of Stochastic Weather Generators (SWG) and Multiple Linear Regression (MLR). MLR was used to develop the relationship between the predictands (Tmax, Tmin and precipitation) and NCEP predictors (large-scale climate data) by yielding regression parameters. These regression parameters were further used for screening and calibration of NCEP predictors. After screening, the NCEP predictors were established on the correlation analysis of large-scale predictors with the local scale predictands. The screening of the GCM’s large-scale predictors is the key to success of the SDSM and yields through regression parameters were established on the correlation matrix and absolute partial correlation [8].
The regression parameters, NCEP predictors and GCM data (CCSM4 and HadCM3) were further used by the SWG scheme to simulate the time series data for the 21st century [13]. The precipitation projections were challenging as compared to the temperature due to its complex nature. The precipitation was modeled by using a conditional sub model and SWG conditioned and applied on the predictors. The conditional sub model was designed for the precipitation [22] as it occurs on each day t or not by using the following equation:
P t = 0 + i = 1 n i u t i + t 1 w t 1
The P t   is the conditional processes for precipitation of each day t; u t i is used to normalized predictor, while i is regression parameter; t 1 and   w t 1 are the conditional possibilities for precipitation occurrence on day t − 1.

2.6.2. LARS-WG Model

LARS-WG is the stochastic weather generator that was applied for the simulation of weather data for both the present and future climatic conditions. The synthetic weather data were generated for three timelines: model calibration (1976–2005), model validation (2006–2020) and scenario generation (2021–2099). For model calibration (1976–2005), the daily observed series data for Tmax, Tmin and precipitation were applied to examine parameters for probability distribution. The observed weather data (Tmax, Tmin and precipitation) were used in LARS-WG to generate time series of erratic length by selecting random values from the stations’ distributions [9]. LARS-WG was based on the semi-empirical distribution that used cumulative probability distribution function to approximate probability of dry and wet series for Tmax and Tmin and dry and wet days for the precipitation [13]. The future climate scenarios were generated for the period (2021–2099) for selected RCPs based on the LARS-WG baseline parameters (1976–2005). In LARS-WG, the differences in current and future climatic changes were incorporated by using bias correction [23], which is the difference of the mean monthly changes of Tmax, Tmin and precipitation. The following equations are used for T i and P i :
T i = T 2021 2099 T 1976 2005
P i = P 2021 2099 P 1976 2005
The T i and P i are to represent long term changes of each month ( 1 i 12 ), for temperature and precipitation, respectively.

2.6.3. SDSM and LARS-WG Performance Evaluations

SDSM and LARS-WG performance evaluations were assessed by using statistical techniques. The statistical errors (Equations (12)–(14)); mean absolute error (MAE), root mean square error (RMSE) and bias were used along with NSE and KGE. The NSE and KGE are hypersensitive to extreme values and inconsiderable to the proportional and additive differences that exist between the observation and model simulations. The equations for NSE and KGI were described in the above section (Equations (4) and (5)). Following are the equations for MAE, RMSE and bias:
MAE = i = 1 n X o X s n
RMSE = 1 n i = 1 n X o X s 2
Bias = i 1 i = n x o n i 1 i = n x s n
X o and X s are the observation and simulated data in the above equations by statistical downscaling, respectively, and n is the number of all samples. Additionally, weight technique is used for performance matrices incorporating all the measurements errors to avoid discrepancies. More weight (0.5) was assigned to measurement errors (MAE, RMSE and bias) while less weight (0.15) was applied to the NSE and KGE due to their oversensitivity to extreme values [6]. Equation (15) is used to rank the SDSM and LARS-WG model based on the weights assigned to the different measurement errors. The performance evaluation of the downscaling methods were summarized by using summation of the weights of each measurement error by using Equation (15).
W i = W MAE MAE i MAE max + W RMSE RMSE i RMSE max + W Bias Bias i Bias max + W KGE KGE i KGE max + W NSE NSE i NSE max
where the i index represents to the downscaling model; W i refers to the overall performance measure.

3. Results and Discussion

3.1. Temperature and Precipitation Trend Analysis for the Baseline Period

The Mann–Kendall trend model was applied on the eleven meteorological stations of the Jhelum River basin. In this test, the Kendall’s tau denotes the strength of the monotonic trend where the value of Kendall’s tau ranges from −1 ≤ τ ≤ 1. The positive Kendall’s tau indicated an increasing trend whereas the negative values denoted decreasing. Based on a 95% significance level, the p value is ≤ = 0.05 and indicated the presence of a significant trend while p value ≥ = 0.05 verified that there was no trend in the time series data. The results found that there were different trends for all stations. The significant trend for Tmax depicted for Station D station have a p value less than 0.05, while the rest of the stations showed no significant trend (Table 4).
The Mann–Kendall trend for Tmin predicted in line trend results of Tmax. The metrological stations D, I and another station, K, have shown a statistical trend. Other stations represented an insignificant trend for the historical time series data (Table 5).
The Mann–Kendall statistics results for precipitation illustrated no significant trend for most of the basin during the historical period. The three stations B, H and K depicted a negative trend for precipitation (Table 6).

3.2. Selection of GCMs

It was difficult to perform statistical downscaling (SDSM and LARS-WG) for all five GCMs for the river basin. Therefore, the selection of GCMs were made by using Equations (4)–(6) and the best GCMs were chosen for further processing based on ranks described in Table 7. The best proved GCMs were CCSM4 and HadCM3, in relationship to the reference data (Tmax, Tmin and precipitation), with 0.71 and 0.68 ranking, respectively (Table 7).

3.3. Screening Predictor Variable of SDSM

The twenty predictors (Table 2) analyzed by absolute partial correlation coefficient (abs P) between the predictors and predictands are summarized in Figure 2. Among the selected thirteen predictors, 500 hPa vorticity (p5_z) depicted the highest relationship for Tmax, Tmin and precipitation while the predictor surface divergence (p_zh) showed the weakest relationship for the predictands. The best predictor for Tmax was surface vorticity (p_z) with highest value (0.59) and minimum relationship (0.22) of predictor with 500 hPa geopotential height (p500). The best predictor for Tmin was 500 hPa vorticity (p5_z) and for precipitation was 850 hPa relative humidity (r850). The weakest relationship of Tmin existed with predictor mean temperature at 2 m (temp) and of precipitation with surface divergence (p_zh). These thirteen predictors were further used in the downscaling process and the selection criteria, based on the previous studies [24].

3.4. Calibration and Validation of SDSM and LARS-WG

The GCMs’ (CCSM4 and HadCM3) data were downscaled by using statistical downscaling methods (SDSM and LARS-WG) during the calibration period (1976–2005). The results of calibrations are plotted against the observed referenced data in Figure 3. The calibration results depicted that SDSM yields more accurate and good quality time series data for Tmax and Tmin against the observations for the calibration period. On the contrary, LARS-WG predicted more accurate precipitation time series data. The precipitation is the complex indicator and LARS-WG proved to be more accurate in recording extreme events and precipitation distribution. The results of the calibration period simulated against RCP-4.5 were more accurate in comparison to the RCP-8.5, as the later simulation were more exaggerated from the former one (Figure 3). The precipitation for the summer months i.e., June, July and August (JJA) were underestimated while for winter months, namely, December, January and February (DJF) were overestimated. The simulated results of Tmax for the calibration period were more accurate as compared to the Tmin which was a bit exaggerated for the winter months (Figure 3).
The statistical downscaling models (SDSM and LARS-WG) were validated for the period of (2006–2020) by using statistical measures MAE, RMSE, bias, KGE and NSE. These measures were weighed as discussed in the Materials and Methods section following Equation (15). The values for the statistical errors were rescaled between 0–1 to avoid negative values for the precipitation. The performance assessment of the models for the validation for each of the climatic parameters are described in Table 8. The results of statistical measures proved SDSM as more efficient (0.67) in comparison to the LARS-WG (0.51) for the validation period for the Jhelum River basin. SDSM was further applied to generate future climate changes for the 21st century due to its better performance. The validation results for RCPs were in line with the calibration period and RCP-8.5 was exaggerated in comparison to the RCP-4.5 (Table 8). The simulated indicator Tmax proved best during the validation and precipitation was depicted the least accurate due to its complex nature (Table 8).

3.5. Future Climate Projections

The projected annual and seasonal changes for Tmax, Tmin and precipitation under both RCPs (RCP-4.5 and RCP-8.5) for the mid-century (2021–2060) and end-century (2061–2099) over the Jhelum River basin are summarized in Table 9. Overall, an annual increase in temperature (Tmax and Tmin) was observed for both RCPs, however, the temperature rise for RCP-8.5 (4.48 °C) was expected to be more significant than RCP-4.5 (2.52 °C) (Table 9). The rate of increase in temperature and precipitation are in line with the previous studies in the basin [25]. As the climate change triggered changes in the river basin by making it hot and wet for the upcoming century, it can trigger other physical changes in the basin by altering the balance of glaciation, deglaciation and flooding.
The 21st century projection for Tmax was about 2.54 °C–4.66 °C and for Tmin was about 2.53 °C–4.52 °C for RCP-4.5 and RCP-8.5 for the whole basin. The seasonal change in projected Tmax and Tmin depicted the highest temperature changes for the winter season (DJF) followed by the pre-monsoon (AMJ) period (Table 9) under both scenarios for mid-century and end-century. This warming of the winter season may trigger the earliest deglaciation period due to a rise in critical temperature by changing the hydrologic cycle. The rise in critical temperature during the winter season may alter the precipitation form from snow to rainfall, consequently shrinking the snow/glacier-covered area in the Jhelum River basin that leads in a reduction of summer river discharge. In addition, the temperature increase for the pre-monsoon period could generate early and rapid snow/glacier melting from a high altitude/elevation of the basin that might cause shrinking of the glacier/snow-covered area of the basin. The temperature increase for RCP-8.5 simulations were more intense than RCP-4.5, which could generate more severe consequences for the cryosphere and hydrosphere of the river basin.
The precipitation projection depicted an increase for both scenarios, but a slight decrease was observed for RCP-8.5 in comparison to the RCP-4.5. The rise in precipitation ranges were between 6–12% over the basin. The rate of increase in precipitation were more considerable for the winter and monsoon season for both scenarios [25]. The minimum increase in precipitation was observed for the pre-monsoon period. The spatial distribution of Tmax, Tmin and precipitation was immensely significant due to the powerful impacts of significant climatic patterns [26,27]. The spatial patterns for precipitation and temperature for the mid- and end-century for both GCMs under climate scenarios are presented in Figure 4 and Figure 5. The spatial distribution of precipitation changes is extremely important because the basin is influenced by two contrasting climatic systems: monsoon during summer and westerlies circulation pattern during winters. The projected rise in winter precipitation is due to the southwest side of the basin that is affected by the westerlies while the rise in precipitation in the northeastern region of the basin is due to the monsoon precipitation during the summer season (Figure 4). The precipitation spatial pattern is heterogenous in the basin and, overall, depicted a rise for both RCPs for the 21st century.
The spatial distribution of temperature projections (Figure 5) depicted that the northeastern part of the basin was more affected from a rise in temperature for both GCMs under climate scenarios. The southeastern part of the basin was noticed to be least impacted for a rise in temperature for 21st century. The transition in class boundaries was observed, especially in the northern part of the basin, which may cause deglaciation due to a rise in temperature in the high altitudinal (4000–6000 masl) areas of the basin. The temperature changes in the northern mountainous region might disturb the balance of the cryosphere and hydrosphere by altering the permanent snow-covered area of the Jhelum River basin. The temperature changes are considered a more significant variable because they trigger changes in precipitation forms and snow/glacial reserves. The projected changes in precipitation and temperature may accelerate the surface runoff during the mid-century and can introduce recurrent droughts in the end-century due to disturbance in the hydrosphere and cryosphere. These changes in temperature and precipitation patterns from the mid- (2021–2060) to the end-century (2061–2099) were more prominent and severe for RCP-8.5 than RCP 4.5.

4. Conclusions and Remarks

This study focused on the climatic changes for the transboundary mountainous region by using CMIP5 GCMs and climate scenarios. The Jhelum River basin is the part of the Himalayas mountains which are affluent in glaciers after polar regions and climate changes can induce tremendous, unexpected variabilities in the basin. GCMs are the most reliable mathematical tools to project climatic changes as it is a global phenomenon and linked through permanent circulations such as atmospheric winds and oceanic currents. The GCM projections for the local/regional basin may cause uncertainties due to the influence of the micro/regional climate. To overcome the uncertainties, five GCMs from CMIP5 project were selected for the study and two GCMs, CCSM4 and HadCM3, were chosen based on the relationship with the reference observed data. The large-scale GCM predictors were downscaled by using two statistical downscaling methods (SDSM and LARS-WG) that can relate large-scale predictors with the local climatic variables (Tmax, Tmin and precipitation).
SDSM was proved to be a more effective and efficient method than LARS-WG in downscaling the climatic variables for the validation period, except for a bit in improved performance of LARS-WG for precipitation due to its wet and dry spells. Both climatic indicators of temperature and precipitation projected a rising trend for the 21st century, but the change was more pronounced for RCP-8.5 than RCP-4.5, which depicted the storyline of the climatic scenario. The seasonal changes depicted that the winter season was more threatened and became warmer and wetter with time which may disturb the existing snow cover. After winter, the pre-monsoon season predicted a rising trend for temperature that may introduce early deglaciation by disturbing the balance of the cryosphere and hydrosphere of the basin.
The study area is a transboundary conflicted region and has a limited number of weather stations that affect the accuracy of the climatic projections. Multiple GCMs can be used to check the most suitable location for the study area but it can make the research more intensive as the focus is on the multiple task of first selecting the GCMs then downscaling by using two comparable methods. The dynamic downscaling can also be used but it is a very intensive and computer-oriented downscaling method difficult for the individual researcher. Two climatic scenarios: one medium stabilization scenario RCP 4.5 and the other, a high emission scenario RCP 8.5, were selected to project climatic changes. The other two RCP-2.5 and RCP-6.0 can be used for future projections.

Author Contributions

S.M. designed the methodology and conceptualization in writing and editing, G.R. worked on the GCMs to select the suitable ones for the study area. M.M. worked on SDSM; S.M and N.A.-A. worked on LARS-WG; K.U. worked on statistical analysis for SDSM and LARS-WG. While M.F.U.M. worked on spatial projections, M.F.U.M. and N.A.-A. helped in data curation, visualization, writing review and editing. N.T.T.L. supervised the research, review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to acknowledge the Pakistan Meteorological Department for providing precipitation and temperature data for the study area. The authors would also like to thank the portal of the Royal Netherlands Meteorological Institute (KNMI) for freely available datasets.

Conflicts of Interest

The authors declare that they have no competing interest for the research.

References

  1. IPCC. The Physical Science Basis; IPCC: Geneva, Switzerland, 2013. [Google Scholar]
  2. Ali, S.; Kiani, R.S.; Reboita, M.S.; Dan, L.; Eum, H.I.; Cho, J.; Dairaku, K.; Khan, F.; Shreshta, M.L. Identifying hotspots cities vulnerable to climate change in Pakistan under CMIP5 climate projections. Int. J. Climatol. 2021, 41, 559–581. [Google Scholar] [CrossRef]
  3. Moazzam, M.F.U.; Vansarochana, A.; Rahman, A.U. Analysis of flood susceptibility and zonation for risk management using frequency ratio model in District Charsadda, Pakistan. Int. J. Environ. Geoinf. 2018, 5, 140–153. [Google Scholar] [CrossRef]
  4. Rahman, G.; Rahman, A.-u.; Ullah, S.; Dawood, M.; Moazzam, M.F.U.; Lee, B.G. Spatio-temporal characteristics of meteorological drought in Khyber Pakhtunkhwa, Pakistan. PLoS ONE 2021, 16, e0249718. [Google Scholar] [CrossRef]
  5. Moazzam, M.F.U.; Rahman, G.; Munawar, S.; Tariq, A.; Safdar, Q.; Lee, B.-G. Trends of Rainfall Variability and Drought Monitoring Using Standardized Precipitation Index in a Scarcely Gauged Basin of Northern Pakistan. Water 2022, 14, 1132. [Google Scholar] [CrossRef]
  6. Fenta Mekonnen, D.; Disse, M. Analyzing the future climate change of Upper Blue Nile River basin using statistical downscaling techniques. Hydrol. Earth Syst. Sci. 2018, 22, 2391–2408. [Google Scholar] [CrossRef] [Green Version]
  7. Saddique, N.; Usman, M.; Bernhofer, C. Simulating the impact of climate change on the hydrological regimes of a sparsely gauged mountainous basin, Northern Pakistan. Water 2019, 11, 2141. [Google Scholar] [CrossRef] [Green Version]
  8. Munawar, S.; Tahir, M.N.; Baig, M.H.A. Twenty-first century hydrologic and climatic changes over the scarcely gauged Jhelum river basin of Himalayan region using SDSM and RCPs. Environ. Sci. Pollut. Res 2022, 29, 11196–11208. [Google Scholar] [CrossRef]
  9. Hansen, J.W.; Mavromatis, T. Correcting low-frequency variability bias in stochastic weather generators. Agric. For. Meteorol. 2001, 109, 297–310. [Google Scholar] [CrossRef]
  10. Dahri, Z.H.; Ludwig, F.; Moors, E.; Ahmad, S.; Ahmad, B.; Shoaib, M.; Ali, I.; Iqbal, M.S.; Pomee, M.S.; Mangrio, A.G. Spatio-temporal evaluation of gridded precipitation products for the high-altitude Indus basin. Int. J. Clim. 2021, 41, 4283–4306. [Google Scholar] [CrossRef]
  11. Araya-Osses, D.; Casanueva, A.; Román-Figueroa, C.; Uribe, J.M.; Paneque, M. Climate change projections of temperature and precipitation in Chile based on statistical downscaling. Clim. Dyn. 2020, 54, 4309–4330. [Google Scholar] [CrossRef]
  12. Hamlet, A.F.; Byun, K.; Robeson, S.M.; Widhalm, M.; Baldwin, M. Impacts of climate change on the state of Indiana: Ensemble future projections based on statistical downscaling. Clim. Change 2020, 163, 1881–1895. [Google Scholar] [CrossRef] [Green Version]
  13. Baghanam, A.H.; Eslahi, M.; Sheikhbabaei, A.; Seifi, A.J. Assessing the impact of climate change over the northwest of Iran: An overview of statistical downscaling methods. Theor. Appl. Climatol. 2020, 141, 1135–1150. [Google Scholar] [CrossRef]
  14. Fowler, H.J.; Blenkinsop, S.; Tebaldi, C. Linking climate change modelling to impacts studies: Recent advances in downscaling techniques for hydrological modelling. Int. J. Climatol. 2007, 27, 1547–1578. [Google Scholar] [CrossRef]
  15. Ali, S.; Eum, H.-I.; Cho, J.; Dan, L.; Khan, F.; Dairaku, K.; Shrestha, M.L.; Hwang, S.; Nasim, W.; Khan, I.A. Assessment of climate extremes in future projections downscaled by multiple statistical downscaling methods over Pakistan. Atmos. Res. 2019, 222, 114–133. [Google Scholar] [CrossRef]
  16. Wilby, R.L.; Dawson, C.W.; Barrow, E.M. SDSM—A decision support tool for the assessment of regional climate change impacts. Environ. Model. Softw. 2002, 17, 145–157. [Google Scholar] [CrossRef]
  17. Tahir, T.; Hashim, A.; Yusof, K. Statistical downscaling of rainfall under transitional climate in Limbang River Basin by using SDSM. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Langkawi, Malaysia, 4–5 December 2017; p. 012037. [Google Scholar]
  18. Saddique, N.; Bernhofer, C.; Kronenberg, R.; Usman, M. Downscaling of CMIP5 models output by using statistical models in a data scarce mountain environment (Mangla Dam Watershed), Northern Pakistan. Asia-Pac. J. Atmos. Sci. 2019, 55, 719–735. [Google Scholar] [CrossRef]
  19. Rather, M.A.; Satish Kumar, J.; Farooq, M.; Rashid, H. Assessing the influence of watershed characteristics on soil erosion susceptibility of Jhelum basin in Kashmir Himalayas. Arab. J. Geosci. 2017, 10, 1–25. [Google Scholar] [CrossRef]
  20. Malik, K.M.; Taylor, P.A.; Szeto, K.; Khan, A.H. Characteristics of central southwest Asian water budgets and their impacts on regional climate. Atmos. Clim. Sci. 2013, 3, 259–268. [Google Scholar] [CrossRef] [Green Version]
  21. Dawood, M. Spatio-statistical analysis of temperature fluctuation using Mann–Kendall and Sen’s slope approach. Clim. Dyn. 2017, 48, 783–797. [Google Scholar]
  22. Choi, K.-S.; Oh, S.-B.; Byun, H.-R.; Kripalani, R.; Kim, D.-W. Possible linkage between East Asian summer drought and North Pacific oscillation. Theor. Appl. Climatol. 2011, 103, 81–93. [Google Scholar] [CrossRef]
  23. Kling, H.; Fuchs, M.; Paulin, M. Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios. J. Hydrol. 2012, 424, 264–277. [Google Scholar] [CrossRef]
  24. Jin, H.; Ju, Q.; Yu, Z.; Hao, J.; Gu, H.; Gu, H.; Li, W. Simulation of snowmelt runoff and sensitivity analysis in the Nyang River Basin, southeastern Qinghai-Tibetan Plateau, China. Nat. Hazards 2019, 99, 931–950. [Google Scholar] [CrossRef]
  25. Munawar, S.; Tahir, M.N.; Baig, M.H.A. Future climate change impacts on runoff of scarcely gauged Jhelum river basin using SDSM and RCPs. J. Water Clim. Change 2021, 12, 2993–3004. [Google Scholar] [CrossRef]
  26. Hassan, Z.; Shamsudin, S.; Harun, S. Application of SDSM and LARS-WG for simulating and downscaling of rainfall and temperature. Theor. Appl. Climatol. 2014, 116, 243–257. [Google Scholar] [CrossRef]
  27. Azmat, M.; Qamar, M.U.; Huggel, C.; Hussain, E. Future climate and cryosphere impacts on the hydrology of a scarcely gauged catchment on the Jhelum river basin, Northern Pakistan. Sci. Total Environ. 2018, 639, 961–976. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Geographical location of the Jhelum River basin.
Figure 1. Geographical location of the Jhelum River basin.
Atmosphere 13 00898 g001
Figure 2. Selected predictors against predictands (Tmax, Tmin and precipitation).
Figure 2. Selected predictors against predictands (Tmax, Tmin and precipitation).
Atmosphere 13 00898 g002
Figure 3. Calibration results of SDSM and LARS-WG models (1976–2005).
Figure 3. Calibration results of SDSM and LARS-WG models (1976–2005).
Atmosphere 13 00898 g003
Figure 4. Spatial projections of precipitation (mm/day): (a) CCSM4 RCP-4.5 2021–2060; (b) CCSM4 RCP-4.5 2061–2099; (c) CCSM4 RCP-8.5 2021–2060; (d) CCSM4 RCP-8.5 2061–2099; (e) HadCM3 RCP-4.5 2021–2060; (f) HadCM3 RCP-4.5 2061–2099; (g) HadCM3 RCP-8.5 2021–2060; (h) HadCM3 RCP-8.5 2061–2099.
Figure 4. Spatial projections of precipitation (mm/day): (a) CCSM4 RCP-4.5 2021–2060; (b) CCSM4 RCP-4.5 2061–2099; (c) CCSM4 RCP-8.5 2021–2060; (d) CCSM4 RCP-8.5 2061–2099; (e) HadCM3 RCP-4.5 2021–2060; (f) HadCM3 RCP-4.5 2061–2099; (g) HadCM3 RCP-8.5 2021–2060; (h) HadCM3 RCP-8.5 2061–2099.
Atmosphere 13 00898 g004
Figure 5. Spatial projections of average temperature (°C): (a) CCSM4 RCP-4.5 2021–2060; (b) CCSM4 RCP-4.5 2061–2099; (c) CCSM4 RCP-8.5 2021–2060; (d) CCSM4 RCP-8.5 2061–2099; (e) HadCM3 RCP-4.5 2021–2060; (f) HadCM3 RCP-4.5 2061–2099; (g) HadCM3 RCP-8.5 2021–2060; (h) HadCM3 RCP-8.5 2061–2099.
Figure 5. Spatial projections of average temperature (°C): (a) CCSM4 RCP-4.5 2021–2060; (b) CCSM4 RCP-4.5 2061–2099; (c) CCSM4 RCP-8.5 2021–2060; (d) CCSM4 RCP-8.5 2061–2099; (e) HadCM3 RCP-4.5 2021–2060; (f) HadCM3 RCP-4.5 2061–2099; (g) HadCM3 RCP-8.5 2021–2060; (h) HadCM3 RCP-8.5 2061–2099.
Atmosphere 13 00898 g005
Table 1. Detail of the Meteorological Stations.
Table 1. Detail of the Meteorological Stations.
Station NameStation IDLat (°N)Long (°E)Elevation (m) amslAnnual Prec (mm)Mean Teemp(°C)
JhelumA33.173.74614125522.2
KotliB33.573.891402142317.5
PlandriC33.7273.711676134616.2
RawlakotD33.8773.682213178012.8
MurreeE33.9273.38845118819.4
GaridopataF34.2273.61702138820.6
MuzaffarabadG34.3873.47996169318.3
BalakotH34.5673.34236313016.18
NaranI34.9173.64222012887
AstoreJ35.174.8221684486.56
PoonchK33.9174.03815151619.4
Table 2. Selected GCMs from CMIP5 used for downscaling techniques.
Table 2. Selected GCMs from CMIP5 used for downscaling techniques.
Modelling CentreGCMResolution
National Center for Atmospheric Research USACCSM40.9° × 1.25°
UK Meteorological Office UKHadCM32.5° × 3.75°
Geophysical Fluid Dynamics Laboratory USAGFDL-CM32° × 2.5°
National Institute for Environmental Studies JapanMRI-CGCM31.12° × 1.12°
Canadian Centre for Climate Modelling and Analysis CanadaCanESM22.79° × 2.8°
Beijing Climate Center, ChinaBCC-CSM1–12.81° × 2.81°
Table 3. Description of NCEP predictor.
Table 3. Description of NCEP predictor.
No.PredictorCodeNo.PredictorCode
1Mean sea level pressuremslp11500 hPa meridional velocityp5_v
2500 hPa relative humidityr50012Surface specific humidityShum
3850 hPa vorticityP8_z13Mean temperature at 2 mtemp
4Surface zonal velocityp_u14Surface airflow strengthp_f
5500 hPa vorticityp5_z15Surface meridional velocityp_v
6Surface vorticityp_z16Surface wind directionp_th
7500 hPa wind directionp5th17Surface divergencep_zh
8850 hPa relative humidityr85018500 hPa airflow strengthp5_f
9Surface zonal velocityp_u19500 hPa zonal velocityp5_u
10850 hPa meridional velocityp8_v20500 hPa geopotential heightp500
Table 4. Mann–Kendall test results for Tmax for Meteorological stations.
Table 4. Mann–Kendall test results for Tmax for Meteorological stations.
Station IDMann–Kendall StatisticsKendall’s TauVariance (S)p Value
(Two Tailed Test)
A380.20009500.2300
B−320.16849500.3145
C540.28429500.0855
D880.46329500.0048
E70.03699490.8456
F160.08429500.6265
G160.08429500.6265
H−44−0.15349500.2845
I240.1879500.0756
J260.0949500.5265
K330.0489500.2651
Table 5. Mann–Kendall test results for Tmin for Meteorological stations.
Table 5. Mann–Kendall test results for Tmin for Meteorological stations.
Station IDMann–Kendall StatisticsKendall’s TauVariance (S)p Value
(Two Tailed Test)
A360.18959500.2561
B440.23169500.1630
C−68−0.35799500.0297
D930.49089490.0028
E360.18959500.2561
F−20−0.10539500.5376
G270.25714080.1982
H−34−0.1978600.2131
I580.28769500.0341
J340.19519500.261
K450.31619500.130
Table 6. Mann–Kendall test results for Precipitation for Meteorological stations.
Table 6. Mann–Kendall test results for Precipitation for Meteorological stations.
Stations IDMann–Kendall StatisticsKendall’s TauVariance (S)p Value
(Two Tailed Test)
A140.07379500.6732
B−20−0.10539500.5376
C340.17899500.2843
D360.18959500.2561
E380.20009500.2300
F540.28429500.0855
G−19−0.18104080.3731
H−24−0.12454500.4786
I280.13459500.2785
J160.09739500.5732
K−13−0.1319500.3761
Table 7. Selection of GCMs.
Table 7. Selection of GCMs.
GCM ModelsPearson’s Correlation Coefficient (r)KGENSERank
TmaxTminPTmaxTminPTmaxTminP
CCSM40.820.930.790.470.490.420.780.890.790.71
GFDL-CM30.490.560.610.390.340.210.510.680.570.48
HadCM30.720.830.770.470.390.320.880.790.910.68
MRI-CGCM30.590.670.510.290.330.270.610.580.640.49
CanESM20.620.570.490.290.340.370.710.670.640.52
BCC-CSM1–10.770.810.670.490.430.330.780.890.810.66
Table 8. Performance of the models during the validation period (2006–2020).
Table 8. Performance of the models during the validation period (2006–2020).
ModelsMAERMSEBiasKGENSE
TmaxTminPTmaxTminPTmaxTminPTmaxTminPTmaxTminP
CCSM4 RCP-4.5
SDSM0.720.930.790.50.690.420.780.890.790.680.530.480.610.720.42
LARS-WG0.490.560.610.90.420.720.510.680.570.480.580.620.620.570.69
CCSM4 RCP-8.5
SDSM0.520.650.670.40.920.510.730.530.650.420.590.610.440.570.59
LARS-WG0.30.490.390.890.840.560.680.620.640.50.460.550.360.380.41
HadCM3 RCP-4.5
SDSM0.790.770.510.410.430.320.610.580.640.490.50.540.590.670.58
LARS-WG0.620.570.490.210.670.710.710.670.640.520.710.650.690.710.65
HadCM3 RCP-8.5
SDSM0.640.810.510.530.810.530.650.570.710.670.640.520.730.690.55
LARS-WG0.650.590.470.760.790.450.290.280.620.510.730.540.630.670.61
Table 9. Projected changes in mean Tmax (°C), Tmin (°C) and total precipitation (%) on annual and seasonal bases under future scenarios (RCP-4.5 and RCP 8.5) for mid-century (2021–2060) and end-century (2061–2099) in comparison to the baseline period (1976–2005) over the Jhelum River basin.
Table 9. Projected changes in mean Tmax (°C), Tmin (°C) and total precipitation (%) on annual and seasonal bases under future scenarios (RCP-4.5 and RCP 8.5) for mid-century (2021–2060) and end-century (2061–2099) in comparison to the baseline period (1976–2005) over the Jhelum River basin.
Season/PeriodRCP-4.5RCP-8.5
CCSM4HadCM3AverageCCSM4HadCM3Average
2021–20602061–20992021–20602061–20992021–20602061–20992021–20602061–20992021–20602061–20992021–20602061–2099
Tmax °C
Annual (J–D)1.362.461.222.621.292.542.344.612.414.712.374.66
Winter (D–F)1.462.591.352.761.42.672.594.82.624.712.64.75
Pre-monsoon (A–J)1.52.541.262.861.382.72.344.712.514.582.424.64
Monsoon (J–S)1.122.11.062.241.092.172.094.342.14.562.094.45
Tmin °C
Annual (J–D)1.352.541.222.521.962.532.434.532.524.512.474.52
Winter (D–F)1.482.761.362.751.422.752.664.652.724.612.694.63
Pre-monsoon (A–J)1.562.541.252.461.402.52.344.592.514.572.424.58
Monsoon (J–S)1.022.341.062.361.042.352.294.362.354.362.324.36
Precipitation (%)
Annual (J–D)6.7712.626.1712.516.4712.567.4110.587.412.57.411.54
Winter (D–F)5.4913.564.8414.95.114.236.5111.687.4514.56.9813.09
Pre-monsoon (A–J)6.5211.455.7810.896.1511.177.539.536.6510.47.099.96
Monsoon (J–S)8.312.877.8911.768.0912.38.2110.558.112.68.111.57
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Munawar, S.; Rahman, G.; Moazzam, M.F.U.; Miandad, M.; Ullah, K.; Al-Ansari, N.; Linh, N.T.T. Future Climate Projections Using SDSM and LARS-WG Downscaling Methods for CMIP5 GCMs over the Transboundary Jhelum River Basin of the Himalayas Region. Atmosphere 2022, 13, 898. https://doi.org/10.3390/atmos13060898

AMA Style

Munawar S, Rahman G, Moazzam MFU, Miandad M, Ullah K, Al-Ansari N, Linh NTT. Future Climate Projections Using SDSM and LARS-WG Downscaling Methods for CMIP5 GCMs over the Transboundary Jhelum River Basin of the Himalayas Region. Atmosphere. 2022; 13(6):898. https://doi.org/10.3390/atmos13060898

Chicago/Turabian Style

Munawar, Saira, Ghani Rahman, Muhammad Farhan Ul Moazzam, Muhammad Miandad, Kashif Ullah, Nadhir Al-Ansari, and Nguyen Thi Thuy Linh. 2022. "Future Climate Projections Using SDSM and LARS-WG Downscaling Methods for CMIP5 GCMs over the Transboundary Jhelum River Basin of the Himalayas Region" Atmosphere 13, no. 6: 898. https://doi.org/10.3390/atmos13060898

APA Style

Munawar, S., Rahman, G., Moazzam, M. F. U., Miandad, M., Ullah, K., Al-Ansari, N., & Linh, N. T. T. (2022). Future Climate Projections Using SDSM and LARS-WG Downscaling Methods for CMIP5 GCMs over the Transboundary Jhelum River Basin of the Himalayas Region. Atmosphere, 13(6), 898. https://doi.org/10.3390/atmos13060898

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