1. Introduction
Most scholastic pieces of work produced by experts and researchers all agreed that greenhouse gases (GHG) are the main culprit for global warming and, as further supported by Molla et al. [
1] and Chang et al. [
2], carbon dioxide (CO
2) accounts for the vast proportion of the emissions. There is a consensus among energy researchers that climate change is being driven by increasing atmospheric greenhouse gases, most notably CO
2. The direct global GHG emissions from the transport sector were estimated to be 8.7 GtCO
2-
eq and accounted for 23% of global energy-related CO
2 emissions.
The reports of Intergovernmental Panel on Climate Change [
3] in 2019 found that 70% of the transport emissions are produced from road vehicles, while aviation, shipping and rail have a share of 12%, 11% and 1%, respectively [
4,
5]. Several gases, collectively labelled GHG, contribute to global warming processes. A decade of rapid development in South Africa’s transport sector has resulted in tremendous environmental pressures, which include energy consumption and consequently carbon emissions.
From the accounts of [
3], contributions of the various gases to global warming are CO
2 (53%), chloro-fluro carbons (CFCs) (24%), methane (CH4) (15%), and nitrous oxide (N
2O) (5%) [
6,
7,
8,
9]. It is observed that the largest contribution to the greenhouse effect is produced by CO
2 emissions. The transportation activities and primarily the combustion of fossil fuels and the resultant carbon emissions are among the activities significantly causing warming. Statistical records in recent times have shown that the ratio of increase in transport sector’s emissions in the developing regions of the world is higher than in the developed regions, specifically, North America and Europe [
10,
11,
12]. With the introduction of additional energy policies, the trend is projected to persist in the developed regions. South Africa is the second largest economy and the biggest energy user in Africa, having the most sophisticated transport system, and remains the largest GHG emitter in the continent [
13].
The CO
2 emissions in South Africa’s transport sector in 2019 was estimated to be around 60 MtCO
2-
eq. Considering the tremendous increase in the vehicle stocks with the volume of both passenger and freight transport, South Africa’s fossil fuel consumption in the transport sector has increased to around 500,000 barrels per day (tb/d) in 2022 and the import reliance will leap by 85% [
14,
15]. These are not only related to energy security and environmental pressures but some prevailing challenges as well. These challenges have drawn nationwide attention as all indicators point to addressing the concerns around CO
2-
eq emissions in the most industrially established provinces in South Africa. The three main economic buoyant provinces, namely, Gauteng, KwaZulu-Natal and Western Cape, are addressed in the paper fitting with the country’s outline plan to reduce energy consumption per unit of gross domestic product (GDP) by 9% and CO
2 emissions per unit of GDP by 8%, in comparison to the available reports of 2015 level. To address the challenges, a series of recommendations have been made both on a national and global scale, as related to fuel consumption limits, energy consumption tax and vehicle emission standards [
13,
16,
17,
18,
19].
However, many of the recommendations fundamentally dealt with energy efficiency and conversions that were deemed not to be associated with the emissions, considering varying parameters [
20,
21]. More factors are deduced to be responsible for South Africa’s transport sector’s emissions [
13]. Hence, it becomes urgently important to develop models for the three major buoyant provinces of South Africa to evaluate and analyse the varying factors responsible for carbon emissions in the transport sector.
South Africa’s GHG and CO
2 emissions relating to environmental science and in the industrial sector have been broadly investigated using an integrated approach, ARDL and wavelet coherence techniques [
22,
23] and in the transport sector covering South Africa and Lesotho [
24]. Except for Oladunni et al. [
13], the previous available studies pay attention to other sectors or GHG emissions being addressed, as combined with no specificity, in most cases neglecting the transport sector and in no way addressing Gauteng, KwaZulu-Natal and Western Cape. The study by Tongwane et al. [
24] examined only road transportation and found that GHG emissions in South Africa are closely related to the number of vehicles on the road and that it is possible to calculate their coefficient ratios. It is largely evident in the analytical methods and empirical findings of Bakker et al. [
25], Capallan-Perez et al. [
26], Alam et al. [
27], Gately et al. [
28] and Diakoilaki et al. [
29] using an integrated approach, bottom-up sectorial approach, ARDL and/or wavelet to perform GHG emissions’ reduction estimates, however, not with the case of applying econometric models.
For the desired level of economic growth, although energy policies for reducing the usage of mobility systems can serve a purpose, they cannot be proven to be sufficient for sustainable development, as the era of COVID-19 made it increasingly evident that a growing need for clean transport systems can drive behavioural changes and support the implementation of new transport technology options. Thus, it becomes imperative to stress the fact that South Africa possesses a vastly developed transport system, having peculiar characteristics in the continent with the account of the three investigated provinces, namely, Gauteng, KwaZulu-Natal and Western Cape. Although, with the use of cutting-edge energy production techniques such as hydrogen-inputs for high performance of heavy-duty gas engines, the contributions of natural gas and renewable energies to having serene residential areas and clean urban cities can be cushioned [
30,
31,
32,
33,
34,
35]. However, reliable estimation undertakings are essential in drawing functional emissions-mitigation methodologies. This research investigates CO
2-
eq emissions in the transportation sector by constructing a panel dataset of the three provinces for the period of five years, which runs from 2016 to 2020 as sourced from public, private, local and international organizations [
36,
37,
38].
At first—Apart from the industrial sector, which consumes more energy globally than any other end-use sector (more than 50%), the transport sector accounts for about 23% of global energy use [
39,
40,
41]. This burning of energy, particularly by combustion of fossils in the transportation sector of South Africa, has resulted to around 13% of CO
2 emissions calling for reduction. If addressed, this will contribute to lowering the global warming intensity.
By looking into several systematic studies on the sectorial emissions in several transportation modes, the current research extensively prioritizes CO
2-
eq emissions with the various impacting forces in the transport sector. Okada [
42] examined CO
2 emissions by factoring in population, pricing, network growth, and metropolitan metro demand. Tongwane et al. [
24] conducted research on road transportation involving two boundary countries, namely South Africa and Lesotho, determining the quantity, kind and carbon content of fuel combusted using similar criteria of [
3]. Cheze et al. [
43] conducted a study on the GHG emissions of air travel and came to some pertinent conclusions about passenger and freight transport. In studies of the transport industry in Thailand and Japan that addressed energy efficiency, some significant mitigation solutions are being derived [
44,
45].
When researching energy with focus on GHG and CO
2 emissions, there are primarily four conventional methodologies that are frequently used. The first is the index decomposition analysis, which breaks down emissions from the transportation sector into four categories: energy intensity, population size, GDP per capita, and carbon intensity. Timilsina and Shrestha [
46] investigated the rise of the transport sector’s CO
2 emissions in Italy and Asia by including fuel-mix and modal shift. Mazzarino [
47] widened the scope of his research to include transport intensity and divided CO
2 emissions in Italy essentially into five transport-type categories. The bottom-up sector-based analysis is the second approach. Using traffic data from the Highway Performance Monitoring System (HPMS) and on-road emissions estimates from the Emissions Database for Global Atmospheric Research Protocols, Gately et al. [
28] examined on-road CO
2 emissions in Massachusetts, the United States (EDGARP). Rio Gonzalez [
48] used input and output techniques created based on bottom-up sector-based to study the structure of CO
2 emissions from land transport in several European Union (UN) nations for the investigation of new energy penetration. The third method, system optimization, has been used extensively to forecast GHG emissions and energy consumption [
49,
50,
51,
52]. Gonzalez et al. [
53] analysed CO
2 emissions from passenger vehicles using panel data approach, as Georgatzi et al. [
54] proceeded to examine the determinants of CO
2 emissions caused by the transport sector in 12 European countries.
The econometric methodology is the fourth method that is highly sensitive. Using panel dataset estimation and econometric methods, Oladunni et al. [
13] investigated GHG emissions and their driving factors in the South African transportation industry. Gonzalez [
55] evaluated the effects of dieselization and induced road transport in Spain as Rentziou [
56] calculated the total number of vehicle miles traveled (VMT) for a specified period for 48 states in the continental United States. For OECD nations, the correlations between gasoline prices, car ownership, income, and gasoline consumption were studied using highly econometric methods [
57,
58]. In addition, the detection of GHG emissions from inland containers was done using econometrically derived time-series data by Dirka and Acciaro [
59].
Due to the increasing pressure from international organizations, it is, therefore, necessary to limit global warming by 1.5 °C [
59,
60]. Shafique et al. [
61] used a panel dataset from 1995 to 2007 accessible for the main Asian economies, including China, South Korea, Singapore, Japan, India, Indonesia, Malaysia, Thailand, and the Philippines, to examine the relationship between transportation, economic growth, and environmental degradation. Of a special case, using LMDI method Wang et al. [
62] studied CO
2 emissions in China’s transport sector. Using a comprehensive decomposition method, Zhou et al. [
63] and Merven et al. [
64] examined the CO
2 emissions in the Chinese transport sector and concluded that the country’s economic expansion is the primary cause of these emissions. It was established that the approach to energy consumption maintains its superiority [
64] in South Africa and [
65] in China, employing bottom-up sector-based analysis on multiple occasions to calculate emissions in the transportation sector. Xu and Lin [
66] optimized the variables affecting CO
2 emissions in the BRICS (Brazil, Russia, India, China, and South Africa) nations by using the system optimization method. Oladunni [
67] applied econometric techniques of a robustly expanded STIRPAT model from an industrial research-based constructed panel dataset to determine the effects of the various factors influencing the GHG emissions in the South African transport sector. The research report covered all the nine provinces of South Africa.
The challenges surrounding energy consumption and GHG-CO2 emissions with the identification of the driving forces continue to be of grave concern to South Africa and its continent, as more study is still needed to catch up with developed countries such as the US, the EU, and China. Previous research papers highlighted some significant flaws, so this work is putting out more improvements to the niche. It is imperative to conduct in-depth research into the three provinces that are the most economically developed, namely Gauteng, KwaZulu-Natal and Western Cape. Hence, to assess the responsiveness of the numerous factors having an impact on CO2 emissions in the transportation sector for the provinces, it is necessary to employ econometric approaches. Thereby, we can determine which important factors have the most and least influence on carbon emissions in the transportation sector.
2. Methodology
Recent energy analyses have shown that fossil fuel burning is the primary source of carbon emissions [
66,
67,
68,
69]. With the rarity of statistical releases, the carbon emissions coefficient method (CECM) is found applicable [
70] in quantifying the quantity of carbon emissions in the three provinces. The CO
2 emissions coefficients, the rate of carbon oxidation reported by the Intergovernmental Panel on Climate Change (IPCC) [
3], as provided by [
71], and the conversion factors from physical units to tons of oil equivalence are all required to calculate carbon emissions in South Africa for the years 2016 to 2020. GHG or carbon emissions from the burning of fossils fuel are generally taken in the standard of CO
2 emissions. Thus, calculated as presented below:
The CO
2-
eqt represents the total carbon emissions in the standard of CO
2 emissions of the equivalence (10
4) in year t. We have
i to denote the eight types of fossil fuels, which are crude oil, motor gasoline, natural gasoline, aviation gasoline, diesel, kerosene, lubricating oil, coal, natural gas, jet fuel and coke. We refer to
ECit as consumption of fuel type of
I (expressed in physical units’ form).
cfi represents the conversion factors from physical unit to the tons of oil equivalent of fuel type
i.
efi denotes the carbon emissions coefficients of fuel type
i (tC/tce). The ratio of
expresses the molecular weight of CO
2 to the atomic weight of carbon (C) and 0
i is the rate of carbon oxidation of fuel type of
i, as shown in Equation (1).
Table 1 presents, in tons, the oil equivalent in their conversion factors from their respective physical units [
3].
To evaluate the effects of various factors on the change in energy consumption over GHG emissions, the IPAT (I = PAT) identity has been functionally applied [
66,
67,
68,
69,
70,
71,
72]. In more general terms, the IPAT is to perform an analysis of the effects of economic activities on the environment. As initially suggested by Ehrlich and Holdren [
73], IPAT possessed a limitation of not permitting hypothesis tests, as the already known values of some specific terms determine the value of the missing term. Hence, to overcome the limitation of IPAT model, the STIRPAT model, known as stochastic impacts by regression on population, affluence, and technology, is established by Dietz and Rosa [
74]. The model of STIRPAT equation is presented as follows:
As
denotes the impacts on the environment as a result of human activities,
represents the population under study. The symbol
stands for affluence taken as per capita GDP. T represents the technological factor while their subscripts suggest that the quantities are differential across units. The constant term of
is fitted to scale the model; b, c and d are the exponential determinants of the parameters
,
and
, respectively where
serves as an error term. The STIRPAT model is formed such that more factors can be incorporated, provided that they are conceptually fitted. Wherein
By its functionality, in empirical investigation, Equation (2) may be converted to logarithmic form:
This model considers seven cogent variables to perform analysis in the transportation sector based on the three major economic provinces of South Africa. The carbon emissions are indicated as CO
2-
eq (in the equivalence of CO
2), which is equivalent to the impacts on the environment, as stated earlier. Passenger vehicles (PV) and economic growth (per capita GDP) are intentionally included to measure social economic activity. Energy intensity (EI), which measures the technological component, is taken into consideration together with energy consumption (EC), urbanization (URB) and transportation infrastructure investments. The CO
2-
eq is the variable that is parametrically dependent on the other assessed variables for the transport sector in this research paper. To foster evaluation and test of hypothesis, all the variables are used in logarithmic form, as follows:
Herein, reports individual effects for the model and β with variable-specific subscripts represent the coefficients to be estimated. As in Equation (2) and denote province and year observed, respectively.
Prior to the model’s estimation, some essential tests were performed. These are Robust Hausman test, panel cointegration test, autocorrelation test, group-wise heteroskedasticity and cross-sectional independence test. The panel cointegration test is used to assess long-run equilibrium relationship among variables, as in the case of Martinez-Zarzoso and Maruotti [
75]. Robust Hausman test is applied to give allowance for group-wise heteroskedasticity, while the autocorrelation and cross-sectional tests exercise independence with the panel dataset, as experimented by Hoechle [
76]. When the case applies, other tests determine estimation methods that are reliable when in violation of the classic hypotheses, as demonstrated by Wooldridge [
77]. The insensitivity of the estimation results from the authenticity, as observed, often emanated from the estimation error of the model parameters and bias of the covariance matrix as stipulated by Kim and Sun [
78].
5. Discussion
The outcomes of the panel assessments reveal many important scenarios for the impacts of the influencing factors and those of carbon emissions. On the surface,
Figure 2 shows a graphical depiction of the supposed normal distributions over the actual datasets for both the dependent variable (CO
2-
eq) and the exogenous variables; their impactful properties, along the period from 2016 to 2020, across the series were determined. For the plot, with regard to CO
2-
eq on the Y-coordinate, a 0.4 graphical unit represents the quantile normal distribution, while a 0.2 graphical unit on the X-coordinate represents 10,000 tons of carbon emissions. It is observed that some linear points for carbon emissions were determined across the series, prior to preceding and exceeding the graphical parametric values. Thus, for all ideal samples, equidistant linearity can be maintained. For EC on the Y-coordinate, a 0.4 graphical unit represents the quantile normal distribution, while a 0.2 graphical unit on the X-coordinate represents the kilogram of oil equivalent per capita. It is observed that some linear points for EC were determined across the series, prior to preceding and exceeding the graphical parametric values. Thus, for all ideal samples, equidistant linearity can be maintained. The case also applies for EI ideal samples’ determination. Except for the differences in the adopted units, which come with variances of units, the parametric values for GDP, PV and TII across the series were graphically determined. However, as the graph depicts for URB, the case of urbanization indicated something extraordinarily sensitive. Some linearity points were maintained where the red sloppy line crosses the parametric values of URB; however, the values, as shown therein, are not concentrated, being rather too slanty. Consequently, the carbon emissions and the impacting exogenous variables (in blue colours), along with their linearity (the red label) in the South African transportation industry, are determined. The red label, which can be taken as constant, as shown in
Figure 2, is applied to improve the predictive ability of the models. The implications in
Figure 2 indicate the ease of variables’ computation to establish a valid relationship among the parametric estimates of the initial regression and the transformed variables.
Furthermore, the blue inclining curves along the two coordinates of Y-X as shown in the graphs of
Figure 2, connote the trends for the parametric values of a given exogenous variable over the determining series. The effects exercised by each of the impacting factors (namely, GDP, EI, EC, URB, TII and PV), and the outcome of CO
2-
eq emissions on air and largely environmental quality are illustrated in
Figure 2. The red linear slopes in the graphical presentations of
Figure 2 are to depict conditional boundaries for which positive increase and or negative decrease could effect changes to the transportation environment. In addition, the proper acceleration of a determined variable to fitting the linear slope can be the required estimates at country’s developmental level across the series.
Figure 3a illustrates the relative changes of CO
2-
eq emissions developed by economic growth.
Figure 3b illustrates the relative changes of CO
2-
eq emissions caused with the results of energy intensity as
Figure 3c indicates the relative changes of CO
2-
eq emissions from the outcome of energy (oil) consumption. The impacted relative differences of CO
2-
eq emissions as forced by urbanization are graphically depicted in
Figure 3d, so as well are the relative differential impacts of transport infrastructure investments on CO
2-
eq emissions are illustrated in
Figure 3e. The effects for the impact of passenger vehicles on carbon emissions in the transport sector of South Africa are also demonstrated in
Figure 3f. The spatial differences depict the closeness value of one variable to the other.
The scattered circled-dotted plots (in blues) are presented as the CO
2-
eq emitted compounds with their relative responses to the respective exogenous impact factors. The preferably Red-falling lines for CO
2-
eq emissions of the transportation industry, or as the case connoted, the Red-rising lines as illustrated in the graphical presentations of
Figure 3 are applied to indicate the directions at which the impact factors are leading. More than any of the investigated variables, the response of passenger vehicles to the estimation models appears to be critical at the three-combined provincial level and for each of the provincial cases, namely, Gauteng, KwaZulu-Natal and the Western Cape. This conforms to the findings of Wang et al. [
98]. Secondly, we found out that the impact of passenger vehicles on carbon emissions decreases from the Gauteng to Western Cape and then to KwaZulu-Natal provinces within the 5-year study period. In the combined case of the three provinces, the results indicated that a 1% increase in passenger vehicles would increase the possibilities of carbon emissions by 0.065%. However, as proved by Habib et al. [
99], the results vary for each of the provincial cases. The notable differences are due to living standard and geographical conditions, as confirmed by He et al. [
100] and Yu and Huang [
101].
As the three studied provinces are among the most urbanized, the impacts of urbanization on CO
2-
eq emissions are not adversely significant because their level of urbanization reaches a high degree of development. Hence, the three provinces become situated places for service industry, as in the works of Li and Li [
102], and even more so in the investigations of Sun et al. [
103]. Urban residents usually have a high social awareness for the merits of low carbon burning, which results in relatively low carbon economic development, as shown by He et al. [
100] and then by Yu and Huang [
101]. As rightly debated, the panel estimation revealed that shifting towards a service economy impacts the overall decrease in the energy intensity and optimization of energy allocations, consequently mitigating carbon emissions in the transport sector of the most urban-wise developed provinces, which tallies to the works put forward by Mulder et al. [
104], followed by Adebayo and Odugbesan [
23].
For the case of the three combined provincial analysis, the GDP is shown to be negative and statistically significant, even at the confidence level of 1%. More so, in the subsamples, GDP was negative and statistically significant for the three separately investigated provinces. Therefore, we confirm the EKC relationship between economic growth and CO
2-
eq emissions in the transport sector. The impact of economic growth on CO
2-
eq emissions in the province of Gauteng is the lowest (−0.0180) among the three provincial cases. This implies that the effects of GDP on carbon emissions are very low in most developed provincial economies. It confirms that the economic growth model of the most developed province possesses the bulkiest environmental efficiency, as suggested in the works of Song and Wang [
105].
The econometric coefficients of energy intensity (EI)—starting from the combined scenario followed by Gauteng, KwaZulu-Natal and Western Cape—are −0.086, 0.549, −0.076 and −0.071, respectively. This shows that energy intensity is one of the critically significant factors that affect CO
2-
eq emissions. It also shows the reduction in energy usage from the less developed to most developed provinces. This conforms to the fact that energy efficiency can be positively correlated with economic buoyancy, as is proved by Pan et al. [
106]. More so, the results of EI, as with CO
2-
eq emissions, possess a direct proportionality constant. The coefficient of transport infrastructural investments for all the investigated cases (−0.038, −0.004, 0.045 and 0.002, respectively) indicate relatively low investments and possess negative impacts on carbon emissions; however, they are not shown to be critical.