So far, only a few studies have evaluated the impact of greenhouse gases emissions on the global and limited area energetics. Furthermore, all of them have concentrated on the increasing of CO2. As new climate projections are now available from a number of climate models under the MPI-ESR-MR experiment, the present study analyses the global and hemispherical energetics under the increase of greenhouse gas forcings that follow Representative Concentration Pathways (RCP26, RCP45, and RCP85). The results have shown a reduction in the LEC intensity as the concentration of greenhouse gases increases, with the RCP85 scenario generating the strongest decrease. For both global and hemispherical domains, zonal kinetic energy is the only energy reservoir which increases in a warmer environment, whereas the conversion between eddy kinetic energy and zonal kinetic energy (CK) is the only energy flux also experiencing an increase. A quantitative analysis of the inner processes involved in the conversion terms shows important changes in the horizontal and vertical eddy-transport of momentum and sensible heat. In the case of CK both vertical and horizontal eddy-transports of momentum play an important role in the increase of zonal kinetic energy for the global domain.
1. Introduction
After Lorenz [1] has derived a set of equations to quantify the energy cycle for the whole atmosphere, many studies have applied the technique to quantify and understand the dynamical processes involved in the energetics of the planet (e.g., [2–9]) and those who are considering just a piece of the atmosphere (e.g., [10–19]). In particular, they differ from each other depending upon the purpose of use and the energetics formalism (space domain, time domain, and mixed space-time domain energetics, see [5]). As discussed in Lorenz [2] and Oort [4], although the absolute value produced by each technique differs from each other, the energy fluxes are qualitatively similar. In general, studies of energetics considering closed domain highlight the maintenance of the general circulation, while studies involving open domain stress the dynamics related to the life cycle of individual atmospheric disturbances.
The increasing of greenhouse gas emissions has amplified the greenhouse efficiency effect by trapping more heat in the mid- and lower troposphere and consequently altering the atmospheric circulation pattern [20]. The consequences of this can be felt through the changes in a broad spatial and temporal-scale of disturbances around both hemispheres as observed in the cases of cyclones and anticyclones (e.g., [21–27]), and tropical cyclones (e.g., [28–30]). For instance, based on a transient simulation run from ECHAM5/OM1 model, [27] revealed that there is a significantly increase in the frequency of the most severe storms over the eastern North Atlantic and North Pacific under an increasing of CO2 relative to the IPCC A1B scenario.
Reference [20] evaluated the atmospheric energetics in the North Hemisphere middle latitude relative to two CO2 concentration conditions (1×CO2 and 2×CO2) from a 10-year-numerical simulation with Canadian Centre for Climate Modelling and Analysis (CCCma) model. The results revealed a reduction in baroclinicity and baroclinic instability for the winter season in the case of doubled CO2 as a consequence of the modified temperature structure (decrease in pole to equator and land/ocean lower troposphere temperature gradients). Furthermore, the “rate of working” of the atmosphere, as characterized by the rate of generation of available potential energy and its conversion to kinetic energy and subsequent dissipation, is reduced in 12% relative to 1×CO2 simulation. The energy cycle was characterized by a decrease in the available potential energies (AZ and AE), with a consequent reduction in the conversion between AZ to AE and AE to KE. The values of zonal kinetic energy (KZ) and the conversion between AE and KE increased, while KE values decreased.
Hernández-Deckers and von Storch [31] based on equilibrium and transient simulations with a coupled atmosphere-ocean GCM showed that the intensity of the Lorenz energy cycle, given by the total conversion rate of potential into kinetic energies, is reduced in a atmosphere with doubled CO2 concentration, corroborating the results of [20]. This reduction is attributed to the decrease in the conversions from potential to kinetic energy of both zonal and eddy modes. Furthermore, they observed an increase (decrease) in the zonal (eddy) kinetic energy reservoir as a response to the increasing of the barotropic (baroclinic) conversion. They found a dual role of the warming pattern in the atmosphere: an increase of AZ in the upper troposphere followed by a weakening of AZ in low levels, with the increase in the static stability being the cause of the former. According to Hernández-Deckers and von Storch [32], based on two equilibrium runs (1×CO2 and 2×CO2) from the coupled atmosphere-ocean ECHAM5/MPI-OM model, the dual response (strengthening in the upper troposphere and weakening below) concerns mainly the transient eddy terms, with the stationary eddies playing a secondary importance.
As seen in the previous works, the atmospheric energetics are broadly modified by significant changes in their composition, as is the case of CO2 increasing. Although the direction of the fluxes evolved in the Lorenz energy cycle is not modified in case of warmed climate drove by the increasing of CO2, the fluxes present changes, which concerns important dynamical processes (e.g., barotropic and baroclinic activities). In this context, the new set of available data relative to the Representative Concentration Pathways is an excellent opportunity to explore the atmospheric energetics under the increasing of greenhouse gases, land use change, and pollution emissions.
The paper is organized as follows: Data and Methodology are discussed in Section 2, and the main results are presented in Section 3, with a discussion of time series of energy and conversion terms and volume-integrated energy components. A follow-up discussion is presented in Section 4, commenting on new insights and perspectives brought by the results.
2. Data and Methodology
In this study, we use a set of four experiments originating from the new Max Planck Institute Earth System Model (MPI-ESM) used in the Coupled Model Intercomparison Project phase 5 (CMIP5). The MPI-ESM consists of the atmosphere model ECHAM6 [34], including a land vegetation module and the ocean model MPIOM [35]. The MPI-ESR-MR model has a spatial resolution equivalent to 1.875°×1.875° on a Gaussian grid, as the model runs with a horizontal resolution of 63 waves triangular spectral truncated and 95 levels in the vertical (T63L95), see [36]. The ocean component of the MPI-ESM, MPIOM, is a free-surface ocean general circulation model formulated on an Arakawa-C grid in the horizontal and a z-grid in the vertical. The MPIOM horizontal resolution translates into a grid ranging from 15 km around Greenland to 185 km near the equator, for more details see [37]. Within the set of runs designed from CMIP5, we have considered three RCPs [38, 39]: (1) one mitigation scenario leading to a very low radiative forcing level of 2.6 W/m2 RCP26 [40], (2) one medium stabilization scenario RCP45 [41], and (3) one representative of a very high emission scenario RCP85 [42]. A historical experiment, representing the control experiment (CTRL), is used here for comparison purposes. The historical experiment covers much of the industrial period from the mid-19th century to near present (1880–2005); however, we used here the period from 1979 to 1999 as the base period. The RCPs experiments continue from the historical run until 2100. The RCPs are characterized in terms of land use cover, greenhouse gas (CO2, CH4, and N2O), and air polluting scenario projections (e.g., SO2 and NOx) [43]. To analyze the impact of the different rate of greenhouse gas emissions on the global and hemispherical LEC, mean daily atmospheric data of zonal and meridional components of the horizontal wind vector (u and v), air temperature (T), vertical component of velocity in pressure coordinates (ω=dp/dt), and geopotential height (z) are used. The atmospheric variables were integrated from the lowest model level, considered to be at the surface (1000 hPa) up to 1 hPa. The impacts are quantified in terms of anomaly produced in the comparison between the RCPs and the control experiments. The period of analysis covers two specific ranges: 1979–1999 and 2079–2099, hereafter called Late 20th century (L20C) and Late 21st century (L21C), respectively.
The simplest representation of the Lorenz energetics is the two-component energy cycle (Figure 1(a)), see [33].
Lorenz energy cycle representation for (a) two-component energy cycle in which available potential energy denoted by the letter A is generated (G), converted (C) to kinetic energy (K), and subsequently dissipated (D). (b) Four-component energy cycle with A and K decomposed into mean and eddy parts. Adapted from [33].
Consider(1a)∂A∂t=-C(A,K)+G,(1b)∂K∂t=C(A,K)-D,
where A is the amount of available potential energy, G its generation rate, C is the energy conversion between A and kinetic energy K, and D is the dissipation rate of K. When C(A,K) is positive, K increases at the expense of A, and vice versa.
The four-component energy cycle is characterized by the subdivision of A and K into their zonal and eddy components and the set of equations can be written as(2a)∂AZ∂t=-CZ-CA+GZ,(2b)∂AE∂t=CA-CE+GE,(2c)∂KZ∂t=CK+CZ-DZ,(2d)∂KE∂t=CE-CK-DE.
In the four-component energy cycle the potential energy is portioned into zonal available potential energy (AZ) and eddy available potential energy (AE). The kinetic energy can be treated similarly as potential energy; that is, the kinetic energy of the flow can be separated into its zonal (KZ) and eddy (KE) parts. In the energy cycle of Lorenz, AZ and AE are produced by generation processes named GZ and GE, while reservoirs of energy are connected by conversion terms CZ, CA, CE, and CK (resp, the conversion from AZ to KZ, AZ to AE, AE to KE, and KE to KZ). The CZ term involves transport of sensible heat by meridional circulation affecting AZ and KZ by altering the variance of zonally averaged temperature and wind. It depends upon the covariance between zonal means of vertical velocity and temperature, and it reflects the growth of KZ at the expense of AZ when there is warm air rising and cold air sinking in the same latitude. The CA involves horizontal and vertical eddy-transports of sensible heat. This term affects AZ and AE by altering the variance of temperature and depends upon the transport of sensible heat along the gradients of zonal averaged temperature. The CE term involves vertical eddy-transports of sensible heat, whose values by unit area are independent of latitude. It affects AE and KE by altering the variance of temperature and wind within latitude circles, and its absolute value indicates to what extent an environment is baroclinic. The term CK involves horizontal and vertical eddy-transports of momentum, and it depends upon the transport of angular momentum along the gradients of angular velocity. The absolute value of CK indicates to what extent an environment is barotropic or supports atmospheric disturbances to intensify via barotropic conversion.
For the case of open domains the four-components energy cycle (Figure 1(b)) must include the boundary and work done processes.
BAZ, BAE, BKZ, and BKE terms, respectively, represent fluxes of AZ, AE, KZ, and KE from the lateral and vertical boundaries. BΦZ and BΦE terms denote the dynamical mechanisms in which they produce or destroy kinetic energy. As described by [44], this is due to the fact that BΦZ, BΦE, CZ, and CE derivatives involve a single term in the form V·∇Φ which represents the appearance of kinetic energy via a cross-isobaric flow towards low pressure. The destruction of kinetic energy occurs when there is cross-isobaric flow towards the high pressure. However, here BΦZ and BΦE are added to the dissipation term.
The last two terms in (3c) and (3d) can be determined in form of residuals, with the incorporation of calculation errors εKZ and εKE in the following form: (4a)RKZ=BΦZ-DZ+ϵKZ,(4b)RKE=BΦE-DE+εKE.
In similar form, the intrinsic errors resulting from the numerical estimations for the other terms in (3a) and (3b) are also included in the estimation of the generation terms:(5a)RGZ=GZ+εAZ,(5b)RGE=GE+εAE.
Detailed mathematical expressions for all components of the energy budget, as described by (3a) and (3d) are presented in the appendix.
Here the Lorenz diagram considers that the conversion terms are always positive, with the arrows indicating the direction of the fluxes. The rates of generation of available potential energy and dissipation of kinetic energy are obtained as residuals by considering that in the long run the gain or loss of potential and kinetic energy is negligible. The four-component energy cycle is used in this study for both global and hemispherical energetics. However, in the first case the boundary fluxes are not included.
3. Results3.1. Time Series of Energy and Conversion Terms
The climatological annual cycles of generation, available potential energy and kinetic energy, with their zonal and eddy components, for the control and RCPs experiments for the whole globe are illustrated in Figure 2. The error bars, indicating the 95% confidence level about the mean, are included in this panel to show the interannual variation experienced by the energy terms. The values of GZ are maximized in May and September followed by two minimum during June and January (Figure 2(a)). The results show diminishing of GZ as greenhouse gases concentration increases. An opposite behavior is observed for GE (Figure 2(b)). In this case GE presents high negative values in the control experiment and declines to minimum negative values for RCP85 experiment. The implications of the fall of GZ and GE (negatively) as the concentration of greenhouse gases increases are that less AZ and AE are created by diabatic processes and consequently affecting the release of KZ and KE in the atmosphere by the conversion terms. The decrease in the generation of available potential energy occurs due to the fact that in a warmed planet the horizontal temperature gradient is reduced, as can be seen from Figure 1 of [20] and Figure 2 of [31].
Climatology of the generation terms GZ (a) and GE (b), and the energy terms AZ (c), AE (d), KZ (e), and KE (f) for a global domain integrated between 1000 and 1 hPa for RCP85, RCP45, RCP26, and control experiments. Error bars indicate 95% confidence intervals about the mean. Units are 105 J/m2 for reservoirs and W/m2 for generation terms.
According to Figure 2(c) the zonal available potential energy is characterized by a well-defined seasonal variation in the case of RCP26 and RCP45, while in the RCP85 experiment the seasonal variations is less visible. The impact of increasing in greenhouse gases is stronger in the RCP45 and RCP85 when compared with the control experiment. The bars of control and RCP26 experiments overlap each other in the major part of the year suggesting a weak impact in comparison with other RCPs. A similar pattern is observed for the eddy available potential energy (Figure 2(d)), in which RCP45 and RCP85 experiments produce the strongest impact. The reduction of the seasonal variations as the concentration of greenhouse gases increases is clear. The results still show that there is decrease of the available potential energies (AZ and AE) due to the drop of the generation terms. As can be seen, the strongest (weakest) weakening in the available potential energy, for both zonal and eddy parts, occurs preferentially in the high (low) range emission experiments. The decreasing in AE as the concentration of greenhouse gases increases is because less AE is destructed or converted into GE. The destruction of GE occurs when colder air is advected to warmer regions and warmer air is advected to colder regions in the same latitude, a process that alters the variance (correlation) of diabatic heat and the temperature field. A brief discussion about the role of this term in produced available potential energy is found in Lorenz [2]; however, a more thorough discussion about this term is limited by the absence of diabatic heating/cooling data from the global circulation models. Figures 2(e) and 2(f) show the values of kinetic energies (KZ and KE) as function of the RCPs and CRTL experiments. The annual cycle of KZ is characterized by a well-defined seasonal variation, with highest (lowest) values from July to September (April to June). One can see that KZ clearly increases, while KE remains approximately unchanged around the year as the concentration of greenhouse gases increases. The increase of KZ occurs in the RCP85 forcing scenario, with minor and almost equally changes in RCP45 and RCP26. As the generation (or destruction) of KZ depends upon the joint effect of the transport of angular momentum along the gradients of angular velocity (CK), transport of sensible heat by meridional circulation (CZ) and the dissipation of kinetic energy by friction, represented by DZ (2c), and an evaluation of the conversion terms are necessary.
The global climatological annual cycles of the vertically-integrated conversion and dissipation terms for the control and RCPs experiments are presented in Figure 3. It can be seen that CA and CK conversions are significantly altered, between the months of June and August, in the RCPs scenarios forcing (Figures 3(b) and 3(d)), with the strongest changes occurring in the RCP85. The results show a decrease in the transfer of energy from AZ to AE as the concentration of greenhouse gases increases, probably because of the reduction of the eddy-transports of sensible heating. This is expected in a warmed environment with less horizontal contrast of temperature. Reference [31] showed that in a warmed planet with doubled CO2 concentration, the meridional temperature gradient decreases due to the anomalous heating in the tropical and polar regions. The changes in CZ and CE are modest and not statistically significant (Figures 3(a) and 3(c)), suggesting a weak relation to the increasing of greenhouse gases concentration. The term CK, which acts around the year as a source of energy for the reservoir of KZ, shows different result. It indicates that in a warmed climate the conversion from KE to KZ increases; that is, there is an increase of eddy-transports of angular momentum along the gradients of angular velocity. Figures 3(e) and 3(f) show an opposite behavior between the dissipation of KZ and KE during the course of the year; that is, RKZ presents the highest values from July to September, when RKE is minimum. The same relation is observed from November to January. These results show that the zonal mean motions are strongly dissipated by friction processes when they are maximized (Figure 2(e)). The same behavior is observed for the eddy motions.
Climatology of the conversion terms CZ (a), CA (b), CE (c), and CK (d) for the whole globe integrated between 1000 and 1 hPa for RCP85, RCP45, RCP26, and control experiments. Error bars indicate 95% confidence intervals about the mean. Units are W/m2.
3.2. Volume-Integrated Energy Components
The mean annual energy cycle (two- and four-component) for the entire globe, southern and northern hemispheres calculated from RCPs and historical experiments are presented, respectively, in Figures 4, 5, and 6. The two-component energy cycle for the global domain (Figure 4(a)) shows an increase in the generation term G as a response to the increasing in greenhouse gases concentration, followed by a decrease in the available potential energy A, conversion term C(A,K) and the dissipation D. However, the total kinetic energy increases. Each of these terms is characterized by the sum of its four-component energy cycle (Figure 4(b)). As can be noted, the increase in G is ~61% (44%, 23%) in the case of RCP85 (RCP45 and RCP26). The others components have the following conditions: decreasing in A of ~28% (16%, 8%), decreasing in C(A,K) of ~4% (2%, 0.2%), increasing in K of ~15% (7%, 4%), and a decreasing in the dissipation term D of ~4% (1%, 0.2%).
Representation of annual mean of spatial domain Lorenz energy cycle for the entire globe in which (a) is representative of a two-component energy cycle and (b) denotes the four-components diagram. Values from up to below in each figure are for RCP85, RCP45, RCP26, and control experiments. Units are 105 J/m2 for reservoirs and W/m2 for conversion, generation, and dissipation terms. In the two-component energy cycle the value of G is given by the sum of its zonal and eddy part (G=GE+GZ), A=AZ+AE, C(A,K)=CZ+CE, K=KZ+KE, and D=RKZ+RKE. Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.
Representation of annual mean of spatial domain Lorenz energy cycle for the Southern Hemisphere in which (a) is representative of a two-component energy cycle and (b) denotes the four-component diagram. The boundaries of available potential energy (BA) and kinetic energy (BK) are included in the two-component energy cycle. Values from up to below in each figure are for RCP85, RCP45, RCP26, and control experiments. Units are 105 J/m2 for reservoirs and W/m2 for conversion, generation, and dissipation terms. Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.
Representation of annual mean of spatial domain Lorenz energy cycle for the Northern Hemisphere in which (a) is representative of a two-component energy cycle and (b) denotes the four-component diagram. The boundaries of available potential energy (BA) and kinetic energy (BK) are included in the two-component energy cycle. Values from up to below in each figure are for RCP85, RCP45, RCP26, and control experiments. Units are 105 J/m2 for reservoirs and W/m2 for conversion, generation, and dissipation terms. Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.
The four-component energy cycle (Figure 4(b)) indicates a general decrease in the values of generation, available potential energy, the conversion terms CZ, CA, and CE, and the dissipation terms as the concentration of greenhouse gases increases. The results show that RCP85 are the experiment with the maximum reduction relative to control experiment. Kinetic energy and the conversion term CK are quantities that increase as the greenhouse gases concentration rises, as can be confirmed from RCPs forcing scenarios. It can be seen that both available potential energies present similar reduction, ~28% for RCP85 and below 10% in the case of RCP26 experiment. The reservoirs of kinetic energy are increased, in which the zonal part presenting the strongest change, that is, an increasing of 21.8% in KZ versus a small increase of ~1% in KE for the RCP85 experiments. This is suggestive of a strengthening of the atmospheric zonal motions. The strength of the baroclinic energy flow, characterized by the conversion from AZ to AE (CA) and AE to KE (CE), is reduced in a condition of warmed climate, as can be seen by the reduction of ~39% in CA and ~1.8% in CE for the RCP85 experiments. In the case of RCP45 and RCP26, the reduction is less strong. The barotropic flow, denoted by the transfer of energy between KE and KZ (CK), is increased showing that the zonal flow grows at the expense of the eddy motion. The conversion term CZ, responsible by the increasing of kinetic energy due to the meridional exchange of mass in the vertical plane, is reduced. These results are attributed to the less available potential energy that is produced in a warmed planet.
The two-component energy cycle for the SH domain displayed in Figure 5(a) shows an increase in the generation term G, followed by a decrease in the available potential energy A, conversion term C(A,K), and the dissipation D. The figure still shows a negative decrease in the flux of A(BA) and an increase in the flux of K(BK). Thus, BA acts as sink of A, and it is less intense for a warmed climate, while BK acts as a source of KZ. As can be observed, the increase in G, in the case of SH domain, is ~34% (24%, 9%) in the case of RCP85 (RCP45 and RCP26). The total available potential energy decreases by ~27% (16%, 8%), followed by a negative decrease of ~5% (3%, 3%) in the BA. The decrease in C(A,K) is ~9% (5%, 3%), with an increasing in K of ~17% (8%, 5%), an increase of ~66% (25%, 8%) in the BK and a decrease in the dissipation term D of ~9% (5%, 3%).
The four-component energy cycle displayed in Figure 5(b) indicates a general decrease in the values of generation, available potential energy, the conversion terms CZ, CA, and CE, and the dissipation terms as the concentration of greenhouse gases increases. Kinetic energy and the conversion term CK are quantities that increase as the greenhouse gases concentration rises, as can be confirmed from RCPs forcing scenarios. In the four-component energy cycle it can be seen that the negative reduction in the boundary terms associated to the available energies is due mainly to the BAZ, as the reduction in BAE is one order magnitude small. The increase in the kinetic energy boundary terms is particularly due to the BKZ term. The energy reservoirs of available energies and kinetic energies, respectively, decrease and increase. For the RCP85 experiment the increase of KE is small (~1%) when compared with increases in KZ (~24%). Both decrease and increase values in the energy reservoirs are more significant in the case of highest greenhouse gases emissions. The baroclinic flow, characterized by positive values of CA and CE, decreases as the planet warms, followed by an increase in the barotropic component, and a reduction in the CZ values. These patterns are similar to the global domain.
Figure 6(a) shows the two-component energy cycle for the NH domain. The results show an increase in the generation term G, a decrease in the available potential energy A, a decrease in BA, a decrease in the conversion term C(A,K), an increase in K, a negative increase in BK, and a decrease in the dissipation term D. In the NH case, BA acts as source of A, and it is less intense for a warmed climate, while BK acts as a sinking of KZ, and it is more intense for a warmed climate. As can be seen, the increase in G, in the case of NH domain, is ~172% (108%, 70%) in the case of RCP85 (RCP45 and RCP26). The total available potential energy has a decreasing of ~33% (20%, 13%), followed by a decrease of ~31% (17%, 10%) in the BA. The decrease in the values of C(A,K) is ~10% (6%, 1%), with an increasing in K of ~11% (6%, 2%), a negative increase of ~64% (42%, 28%) in the BK, and a decrease in the dissipation term D of ~10% (3%, increase of 1%).
In the four-component energy cycle for the NH (Figure 6(b)), observed is a pattern similar to the global and SH domain, that is, reduction in the values of available potential energies, followed by a respective small and a significant increase in KE and KZ, reduction in the generation terms, reduction in the baroclinic flux, increase in CK, and reduction in the dissipation terms. The results also show a reduction in the boundary terms associated to the available energies is mainly due to the BAZ, as the reduction in BAE is inexistent. The negative increase in the kinetic energy boundary terms is particularly strong in both terms.
As can be seen, the reductions and increases in the energetics still depend on the domain evaluated. For example, in the case of RCP85 projection the increase of KZ for the global domain is ~21.8% relative to the control experiment, while KZ in the SH and NH increases about 24.6% and 18.0%, respectively. The same pattern of increase in KZ (and the remaining reservoirs) is followed by the others RCPs projection, though with less intensity.
For all analyzed domains the intensity of the LEC, measured by the conversion between potential and kinetic energy in both zonal and eddy modes, is reduced in a planet with more greenhouse gases, as noted from the values produced by the RCPs experiments. The largest reductions occur for the NH (10.35%), when compared with global (4.68%) and SH (9.57%) domains. The strongest increase is seen to occur in the total kinetic energy for the SH (~17%), followed by an increase of 41.02% in CK. This result indicates that the increase in KZ is related to the increase in CK, as KE acts as a source of KZ. However, another form to attribute the increase of KZ due to CK is to quantifying what dynamic process was the most important to the increase of KZ. This can be achieved by evaluating or quantifying the terms involved in the calculation of CK, see equation in the appendix. The first three terms of CK are related to the horizontal transport of momentum, while the last two terms deal with the vertical transport of momentum. The quantification of their changes for the RCPs experiments can bring us important characteristics about the atmospheric behavior in a climate change perspective.
For the SH (NH, Global) domain, the increase in CK, in the RCP85, is related to an increase of 47.8% (5.2%, 38.3%) in the horizontal eddy-transport of momentum and an increase of 23.9% (24.6%, 21.0%) in the vertical eddy-transport of momentum relative to control experiment. For the experiments RCP26 and RCP45 the increase in the horizontal (vertical) eddy-transport of the SH is 22.8% (8.4%) and 29.7% (10.1%), respectively, showing that in the SH the increase in greenhouse gases produces an increase in both horizontal and eddy transport of momentum, with the former being more strong. In the case of the global domain, for RCP26 and RCP45, the increase of horizontal (vertical) eddy-transport of momentum is 18.7% (6.1%) and 25.0% (9.6%), respectively, and for the NH domain, it corresponds to −1.5% (2.1%) and 2.2% (9.5%), respectively.
Calculations of changes for individual terms of CA indicate that the reduction in this component of the baroclinic conversion is due both horizontal and vertical eddy-transport of sensible heat. For the RCP85 experiment (RCP45 and RCP26) in the global domain, the reduction in the horizontal eddy-transport of sensible heat is 38.7% (22.9% and 12.4%), while for the vertical eddy-transport of sensible heat, the change is of 77.7% (55.5% and 33.3%). For the SH and NH the reduction in the horizontal eddy-transport of sensible heat is of 35.1% (20.7% and 36.4%) and 39.2% (24,0% and 14.9%), respectively, for the RCP85 (RCP45 and RCP26). The reduction in the vertical eddy-transport of sensible heat for the SH and NH is 22.2% (33.3% and 55.5%) and 70% (52% and 38%), respectively, for the RCP85 (RCP45 and RCP26). The reduction in the values of CE term is below 3% relatively to the control experiment for all RCPs.
4. Concluding Remarks
In this work a global and hemispherical energetic study, using the traditional Lorenz two- and four-component energy cycle, is carried out to access the impact of the radiative forcing relative to the greenhouse gases emission on the global and hemispherical LEC. The energetic analysis is based on the RCP85, RCP45, and RCP26 forcing experiments and a control simulation from the MPI-ESM-MR ocean-atmosphere coupled model. The results show a decrease in most of the energy terms from the Lorenz energetics for the global domain, with few exceptions, suggesting a weakness of the LEC in the L21C period, as previous works have shown. However, here we have presented the changes in the Lorenz energetics related to the new RCPs scenarios, which are the results of the increase in greenhouse gases. In the previous literature, the energetics are quantified and evaluated under a condition of CO2 increasing only. In this context, we have evaluated and found that the increase in the barotropic component of the energy cycle, which is the main term responsible to the increase in the zonal component of the kinetic energy, is related to the increase in the horizontal transport of momentum, with a less increase in the vertical component. The reduction of the values in the energy cycle is evident and dependent of the RCP experiment, with the strongest reduction occurring in the highest radiative forcing experiment, RCP85. The results also show an increase in CK, KZ, and KE components of the LEC. Similar reduction in the intensity of the LEC is observed for the SH and NH, with RCP85 emission scenario influencing the most intense decrease in the energy terms. In Particular, for this scenario the energetic projection for the SH indicates an increase in KZ of 24.6%, which is higher than the global (21.8%) and NH (18.0%) domains. The anomalous increase of KZ in the SH (NH and Globe) for the L21C period is attributed to an increase of 47.8% (5.2%, 38.3%) and 23.9% (21.0%, 24.6%), respectively, in the horizontal and vertical eddy-transport of momentum in RCP85. The increase in these transports are less intense for RCP45 and RCP26 mitigated scenarios. The reduction in the horizontal eddy-transport of sensible heat generates a decrease in the conversion between AZ and AE, in the RCPs scenarios. These results contribute to the findings of [20, 31] who have shown that CK increases in a warmed climate; however, our results indicates that the increase in CK is mainly due the increase in the horizontal components of the momentum transport. The way and reason the increase of greenhouse gases emissions impact the transport of momentum and heat will be addressed in another work.
Although the results presented here are only based on one CMIP5 model and they do not probably represent a unique picture of the future behavior of energy cycle, they suggest important changes in the energetics of the atmosphere under a warmed climate, which can motivate further investigation on the atmosphere global circulation.
AppendixLorenz Energetics
Equations of Lorenz in spatial domain.
Mathematical expression of components in the energy balance (3a)–(3d).
The zonal mean of a variable X between longitudes λ1 and λ2 is given by
(A.1)[X]λ=1λ2-λ1∫p1p2Xdλ.
The eddy component of this variable is
(A.2)(X)λ=X-[X]λ.
The mean of X over an area bounded by longitudes λ1 and λ2 and latitudes ϕ1 and ϕ2 is
(A.3)[X]λϕ=1λ2-λ11sinϕ2-sinϕ1∫ϕ1ϕ2∫λ1λ2Xcosϕdλdϕ.
Define the quantity as
(A.4)([X]λ)ϕ=[X]λ-[X]λϕ.
The four energy forms in the SDLEC are
(A.5)AZ=∫p1p2[([T]λ)ϕ2]λϕ2[σ]λϕdp,AE=∫p1p2[(T)λ2]λϕ2[σ]λϕdp,KZ=∫p1p2[[u]λ2+[v]λ2]λϕ2gdp,KE=∫p1p2[(u)λ2+(v)λ2]λϕ2gdp,
where p1 and p2 are, respectively, the upper and lower pressure boundaries, T is the temperature, g the magnitude of the acceleration of gravity, u and v are the eastward and northward components of the wind, respectively, and σ is the static stability parameter which is given by
(A.6)[σ]λϕ=[gTcp-gpR∂T∂p]λϕ,CZ=∫p1p2-[([ω]λ)ϕ([T]λ)ϕ]λϕRgpdp,CE=∫p1p2-[(ω)λ(T)λ]λϕRgpdp,CA=∫p1p2-(×∂∂p(([T]λ)ϕpR/c_p[σ]λϕ)]λϕ[(v)λ(T)λ2[σ]λϕr∂([T]λ)∂ϕ]λϕ000000000+[(([T]λ)ϕpR/c_p[σ]λϕ)(ω)λ(T)λpR/cp000000000000×∂∂p(([T]λ)ϕpR/cp[σ]λϕ)]λϕ)dp,CK=∫p1p21g[(v)λ(u)λcosϕr∂∂ϕ([u]λcosϕ)]λϕdp+∫p1p21g[(v)λ2r∂[v]λ∂ϕ]λϕdp+∫p1p21g[[v]λ(u)λ2tanϕr]λϕdp-∫p1p21g[(ω)λ(u)λ∂[u]λ∂p]λϕdp+∫p1p21g[(ω)λ(v)λ∂[v]λ∂p]λϕdp,
where r denotes the mean radius of the Earth.
The generation of APE terms and kinetic energy dissipation terms are
(A.7)GZ=∫p1p2[([q]λ)ϕ([T]λ)ϕ]λϕcp[σ]λϕdp,GE=∫p1p2[(q)λ[T]λ]λϕcp[σ]λϕdp,
where q denotes the diabatic processes associated with the generation of APE.
The energy transport integrals are
(A.8)BAZ=c1∫p1p2∫λ1λ212[σ]λϕ×(2([T]λ)ϕ(T)λu+([T]λ)ϕ2u)λ1λ2dϕdp+c2∫p1p212[σ]λϕ00×(2[(v)λ(T)λ]λ([T]λ)ϕcosϕ0000000000000+([T]λ)ϕ2[v]λcosϕ)ϕ1ϕ2dp-12[σ]λϕ(([T]λ)ϕ2+[[ω]λ([T]λ)ϕ2]λϕ[2(ω)λ(T)λ]λ([T]λ)ϕ0000000000+[[ω]λ([T]λ)ϕ2]λϕ)p1p2,BAE=c1∫p1p2∫ϕ1ϕ212[σ]λϕ[u(T)λ2]λ1λ2dϕdp+c2∫p1p212[σ]λϕ([(T)λ2]λcosϕ)ϕ1ϕ2dp-([ω(T)λ2]λϕ2[σ]λσ)p1p2,BKZ=c1∫p1p2∫ϕ1ϕ212g(u[u2+v2-(u)λ2-(v)λ2])λ1λ2dϕdp+c2∫p1p212g([vcosϕ[u2+v2-(u)λ2]]λ)ϕ1ϕ2dp-(12g[ω[u2+v2-(u)λ2-(v)λ2]]λϕ)p1p2,BKE=c1∫p1p2∫ϕ1ϕ212g(u[(u)λ2+(v)λ2])λ1λ2dϕdp+c2∫p1p212g([vcosϕ[(u)λ2+((v)λ2)]]λ)ϕ1ϕ2dp-(12g[ω[(u)λ2+(v)λ2]]λϕ)p1p2,
where c1=-[r(λ2-λ1)(sinϕ2-ϕ1)]-1 and c2=-[r(sinϕ2-ϕ1)]-1.
Finally the integrals for BΦZ and BΦE are
(A.9)BΦZ=c1∫p1p2∫ϕ1ϕ21g([v]λ([Φ]λ)ϕ)λ1λ2dϕdp+c2∫p1p21g(cosϕ[v]λ([Φ]λ)ϕ)ϕ1ϕ2dp-1g([([ω]λ)ϕ([Φ]λ)ϕ]λϕ)p1p2,BΦE=c1∫p1p2∫ϕ1ϕ21g((u)λ(Φ)λ)λ1λ2dϕdp+c2∫p1p21g([(v)λ(Φ)λ]cosϕ)ϕ1ϕ2-1g([(ω)λ(Φ)λ]λϕ)p1p2,
where Φ(=gz) is the geopotential.
Acknowledgments
José Augusto P. Veiga and Tercio Ambrizzi were supported by CNPq (150356/2011-0). Tercio Ambrizzi also has the support from the FAPESP (08/58101-9), CLARIS LPB, and INCLINE/USP.
LorenzE. N.Available potential energy19557157167LorenzE. N.1967World Meteorological OrganizationWiin-NielsenA.On the intensity of the general circulation of the atmosphere19686431132310.1029/RG006i004p00559OortA. H.On estimates of the atmospheric energy cycle19649211483493PeixotoJ. P.OortA. H.The annual distribution of atmospheric energy on a planetary scale1974201521492159KungE. C.TanakaH.Energetics analysis of the global circulation during the special observation periods of FGGE19834011257525922-s2.0-002096833110.1175/1520-0469UlbrichU.SpethP.The global energy cycle of stationary and transient atmospheric waves: results from ECMWF analyses1991453-41251382-s2.0-000992448010.1007/BF01029650LiL.IngersollA. P.JiangX.Lorenz energy cycle f the global atmosphere based on reanalysis datasets200734L1681310.1029/2007GL029985MarquesC. A. F.RochaA.Corte-RealJ.Global diagnostic energetics of five state-of-the-art climate models2011369-10176717942-s2.0-7995457473310.1007/s00382-010-0828-9LauN.-C.The structure and energetics of transient disturbances in the northern hemisphere wintertime circulation19793669829952-s2.0-0018721523DepradineC. A.Energetics of large-scale motion in the tropics during GATE at 250 mb1980108788689510.1175/520-0493MichaelidesS. C.Limited area energetics of Genoa cyclogenesis1987115113262-s2.0-002349908810.1175/1520-0493MichaelidesS. C.A spatial and temporal energetics analysis of a baroclinic disturbance in the Mediterranean19921207122412432-s2.0-0027009210BakerW. E.BrinY.A comparison of observed and forecast energetics over North America19851114686416632-s2.0-002221003710.1002/qj.49711146819KrishnamurtiT. N.SinhaM. C.JhaB.MohantyU. C.A study of South Asian monsoon energetics19985515253025482-s2.0-003245324810.1175/1520-0469(1998)055<2530:ASOSAM>2.0.CO;2 KrishnamurtiT. N.PattnaikS.StefanovaL.The hurricane intensity issue200511371886191210.1175/MWR2954.1VeigaJ. A. P.PezzaA. B.SimmondsI.Silva DiasP. L.An analysis of the environmental energetics associated with the transition of the first South Atlantic hurricane20083515L158062-s2.0-5374909546910.1029/2008GL034511PezzaA. B.VeigaJ. A. P.SimmondsI.KeayK.MesquitaM. D. S.Environmental energetics of an exceptional high-latitude storm201011139452-s2.0-7795445211010.1002/asl.253PezzaA. B.GardeA. L.VeigaJ. A. P.SimmondsI.Large scale features and energetics of the hybrid subtropical low “Duck” over the Tasman Sea201310.1007/s00382-013-1688-xBoerG. J.Some dynamical consequences of greenhouse gas warming19953347317512-s2.0-2134443250310.1080/07055900.1995.9649551AgeeE. M.Trends in cyclone and anticyclone frequency and comparison with periods of warming and cooling over the Northern Hemisphere19914226326710.175/1520-0442McCabeG. J.ClarkM. P.SerrezeM. C.Trends in Northern Hemisphere surface cyclone frequency and intensity20011412276327682-s2.0-003587465910.1175/1520-0442SerrezeM. C.CarseF.BarryR. G.RogersJ. C.Icelandic low cyclone activity: climatological features, linkages with the NAO, and relationships with recent changes in the Northern Hemisphere circulation19971034534642-s2.0-694425669310.1175/1520-0442PezzaA. B.AmbrizziT.Variability of Southern Hmisphere cyclone and anticyclone behavior: further analysis20031671075108310.1175/1520-0442(2003)016<1075:VOSHCA>2.0.CO;2PezzaA. B.SimmondsI.RenwickJ. A.Southern Hemisphere cyclones and anticyclones: recent trends and links with decadal variability in the Pacific Ocean20072711140314192-s2.0-3454859160310.1002/joc.1477BlackM. T.PezzaA. B.A universal, broad-environment energy conversion signature of explosive cyclones201340245245710.1002/grl.50114UlbrichU.LeckebuschG. C.GriegerJ.Are greenhouse gas signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?2013221616810.1127/0941-2948/2013/0420SugiM.YoshimuraJ.Decreasing trend of tropical cyclone frequency in 228-year high-resolution AGCM simulations20123919L1980510.1029/2012GL053360OouchiK.YoshimuraJ.YoshimuraH.MizutaR.KusunokiS.NodaA.Tropical cyclone climatology in a global-warming climate as simulated in a 20 km-mesh global atmospheric model: frequency and wind intensity analyses20068422592762-s2.0-3374519523610.2151/jmsj.84.259HoyosC. D.AgudeloP. A.WebsterP. J.CurryJ. A.Deconvolution of the factors contributing to the increase in global hurricane intensity2006312577094972-s2.0-3364569357910.1126/science.1123560Hernández-DeckersD.von StorchJ.-S.Energetics responses to increases in greenhouse gas concentration20102314387438872-s2.0-7795441936210.1175/2010JCLI3176.1Hernández-DeckersD.von StorchJ.-S.The energetics response to a warmer climate: relative contributions from the transient and stationary eddies20112110512010.5194/esd-2-105-2011BoerG. J.LambertS.The energy cycle in atmospheric models20083043713902-s2.0-3834909141910.1007/s00382-007-0303-4StevensB. M.GiorgettaM.EschM.Atmospheric component of the MPI-M earth system model: ECHAM620125214617210.1002/jame.20015CruegerT.HoheneggerC.MayW.Tropical precipitation and convection changes in the Max Planck institute earth system model (MPI-ESM) in response to CO2 forcing20135859710.1002/jame.20012GiorgettaM. A.JungclausJ.ReickC. H.Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the coupled model intercomparison project phase 5201310.1002/jame.20038JungclausJ. H.FischerH. N.HaakH.Characteristics of the ocean simulations in the Max Planck institute ocean model (MPIOM) the ocean component of the MPI-earth system model20135242244610.1002/jame.20023van VuurenD. P.EdmondsJ.KainumaM.RiahiK.ThomsonA.HibbardK.HurttG. C.KramT.KreyV.LamarqueJ.-F.MasuiT.MeinshausenM.NakicenovicN.SmithS. J.RoseS. K.The representative concentration pathways: an overview201110915312-s2.0-8005390302410.1007/s10584-011-0148-zTaylorK. E.StoufferR. J.MeehlG. A.An overview of CMIP5 and the experiment design20129348549810.1175/BAMS-D-11-00094.1van VuurenD. P.StehfestE.den ElzenM. G. J.KramT.van VlietJ.DeetmanS.IsaacM.GoldewijkK. K.HofA.BeltranA. M.OostenrijkR.van RuijvenB.RCP2.6: exploring the possibility to keep global mean temperature increase below 2°C20111091951162-s2.0-8005389265910.1007/s10584-011-0152-3ThomsonA. M.CalvinK. V.SmithS. J.KyleG. P.VolkeA.PatelP.Delgado-AriasS.Bond-LambertyB.WiseM. A.ClarkeL. E.EdmondsJ. A.RCP4.5: a pathway for stabilization of radiative forcing by 21002011109177942-s2.0-8005391261910.1007/s10584-011-0151-4RiahiK.RaoS.KreyV.ChoC.ChirkovV.FischerG.KindermannG.NakicenovicN.RafajP.RCP 8.5—a scenario of comparatively high greenhouse gas emissions2011109133572-s2.0-8005390792910.1007/s10584-011-0149-yMossR. H.KathyA.EdmondsJ. A.HibbardK. A.ManningM. R.RoseS. K.Van VuurenD. P.CarterT. R.EmoriS.KainumaM.KramT.MeehlG. A.MitchellJ. F. B.NakicenovicN.RiahiK.SmithS. J.StoufferR. J.ThomsonA. M.WeyantJ. P.WilbanksT. J.The next generation of scenarios for climate change research and assessment201046372827477562-s2.0-7674909633810.1038/nature08823MuenchH. S.On the dynamics of the wintertime stratosphere circulation196522340360