This study analyzes nonlinear behavior links with atmospheric teleconnections between hydrologic variables and climate indices using statistical models during warm season (June to September) over the Korean Peninsula (KP). The ocean-related major climate factor, which is the El Niño-Southern Oscillation (ENSO) was used to analyze the atmospheric teleconnections by principal component analysis (PCA) and a singular spectrum analysis (SSA). The nonlinear lag time correlations between climate indices and hydrologic variables are calculated by Mutual Information (MI) technique. The nonlinear correlation coefficients (CCs) by MI were higher than linear CCs, and ENSO shows a few months of lag time correlation. The warm season hydrologic variables in KP shows a significant increasing tendency during the warm pool (WP), and the cold tongue (CT) El Niño decaying years shows a significant decreasing tendency, while the La Niña year shows slightly above normal conditions, respectively. A better understanding of the relationship between climate indices and streamflow, and their local impacts can help to prepare for the river discharge management by water managers and scientists. Furthermore, these results provide useful data for policy makers and end-users to support long-range water resources prediction and water-related policy.
The changing of global climate characteristics has induced more frequent extreme climate events, which are closely related to the seasonal variations of regional water resource variabilities (Wang et al., 2000; Pizarro and Lall, 2002; Yoon, 2015). Many climate scientists believe that extreme climate characteristics have changed with the increase in the Earth's temperature, and the frequency of occurrence and intensity of extreme climactic events is expected to gradually increase (IPCC, 2007). In particular, seasonal variations of hydro-meteorological variables are closely related to the characteristic changes in the global climate. For the efficient management and prediction of water resources, it is necessary to investigate the relationship between climate factors and hydrological indicators (Kim et al., 2012; Yoon et al., 2013).
The effect of local hydro-meteorological long-term changes do not occur independently and their correlations appear to be based on large-scale circulation pattern changes. It can be predicted a teleconnection based abnormal sea surface temperature (SST), and its nonlinear climate systems behavior in ocean-related climate indices such as the El Niño-Southern Oscillation (ENSO) in the Tropical Pacific Ocean (TPO) region (Piechota and Dracup, 1996; Piechota et al., 1998; Klein et al., 1999). Based on many studies on ENSO, which is major sources of atmospheric large-scale environmental changes, it is well known that they are closely correlated to seasonal variations of local hydrological pattern changes such as precipitation and streamflow (McPhaden et al., 2006).
The different phases of ENSO, which manifests as abnormally warm SST conditions in the central Pacific Ocean region, and its physical mechanism have been actively examined by the TPO coastal region countries scientists. These phases have been called El Niño Modoki, Central Pacific (CP) El Niño, and warm-pool (WP) El Niño etc. Several studies have shown this abnormal SST conditions will occur more frequently due to current conditions of global warming and changing climate. As a representative study, Feng et al. (2010) revealed that the El Niño Modoki phenomenon is affecting areas across of the world through significant changes of the large-scale climate and hydrological circulation patterns. The conventional climate indices Niño3 and Niño4 in the TPO region were used to identify the two types of ENSO (Ashok and Yamagata, 2009; Kug et al., 2009; Yeh et al., 2009). However, Ren and Jin (2011) argued for a distinction in the different types of El Niño, cold-tongue (CT) El Niño and WP El Niño, by using the Niño3 and Niño4 SST observation data.
In this study, we analyzed the influence of the long-term hydro-meteorological variability over the East Asia (EA) and the Korean Peninsula (KP) by looking at the conventional strong CT El Niño, different types of WP El Niño, and strong La Niña events as defined by Ren and Jin (2011). Specifically, this study had three main objectives: (1) To analyze the significant large-scale pattern changes and long-term hydrological variability in KP region by the ocean-related abnormal climate phenomenon in accordance with the ENSO; (2) to investigate linkages of the atmospheric teleconnections using a statistical correlation analysis between climatic indices and hydro-meteorological variables over KP; (3) to perform a diagnostic study on the possibility of seasonal prediction for water resources management using climate indices.
The five major river basins (Han River basin, Nakdong River basin, Geum
River basin, Seomjin River basin, and Youngsan River basin) and their
sub-watersheds are located in the Korean peninsula (126
ENSO data used in this research was obtained from the National Oceanic and
Atmospheric Administration (NOAA), which provides monthly ENSO data from the
four sites (Niño1
The 1966 to 2012 precipitation and streamflow data was obtained from WAMIS in Korea. The average precipitation was calculated by the Thiessen polygon network method from the 125 precipitation gauge stations in the five major river basins of the KP. The streamflow data was calculated using the Precipitation-Runoff Modeling System (PRMS), a widely applied semi-distributed long-term runoff simulation model developed by the US Geological Survey (Leavesley et al., 1983; Dressler et al., 2006). By using PRMS to calculate the natural streamflow, the effect of artificial dams was eliminated from the streamflow data.
ENSO phenomenon in the TPO region induced large-scale atmospheric climate pattern changes. To understand their local hydro-environmental impacts in the KP, it is necessary to analyze the linear and non-linear correlations and its lag time correlations between climate indices and hydro-meteorological variables. In general, abnormal SST in the TPO region has triggering effects for the troposphere temperature rise. A teleconnection, enabled by the jet stream, connects the abnormal SST variation in these subtropical regions to hydrological variations worldwide.
To explore this, we performed a teleconnection based time series analysis and analyzed its hydrological local impacts in the KP region, shown in Fig. 1. First, we analyzed stationary processes by the logarithmic transformation in each time series data set, such as ocean-related climate indices and hydrological data. Second, we analyzed the low frequency for the time series data sets by PCA using a SSA. Third, the joint probability distribution density function was derived from the SSA results using a kernel density method. The linear and non-linear lag time correlation of the joint probability distribution density function, derived using the MI technique, was used to analyze the atmospheric teleconnections. Finally, we assessed the impact of this teleconnection on hydro-meteorological variability and the resulting implications for local water resource management.
The emerging issue of the new type of El Niño phenomenon has caused
significant changes in the hydro-meteorological patterns throughout the
world. However, quantitative studies on features of CT/WP El Niño were
found to be relatively insufficient for the Korean watersheds. The effect of
different phases of ENSO on the summer precipitation and streamflow
characteristics was analyzed by adopting the classification criteria
proposed by Ren and Jin (2011) for hydrologic variables from the five major
river basins in Korea and their 113 sub-basins. We calculated these
variables and divided into the two different indices (
Singular spectrum analysis (SSA) is a non-parametric spectral estimation
method that reduces dispersion by changing the coordinates in time-series in
derived techniques from principal component analysis (PCA), which enables
the extraction of information from noisy time series. By removing the
non-harmonic components from the original time series data, one can more
easily understand its long-term frequency and trend (Moon and Lall, 1996).
SSA embeds a time series data
Flow chart of teleconnection based seasonal prediction study between climate indices and hydro-meteorological variables.
To separate the non-harmonic components, the size of the eigenvalues is
defined by the orthogonal process. The orthogonal function is calculated
between
Mutual information (MI) is a quantitative measurement of how much one random
variable (
If we have two kinds of time series data sets such as (
In this study, we performed a linear regression (LR) approach with the estimated average MI value and a nonlinear regression using the method of Kernel density estimation (KDE). The 90 % confidence bounds were estimated using 1000 bootstrap resampling replications, allowing for a more accurate calculation of the confidence bounds given the limited data.
A comparison of Niño3 and Niño4 anomalies using a scatter plot
revealed a correlation of 0.77 (
Composite sea surface temperature anomalies (SSTA) in the
TPO region in developing phases of two types of El Niño and La Niña
during December–January.
Through this classification research, we confirmed that the WP type El Niño is occurring more frequently than the CT type El Niño in the Tropical Pacific Ocean in recent years. From the changes in frequency of WP to CT, we can see that using the subdivision proposed by Ren and Jin is necessary for more accurate analysis of the impacts of these major climatic events.
Figure 2 shows the composite sea surface temperature anomalies of the mean SST in the TPO region in developing phases of two types of El Niño and La Niña during December–January. The SST data was obtained from the National Centers for Environmental Prediction (NCEP) and the National Center for Atmospheric Research (NCAR) reanalysis version 2 and the 30-year climatology for the period of 1981–2010. In the case of CT El Niño years, the above-normal anomalies occurred in the central Pacific and eastern Pacific regions while below-normal anomalies occurred throughout the western Pacific region (Fig. 2a). This warm anomaly grew stronger and extended to the west during the fall. In the winter, a positive SST anomaly pattern strengthened, but as spring arrived the CT El Niño weakened and disappeared. However, the case of WP El Niño pattern shows a slight yet clear difference when compared to the CT El Niño. In the WP El Niño years, above-normal SST anomalies clearly appeared in the central Pacific region and a slightly above-normal SST anomaly was observed in the eastern Pacific regions (Fig. 2b). This central Pacific warm anomaly grew stronger and extended to the east during the fall. In the wintertime, a positive SST anomaly pattern strengthened, but as spring arrived, the WP El Niño weakened and disappeared. In the case of La Niña years, a cold SST anomaly pattern was located in the eastern Pacific and central Pacific Ocean (Fig. 2c). This cold anomaly grew stronger and extended to the west during the fall. In the wintertime a negative SST anomaly pattern strengthened but, as spring arrived, the La Niña weakened and disappeared. Many recent studies (Kug et al., 2009; Feng et al., 2010) show these different SST anomaly patterns may affect air circulation and they are findings significant changes in different regions in seasonal precipitation over the Pacific region countries.
Figure 3 shows the composite anomalies (1981–2010 climatology) of Global
Precipitation Climatology Center (GPCC) precipitation from June to September
during strong ENSO events over the North East Asia domain. We extracted
monthly precipitation for 1901–2013 from the GPCC full quality control data
reanalysis version 6.0 at
Composite anomalies (1981–2010 climatology) of GPCC Full V6 combined precipitation (JJAS: June to September) during ENSO events over the North East Asia domain.
Data for two different CT/WP El Niño years, and La Niña years allowed for a characteristic analysis of GPCC rainfall over the Korean Peninsula. Figure 3a to c shows the regional scale variability of Korean summer (June–September) precipitation during CT/WP El Niño and La Niña years. During CT El Niño years (Fig. 3a), the GPCC precipitation shows anomalously dry conditions (10–50 mm lower than long-term normal) on the KP, while WP El Niño years show anomalously wet conditions (20–50 mm higher than long-term normal) (Fig. 3b). As well, La Niña years showed slightly wet conditions: 0–10 mm higher than the long-term normal. While the results of this study rely on short observations and relatively limited data set, these findings are important to seasonal prediction for streamflow in basin scale and establishing adaptive solutions to changing climate related issues.
Figure 4 shows nonlinear and linear correlation coefficients (CCs) with their
90 % confidence bounds between climate index and hydro-meteorological
variable in the KP, using KDE and LR approaches respectively. As described
in the methodology, we obtained the non-linear behaviors and their
correlations between climate indices and hydro-meteorological variable from
a MI analysis. The linear and non-linear correlation of the
Nonlinear and linear CCs with their 90 % confidence bounds between climate index and hydro-meteorological variable in the Korean Peninsula, using KDE and LR approaches, respectively. The confidence bounds are given by 5 and 95 % quantiles of 1000 bootstrap replications.
Composite anomalies of annual and JJAS (June to September) rainfall during CT/WP El Niño and La Niña years. The hatched polygons indicate statistically significant changes in annual rainfall based on the 10 % significance level.
Same as Fig. 5, but for streamflow.
Scatter plot of percentage change anomaly (departures from
the 1971–2000 normal) changes between precipitation and streamflow in
different climate conditions.
Figures 5 and 6 show composite anomalies of annual and JJAS (June to September)
rainfall and river discharge during CT/WP El Niño and La Niña years
in the basins of the KP. The hatched polygons indicate statistically
significant changes in annual rainfall based on the 10 % significance
level by Student's
During CT El Niño events (Fig. 5a, d), the annual/JJAS rainfall on the
KP was
As a result of this composite analysis, a substantial increasing tendency of annual/JJAS rainfall and streamflow was found during WP El Niño events, while a significant decreasing trend was found during CT El Niño events, as compared to long-term normal years. In addition, the annual/JJAS rainfall and streamflow during La Niña event shows a slight tendency above normal.
Figure 7 shows a scatter plot of percentage change anomaly (departures from
the 1971–2000 normal) changes between precipitation and streamflow during
annual and JJAS season in different climate conditions of the ENSO events.
ENSO conditions clearly show that the differences during the WP El Niño
year (
Throughout the results, different phase of ENSO events clearly show influence on hydrologic seasonal variability in the KP and its sub-watersheds. This shows a linkage of possible teleconnections and characteristic changes between major climactic events and hydrological indicators.
The mechanism of major climate phenomena such as El Niño and La Niña has not yet been fully understood. It is currently impossible to completely simulate these large-scale climate phenomena using the physical climate model, which would require perfect knowledge of the global environment and its weather conditions. Furthermore, analyzing climatic phenomena through the physical model is only possible in limited geographical areas, as predicting involves considerable difficulties. With these considerations, this study analyzed non-linear behavior links of the atmospheric teleconnections between hydrologic variable and climate indices using statistical models. One disadvantage of the statistical approach is that the limited number of observation reduces the significance level of the results. This is why the statistical model is utilized to complement to the prediction result of the physical model, resulting in more significant results.
Results from composite analysis in JJAS seasonal rainfall and river discharge during the CT/WP El Niño and the La Niña events in the KP and its sub-watershed show that the WP El Niño years show a substantially increasing tendency than normal years. The CT El Niño decaying years show a significantly decreasing tendency compared to the normal years. In addition, the La Niña years shows a slightly increasing tendency and normal condition. Also the MI based nonlinear and linear teleconnections using KDE and LR between climate indices and hydrologic variables show that the nonlinear CCs were higher than linear correlations, and the ENSO has a few months lag time correlation on rainfall and streamflow anomalies over the KP region.
This result shows a methodology of the teleconnection based long-range water resources seasonal prediction for reducing a changing climate uncertainty, when abnormal SSTA (such as, El Niño, La Niña) phenomenon occurs in the TPO regions. Furthermore, these results are useful water managers and end-users to support long-range water resources prediction and water-related policy making. A better understanding on hydrologic natural and anthropogenic drivers of this relationship between ENSO and streamflow could help to scientists to prepare river discharge control in the major river basins in KP. In addition, this research will facilitate well-informed decision-making in minimizing climate-related disasters such as flood and drought on agriculture sector, and supports optimal development of strategies ensuring use of water resources in best management under a changing climate.
The data set of precipitation and streamflow over Korean Peninsula are
respectively available at the Korean government's Water Resources Management
Information System (WAMIS),
This research was supported by the APEC Climate Center and we appreciate the reviewers for their valuable comments of this article.