Introduction

The threat of global warming has raised the level of awareness throughout the world regarding the need to minimize the critical situation confronting all societies (Hasanov et al. 2021). Human activities on the planet’s surface are the sole contributor to global warming, resulting in environmental destruction due to global warming (Kirikkaleli & Adebayo 2020). With the advent of global warming, countries have been charged with finding individual and collaborative ways of thinking and working towards mitigating global warming. Climate change is a global problem that has strengthened the international and domestic consciousness to identify ways of mitigating the growing trend (Olanrewaju et al. 2021; Adebayo 2021a). Emissions from diverse energy sources, particularly fossil fuels and other non-renewable sources of energy, are dispersed as pollutants into the air. These are responsible for adversely affecting both the climate and the welfare of the people. Not only are the pollutants released into the environment but they also have connections to bodies of water and wetlands, which can damage or poison marine life (Kirikkaleli & Adebayo2021). The pollution of both water bodies and air has a detrimental impact on society by impacting the populace’s living conditions, health, and nutrition (Zhang et al. 2021; Udemba et al. 2021).

According to various measures, growth in the economy has been identified to have a disastrous impact as a result of pollution. Various economic practices, both directed at and based on economic growth, contribute to pollutants’ emissions (Ayobamiji & Kalmaz 2020; Udemba 2020). Such practices from multiple sectors (petroleum sector, manufacturing, oil extraction, agriculture) of the economy that cause GDP growth also trigger pollution (Umar et al. 2020). This pollution from various sources and economic sectors negatively impacts human well-being through various types of diseases such as cancer and heart disease (Adedoyin et al. 2020; Oluwajana et al. 2021; He et al. 2021). Manufacturing practices such as the use of heavy duty equipment with the potential to burn huge amounts of fossil fuels, production processes, the delivery of products from the point of origin to the last or final customer with vehicles releasing CO2 emissions via exhaust pipes, waste dumping in bodies of water, and electricity generation using other fossil fuels and coal all lead to environmental degradation.

Since 2016, South Korea’s total energy use has remained stable after a period of constant growth (2.7%/year between 2000 and 2016). South Korea’s economic performance’s driving forces are oil and coal, which account for 37% and 26% of energy needs, respectively. This shows the reliance of the South Korean economy on oil and coal as an energy source. In 2019, South Korea consumed 2.5 million barrels per day of petroleum and other liquids, ranking it the 8th largest consumer globally (EIA, 2020). This implies that, as a consequence of the rising demand for oil and coal in South Korea, CO2 emissions will certainly rise. Existing policies are projected to result in an emissions level of 665 to 743 MtCO2e/year in 2030, subject to the ongoing impact of the COVID-19 crisis, without pollutions from LULUCFFootnote 1 (EIA, 2020). In its NDC, South Korea is officially committed to achieving 539 MtCO2e/year except for LULUCF. To achieve this target, South Korea will have to improve its climate policies significantly, even more so if the government intends to achieve carbon neutrality by 2050, such as by amending and enhancing its 2030 NDC to be compliant with the Paris Accord. South Korea is ranked among the world’s top five importers of coal, liquefied natural gas (LNG), and total petroleum liquids. According to the EIA (2020), in 2019, the total primary energy consumption of South Korea consisted of coal (28%), petroleum and other liquids (43%) natural gas (16%), renewable sources (3%), and nuclear (10%), as depicted in Fig. 1. Based on the fascinating energy mix of South Korea, the present research aims to examine the effects of the positive development and growth of South Korea, with an emphasis on the energy consequences of South Korea’s economic performance via urbanization, gross capital formation, and use of energy. This research, therefore, explores the effects of CO2 pollution, urbanization, gross capital development, and use of energy on South Korea’s economic performance.

Fig. 1
figure 1

South Korea total primary energy consumption by fuel 2019

It is important to recognize that, as a foremost economy in terms of growth, it is necessary to explore the economy and render appropriate economic sustainability suggestions based on analytical results. Based on the results, a detailed analysis of the South Korean economy’s sustainability will allow us to develop sustainable policies to answer questions including (a) Can South Korea diversify its policies regarding its energy mix by embracing renewable energy to boost its green economy? (b) Can South Korea introduce innovative pollution mitigating measures without weakening sustainable economic development? It is also essential to remember that, as the South Korean economy grows, the position of the nation in the international rankings of CO2 is very vulnerable and the levels are rising at the same rate as those of other nations considered to be the major global emitters such as India, China, US, Japan, and Russia. In South Korea and Japan, which have the same economic characteristics and demographic composition, there is a disparity between the trend of output growth and the level of pollution. In both countries, this pattern has contributed to the emissions of CO2. That being said, considering the rapid population growth rate, significant efforts have been made to mitigate the detrimental consequences of global warming without impacting GDP growth. This is the inspiration for the researchers to investigate the variables illustrated in this report, as we strive to utilize the implications of this research to make policy recommendations for government administrators and stakeholders.

Given this progress, despite the optimistic and substantial growth of the South Korean economy, there has been limited emphasis on examining this trend’s significance. Nevertheless, this research intends to examine the economic performance in South Korea amidst CO2 emissions. The current research is distinct from the existing studies because it accounts for other economic growth determinants, including energy usage, CO2 emissions, gross capital formation, and urbanization. In addition to the autoregressive distributed lag (ARDL), fully modified ordinary least squares (FMOLS), dynamic ordinary least squares (DOLS), and Gradual shift tests, the study employs the wavelet coherence test. The advantage of the wavelet coherence test is that it can capture both correlation/causality between series at different frequencies and time periods. This report expands/complements the discussion on the South Korean economy on the growth-energy and pollution nexus and expands on India’s research by Udemba et al. (2021). The research is inspired by the Sustainable Development Goals (SDGs-7, 8, 12, and 13) and discusses specific energy use concerns (SDG-7) with a particular emphasis on green and sustainable use of energy (SDGs 7 and 12) to meet the 2020 Agenda. This is to avoid problems associated with economic growth (SDGs-8) and climate change (SDGs-13). The present research is considered particularly timely and deserving of inquiry, especially in the current era in which responsible energy use and environmental protection are increasingly being targeted.

The concluding part of this report is planned in the following way: a short synopsis of the previous studies and theoretical framework is presented in the “Literature review” section. The “Data and methodology” section covers the data and the methodologies, whereas the “Findings and discussion” section presents the analytical results. The “Conclusion and policy directions” section presents the conclusion and policy directions.

Literature review

This part is divided into two, namely, (i) empirical review-which discussed the previous studies conducted regarding the relationship between economic growth, gross capital formation, CO2 emissions, urbanization, and energy consumption, and (ii) the theoretical framework-which discuss the environmental Kuznets curve theory.

Empirical review

An overview of the relevant literature on this topic will be discussed by reviewing the connections observed between the dependent variable (GDP) and its regressors (CO2 emissions, urbanization, gross capital formation, and energy use). There is no consensus in the literature on the relation between GDP and these regressors as a result of the mixed outcomes, which has led to an increase in interest in this subject matter. The study of Teng et al. (2020) found that GDP increased CO2 emissions for ten different OECD economies during the period between 1985 and 2018. However, Ayobamiji and Kalmaz (2020) employed the wavelet technique to capture the time-frequency dependency between CO2 and real output, which is consistent with the results of Teng et al. (2020). The study of Ahmed et al. (2020) revealed that GDP exerts a positive impact on CO2 emissions in G7 economies. Aye and Edoja (2017) found a negative link between GDP and CO2 emission in 31 developing countries. Salahuddin et al. (2018) showed no association between CO2 and real output. In Kuwait, Wasti and Zaidi (2020) concluded that energy consumption and CO2 emission accelerate GDP. Chontanawat (2020) and Gorus and Aydin (2019) suggested that there was no causal association between CO2 and real GDP in ASEAN economies for the period from 1971 to 2015. Kirikkaleli and Adebayo (2020), Wang et al. (2019), Kirikkaleli et al. (2020), Aydoğan and Vardar (2020), and Jafari et al. (2015) revealed a one-way casual interaction from GDP to CO2 emissions. However, while Gao and Zhang’s (2021) study showed that there is a unidirectional causal link from CO2 emission to GDP, a bidirectional causal link between CO2 emission and GDP was revealed by Wu et al. (2018). The study of Ahmed et al. (2019) revealed a positive relationship between GDP and CO2 emissions. Bouznit and Pablo-Romero (2016) examined the interconnection between CO2 emissions and GDP in Algeria from 1970 to 2010 utilizing the ARDL approach. The results showed that there was a positive association between CO2 emission and GDP. The study of Odugbesan and Adebayo (2020) in South Africa also revealed a positive interaction between CO2 emissions and GDP. In Japan, the study of Ahmed et al. (2021) and Adebayo (2021b) established a positive association between GDP growth and environmental degradation. Al-Mulali (2011) established a positive connection between CO2 emissions and GDP in the MENA region. Furthermore, they found evidence of a two-way causal link between CO2 emissions and GDP. Awosusi et al. (2020) employed panel data from 1980 to 2018 for the MINT economies. The results showed that there is no significant link between CO2 emission and GDP. However, economic growth Granger causes CO2 emissions. Adebayo (2020) also employed the ARDL and wavelet coherence methods to examine the long-run and causal relationship between CO2 emissions and GDP in Mexico. The results showed a positive link between these variables. In terms of causality, they identified a two-way interaction between CO2 emissions and GDP. Zhang et al.’s (2021) study revealed a different causal interconnection between CO2 emissions and GDP in Malaysia, namely a one-way link from GDP to CO2 emissions.

The study of Khobai and Le Roux (2017) established a bidirectional link between GDP and energy usage. Muhammad (2019) examined the link between energy usage and GDP in the MENA economies from 2001 to 2017, suggesting a negative linkage between energy use and GDP. Shahbaz et al. (2018) established a positive interaction between energy consumption and GDP in the top ten energy-consuming economies utilizing the quantile-on-quantile (QQ) approach for the period 1960Q1 to 2015Q4. This is consistent with the study conducted by Magazzino (2018) in Italy. Mutascu (2016) explored the causal association between economic growth and energy use in the G7 nations. The author revealed a bidirectional link between economic growth and energy use in the USA, Canada, and Japan, but no causality was evident in the UK and Italy. Yang and Zhao (2014) examined the association between energy consumption and economic growth from 1970 to 2008 in India, utilizing the Granger causality test and DAG. The author revealed that there was a unidirectional linkage from energy consumption to GDP. Faisal et al. (2016) utilized the TY causality to examine the link between GDP and Russia’s energy consumption. The authors revealed no causal link between these two variables. Ha and Ngoc (2020) also employed Toda-Yamamoto causality on data covering the period from 1971 to 2017 in Vietnam. The authors revealed a two-way causal link between GDP and energy use. Baz et al. (2019) confirmed a positive shock moving from energy consumption to GDP in Pakistan from 1971 to 2014. Rahman et al. (2020) revealed that energy consumption positively affected China’s GDP, covering the period from 1981–2016.

The study of Nathaniel and Bekun (2020) examined the association between urbanization and GDP in Nigeria covering the period from 1971 to 2014 by employing the Bayer and Hanck cointegration test, ARDL, FMOLS, DOLS, CCR, and VECM Granger causality. They found that urbanization negatively inhibits GDP, and there is a bidirectional link between urbanization and GDP. Nguyen and Nguyen (2018) found that urbanization positively affects GDP in the ASEAN countries. Ali et al. (2020a) examined the association between urbanization and GDP using the Maki cointegration test, FMOLS, DOLS, CCR, and VECM Granger causality covering the period from 1971 to 2014. The authors found that urbanization hinders GDP in Nigeria, and there is unidirectional causality from urbanization to GDP. Zheng and Walsh (2019) concluded that urbanization is a major contributor to GDP in China. Yang et al. (2017) found a positive association between GDP and urbanization.

Numerous studies have been conducted in terms of the linkage between gross capital formation and GDP growth, although their findings are mixed. For instance, Topcu et al. (2020) explored the interaction between gross capital accumulation and GDP by using the panel vector (PVAR) for the period between 1980 and 2018 for 124 economies. The author concluded that the impact of gross capital formation differs based on the country’s income level. Etokakpan et al. (2020) examined the association between gross capital accumulation and GDP in Malaysia covering the period 1980–2014, employing the Bayer and Hanck cointegration tests, ARDL, and Granger causality. The authors concluded that an increase in gross capital formation would increase GDP. Kong et al. (2020) employed recent panel techniques to examine the relationship between gross capital formation and GDP for 39 African economies. The authors established a positive link between gross capital formation and GDP. Furthermore, a bidirectional causal link was also evident between these two variables. Boamah et al. (2018) also found similar results for 18 Asian nations by employing panel data covering the period from 1990 to 2017. Table 1 presents a synopsis of related studies.

Table 1 Synopsis of studies

Theoretical framework

This study’s theoretical work is based on the EKC which was built on the Kuznets curve of Kuznets (1955), which was centered on income inequality. This theory explains the increasing trend of inequality and income per capita. There is a turning point along the curve that shows where the farmers’ per capita income, who exit the farming practices to take up white collar jobs in urban regions is increasing, which closes the large gap between rich and poor. Environmental economists such as (Panayotou 1997; Grossman and Krueger 1991) improved this theory by examining the association between economic growth and environmental quality. The effect of GDP growth on the quality of any economy’s environment arises in 3 phases—scale effect, structural effect, and composite effects. In the first phase, environmental degradation is experienced but reaching a point (turning point), the environmental quality begins to improve due to development in innovations and increasing environmental consciousness. This first phase is termed the scale effects. This phase is related to developing nations because non-renewable energy sources promote their economic and production activities. The structural and composite effects are regarding as the turning point. This is associated with developed countries, where most of their economic activities are service and technology driven.

Data and methodology

Data

The present research explores the impact of CO2 emissions (CO2) on economic growth (GDP) and also the role of gross capital formation (GCF), energy use (EC), and urbanization (URB) in South Korea using data spanning between 1965 and 2018. In the case of South Korea, the current research was conceived with the perspective of examining the connections between GDP growth, CO2 pollution, urbanization, and energy use. The empirical modeling is based on the ARDL technique. This analysis is based on the study by Udemba (2020) and Nathaniel et al. (2020) by adjusting for further catalysts of growth that have been overlooked in the literature, including growth theory caused by the urban population. As shown in the model of Solow growth regarding capital and labor contribution. For South Korean cases that have the same economic characteristics, urban populations are included in our sample scenario. The parameters utilized are transmuted into a logarithm. This was conducted to ensure data is normally distributed (Rjoub et al. 2021; Kirikkaleli et al. 2020). Table 2 illustrates the data source, measurement, and unit of measurement. Also, the flow of analysis is depicted in Fig. 2. The study economic function and econometric model are depicted in Eqs. (1) and (2):

$$ {\mathrm{GDP}}_t=f\left({\mathrm{CO}}_{2t},{\mathrm{URB}}_t,{\mathrm{EC}}_t\ {\mathrm{GCF}}_t\right) $$
(1)
$$ {\mathrm{GDP}}_t={\vartheta}_0+{\vartheta}_1{\mathrm{CO}}_{2t}+{\vartheta}_2{\mathrm{URB}}_t+{\vartheta}_3{\mathrm{EC}}_t+{\vartheta}_3{\mathrm{GCF}}_t+{\varepsilon}_t $$
(2)
Table 2 Variable units and sources
Fig. 2
figure 2

Analysis flowchart

In Eq. (1), GDP, CO2, GCF, EC, and URB represent economic growth, CO2 emissions, gross capital formation, energy consumption, and urbanization.

Methodology

Analysis of correlation is used to verify the association of two time series data. The correlation can be defined as follows:

$$ \mathrm{Corr}\ \left(\mathrm{X},\mathrm{Y}\right)=\frac{\mathrm{Cov}\ \left(X,Y\right)}{\sqrt{\mathrm{Var}(X)\mathrm{Var}(Y)}} $$
(3)

where the covariance between the two time series (X,Y) is denoted as Cov(X,Y) while Var (X) and Var(Y) denote the value of the two time series of X and Y, respectively.

To achieve this purpose, ADF, PP, and Zivot and Andrews (ZA) unit root tests were employed to establish the order of integration. However, the ZA can capture both the stationarity property and structural break. This study employed the ARDL approach because this technique accommodates a limited number of observations (Ayobamiji and Kalmaz 2020; Kirikkaleli et al. 2018; Adebayo and Kirikkaleli 2021). Its model is suitable for a model with different lags and mixed order of integration. It is also beneficial because of the attributes of revealing coefficients in the short and long run simultaneously and solves the problem of autocorrelation. This makes the credibility of the formulated policy of this study to be effective. ARDL modeling is defined in the equation below:

$$ \varDelta \mathrm{In}{\mathrm{GDP}}_t={\alpha}_0+\sum \limits_{i=1}^t{\alpha}_1\varDelta \mathrm{In}{\mathrm{GDP}}_{t-i}+\sum \limits_{i=1}^t{\alpha}_2\varDelta \mathrm{In}\mathrm{C}{0_2}_{t-i}+\sum \limits_{i=1}^t{\alpha}_3\varDelta \mathrm{In}{\mathrm{EN}}_{t-i}+\sum \limits_{i=1}^t{\alpha}_4\varDelta \mathrm{In}{\mathrm{GCF}}_{t-i}+\sum \limits_{i=1}^t{\alpha}_5\varDelta \mathrm{In}{\mathrm{URB}}_{t-i}+{\beta}_1\mathrm{In}{\mathrm{GDP}}_{t-1}+{\beta}_2\mathrm{In}{{\mathrm{CO}}_2}_{t-1}+\kern0.5em {\beta}_3\mathrm{In}{\mathrm{EN}}_{t-1}+{\beta}_4\mathrm{In}{\mathrm{GCF}}_{t-1}+{\beta}_5{\mathrm{InURB}}_{\mathrm{t}-1}+\rho {\mathrm{ECT}}_{t-i}+{\varepsilon}_t $$
(4)

where αi=5 and βi=5 are the long- and short-run parameters respectively, ρ denotes parameter for ECTt − i, εt and Δ denote the error term and first difference respectively, and ECTt − i represents the error correction term, which is the adjusted speed to long-run balance from short-run shock. The ARDL hypotheses are written below. The null hypothesis reiterates that there is no cointegration presence in the model, while the alternate hypothesis affirms a contradictory view, which is the presence of cointegration.

$$ {H}_{\mathrm{o}}:{\alpha}_0={\alpha}_1={\alpha}_2={\alpha}_3={\alpha}_4={\alpha}_5 $$
(5)
$$ {H}_1:{\alpha}_0\ne {\alpha}_1\ne {\alpha}_2\ne {\alpha}_3\ne {\alpha}_4\ne {\alpha}_5 $$
(6)

The testing procedure in the ARDL model is by comparing the F or T statistics calculated with the critical bound (lower and upper bound). Normality test, heteroscedasticity test, Ramsey RESET, and serial correlation are the diagnostic tests undertaken to examine for BLUE. Model stability was checked employing the cumulative sum of recursive residuals (CUSUM) and cumulative sum of squares of recursive residuals (CUSUM of squares). Furthermore, the long-run coefficients of the ARDL model were verified using the FMOLS and DOLS tests.

Wavelets coherence is employed to detect the time-frequency dependence of energy consumption, CO2 emissions, urbanization, and gross capital formation on economic growth. Time-frequency dependence puts into account the changes over time and how the relationship varies from one frequency to another becomes essential and strategic in the formulation of policies (Beton Kalmaz and Adebayo 2020; Mutascu 2018; Alola and Kirikkaleli 2019; Umar et al. 2020; Alola and Kirikkaleli 2020). The Morlet wavelet function was employed since it brings balance between phase and amplitude. Morlet wavelet function is defined as follows:

$$ \mathcal{w}\left(\mathcal{n}\right)={\pi}^{-\frac{1}{4}}{e}^{-i\mathcal{wn}}{e}^{-\frac{1}{2}{\mathcal{n}}^2} $$
(7)

Note: non-dimensional frequency was used by \( \mathcal{w} \); i denotes \( \sqrt{-1} \) p(\( \mathcal{n} \)). Using the time and space, with \( \mathcal{n} \) = 0,1, 2, 3……N-1, the time series continuous wavelet transformation (CWT) is defined as:

$$ {\mathcal{w}}_{k,f}\left(\mathcal{n}\ \right)=\frac{1}{\sqrt{h}}\mathcal{w}\left(\frac{n-k}{f}\right),\kern0.75em k,f\in \mathbb{R},f\mathbf{\ne}0 $$
(8)

where k and f symbolize time and frequency, respectively. CWT helps the cross-wavelet analysis to interrelate between two variables (Kirikkaleli 2019). The CWT equation is written below as follows:

$$ {\mathfrak{w}}_p\left(k,f\right)={\int}_{-\infty}^{\infty }p(n)\frac{1}{\sqrt{f}}\mathcal{w}\left(\frac{\overline{n-k}}{f}\right) dn $$
(9)

The local variance was revealed using the wavelet power spectrum (WPS). The equation defining the WPS is as follows:

$$ {\mathrm{WPS}}_p\left(k,f\right)={\left|{W}_p\left(k,f\right)\right|}^2 $$
(10)

To examine the co-movement between two time series, the wavelet coherence approach (WTC) was used, which is defined in the equation below:

$$ {R}^2\left(k,f\right)=\frac{{\left|S\left({f}^{-1}{W}_{pj}\left(k,f\right)\right)\right|}^2}{S\left({f}^{-1}{\left|{W}_p\left(k,f\right)\right|}^2\right)S\left({f}^{-1}{\left|{W}_j\left(k,f\right)\right|}^2\right)} $$
(11)

where the smoothing operator to both time and scale with 0 ≤ R2(k,f) ≤ 1 is denoted as S. WTC can also detect the phase difference ϕpq of the two time series; it is defined in this form:

$$ {\phi}_{pq}\left(k,f\right)={\tan}^{-1}\left(\frac{L\left\{S\left({f}^{-1}{W}_{pj}\left(k,f\right)\right)\right\}}{O\left\{S\left({f}^{-1}{W}_{pj}\left(k,f\right)\right)\right\}}\right) $$
(12)

where L denotes an imaginary operator while O stands for a real part operator.

Apart from the wavelet coherence approach, the gradual shift causality test developed by Nazlioglu et al. (2016) was utilized to establish the direction of causation between two variables. Nazlioglu et al. (2016) employed the Toda and Yamamoto (1995) and Fourier approximation, which captures the structural changes during the period of coverage (Kirikkaleli and Gokmenoglu 2020; Gokmenoglu et al. 2019). This technique helps to overcome the inaccuracies and inconsistencies associated with the VAR model. Using the modified VAR model stated in the equation below:

$$ {y}_t=\upsigma \left(\mathrm{t}\right)+\kern0.5em {\beta}_1{y}_{t-1}+\dots +{\beta}_{p+d\max }{y}_{t-\left(p+d\max \right)}+{\varepsilon}_t $$
(13)

where yt symbolizes variable used; σ symbolizes intercept; β symbolizes coefficient matrices; ε symbolizes the error term; and t symbolizes time function. The Fourier approximation with cumulative frequencies is defined as:

$$ \upsigma \left(\mathrm{t}\right)={\sigma}_0+\sum \limits_{k=1}^n{\gamma}_{1k}\sin \left(\frac{2\pi kt}{T}\right)+\sum \limits_{k=1}^n{\gamma}_{2k}\cos \left(\frac{2\pi kt}{T}\right) $$
(14)

where γ2k and γ1k measures the displacement and frequency amplitude respectively and the number of frequencies is denoted as n. Fourier Toda-Yamamoto causality with cumulative frequencies (CF) is defined as follows in:

$$ {y}_t={\sigma}_0+\sum \limits_{k=1}^n{\gamma}_{1k}\sin \left(\frac{2\pi kt}{T}\right)+\sum \limits_{k=1}^n{\gamma}_{2k}\cos \left(\frac{2\pi kt}{T}\right)+{\beta}_1{y}_{t-1}+\dots +{\beta}_{p+d\max }{y}_{t-\left(p+d\max \right)}+{\varepsilon}_t $$
(15)

where approximation frequency is symbolized as k. For the Fourier Toda-Yamamoto causality with single frequencies, single-frequency components are defined in Eq. (11) as follows:

$$ \upsigma \left(\mathrm{t}\right)={\sigma}_0+{\gamma}_1\sin \left(\frac{2\pi kt}{T}\right)+{\gamma}_2\cos \left(\frac{2\pi kt}{T}\right) $$
(16)

The Fourier Toda-Yamamoto causality with single frequencies (SF) is defined as follows:

$$ {y}_t={\sigma}_0+{\gamma}_1\sin \left(\frac{2\pi kt}{T}\right)+{\gamma}_2\cos \left(\frac{2\pi kt}{T}\right)+{\beta}_1{y}_{t-1}+\dots +{\beta}_{p+d}{y}_{t-\left(p+d\right)}+{\varepsilon}_t $$
(17)

Findings and discussion

The present paper aims to examine the connection between CO2 emissions (CO2) and economic growth (GDP) as well as the role of urbanization (URB), gross capital formation (GCF), and energy usage (EC) in South Korea between 1965 and 2018. Figure 3 (correlation box) depicts the correlation among the parameters. In Fig. 3, the blue color illustrates a positive correlation, while the red color illustrates a negative correlation between variables. The correlation outcomes show that all the variables have a strong correlation with each other with the exemption of GCF with a weak correlation with other indicators. Consequently, the elementary summary statistical characteristics that report the measure of central tendencies and dispersion outlined in Table 3 show that urbanization reveals the highest average followed by economic growth, energy use, and the least gross capital formation. All series shows negative skewness with light tail, and the kurtosis revealed that data are normally distributed with the exemption of GCF. Furthermore, the study utilized the aforementioned methods to test the stationarity features of the data. The outcomes show that the parameters are integrated in mixed order, i.e., I(I) and I(0), as depicted in Tables 4 and 5.

Fig. 3
figure 3

Correlation between GDP, URB, EC, GCF, and CO2

Table 3 Descriptive statistics
Table 4 Traditional unit root tests
Table 5 ZA unit root test

After the stationarity characteristics of the series are affirmed, we proceed to estimate the ARDL framework, which is reported in Table 6. The outcomes are presented in such a manner that the existence of a long-run association between the parameters is verified if the estimated value of the F-test is greater than the values of both limits (lower and upper bound). In this scenario, we fail to accept the null hypothesis of no cointegration. Nonetheless, if the F-statistics is less than the lower bound critical value, the alternative explanation of the presence of cointegration is dismissed, although, if the F-stats breakdown between the two thresholds critical values, the finding is called inconclusive.

Table 6 Bound test

We proceed to estimate the long-run and short-run association between variables of concern, and findings are depicted in Table 7. Appropriate lag selection is essential when applying the ARDL. Thus, we utilized the AIC criteria proposed by Akaike (1987). As stated by Udemba et al. (2021) and Zhang e al. (2021), the AIC is preferred for selection lag due to its superior characteristics. The model goodness of fit is depicted by the R2 (0.99) and Adj R2 (0.98), respectively. The outcomes of the R2 and Adj R2 illustrate that 99% and 98% variation in GDP can be explained by urbanization, gross capital formation, CO2 emissions, and energy consumption, and the remaining percentage can be attributed to error. The speed of adjustment is seen to facilitate long-term convergence between the parameters with a significant and negative error correction model coefficient (ECM). The outcome of the ECT is 0.22, which illustrates evidence of cointegration among the parameters, and this signifies the capability of the model to witness 22% speed of adjustment to verify the alignment to equilibrium in the long run on GDP due to the effect of the regressors (URB, EC, CO2, and GCF). The outcomes of the ARDL illustrate the different linkage between GDP and the regressors (URB, EC, CO2, and GCF). For instance, there is proof of a negative linkage between GDP and CO2 emissions in the short run. This illustrates that an increase in CO2 emissions will mitigate GDP, which infers that an increase in CO2 by 1% in the short run is accompanied by a 0.28% decrease in GDP when other indicators are held constant. This result concurs with the findings of Lee (2013) for G-20 nations and Udemba et al. (2021) for India. Furthermore, urbanization and energy consumption exert a positive impact on economic growth in the short run, which illustrates that 1% in urbanization and energy consumption will increase economic growth by 1.92% and 0.23%, respectively. This is a strong indicator that the intensity of energy and the populace are critical indicators in the growth of the South Korean economy. The industrial sector uses a momentous quantity of energy as one of the main contributors and generators of the South Korean economy, which in turn affects income positively. This outcome aligns with the results of Zhang et al. (2021) for Malaysia and Kirikkaleli and Demet (2020) for Turkey. Moreover, there is no evidence of interconnection between gross capital formation and GDP in the short run, which aligns with the study of Olanrewaju et al. (2021) for Thailand and Zhang et al. (2021) for Malaysia.

Table 7 ARDL long-run and short-run results

In the long run, there is evidence of a long-run association between GDP and CO2 emissions, which illustrates that a 0.15% increase in GDP is linked with a 1% upsurge in CO2 emissions. This is not unexpected given that the South Korean economy is primarily an investment-oriented and manufacturing economy that relies heavily on the utilization of energy; nonetheless, a positive side of this is the potential to minimize CO2 by shifting the energy mix to include more renewable options such as wind and solar energy (renewables). Furthermore, urbanization influences economic growth positively in South Korea, which is in harmony with a positive linkage with GDP growth and urbanization. We see that a 1.92% increase in GDP growth is due to a 1% increase in the urban population. We make the claim based on the empirical revelation that the growing population in South Korea is productive to her economic trajectory. Nevertheless, there is a necessity for cautiousness on the part of policymakers to match urban infrastructure and amenities in the rural area. This is to avoid the rush to urban cities, given that most government officials develop urban areas more than rural areas. Otherwise, the urban infrastructure might be overwhelmed and might impede economic growth in the long run. We noticed that energy use improves economic growth as we noticed that a 0.23% increase in economic growth is due to an increase in energy use by a magnitude of 1%. This outcome gives credence to the energy-induce growth hypothesis, which complies with the study of Nathaniel et al. (2020) for Nigeria and also in Pakistan by Shahbaz et al. (2012). This outcome implies that the South Korean economy is energy driven and cannot embark on conservative energy strategies as such action will compromise economic growth.

Furthermore, various post-estimation tests are conducted. The outcomes of the normality, serial correlation, heteroscedasticity, and Ramsey tests show that the model is well specified, and there is no serial correlation. Moreover, the outcomes of the CUSUM and CUSUMSQ depict in Figure 4a, and 4b, respectively, which exemplifies that the model is stable.

Fig. 4
figure 4

a CUSUM. b CUSUM of squares

To confirm the ARDL long-run estimations’ outcomes, the current study utilized the FMOLS and DOLS, which are portrayed in Table 8. The outcomes show that CO2 emissions. Energy consumption and urbanization enhance economic growth in the long run, while there is no significant interconnection between gross capital formation and economic growth.

Table 8 FMOLS and DOLS outcomes

The current research further utilizes the wavelet coherence (WTC) test to catch the causality and correlation between economic growth and the regressors (CO2, EC, URB, and GCF). This method is crafted from physics to obtain formerly undetected information. Therefore, the research investigates the connection in the short, medium, and long run between GDP and its regressors. The cone of influence (COI) is the white cone where discussion is carried out in the WTC. The thick black contour illustrates a level of significance based on simulations of Monte Carlo. Figure 5a–5d 0–4, 4–8, and 8–16 illustrate short, medium, and long term, respectively. Moreover, the vertical and horizontal axis in the figure depicts frequency and time, respectively. The blue and yellow colors depict low and high dependence between the series. In-phase and out-of-phase connections are depicted by rightward and leftward arrows, respectively. Moreover, the rightward-down (leftward-up) illustrates that the first variable lead (cause) the second parameter while the rightward-up (leftward-down) depicts that second parameter lead (cause) the first parameter. Figure 5a illustrates the WTC between GCF and GDP between 1965 and 2019. At different frequencies between 1970 and 2016, the arrows are rightward-up, showing the positive interconnection between the series with GCF leading. Figure 5 b depicts the WTC between EN and GDP between 1965 and 2018 in South Korea. The majority of the arrows are rightward-up, which illustrates a strong connection (dependency) at different frequencies with energy use leading. Figure 5c illustrates the WTC between URB and GDP in South Korea between 1965 and 2019. At different scales, from 1970 to 2016, the bulk of arrows are facing rightward-up, which show positive association (dependency) at different frequencies. Furthermore, the rightward-up arrows show that URB lead (cause) GDP. Figure 5 d illustrates the WTC between URB and GDP in South Korea between 1965 and 2019. From 1995 to 2003, the bulk of arrows are facing rightward-up at high and medium frequencies, which shows a positive association (dependency) between GDP and GCF. Furthermore, the rightward-up arrows show that GCF lead (cause) GDP. The outcomes from the wavelet coherence test comply with the results of DOLS, FMOLS, and ARDL.

Fig. 5
figure 5

a WTC between GDP and CO2 emissions. b WTC between GDP and energy consumption. c WTC between GDP and urbanization. d WTC between GDP and gross capital formation

To capture the causal association among the indicators, the gradual shift test is applied, and the results of this technique are presented in Table 9. The outcomes of the gradual shift causality revealed: (i) a one-way causality from energy consumption to GDP at a significance level of 1%. This suggests that South Korea’s economic growth is not immune from development caused by fossil fuels. This is because South Korea’s economy is heavily reliant on production and manufacturing operations. This outcome corresponds to the findings of Ramakrishna and Rena (2013) and Udemba et al. (2021); (ii) a unidirectional causality from CO2 emissions to GDP growth at a 1% level of significance, which implies that CO2 emissions are a strong predictor of GDP growth in South Korea. In this regard, South Korean policymakers should develop policies in line with the nation’s energy portfolio diversification. This finding is consistent with the results of Zhang et al. (2021) for Malaysia, Adebayo (2020) for Indonesia, and Awosusi et al. (2020) for the MINT economies; (iii) a one-way causality from GDP to urbanization at the significance level of 1%, which suggests that GDP can predict URB. This finding has consequences for South Korea’s economic growth in terms of rapid urbanization, where investment will be attached mainly to the vulnerable population with the promise of improvement in welfare. This outcome concurs with the findings of Awosusi et al. (2020), Udemba et al. (2021), and Zhang et al. (2021).

Table 9 Gradual shift causality test

Conclusion and policy directions

The current study adds to the previously existing literature by assessing the linkage between economic growth, CO2 emissions, energy usage, urbanization, and gross capital formation in South Korea using yearly data stretching between 1965 and 2019. To accomplish the stated objectives, the ARDL bounds test, the gradual shift causality test, and the novel wavelet coherence test are utilized. The outcomes show a mix (significant and insignificant) of associations between economic growth and the regressors. The outcome of the bounds test reveals that all the indicators have a long-run interconnection. Furthermore, the outcomes of the ARDL long-run and short-run estimations show that energy usage, urbanization, and CO2 emissions enhance the economic performance of South Korea, while gross capital formation exerts an insignificant impact on the economic performance of South Korea. Furthermore, we applied the novel wavelet test to capture the correlation and causal association between economic growth and the regressors. The wavelet analysis findings revealed a positive connection between economic growth and the regressors with the exemption of gross capital formation, which has a weak interconnection with economic performance. Furthermore, the wavelet coherence test outcomes provide further support for the ARDL, FMOSL, and DOLS tests. The gradual shift causality test outcomes provide intuition and credibility to the linkage among economic growth and urbanization, energy usage, gross capital formation, and CO2 emissions.

This research’s outcomes have aided us in embracing the promotion of energy intensity diversification of South Korea. This could be achieved by implementing a more ambitious green energy initiative that will maintain the nation’s economic momentum. The design and execution of successful policies to regulate South Korean energy and manufacturing sector practices will improve its sustainable growth. This will continue to regulate the CO2 pollution levels in the nation if the government sets emission restrictions on firms and factories who are emitting CO2 emissions. The threat of punitive action or high taxes on infringers of this policy will deter environmental pollution. Also, energy usage should be embraced by incorporating sustainable (renewable) energy sources, including hydropower, oceanic, and wind energy sources. Furthermore, South Korea should be careful when formulating policies that will stimulate economic growth at the expense of environmental degradation. Implementing the aforementioned policies will help to maintain sustainable economic development and South Korea’s proven environmental performance. This study’s outcome could also have a positive impact on neighboring nations who are willing to take the steps suggested in this paper to strengthen their sustainable growth. In conclusion, this study has examined the nexus between energy use, urbanization, energy consumption, gross capital formation, and economic growth in South Korea using recent time series data. Further studies can be conducted for other emerging nations while considering asymmetric in the econometrics modeling or the use of micro-disaggregated data. Furthermore, other studies can account for other drivers of growth that have not been explored in this study.