Chaos, Fractionality, Nonlinear Contagion, and Causality Dynamics of the Metaverse, Energy Consumption, and Environmental Pollution: Markov-Switching Generalized Autoregressive Conditional Heteroskedasticity Copula and Causality Methods

: Metaverse (MV) technology introduces new tools for users each day. MV companies have a significant share in the total stock markets today, and their size is increasing. However, MV technologies are questioned as to whether they contribute to environmental pollution with their increasing energy consumption (EC). This study explores complex nonlinear contagion with tail dependence and causality between MV stocks, EC, and environmental pollution proxied with carbon dioxide emissions (CO 2 ) with a decade-long daily dataset covering 18 May 2012–16 March 2023. The Mandelbrot–Wallis and Lo’s rescaled range (R/S) tests confirm long-term dependence and fractionality, and the largest Lyapunov exponents, Shannon and Havrda, Charv â t, and Tsallis (HCT) entropy tests followed by the Kolmogorov–Sinai (KS) complexity measure confirm chaos, entropy, and complexity. The Brock, Dechert, and Scheinkman (BDS) test of independence test confirms nonlinearity, and White‘s test of heteroskedasticity of nonlinear forms and Engle’s autoregressive conditional heteroskedasticity test confirm heteroskedasticity, in addition to fractionality and chaos. In modeling, the marginal distributions are modeled with Markov-Switching Generalized Autore-gressive Conditional Heteroskedasticity Copula (MS-GARCH–Copula) processes with two regimes for low and high volatility and asymmetric tail dependence between MV, EC, and CO 2 in all regimes. The findings indicate relatively higher contagion with larger copula parameters in high-volatility regimes. Nonlinear causality is modeled under regime-switching heteroskedasticity, and the results indicate unidirectional causality from MV to EC, from MV to CO 2, and from EC to CO 2 , in addition to bidirectional causality among MV and EC, which amplifies the effects on air pollution. The findings of this paper offer vital insights into the MV, EC, and CO 2 nexus under chaos, fractionality, and nonlinearity. Important policy recommendations are generated.


Introduction
Natural environmental systems can exhibit chaotic behavior due to the complexity of interactions between various factors [1].Human activities, including those related to the development of digital technologies, can introduce disturbances that affect the balance of these systems.Still, this topic has found very limited applications in recent environmental studies.Several studies are currently exploring the potentially chaotic nature of environmental pollution.In 1963, Lorenz introduced a valuable tool to depict the dynamic characteristics of such systems [2].The atmospheric system exhibits a complex and nonlinear structure, encompassing both stochastic and deterministic elements.Chaotic structures in the average ozone concentration within the atmospheric system are empirically shown in the literature [3,4].Further, [5] demonstrated nonlinearity in NO 2 and CO using the Volterra-Wiener-Korenberg (VWK) method.Meanwhile, [6] uncovered evidence of chaotic structures in NO, CO, SO 2 , PM10, and NO 2 with the utilization of the Lyapunov exponent.In addition, [7] detected chaotic characteristics in the air pollution index (API) in Lanzhou, northwest China; [8] explored chaotic dynamics of air pollution in India, China, and Türkiye; and [9] examined chaotic dynamics of air temperature with Lie group algebra.While these studies explored the factors affecting environmental pollution and the effects of metaverse and energy consumption, their chaotic structure with longterm dependence leading to nonlinearity in contagion and causal relations have not been explored in the literature.
The metaverse (MV) has an important potential to transform our daily lives in the near future, and MV technologies will also bring new business opportunities.In fact, by connecting virtual and physical spaces with virtual reality (VR), augmented reality (AR), and internet technologies (IT), the MV is expected to host a growing number of activities of people, which include working, shopping, attending schools, games, and entertainment.According to a survey conducted by Gartner, a leading information technology research and consultancy company, it is projected that by 2026, 25% of the global population will spend at least one hour per day in the MV, engaging in activities such as work, shopping, education, socializing, and entertainment [10].
Digital technologies, especially the MV, could have significant impacts on the chaotic nature of environmental pollution.It is very important to study the chaotic nature of both digital technology and environmental pollution under the influences of human behavior and new digital technologies.The chaotic structure relationship between energy prices, the metaverse, and environmental pollution is not inherently direct, but one can explore potential connections and interactions among these concepts.It is important to note that the metaverse is a concept that primarily refers to a virtual shared space created by the convergence of physical and virtual reality, including the sum of all virtual worlds, augmented reality, and the internet.The infrastructure required for the metaverse, including data centers, can contribute to environmental concerns if not managed sustainably.Energy consumption and electronic waste are potential issues.The development and maintenance of the infrastructure supporting the metaverse, such as servers and data centers, can contribute to pollution.
The massive energy requirement of MV technology and its significant environmental implications have been long foregone [11].A projection from Global Bank Citi shows an even more serious situation: the number of MV users could reach up to 5 billion by 2030, and the MV economy could grow to a value of USD 8 to USD 13 trillion [12].
Establishing the MV requires both physical infrastructure and a reliable supply of electricity.Essential components of the infrastructure include hardware devices like VR headsets, which are crucial for delivering immersive virtual experiences to users.Moreover, upgrades to communication networks and data centers are necessary to enhance bandwidth, computing power, and speed, enabling the accommodation of a larger number of users and meeting their increased online activity requirements [13].
The impacts of the MV on energy consumption (EC) and environmental pollution have raised many concerns.As the MV expands and becomes more integrated into our daily lives, it will have various effects on energy and the environment.Regarding the energy problem, the operations of the MV necessitate a significant amount of energy required for servers, data centers, and devices such as the 3D, VR, and AR equipment involved in creating and maintaining the virtual environment [14].Since no data exist specifically for the MV's energy consumption in financial transactions, Bitcoin could be used instead for such analysis.It was shown that the mining of a single Bitcoin in 2021 led to the level of electricity consumption of a typical household over 12 years [15].The yearly Bitcoin mining in the world sums to 200 TWh terawatt-hours, equivalent to the energy consumption of Thailand in a year [16].For 2019, the energy consumption of Bitcoin mining was estimated to be significantly higher than the total energy consumption of Austria [17].Further, Ref. [18] proposes a method to calculate the energy consumption of cryptocurrencies, and Digiconomist estimates Bitcoin's peak yearly consumption as 204 TWh, which is expected to result in 95 million tons of CO 2 emissions per year [19].At the same time, artificial intelligence (AI) technology has the potential to generate and populate virtual environments within the MV, leading to an extensive amount of energy demand.Some examples have been given for recent popular AI systems such as the OpenAI's DALL-E.In those cases, these technologies are shown to require substantial computational power for both training the models and operating the environment.As a result, their environmental effects should be seriously considered.A recent study points out the inefficiency of the electricity grid and distribution for the proposed MV; however, the MV could also have a positive effect in the future if the technology is used to increase distributional efficiency with an effective matching of energy supply and demand [20].
Though this is the case, the investigation of the MV, energy consumption, and carbon dioxide (CO 2 ) emissions nexus is of vital importance to be measured empirically.To address these environmental concerns, it is crucial to implement sustainable practices and policies within the MV ecosystem.The high energy demand of the MV leads to increased greenhouse gas (GHG) emissions [14], which includes CO 2 emissions, a major contributor to the sera effect.The interconnectedness of the MV and the activities conducted are criticized for creating a significant carbon footprint [11].Carbon footprint encompasses the emissions associated with the infrastructure, data storage, data transmission, and user activities [14].In addition to the environmental effects of EC, further effects can be attributed to electronic waste generation (e-waste) and carbon the footprint associated with MV technologies and activities.
Although the impact of the metaverse on high energy demand and environmental pollution is emphasized by the above studies, these studies did not analyze chaotic structure.Existing explanations for analyzing the chaotic nature of environmental pollution remain incomplete without a comprehensive investigation of chaotic structures using the largest Lyapunov (λ) and Shannon entropy tests and without examining the long-term dependence and fractionality in environmental pollution.As we emphasized above, the chaotic structure is an important phenomenon of these series, and the analyzed variables have complex associations that could be altered even after a small shock is applied to their dynamics, leading to changes in their structures.To our knowledge, no empirical analyses exist for the assessment of contagion, tail dependence and causality dynamics, or possibly chaotic and fractional dynamics coupled with nonlinear relations between the MV, EC, and CO 2 .
With a dataset that consists of daily observations for CO 2 emissions, energy consumption by the major cryptocurrency, and the MV stock index for a period covering 18 May 2012-16 March 2023, this study explores chaotic structure by using Lyapunov exponent tests [21] and Shannon entropy [22], long-term dependence and fractionality with rescaled range (R/S) and Lo's modified R/S [23] tests, Geweke and Porter-Hudak's and Robinson and Henry's tests [24,25] for long memory, and estimation of fractional parameters.Afterward, the series and the relations among these series are modeled with nonlinear models and nonlinear causality tests aiming at capturing contagion, tail dependence, and causality.
Evidence of chaos is investigated by using the largest Lyapunov exponent (λ) and Shannon entropy [21,22].λ is used to determine the sensitivity to primal conditions and identify the presence of chaotic motion.Hence, the detection of chaos in our study is conducted following a set of studies with this respect [21,[26][27][28][29][30][31].Further, the series is analyzed with the HCT entropy measure first developed by Havrda and Charvât [32] and afterward by Tsallis [33].In addition, the series is analyzed with the Kolmogorov-Sinai (KS) complexity measure to examine different aspects of the dynamics for their dependence on initial conditions, chaos, complexity, and unpredictability [34].In the second stage, longterm dependence is tested with the Malderbrot and Wallis R/S test [35] and Lo's modified R/S test [23].Long-memory and fractional parameters are further tested and estimated with Geweke and Porter-Hudak's [24] and Robinson and Henry's methods [25].The series are examined for nonlinearity with the BDS test introduced by [36,37].Heteroskedasticity is tested with two methods.The first is White's Lagrange multiplier (LM) test [38], which is a robust test used to examine the existence of nonlinear forms of heteroskedasticity.The second is Engle's ARCH-LM test, which follows the LM testing procedure to examine the non-existence of autoregressive conditional heteroskedasticity (ARCH) processes [39].
If the methods at the previous stages point at a chaotic structure, complexity, nonlinearity and fractionality, we use the Markov-Switching Generalized Autoregressive Conditional Heteroskedasticity Copula Causality (MS-GARCH-Copula Causality) method, a novel hybrid method benefiting from Markov-Switching GARCH models in marginal distributions, copula modeling for joint distributions for regime-dependent tail dependence and contagion, and regime-dependent causal inference as an extension of Kim's GARCH-Copula Causality to regime-dependent nonlinear causality tests [40].The approach has found recent applications for a nonlinear inquiry of regime-dependent contagion, tail dependence, and causality dynamics [41,42].As it will be discussed, the MS-GARCH-Copula Causality method allows for the determination of ARCH and GARCH effects in different regimes in addition to persistence, contagion, and Granger causality in the presence of highly complex dynamics.
This manuscript consists of five sections.Section 1 covers the introduction and literature.Section 2 focuses on the methodology, Section 3 provides the dataset and empirical results, and Section 4 includes the conclusion and policy implications.

Methodology
This study uses the Markov-Switching (MS) Autoregressive Moving Average (ARMA) Generalized Autoregressive Conditional Heteroskedasticity (GARCH) Copula Causality, hence forth, the MS-GARCH-Copula Causality method [40].This new method generalizes Kim's [43] causality tests to copula and GARCH with Markov regime switching.A time series r t is assumed to follow a nonlinear conditional mean, as in Equation (1), and nonlinear conditional variance, as in Equation (2) [44,45], where In Equation ( 2), non-negativity constrained parameters ϕ, α, κ, β > 0 allow for strictly positive variance processes similar to the single-regime GARCH model of [46].Switching between distinct regimes is governed with a process of Markov-switching as As shown by [44,45], there are two methods to describe σ t,(st) and ε 2 t−1,(st) , i.e., Henneke et al.'s method or Francq and Zakoian's to obtain stationarity [47,48].In this study, in the estimation of the model, we follow the second method as used in a set of studies [40,42].
The model given in Equations ( 1)-( 5) assumes MS-GARCH processes for marginal distributions.Their joint distributions are modeled with regime-dependent copulae based on Clayton, Student's t, and Gumbel copulae.With the use of Sklar's extended theorem, single-regime copula modeling could be extended to allow nonlinear and regime-dependence; therefore, dependence characteristics vary for each observation [40].The first copula is the Student's t copula, which is written as which allows for the modeling of a symmetric tail dependence by assuming λ LT = λ UT , where the copula parameters for the lower and upper tails are Further, the last two copula functions are asymmetric in nature.Assuming a copula function of the form below the Gumbel copula is an asymmetric copula function with and θ ∈ {1, +∞} is the positively defined power coefficient.If instead, the following function is maintained one obtains the second asymmetric copula function, the Clayton copula, for which the power term is defined as θ ∈ {0, +∞}; for the upper tail coefficient, λ U = 0, and for he lower tail, λ L = 2 −θ −1 .The log-likelihood (LL) function is derived for maximization The regime-switching probability of each regime is conditional on γ t−1 , the available information set at period t−1, so that P(s t = j|γ t−1 ) and P(s t = j|γ t ).Following [49], Kim's filter is applied to recover the probabilities as and The transition and the smoothed transition probabilities are attained from and P s t=j γ t = P s t=j γ t ∑ n i=1 where the dependence follows a time-varying process The copula-based Granger causality is a newly suggested method to investigate causality between analyzed series [43], where Kim generalizes Granger causality [50] to copulae.In this paper, the approach is extended to Markov regime switching both in the marginal and joint distributions.
Following the method, marginals are modeled by MS-ARMA-GARCH processes, and joint distributions with copula are extended to nonlinear copula Granger causality [40].The method generalizes GARCH-Copula Causality to regime-switching GARCH-Copula Causality, which helps the researcher as it is an effective method to provide causality direction, if it exists, within each regime.Other nonlinear extensions of copula to causality, readers are referred to [51,52].
Let A = a t , B = b t , and assume A ↛ B denotes no causality from A to B and A → B reflects the opposite: causality in the direction from A to B. Hence, the null hypothesis of no Granger causality is tested using Equation ( 16) below where b n t = (b t , . . . ,b t−n+1 ) and a m t = (a t , . . . ,a t−m+1 ) show past information regarding B and A with orders n and m, and f is the conditional probability of density equality, which reflects equality of prediction performance of the bivariate setting to the univariate case.Hence, to test the null hypothesis of no Granger causality from A to B, H 0 : A↛B, against the alternative of Granger-causality (GC) running from A to B, H 1 : A → B, the likelihood ratio test statistic is obtained It is known that copula at high dimensions could be recursively reduced towards a set in low dimensions.The resulting representation eases implementation in addition to efficiency during the estimation of GC with the utilization of conditional copula.Hence, the GC representation could be restated as and in Equation ( 18), h represents the conditional joint density of (A, B) and f and g denote the marginal density of B and A. Hence, the conditional joint density is where c is the copula density function, u = F(b t+1 |b n t ) and v = G(a m t |b n t ), and F and G represent the conditional marginal distributions of B and A. By substituting Equation (19) into Equation (18), copula-based GC causality from A to B (CGC A→B ) is attained from The suggested approach integrates nonlinearity in light of regime dependence.The copula-based methods are shown to have good properties for cases with variables and residuals that deviate from the Gaussian distribution assumption [53].Kim et al. suggested the statistical test and statistical properties of copula nonlinear causality [43].Their method leads to copula directional dependence and copula nonlinear Granger causality tests [43].The MS-GARCH-Copula Causality method simultaneously helps in the determination of the regime-dependent direction of causality, and it permits the researcher to study nonlinear contagion effects and the conditional tail dependence between the analyzed variables [40].

Data
This paper uses a dataset that consists of daily observations for CO 2 emissions, energy consumption, and metaverse stock index variables for a sample covering 18 May 2012-16 March 2023.CO 2 emissions are denoted as CO 2 and represent daily atmospheric CO 2 emissions PPM at the global level obtained from the CO 2 Earth Database.The CO 2 Earth Database utilizes independent CO 2 monitoring from various observatories at different geographies at different latitudes from four core observatories with two independent readings based on NOAA and Scripps methods for confirmatory purposes.
The metaverse stock index (MV) is produced from Roundhill's metaverse index.A detailed explanation of the data selection is below.The unavailability of MV data or its availability affected the quantitative research in the empirical literature, especially those studies aiming at the utilization of long time series datasets.As a result, the applications are limited.The same also holds for EC data resulting specifically from the MV.For the former, we focused on MV company stock market capitalization, and we assumed the MV stock market fluctuations under the influence of the developments regarding the MV.In the Eikon Reuters Database, two main MV stock market indices are Roundhill's and Solactive's MV indices.Both indices are calculated from metaverse companies' stock prices weighted by the relative size of each MV company's market capitalization in U.S. dollars.Roundhill includes 52 MV companies and Solactive includes 30.With a larger coverage, the former also is a shorter time series.To overcome this problem, we focused on Roundhill's index; however, we reduced the number of MV companies to 36, still covering more companies than Solactive's, while achieving a decade-long dataset.This longer MV index was checked for distributional characteristics by comparison to Roundhill's index, and both series closely resembled each other.Pearson's rho is 0.999, and independence tests indicated no difference in the variance processes in addition to similar skewness and kurtosis statistics.
To overcome the unavailability of energy consumption data in daily frequency, following the literature, in Section 1, which links cryptocurrencies and the MV, we based our analysis on the energy consumption (EC) of Bitcoin derived from the Cambridge Bitcoin Energy Consumption Index, published by the University of Cambridge, Centre for Alternative Finance.EC data is given in terawatt-hours (Twh).All variables are subject to natural logarithms and are first differenced.Therefore, the series analyzed represents daily % changes in the analysis.
The descriptive statistics for the variables are reported in Table 1.The JB test results indicate the non-normality of all series.Excess kurtosis exists for all series, which indicates a leptokurtic distribution with heavy tails, generally considered a sign of heteroskedasticity.

Results
The results were obtained in three stages: descriptive results, cointegration results, and causality results.Three stages are given as follows. 1.
Long-term dependence and fractionality.Fractionality and long-term dependence are examined for the dataset with Malderbrot and Wallis' R/S and Lo's modified R/S tests.Long-memory and fractional parameters are estimated, and further tests for statistical significance are conducted with the Geweke and Porter-Hudak [24] and Robinson-Henry methods [25].

2.
Determination of chaos, entropy, and complexity.To determine chaotic dynamics, the largest Lyapunov (λ) [21] exponents are calculated.The λ results provide insight into the average exponential rate of divergence or convergence of nearby orbits in the phase space.A positive Lyapunov exponent (λ > 0) indicates chaos, signifying the divergence of nearby trajectories.Conversely, when λ ≤ 0, the time series corresponds to regular motion.A λ value equal to zero indicates a bifurcation.A λ value greater than 1 signifies a deterministic chaotic process.If λ ≤ 0 and λ is greater than 1, we do not explore cointegration and Granger causality tests.However, when λ is less than 1 and falls between 0 and 1 (0 < λ < 1) and the Shannon entropy yields a positive value, it suggests uncertainty or a random process.Entropy is evaluated with Shannon entropy [22] and by Havrda-Charvât-Tsallis (HCT) entropy measures [32,33].Kolmogorov-Sinai complexity (KS) is evaluated following the method of [34].The measures above not only give information on the existence of chaos but also the degree of randomness and complexity.Following the steps above, we delve into the nonlinearities and the stochastic process of the variables.

3.
Nonlinearity testing.The BDS test [36,37] is used in various dimensions to determine the nonlinear structure.

4.
Unit root and stationarity testing.Following BDS tests, the unit root and stationary are tested with linear and nonlinear methods.To this end, traditional linear Augmented Dickey-Fuller (ADF) [54]  Modeling series and their relationships are modeled for their marginal and joint distributions with Markov-Switching Generalized Autoregressive Conditional Heteroskedasticity Copula (MS-GARCH-Copula) modeling.The method allows for the modeling of nonlinear and heteroskedastic series for their tail dependence and contagion dynamics.6.
Causality testing.The MS-GARCH-Copula Causality test is an extension of linear Granger causality testing to nonlinear and heteroskedastic series.The method benefits from GARCH-Copula and causality tests.7.
Comparative analysis and robustness.After evaluating the models with diagnostics tests and following nonlinear causality testing in the previous stage, the directions of causality are determined under distinct regimes.At this stage, the results are compared to single-regime GARCH-Copula Causality test results.The findings are examined for possible deviations and differences from the nonlinear methods.

Long-Term Dependence and Fractionality
According to the descriptive tests in Table 1, the series were suspected to have heteroscedasticity in the daily percentage change series.As to be seen in the next section, the White and ARCH-LM tests, the former testing nonlinear forms of heteroskedasticity, and the latter testing ARCH-type effects, confirm these findings.The long memory of the series was tested for return series with Mandelbrot and Wallis R/S and Lo's modified R/S tests.Results are given in Table 2.The R/S test results favor the existence of long-term dependence for all level series.The fractional differentiation parameters are calculated to be 0.990 to 0.995 for all level series.After the first differencing, both R/S tests reported no long-term dependence for EC and the MV.For CO 2 , in contrast, the results indicate long-term dependence.While this is the case, the d parameter is estimated as 0.0525, which is statistically insignificant.For the MV, the d parameter is statistically insignificant; however, for EC, though significant, the parameter is estimated as −0.0395, which is very close to zero, for the first-differenced series.Following the literature, squared first-difference series are assumed to measure the volatility of the series.It should be noted that the results in Table 1 also indicated heteroskedasticity, which will be confirmed by the White and ARCH tests in the next section.Due to the purpose of volatility modeling in this study, if the squared first-differenced series are evaluated, the R/S tests result in the acceptance of long-term dependence and the d-parameters are statistically significant for all volatility variables.Hence, for daily percentage change series, long memory cannot be rejected for CO 2 series, and squared series are taken as proxies for daily volatility and the findings favor long memory and fractionality at the 1% significance level.Interestingly, the d parameter for CO 2 exhibits the largest estimate, 0.84, indicating strong persistence and long-term dependence.For the EV and MV volatility series, the d parameter estimates are 0.475 and 0.263, which are in the range of 0 < d < 0.5, indicating a relatively shorter dependence structure compared with the CO 2 series analyzed.

Chaos and Entropy Tests
After the non-normality and fractionality tests, the analysis concentrated on the chaotic structure.The results for the maximum Lyapunov exponents (λ) calculated for the firstdifferenced series are given in Table 3.The findings indicate that all variables are subject to chaotic behavior, given that λ estimates are bounded between 0 < λ < 1 in all cases.Given that λ > 0, all variables exhibit a divergence following the initial conditions.We noted a clear indication of chaos for all variables, especially for the environment variable for which the positive Lyapunov estimates take values from 0.111 to 0.135.According to the results, we conclude that all variables possess chaotic dynamics.The existence of chaos could result from the noise present, delay time, and the size of the sample [57].The tests conducted for level series also confirmed that the series followed chaotic processes.To save space, the results for level series are not reported, but they are available upon request.
The use of statistical techniques to determine the complexity in nonlinear systems is effectively conducted with entropy [28,58,59].Entropy is a measure of the mechanism's rate of information generation, and to examine chaotic dynamics with entropy, the Shannon entropy results are reported in Table 4.The Shannon entropy results indicate that the degree of uncertainty is highly significant given that all are distinct from zero.The Shannon entropy serves as a randomness or uncertainty meter.According to the results, all variables exhibit completely random or uncertain behavior.For interpretability, Shannon entropy is transformed to the range of [0, 1] by utilizing the minimum and maximum bits.If equal to 1, the results indicate complete randomness and the series could be predicted, while equality to zero indicates that the incident was entirely predictable.A score approaching one indicates a larger departure from the long-term equilibrium.Furthermore, entropy can be understood as a measurement of the information distortion reflected by EC, the MV, and CO 2 .Since the values of Shannon entropy are significantly larger than zero, the results confirm the significant degree of uncertainty in the MV, EC, and CO 2 .The series is further evaluated with the HCT divergence measure, first developed by Havrda and Charvât [32] and afterward by Tsallis [33] for confirmatory reasons.The results of the HCT measure confirm the results obtained by Shannon entropy.
Lastly, we used the KS complexity measure to evaluate the three series analyzed for complexity.In the calculation of the KS complexity measure, we utilized the method by Kaspar and Schuster [34].The KS complexity measure is related to the complexity or unpredictability of a system, and higher values could be interpreted as an indication of unpredictability due to the chaotic characteristics of the series.Compared with the Shannon and HCT measures of entropy, the KS complexity measure focuses on specialized aspects of the series.Both the Shannon entropy and KS measures aim to capture aspects of information and uncertainty; however, they are applied in different contexts.Shannon entropy concentrates on the probability distributions and information and KS entropy concentrates on the behavior of dynamic systems and the predictability of trajectories.Hence, Shannon entropy provides a measure of the information content of a source or distribution, whereas KS entropy is more concerned with the dynamical aspects of systems and their sensitivity to initial conditions.Overall, the results of the KS measure indicate the confirmation of chaotic dynamics and randomness in addition to complexity, and the results obtained from the Shannon, HCT, and KS measures indicate randomness and chaotic behavior in the MV, EC, and CO 2 series analyzed.
The results for energy consumption, air pollution, and the metaverse show deviations in the series from their long-term equilibrium and the existence of a notable and discernible increase in uncertainty in these variables coupled with a decline in the degree of information being carried.Overall, the results indicate uncertainty and randomness in all series analyzed.

Nonlinearity Test Results
The BDS test is a test that benefits from correlation dimensions to examine independence under the null to be tested against the alternative of nonlinearity and dependence.The results are reported in Table 5.The results confirmed that all series are nonlinear for all dimensions at the conventional levels of statistical significance.The ARCH-LM tests indicated ARCH-type heteroskedasticity in all series.The White test was conducted with cross-terms.The test is known to be robust against various forms of heteroskedasticity, especially of nonlinear forms.The White tests led to the conclusion of heteroskedasticity in each series.
The second section of Table 6 presents unit root tests for linear and nonlinear series.The augmented Dickey-Fuller (ADF) linear unit root test indicates that all series are I(1) processes, and to avoid spurious results in econometric models, one needs to utilize firstdifferenced series.The Kwiatkowski-Phillips-Schmidt-Shin (KPSS) stationarity test is robust to breaks and nonlinearity.Kapetanios-Shin-Snell (KSS) tests the null of the unit root against STAR-type nonlinear processes.Both the KPSS and KSS tests result in the conclusion that all series are integrated of order one, and the log-first-differenced series are utilized in the econometric modeling.
If an overall look is presented, the R/S test results favored long-term dependence that became most evident for level series, for which the estimated d parameters were very close to 1.In the unit root tests above, d = 1 by nature of these tests for level series.One could assume that unit root testing is to be in line with the d parameter estimates since the tests assume d = 1 for level variables during the testing process.In the next section, the analyses are conducted by maintaining this strategy of assuming d = 1 for all series following the unit root test and d parameter estimation results for the level series.

MS-GARCH-Copula Estimation Results
Overall, the findings so far indicated chaotic dynamics and nonlinear behavior with heteroskedasticity.In this section, the MS-GARCH-Copula models are estimated, and the results are reported in Table 7.Before the interpretation of the parameter estimates for contagion, a set of measures should be evaluated to define the characteristics of each regime to be analyzed.The regime transition probabilities are reported for this purpose.P(st|st−1) is the conditional probability of the state at period t conditional on the state at period t -1, and in the table, P(1|1) stands for P(st = 1|st−1 = 1).The persistence of each regime is dictated and depends on the regime-switching probability estimates.Similarly, P(2|2) designates P(st = 2|st−1 = 2).Both P(1|1) and P(2|2) are statistically significant at the 5% significance level.For all MS-GARCH marginal models modeling CO 2 , EC, and MV series, the transition probability estimates are no less than 0.7, suggesting a strong level of persistence in each regime; however, P(1|1) > P(2|2), and the first regime is relatively more persistent for all variables except CO 2 .Therefore, for EC and the MV, the second regime is more persistent, while it is the opposite for CO 2 .However, it should be remembered that both regimes have a strong level of persistency.Since the standard error estimates are lower in Regime 1 relative to Regime 2 for all variables, Regime 1 is characterized as the low-variance or, in other words, low-volatility regime.Hence, Regime 2 indicates a relatively higher volatility regime for the MV, EC, and CO 2 .
If the parameter estimates are evaluated, we observe distinct ARCH and GARCH parameter estimates in each regime, confirming asymmetry in volatility dynamics.In addition, the stability criterion, ARCH + GARCH < 1 is satisfied for each regime of the models for the MV, EC, and CO 2 .Finally, the remaining unmodelled volatility is tested with ARCH-LM tests for each model, and we report the p-values in the last column.At a 5% significance level, the null hypothesis of no ARCH-type heteroskedasticity cannot be rejected, leading to findings suggesting that ARCH-type heteroskedasticity dynamics are effectively captured.The last column also reports various statistics including LogL, RMSE, and MAE, suggesting satisfactory results for the fit of the models.For contagion and tail dependence, the copula parameter estimates are given in the last two rows.The copula parameters are measures of tail dependence and are important measures of contagion.The copula parameter estimates in Regimes 1 and 2 are 0.527 and 0.5665 for the CO 2 -EC pair, pointing at the strong degree of positive dependence and contagion between CO 2 emission concentration and energy consumption in both regimes, though relatively larger in the second regime, that is, the high-volatility regime.The CO 2 -MV pair leads to an even larger degree of positive tail dependence of 0.7991 and 0.7614 in Regimes 1 and 2. For the EC-MV pair, the copula estimates are 0.56 in the low-volatility regime and 0.59 in the high-volatility regime.The results indicate a strong degree of contagion among all variables with close but varying magnitudes for regimes and provide support for positive associations.The significance of the coefficients and their sizes suggest a moderate level of tail dependence amid the metaverse and energy consumption and moderate to high levels of dependence between CO 2 emission concentration and the metaverse.

Nonlinear Causality Test Results with Regime-Switching and Copulae
MS-GARCH-Copula models are generalized to MS-GARCH-Copula Causality tests following Section 2. Test results are given in Table 8 below for the regime-dependent causality tests that were conducted.The tests aim to examine and determine the direction of causality in each regime.Causality test results also provide a tool to produce policy recommendations.
The causality between the MV and CO 2 emission concentration test results is given in the first row for each regime.The test results indicate unidirectional causality from the MV to CO 2 in both the low-and high-volatility regimes.Since the causality in the opposite direction does not hold, the direction of causality is confirmed to be regimeindifferent.Causal tests given in the second row evaluate the causal links between EC and CO 2 emission concentration.The results confirm unidirectional causality in both regimes from CO 2 to EC.The last row reports the causality tests between the MV and EC, where the findings favored unidirectional causality from MV to EC in Regime 1.However, in the high-volatility regimes denoted as Regime 2, bidirectional causality could not be rejected.
The results provide evidence of feedback effects in Regime 2 only, and no feedback effect exists from EC to the MV in the low-volatility regimes.Overall, this investigation suggests asymmetry in causality dynamics with varying test statistics and even differentiated results depending on the regime.We aim to examine whether the findings obtained by asymmetric causality testing provided deviations from the non-regime-switching methods because of omitting the nonlinearity of the Markov-switching form.Due to the highly volatile nature of the series, linear Granger causality tests based on vector autoregressive (VAR) models would be inappropriate.Therefore, as a baseline model, we provided single-regime GARCH-Copula-Causality test results [43].The results are given in Table 9.According to the results reported in the first row of Table 9, unidirectional causality is confirmed from the MV to CO 2 , similar to the Markov-switching version reported in Table 8.The second row of Table 9 reports the linear causality test results between CO 2 and EC and, accordingly, no causality exists among them.In contrast, our results in Table 8 suggested no causal link in Regime 1 only, but unidirectional causality from EC to CO 2 in Regime 2, the high-volatility regime.The last column of Table 9 is in line with Table 8's last column, confirming the rejection of non-causality and acceptance of causality between the MV and EC.A summary of comparative results is reported in Table 10.

Discussion with Comparisons to the Existent Literature
The results of the linear approach should be evaluated with caution since all variables exhibit nonlinear behavior, and omitting nonlinearity could lead to inefficiencies possibly with deviated directions causality.
Overall, the results indicate the effects of the metaverse on EC and CO 2 .Here, we aim to compare such findings with the existing literature.However, the number of studies is limited.One reason is data availability.In terms of CO 2 and the MV, we found no studies with similar empirical investigations as our paper.In the neighboring literature, Ferreira et al. show positive and mid-to-high level dependence relationships for fintech, Bitcoin, and CO 2 emissions [60].The effectiveness of nonlinear causality modeling through copulae with GARCH marginals is shown to be efficient for nonlinear relationships and series [43].Compared with Kim et al.'s study, this study extends GARCH-copula to MStype regime switching, and the findings confirm the nonlinear effects of the MV on EC and CO 2 .Such effects are also shown to exist for different financial markets including the metaverse [41].Non-fungible tokens (NFTs) are considered by some as a factor to integrate into the metaverse in the near future.Among the empirical research, few research papers point to NFT market volatility and transmissions to other markets.The dependence between NFT markets, its three sub-markets, and commodity and equity markets including gold and oil prices is shown with quantile cross-spectral coherency and quantile regression methods [61].The NFT market is shown to be in relation to geopolitical risks and equity market volatility [61].
Our paper confirms the dependence between the MV and energy consumption in addition to a strong level of tail dependence with CO 2 emission concentration, a measure of air pollution, under distinct regimes.With a different method and data, Kshetri and Dwivedi suggest the MV could have carbon footprints, which could have environmental implications due to their energy requirements [14].Another study is Zhao and You's paper, which highlighted expectations for the consumer-side MV to contribute to environmental degradation [62].These explanations focus on the technological adaptation level of the MV, and our study has a clear difference in terms of method and theory.While these studies focus on the level of MV technology and its effects on energy consumption, the unavailability of energy consumption directly from MV technology leads to an impossibility of empirical analysis.However, our paper provides a different aspect of the MV.We argue that MV companies have reached a considerable size in the stock market, so the fluctuations in MV stocks could have effects or associations between energy consumption and CO 2 emissions, and the results indicate such findings.In fact, the stock market size of the metaverse and its volatility have positive effects on both EC and CO 2 ; however, these results are evident as tail dependence, confirming contagion at lower and upper tails.Overall, there is a lack of empirical research on the nexus between the metaverse, EC, and CO 2 emission concentration, and an important reason is the availability of data, despite existing projections in the literature.
A limited amount of the literature follows a similar approach in the neighbor literature on cryptocurrencies.Syuhada et al. examine Mana and Theta metaverse cryptocurrencies using connectedness measures, and they highlight connectedness with proposed aggregate Value at Risk (VaR) and Expected Shortfall (ES) measures and copulae modeling [63].Nonlinearity modeling is also shown to improve efficiency in dual rotor systems [64].
By putting forth nonlinear and heteroskedastic dynamics in the MV-EC-CO 2 nexus, our research provides an early empirical analysis with the hope of future work that utilizes new approaches after the availability of metaverse data at an aggregate level.Other than the data issue, this paper also contributes in terms of showing chaotic dynamics of the data analyzed.In the case of the dataset analyzed, this paper confirms that causality dynamics could be modeled with MS-GARCH-Copula-Causality, which allows for nonlinearity and regime-dependent contagion and causality relations.

Conclusions and Policy Recommendations
The assessment of the metaverse and its links to energy and the environment is of crucial importance.This paper aimed to fill the research gap regarding the examination of the metaverse, energy consumption, and CO 2 emission concentrations.The metaverse index utilized in this study comprises a large set of high-tech companies, and it is argued that the dataset characterized with nonlinearity, chaos, and fractionality in certain cases could be modeled with MS-GARCH-Copula and Causality methods for contagion and causality dynamics.Due to the characteristics determined with relevant tests for the data, such a setting is crucial to examine the tail-dependence and associations among metaverse stocks that are under the influence of metaverse technology fluctuations and shocks, possibly with contagion effects on energy consumption and global air pollution at the upper and lower tails of the dataset.
The nonlinearity and fractionality tests favored fractionality, nonlinearity, and longterm dependence in all series analyzed, which became especially evident for the volatility of all series analyzed.Further, chaotic structure, complexity and entropy could not be rejected in CO 2 emission concentrations, the metaverse, or energy consumption.These effects became more vital at the upper and lower tails of the distributions of data, characterized by two distinct regimes in this study, noted as the low-and high-volatility regimes, respectively.Further, fractionality was most evident in the CO 2 emissions concentrations.
Following the determination of chaotic and nonlinear processes in addition to fractionality and long-term dependence, the series were modeled with MS-GARCH-Copula and MS-GARCH-Copula-Causality methods in the modeling stage.The results and conclusions obtained from the MS-GARCH-Copula Causality method are summarized in Figure 1, followed by the results gathered from the GARCH-Copula Causality, which is the single-regime approach.Compared with the former, the latter omits nonlinearity; however, this model deviates from traditional Granger causality tests since it specifically focuses on controlling heteroskedasticity in the series analyzed.
The nonlinear MS-GARCH-Copula and Causality methods provided important insights into two distinct regimes characterized as low-and high-volatility regimes, and according to our findings, tail dependence and contagion existed in both regimes.The level of dependence measured with copula parameter estimates ranged between 0.53 and 0.79 depending on the regime, which pointed to a moderate to high level of positive association between each series, confirming the contagion and tail-dependence between the metaverse, energy consumption, and CO 2 .
After the determination of tail dependence with copula parameters, the direction of causality was determined with the utilized nonlinear approaches.The findings indicated unidirectional causal links from the metaverse to energy consumption, from the metaverse to CO 2 , and from CO 2 to energy consumption in both the low-and high-volatility regimes.Further, unidirectional causality existed from the metaverse to CO 2 , while the relationship became bidirectional in Regime 2, indicating a cyclical pattern of causality with feedback effects in the high-volatility regime.It is evident that during more turbulent periods, the connection between the metaverse and energy consumption is subject to a cyclical self-feeding pattern, and caution is needed for their environmental effects.Hence, in the future, unless the energy dependence of the metaverse can be solved, the environmental effects could be more harmful.Overall, a not-so-favorable result is a positive association between metaverse stocks and EC and CO 2 emissions no matter the regime type.Given the strong level of association, the upward trend in the metaverse necessitates urgent policies to encourage energy-efficient and cleaner technological advancements to reduce the footprint of such technologies.Lastly, the findings were compared with the single-regime GARCH-Copula Causality results for comparative and robustness purposes.Ignoring nonlinearity in the series could lead to deviations in the direction of causality, therefore resulting in inefficiencies in policy applications aiming at controlling the metaverse's effects on energy use and the environmental degradation.
One could argue that many dominant factors that exist, which affect CO 2 emissions and its concentration, the metaverse, and Bitcoin's energy consumption could have limited effects comparatively.This study does not have such a proposal; however, it is a pivotal study that directs attention to recent metaverse technologies in addition to Bitcoin's energy consumption.Further, as discussed in Sections 1 and 2, the energy consumption of such technologies and their environmental effects are shown to amount to those of a set of countries.This paper suggests that close attention is needed for such new technologies, and the findings confirm their effects are becoming evident on the air pollution concentration.
The findings pointed to fractionality, chaos, and nonlinearity in all series.The modeling of the metaverse, Bitcoin's energy consumption, and air pollution with daily datasets should aim at integrating these characteristics of data in future studies.In this respect, this study confirmed tail dependence and contagion with asymmetric copula parameters in addition to causality among the metaverse, energy consumption, and CO 2 emissions.Hence, this study is pivotal in highlighting the implications of the metaverse and Bitcoin's energy consumption and their impacts on global CO 2 emissions concentrations.
Two sets of policies could be generated.First, on the metaverse company side, to mitigate CO 2 emissions, decision-makers should aim to reduce the carbon footprint of the metaverse drastically in manufacturing.However, there is a policy gap for such a commitment, and economic incentives could be generated by a set of policies including carbon tax and carbon subsidies and campaigns for consumer awareness.Second, on the national level, policies should aim to improve energy efficiency, renewable energy investments, and carbon capture technologies to support the mitigation of energy consumption and the environment effects of future metaverse technologies.Note.*** indicates statistical significance at the 1% significance level.The findings of our study, combined with the existing literature, have enabled us to generate a comprehensive set of policy recommendations.The trajectory of the metaverse phenomenon will play a pivotal role in mitigating its environmental impacts.Analogous to the pervasive rise of the internet, e-commerce, and digital devices over the past two to three decades, the metaverse harbors the potential for substantial environmental repercussions unless its energy policies undergo careful reevaluation.As the metaverse gains broader societal accessibility, akin to the ubiquitous presence of the internet and smart devices in contemporary lives, the escalating energy demands associated with metaverse activities will inevitably lead to heightened emissions.To curtail the pace and magnitude of these adverse effects, implementation of intelligent automation technologies to enhance energy efficiency, prioritizing the adoption of renewable energy sources to reduce reliance on fossil fuels and mitigate GHG, e-waste management, and encouragement of recycling in metaverse companies and the metaverse community is of vital importance.In addition, emphasizing the use of eco-friendly materials and efficient recycling methods is important, including recovering resources from used equipment, for example, batteries.Analogous to the pervasive rise of the internet, e-commerce, and digital devices over the past two to three decades, the metaverse harbors the potential for substantial environmental repercussions unless its energy policies undergo careful reevaluation.As the metaverse gains broader societal accessibility, akin to the ubiquitous presence of the internet and smart devices in contemporary lives, the escalating energy demands associated with metaverse activities will inevitably lead to heightened emissions.To curtail the pace and magnitude of these adverse effects, implementation of intelligent automation technologies to enhance energy efficiency, prioritizing the adoption of renewable energy sources to reduce reliance on fossil fuels and mitigate GHG, e-waste management, and encouragement of recycling in metaverse companies and the metaverse community is of vital importance.In addition, emphasizing the use of eco-friendly materials and efficient recycling methods is important, including recovering resources from used equipment, for example, batteries.

Figure 1 .
Figure 1.A summary of the conclusions and empirical results.

Figure 1 .
Figure 1.A summary of the conclusions and empirical results.Note.*** indicates statistical significance at the 1% significance level.The findings of our study, combined with the existing literature, have enabled us to generate a comprehensive set of policy recommendations.The trajectory of the metaverse phenomenon will play a pivotal role in mitigating its environmental impacts.Analogous to the pervasive rise of the internet, e-commerce, and digital devices over the past two to three decades, the metaverse harbors the potential for substantial environmental repercussions unless its energy policies undergo careful reevaluation.As the metaverse gains broader societal accessibility, akin to the ubiquitous presence of the internet and smart devices in contemporary lives, the escalating energy demands associated with metaverse activities will inevitably lead to heightened emissions.To curtail the pace and magnitude of these adverse effects, implementation of intelligent automation technologies to enhance energy efficiency, prioritizing the adoption of renewable energy sources to reduce reliance on fossil fuels and mitigate GHG, e-waste management, and encouragement of recycling in metaverse companies and the metaverse community is of vital importance.In addition, emphasizing the use of eco-friendly materials and efficient recycling methods is important, including recovering resources from used equipment, for example, batteries.

Table 1 .
Descriptive Statistics for log-first-differenced series.
Note.For the Jarque-Bera test of normality, JB is the Chi-square test statistic, and the p-values are in brackets.Mean is arithmetic mean, med. is median, mx. and mn. are maximum and minimum, and skw.and kr. are skewness and kurtosis.
the Kwiatowski et al.(KPSS)[56] test of stationarity is used.While the first assumes a linear unit root, the second tests the unit root against smooth transition autoregressive (STAR) type nonlinearity.Lastly, KPSS is a nonparametric test known to perform well under series with structural breaks and nonlinearity.5.
is followed by a nonlinear Kapetanios et al. (KSS) [55] unit root, and

Table 2 .
R/S test results and fractional difference parameter estimates.

Table 4 .
Entropy and complexity measures.
BDS test statistics are statistically significant evidence of nonlinearity for all series at all dimensions.