A Global and Hemispherical Analysis of the Lorenz Energetics Based on the Representative Concentration Pathways Used in CMIP 5

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 CO 2 . 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 (C K ) 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 C K both vertical and horizontal eddy-transports of momentum play an important role in the increase of zonal kinetic energy for the global domain.


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][3][4][5][6][7][8][9]) and those who are considering just a piece of the atmosphere (e.g., [10][11][12][13][14][15][16][17][18][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][22][23][24][25][26][27]), and tropical cyclones (e.g., [28][29][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 CO 2 relative to the IPCC A1B scenario.
Reference [20] evaluated the atmospheric energetics in the North Hemisphere middle latitude relative to two CO 2 concentration conditions (1 × CO 2 and 2 × CO 2 ) 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 2 Advances in Meteorology for the winter season in the case of doubled CO 2 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 × CO 2 simulation.The energy cycle was characterized by a decrease in the available potential energies (  and   ), with a consequent reduction in the conversion between   to   and   to   .The values of zonal kinetic energy (  ) and the conversion between   and   increased, while   values decreased.
Hernández-Deckers and von Storch [31] based on equilibrium and transient simulations with a coupled atmosphereocean 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 CO 2 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   in the upper troposphere followed by a weakening of   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 × CO 2 and 2 × CO 2 ) 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 CO 2 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 CO 2 , 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.

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 freesurface 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/m 2 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 (CO 2 , CH 4 , and N 2 O), and air polluting scenario projections (e.g., SO 2 and NO  ) [43].To analyze the impact of the different rate of greenhouse gas emissions on the global and hemispherical , mean daily atmospheric data of zonal and meridional components of the horizontal wind vector ( and V), air temperature (), vertical component of velocity in pressure coordinates ( = /), and geopotential height () 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]. Consider where  is the amount of available potential energy,  its generation rate,  is the energy conversion between  and kinetic energy , and  is the dissipation rate of .When (, ) is positive,  increases at the expense of , and vice versa.
Figure 1: Lorenz energy cycle representation for (a) two-component energy cycle in which available potential energy denoted by the letter  is generated (), converted () to kinetic energy (), and subsequently dissipated ().(b) Four-component energy cycle with  and  decomposed into mean and eddy parts.Adapted from [33].
The four-component energy cycle is characterized by the subdivision of  and  into their zonal and eddy components and the set of equations can be written as In the four-component energy cycle the potential energy is portioned into zonal available potential energy (  ) and eddy available potential energy (  ).The kinetic energy can be treated similarly as potential energy; that is, the kinetic energy of the flow can be separated into its zonal (  ) and eddy (  ) parts.In the energy cycle of Lorenz,   and   are produced by generation processes named   and   , while reservoirs of energy are connected by conversion terms   ,   ,   , and   (resp, the conversion from   to   ,   to   ,   to   , and   to   ).The   term involves transport of sensible heat by meridional circulation affecting   and   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   at the expense of   when there is warm air rising and cold air sinking in the same latitude.The   involves horizontal and vertical eddy-transports of sensible heat.This term affects   and   by altering the variance of temperature and depends upon the transport of sensible heat along the gradients of zonal averaged temperature.The   term involves vertical eddy-transports of sensible heat, whose values by unit area are independent of latitude.It affects   Consider ,   ,   , and   terms, respectively, represent fluxes of   ,   ,   , and   from the lateral and vertical boundaries.Φ and Φ terms denote the dynamical mechanisms in which they produce or destroy kinetic energy.As described by [44], this is due to the fact that Φ, Φ,   , and   derivatives involve a single term in the form ⋅∇Φ which represents the appearance of kinetic energy via a crossisobaric flow towards low pressure.The destruction of kinetic energy occurs when there is cross-isobaric flow towards the high pressure.However, here Φ and Φ 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   and   in the following form: 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: 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.

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   are maximized in May and September followed by two minimum during June and January (Figure 2(a)).The results show diminishing of   as greenhouse gases concentration increases.An opposite behavior is observed for   (Figure 2(b)).In this case   presents high negative values in the control experiment and declines to minimum negative values for RCP85 experiment.The implications of the fall of   and   (negatively) as the concentration of greenhouse gases increases are that less   and   are created by diabatic processes and consequently affecting the release of   and   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].
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 (  and   ) 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   as the concentration of greenhouse gases increases is because less   is destructed or converted into   .The destruction of   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 (  and   ) as function of the RCPs and CRTL experiments.The annual cycle of   is characterized by a well-defined seasonal variation, with highest (lowest) values from July to September (April to June).One can see that   clearly increases, while   remains approximately unchanged around the year as the concentration of greenhouse gases increases.The increase of   occurs in the RCP85 forcing scenario, with minor and almost equally changes in RCP45 and RCP26.As the generation (or destruction) of   depends upon the joint effect of the transport of angular momentum along the gradients of angular velocity (  ), transport of sensible heat by meridional circulation (  ) and the dissipation of kinetic energy by friction, represented by   (2c), and an evaluation of the conversion terms are necessary.
The global climatological annual cycles of the verticallyintegrated conversion and dissipation terms for the control and RCPs experiments are presented in Figure 3.It can be seen that   and   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   to   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 CO 2 concentration, the meridional temperature gradient decreases due to the anomalous heating in the tropical and polar regions.The changes in   and   are modest and not statistically significant (Figures 3(a       that in a warmed climate the conversion from   to   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   and   during the course of the year; that is,   presents the highest values from July to September, when   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.4(a)) shows an increase in the generation term  as a response to the increasing in greenhouse gases concentration, followed by a decrease in the available potential energy , conversion term (, ) and the dissipation .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  is ∼61% (44%, 23%) in the case of RCP85 (RCP45 and RCP26).The others components have the following conditions: decreasing in  of ∼28% (16%, 8%), decreasing in (, ) of ∼4% (2%, 0.2%), increasing in  of ∼15% (7%, 4%), and a decreasing in the dissipation term  of ∼4% (1%, 0.2%).

Volume-Integrated Energy
The four-component energy cycle (Figure 4(b)) indicates a general decrease in the values of generation, available potential energy, the conversion terms   ,   , and   , 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   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   versus a small increase of ∼1% in   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   to   (  ) and   to   (  ), is reduced in a condition of warmed climate, as can be seen by the reduction of ∼39% in   and ∼1.8% in   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   and   (  ), is increased showing that the zonal flow grows at the expense of the eddy motion.The conversion term   , 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  domain displayed in Figure 5(a) shows an increase in the generation term , followed by a decrease in the available potential energy , conversion term (, ), and the dissipation .The figure still shows a negative decrease in the flux of (  ) and an increase in the flux of (  ).Thus,   acts as sink of , and it is less intense for a warmed climate, while  acts as a source of   .As can be observed, the increase in , in the case of  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   .The decrease in (, ) is ∼9% (5%, 3%), with an increasing in  of ∼17% (8%, 5%), an increase of ∼66% (25%, 8%) in the   and a decrease in the dissipation term  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   ,   , and   , and the dissipation terms as the concentration of greenhouse gases increases.Kinetic energy and the conversion term   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   , as the reduction in   is one order magnitude small.The increase in the kinetic energy boundary terms is particularly due to the   term.The energy reservoirs of available energies and kinetic energies, respectively, decrease and increase.For the RCP85 experiment the increase of   is small (∼1%) when compared with increases in   (∼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   and   , decreases as the planet warms, followed by an increase in the barotropic component, and a reduction in the   values.These patterns are similar to the global domain.
Figure 6(a) shows the two-component energy cycle for the  domain.The results show an increase in the generation term , a decrease in the available potential energy , a decrease in   , a decrease in the conversion term (, ), an increase in , a negative increase in   , and a decrease in the dissipation term .In the  case,   acts as source of , and it is less intense for a warmed climate, while   acts as a sinking of   , and it is more intense for a warmed climate.As can be seen, the increase in , in the case of  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   .The decrease in the values of (, ) is ∼10% (6%, 1%), with an increasing in  of ∼11% (6%, 2%), a negative increase of ∼64% (42%, 28%) in the   , and a decrease in the dissipation term  of ∼10% (3%, increase of 1%).
In the four-component energy cycle for the  (Figure 6(b)), observed is a pattern similar to the global and  domain, that is, reduction in the values of available potential energies, followed by a respective small and a significant increase in   and   , reduction in the generation terms, reduction in the baroclinic flux, increase in   , 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   , as the reduction in   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   for the global domain is ∼21.8% relative to the control experiment, while   in the  and  increases about 24.6% and 18.0%, respectively.The same pattern of increase in   (and the remaining reservoirs) is followed by the others RCPs projection, though with less intensity.
For all analyzed domains the intensity of the , 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  (10.35%), when compared with global (4.68%) and  (9.57%) domains.The strongest increase is seen to occur in the total kinetic energy for the  (∼17%), followed by an increase of 41.02% in   .This result indicates that the increase in   is related to the increase in   , as   acts as a source of   .However, another form to attribute the increase of   due to   is to quantifying what dynamic process was the most important to the increase of   .This can be achieved by evaluating or quantifying the terms involved in the calculation of   , see equation in the appendix.The first three terms of   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  (, Global) domain, the increase in   , 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 eddytransport of momentum relative to control experiment.For the experiments RCP26 and RCP45 the increase in the horizontal (vertical) eddy-transport of the  is 22.8% (8.4%) and 29.7% (10.1%), respectively, showing that in the  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) eddytransport of momentum is 18.7% (6.1%) and 25.0% (9.6%), respectively, and for the  domain, it corresponds to −1.5% (2.1%) and 2.2% (9.5%), respectively.
Calculations of changes for individual terms of   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 Advances in Meteorology change is of 77.7% (55.5% and 33.3%).For the  and  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  and  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   term is below 3% relatively to the control experiment for all RCPs.

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 .The energetic analysis is based on the RCP85, RCP45, and RCP26 forcing experiments and a control simulation from the MPI-ESM-MR oceanatmosphere 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  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 CO 2 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  experiment, with the strongest reduction occurring in the highest radiative forcing experiment, RCP85.The results also show an increase in   ,   , and   components of the .Similar reduction in the intensity of the  is observed for the  and , with RCP85 emission scenario influencing the most intense decrease in the energy terms.In Particular, for this scenario the energetic projection for the  indicates an increase in   of 24.6%, which is higher than the global (21.8%) and  (18.0%) domains.The anomalous increase of   in the  ( 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   and   , in the RCPs scenarios.These results contribute to the findings of [20,31] who have shown that   increases in a warmed climate; however, our results indicates that the increase in   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.northward components of the wind, respectively, and  is the static stability parameter which is given by and   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   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   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.
) and 3(c)), suggesting a weak relation to the increasing of greenhouse gases concentration.The term   , which acts around the year as a source of energy for the reservoir of   , shows different result.It indicates

Figure 2 :
Figure 2: Climatology of the generation terms   (a) and   (b), and the energy terms   (c),   (d),   (e), and   (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 10 5 J/m 2 for reservoirs and W/m 2 for generation terms.

Figure 4 :
Figure 4: Representation of annual mean of spatial domain Lorenz energy cycle for the entire globe in which (a) is representative of a twocomponent 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 10 5 J/m 2 for reservoirs and W/m 2 for conversion, generation, and dissipation terms.In the two-component energy cycle the value of  is given by the sum of its zonal and eddy part ( =   +   ),  =   +   , (, ) =   +   ,  =   +   , and  =   +   .Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.

Figure 5 :
Figure 5: 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 (  ) and kinetic energy (  ) 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 10 5 J/m 2 for reservoirs and W/m 2 for conversion, generation, and dissipation terms.Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.

Figure 6 :
Figure 6: 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 (  ) and kinetic energy (  ) 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 10 5 J/m 2 for reservoirs and W/m 2 for conversion, generation, and dissipation terms.Arrows indicate the direction corresponding to positive values; negative values imply opposite flux direction.