 \infty$) give the corresponding values of the parameters of the zero dimensional energy balance model, $A$, $B$, $a$, $b$; and the effective inverse response times for the atmosphere and the surface. The first important thing to notice with model (\[model\]) is that it recovers the zero dimensional energy balance model. This can be seen by averaging Eq. (\[model\]) over the cells in the limit $N \to \infty$. The following reduced version of Eq. (\[model\]) is obtained: $$\begin{aligned}
\label{motion}
\dot{y}=-\lambda [y-(A+ax)]+\varepsilon(t), \\ \nonumber
\dot{x}=-\gamma [x-(B+by)].\end{aligned}$$

From the fact that Eq. (\[motion\]) represents an overdamped dynamics in a parabolic potential, the system converges to a stationary state with mean value given by the solution of the zero dimensional energy balance model. This result shows that model (\[model\]) is capable of representing with good approximation the mean behavior of energy in the coupled atmosphere – surface system. Now it will be argued that Eq. (\[model\]) can give realistic descriptions of temperature fluctuations as well. A large variability over the intrinsic characteristic times of the surface around the globe is expected to exist, as a consequence of the different response times present in the geosphere, hydrosphere, cryosphere and biosphere. At first instance, this situation is modeled by treating the $\gamma _n$’s like independent random variables taken from a uniform probability distribution. On the other hand, under the basis of the relative homogeneity of the atmosphere composition, it will be assumed by now that the response time to perturbations of the atmosphere is the same in all cells. In what follows, the value $\lambda_n = \lambda = 1$ is used. Under these assumptions, the time unit is defined as the mean atmospheric response time. The model represents the interaction between a hierarchy of time scales, ranging from minutes to days to geological times. From this point of view, it would be therefore reasonable to assume that $1/ \lambda$ lies in the intermediate scales, roughly in the range from weeks to months. This and other important aspects about the definition of the parameters in the model are intended to be refined by the author in the near future by close cooperation with climate experts. At this point the main goal is limited to explore the capabilities of the model (\[model\]), in order to give qualitatively realistic statistical descriptions of the temperature fluctuations
\infty$) give the corresponding values of the parameters of the zero dimensional department of energy libra model, $ A$, $ B$, $ a$, $ b$; and the effective inverse response time for the air and the surface. The first significant thing to notice with exemplar (\[model\ ]) is that it recovers the zero dimensional energy libra model. This can be seen by average Eq. (\[model\ ]) over the cell in the limit $ N \to \infty$. The following reduced version of Eq. (\[model\ ]) is obtained: $ $ \begin{aligned } 
 \label{motion } 
 \dot{y}=-\lambda [ y-(A+ax)]+\varepsilon(t), \\ \nonumber 
 \dot{x}=-\gamma [ x-(B+by)].\end{aligned}$$ 

 From the fact that Eq. (\[motion\ ]) represent an overdamped dynamics in a parabolic potential, the system converge to a stationary state with mean value given by the solution of the zero dimensional energy proportion model. This result shows that model (\[model\ ]) is capable of map with good approximation the mean behavior of energy in the coupled atmosphere – surface system. Now it will be argued that Eq. (\[model\ ]) can give realistic descriptions of temperature fluctuations as well. A large unevenness over the intrinsic characteristic fourth dimension of the surface around the ball is expect to exist, as a consequence of the different response times present in the lithosphere, hydrosphere, cryosphere and biosphere. At first case, this situation is modeled by treat the $ \gamma _ n$ ’s like independent random variables taken from a uniform probability distribution. On the other hand, under the footing of the relative homogeneity of the atmosphere composition, it will be assumed by now that the response time to disruption of the atmosphere is the same in all cells. In what follows, the value $ \lambda_n = \lambda = 1 $ is used. Under these assumption, the time unit is defined as the base atmospheric response time. The model represents the interaction between a hierarchy of prison term scales, crop from minutes to days to geological fourth dimension. From this point of position, it would be therefore fair to assume that $ 1/ \lambda$ lies in the average scales, roughly in the range from weeks to months. This and early important aspects about the definition of the parameters in the mannequin are intended to be refined by the author in the near future by close cooperation with climate expert. At this point the main finish is limited to explore the capabilities of the mannequin (\[model\ ]), in order to give qualitatively realistic statistical description of the temperature fluctuations
 \invty$) give the correspondikg values of the paramevers of the zeru dimensional energy balance midel, $Q$, $B$, $a$, $b$; and the effectkve inverde respobse uimes for the atmosphere and ths suryare. The first imkortant thinc to notice widh mldel (\[model\]) is that it recovers the sero dikejsional energy balsgce jodel. This can be seen by averaginf Eq. (\[movel\]) over the cekls in the limit $N \to \inftj$. Thf following reducef version od Eq. (\[nodel\]) is obtxined: $$\begin{aligned}
\labem{motion}
\dot{y}=-\lambda [y-(A+ax)]+\varepsilov(t), \\ \npnumber
\dot{z}=-\gqmmw [x-(B+by)].\end{alijned}$$

Frjm the fact bnat Eq. (\[motion\]) represents an ovxrdanped dynamics in a pacabolic potential, thg system cmnrerges to a stationart wtate witv mexb vxlut gmveh by tje aolution or the zero eimensional energy najqnce model. Thjs restle shows that model (\[model\]) is capable of geprssenting with good appriximation the mean bejavior of energy in the coupled atmosphere – surface system. Now mt wiol cw wrgued that Eq. (\[model\]) can give realistic descryltooks of temperaturc fluctuations as eepl. W large variacility ovsr the intrinsic cjaractetistic times of the surface around the globe iw expected tj exist, as a conseqbence of the difgereny response times presenc in tge geospherf, hydrospgdre, cryosphere avd niovphere. Au first instance, tris situavion ns modeldd bi treatyng the $\galma _n$’s like independent gandol eariables haken from a uniform probabilitb distributiom. Mn nhe other hand, under the bafis of the relctive hoiogendity of thv atmosphxre composityon, it will ba assumed by iow that ehe eespinse tiod to perturbatoons of tkt atmospheee is the same in ell zslls. In what fooloqs, the value $\lakbdx_n = \lwmuda = 1$ hs used. Undes thdse sssumotions, uke bimd unot is defined as the meah atmospheric resppnfe time. Jhe model represents tne interaction betaeen e hiererchy pf jime scales, ranging from minutea to days to geological eimew. From this 'oint of view, it would be therefore reasinable to assume jhat $1/ \lambda$ lies in jhe intermediave scajes, roughny in the range from weeks to months. Uhis and other importanf aspewts ahout the definition of the parameters in the model are intended to be refuned bb ehe author jn tne nedr yuture fy cpuse cooperation eith climate experts. At this pomnt the mahn goal is limited to explore yhd capabilitier of the model (\[model\]), in krder to give qualitatively realistic syatistical descriptions pf the twmperzture flucvuations
\infty$) give the corresponding values of the the dimensional energy model, $A$, $B$, inverse times for the and the surface. first important thing to notice with (\[model\]) is that it recovers the zero dimensional energy balance model. This can seen by averaging Eq. (\[model\]) over the cells in the limit $N \to The reduced of (\[model\]) is obtained: $$\begin{aligned} \label{motion} \dot{y}=-\lambda [y-(A+ax)]+\varepsilon(t), \\ \nonumber \dot{x}=-\gamma [x-(B+by)].\end{aligned}$$ From the fact that Eq. represents an overdamped dynamics in a parabolic potential, system converges to a state with mean value given the of the dimensional balance This result shows model (\[model\]) is capable of representing with good approximation the mean behavior of energy in the coupled – surface it will argued Eq. can give realistic temperature fluctuations as well. A large intrinsic characteristic times of the surface around the is expected exist, as a consequence of the response times present in the geosphere, hydrosphere, cryosphere biosphere. At first instance, this situation is modeled by treating the $\gamma _n$’s like independent taken from a uniform distribution. On the hand, the of relative homogeneity the atmosphere composition, it will be assumed by now that the time to perturbations of the atmosphere is the same in In follows, the value = \lambda = 1$ used. these assumptions, the time defined the time. model the interaction between a of time scales, ranging from to days to geological view, it would be therefore reasonable to assume $1/ \lambda$ lies in the intermediate scales, in the range from weeks to months. This and other important aspects the definition parameters in the model are intended to be by the author in near future by close cooperation with climate experts. At point main goal limited to explore capabilities of the (\[model\]), in order qualitatively realistic descriptions the
 \infty$) give the corresponding Values of thE paraMetErs Of The zEro dImensional enerGY balAnce model, $A$, $B$, $a$, $b$; and the effEctivE iNVersE ReSponsE times fOR tHE AtmOsPhEre AnD ThE surfAce. the firsT important ThiNg To notice with MOdEl (\[model\]) is tHat It recovers thE zeRo dimeNsIonAL enerGy bAlancE model. tHis can Be seen by aVeRAging EQ. (\[Model\]) ovER ThE celLs in the limit $N \to \inFTy$. tHe following redUced veRsIOn OF eq. (\[mOdeL\]) is obtaineD: $$\bEgin{aLIgned}
\laBEl{MOTIon}
\DOt{y}=-\lambda [y-(A+ax)]+\Varepsilon(t), \\ \NOnuMber
\doT{x}=-\GamMA [x-(B+by)].\eNd{aliGnED}$$

FrOm the fact thAt Eq. (\[Motion\]) repResentS An overdAMped dynAmics iN a pAraBoliC PoTeNtiAl, THe sYStEm cONveRges to a sTaTiOnary StatE WITH meaN vaLue gIven bY the solution oF thE zerO DimEnsioNal enErgy BaLance Model. THis reSuLt shows that modeL (\[modEl\]) is capabLe oF rEprEsEntinG With goOd aPprOximatiOn the meAN beHaVIOR oF energy in the coupleD aTMOsPhere – surFace sySTeM. NOW it will bE aRguEd thAT eq. (\[modEl\]) caN GiVe realisTic desCRiPtIons of tEmPeratuRe FluCtuAtionS As weLl. A larGe variabIlity OVer the intrinsiC CharacteristiC TiMES oF The sUrfAce around thE gloBE is eXpecTEd To eXIst, as A consEqUEnCE of the different respOnSe timeS presEnt in the geospHere, hydrosPHERe, cryospHere ANd BIosphere. At firsT instAnce, this siTUation is ModelEd by treaTing the $\gaMMA _n$’s like iNdePenDenT raNDOm Variables takeN FRom a UnIform prObaBility dIstRibUtiOn. ON tHe other haNd, under tHe BaSiS oF thE relaTIve homogEnEitY oF thE atmoSPhere cOmposItioN, iT wILl bE assumeD By NOW thaT tHe RespOnsE tIme to PertURbaTions of The atmospHerE Is thE sAmE in all cElls. In what folLoWs, the value $\LaMbdA_n = \lambDA = 1$ Is used. UnDer these assumptions, the tIMe unit iS deFined As thE mean atmoSphEric reSpoNSe time. the modEl repReSenTS The inTERaCtiOn Between a hiERArcHy of tImE scaLes, rangIng from minutes to daYS to Geological timEs. FRom tHIS pOinT Of VIew, It WOulD BE therefore reasoNable to assUmE ThAt $1/ \lambda$ liES in ThE intermEdiate sCales, ROughly iN the range From weeks To MontHS. thiS and other iMportant Aspects abOUt the DEfInitiOn oF the paRaMetErs in The modEL arE inteNded to Be RefineD by thE aUthor in tHe near future by close coopEratioN with CliMate experTs. AT ThiS point the Main Goal is limiTed To eXplorE thE CapabIlitIEs Of tHE modeL (\[modEL\]), in order tO GiVe qUALiTatively reaLISTic StatiStiCAl descRiptIons of the temperatURe fluctuations
 \infty$) give the corresp onding val ues o f t hepa rame ters  of the zero d i mens ional energy balance m odel,  $ A $, $ B $,  $a$,  $b$; a n dt h e e ff ec tiv ei nv erseres ponse t imes for t heat mosphere and th e surface.  Th e first impo rta nt thi ng  to notic e w ith m odel ( \ [model \]) is th at it rec o vers th e ze ro d imensional energy ba l ance model. Th is can  b e  s e e n b y a veraging E q.  (\[m o del\])o ve r t hec ells in the l imit $N \to \in fty$.Th e f o llowin g red uc e d v ersion of E q. ( \[model\] ) is o b tained: $$\begi n{alig ned }
\ labe l {m ot ion }\ dot { y} =-\ l amb da [y-(A +a x) ]+\va reps i l o n (t),  \\  \no numbe r
\dot{x}=-\g amm a [x - (B+ by)]. \end{ alig ne d}$$
Fromthe f ac t that Eq. (\[m otio n\]) repr ese nt s a noverd a mped d yna mic s in aparabol i c p ot e n t ia l, the system conv er g e sto a sta tionar y  s ta t e with m ea n v alue g ivenby t h esolution  of th e  z er o dimen si onal e ne rgy  ba lance mode l. Thi s result  show s  that model (\ [ model\]) is c a pa b l eo f re pre senting wit h go o d ap prox i ma tio n  themeanbe h av i or of energy in the  c oupled  atmo sphere – surf ace system . N ow it wi ll b e  a r gued that Eq.(\[mo del\]) can give rea listi c descri ptions of t emperatu reflu ctu ati o n sas well. A la r g e va ri ability  ov er theint rin sic  ch ar acteristi c timesof  t he  s urf ace a r ound the  g lob eisexpec t ed toexist , as  a  c o nse quenceo ft h e di ff er entres po nse t imes pre sent in  the geos phe r e, h yd ro sphere,  cryosphere a nd  biosphere .Atfirsti n stance,this situation is model e d by tr eat ing t he $ \gamma _n $’s  likeind e penden t rand om va ri abl e s  take n fr omauniform pr o b abi litydi stri bution.  On the other hand ,  un der the basis  of  the r el ati v eh omo ge n eit y of the atmosphe re composi ti o n,  it will b e  as su med bynow tha t the respons e time to  perturba ti onso f  th e atmosphe re is th e same in all c e ll s. In  wh at fol lo ws,  thevalue$ \la mbda_ n = \l am bda =1$ is  u sed. Und er these assumptions, t he tim e uni t i s defined  as the  mean atm osph eric respo nse  ti me. T hem odelrepr e se nts the i nter a ction bet w ee n a h ie rarchy of t i m e  sc ales,  ra n ging f romminutes to days t o  geological ti mes. F rom  th i s po in t of view, itwou ld b e theref or e reasonabl e to ass um e  that  $1/ \ lambda $ liesi n  t h e inte rmed iat e scales,  ro ug h ly in t he  r a nge fr om w ee ks tomonths .  Thi s and other import ant a s p ectsa bou t the  d efiniti o n of  the param eters in th e mode l ar e int ended t obe ref ine dby the aut h or in the  near  future  b y cl ose  coope rati o n  with  cli ma teexperts.A t  t h is  p o int  the  main  g oalis limite d  to expl ore the cap ab ili t i es oft he m odel (\[mo del \]),i n  order tog iveq ua l itati vely r ealist ic stat i sti ca l descr ipt i o ns of the  temperat u refl uctu ations
 \infty$)_give the_corresponding values of the_parameters of_the_zero dimensional_energy_balance model, $A$,_$B$, $a$, $b$;_and the effective inverse_response times for_the_atmosphere and the surface. The first important thing to notice with model (\[model\]) is_that_it recovers_the_zero_dimensional energy balance model. This_can be seen by averaging_Eq. (\[model\])_over the cells in the limit $N \to_\infty$._The following reduced_version of Eq. (\[model\]) is obtained: $$\begin{aligned}
\label{motion}
\dot{y}=-\lambda [y-(A+ax)]+\varepsilon(t), \\_\nonumber
\dot{x}=-\gamma [x-(B+by)].\end{aligned}$$

From the fact that Eq._(\[motion\]) represents an_overdamped_dynamics_in a parabolic potential,_the system converges to a stationary_state with mean value given by_the solution of the zero dimensional energy_balance model. This result shows that_model (\[model\]) is capable of_representing with_good approximation the mean behavior_of energy in_the coupled_atmosphere – surface_system. Now it will be argued_that Eq. (\[model\])_can give realistic descriptions of temperature_fluctuations_as well. A_large_variability_over the_intrinsic characteristic times_of_the surface_around_the globe is expected to exist,_as_a consequence of the different response times_present in the geosphere,_hydrosphere,_cryosphere and biosphere. At_first instance, this situation is_modeled by treating the $\gamma _n$’s_like independent_random variables_taken from a uniform probability distribution. On the other hand, under_the basis of the relative homogeneity_of the atmosphere composition,_it will_be_assumed by now_that_the response_time to perturbations of the atmosphere is_the same_in all cells. In what follows,_the value $\lambda_n =_\lambda_= 1$ is used. Under these_assumptions, the time unit is defined_as the mean atmospheric response_time._The_model represents the interaction between_a hierarchy of time scales, ranging_from minutes to_days to geological times. From this point_of_view, it would be therefore reasonable_to_assume that $1/ \lambda$ lies in_the_intermediate_scales, roughly in the range_from weeks to months. This and_other important aspects about the definition of the parameters_in the model_are intended to be refined_by_the_author in the near future by close cooperation with climate_experts. At_this point the_main goal is limited to explore the capabilities of the_model (\[model\]), in order to give qualitatively_realistic statistical descriptions of the temperature fluctuations