Global and Regional Temperature Projections Using the Fractional Energy Balance Equation

Author(s):  
Roman Procyk ◽  
Shaun Lovejoy ◽  
Raphaël Hébert ◽  
Lenin Del Rio Amador

<p>We present the Fractional Energy Balance Equation (FEBE): a generalization of the standard EBE.  The key FEBE novelty is the assumption of a hierarchy of energy storage mechanisms: scaling energy storage.  Mathematically the storage term is of fractional rather than integer order.  The special half-order case (HEBE) can be classically derived from the continuum mechanics heat equation used by Budyko and Sellers simply by introducing a vertical coordinate and using the correct conductive-radiative surface boundary conditions (the FEBE is a mild extension).</p><div> <p>We use the FEBE to determine the temperature response to both historical forcings and to future scenarios.  Using historical data, we estimate the 2 FEBE parameters: its scaling exponent (<em>H</em> = 0.38±0.05; <em>H</em> = 1 is the standard EBE) and relaxation time (4.7±2.3 years, comparable to box model relaxation times). We also introduce two forcing parameters: an aerosol re-calibration factor, to account for their large uncertainty, and a volcanic intermittency exponent so that the intermittency volcanic signal can be linearly related to the temperature. The high frequency FEBE regime not only allows for modelling responses to volcanic forcings but also the response to internal white noise forcings: a theoretically motivated error model (approximated by a fractional Gaussian noise). The low frequency part uses historical data and long memory for climate projections, constraining both equilibrium climate sensitivity and historical aerosol forcings. <span>Parameters are estimated in a Bayesian framework using 5 global observational temperature series, and an error model which is a theoretical consequence of the FEBE forced by a Gaussian white noise.</span></p> <p>Using the CMIP5 Representative Concentration Pathways (RCPs) and CMIP6 Shared Socioeconomic Pathways (SSPs) scenario, the FEBE projections to 2100 are shown alongside the CMIP5 MME. The Equilibrium Climate Sensitivity = 2.0±0.4 <sup>o</sup>C/CO<sub>2</sub> doubling implies slightly lower temperature increases.   However, the FEBE’s 90% confidence intervals are about half the CMIP5 size so that the new projections lie within the corresponding CMIP5 MME uncertainties so that both approaches fully agree.   The mutually agreement of qualitatively different approaches, gives strong support to both.  We also compare both generations of General Circulation Models (GCMs) outputs from CMIP5/6 alongside with the projections produced by the FEBE model which are entirely independent from GCMs, contributing to our understanding of forced climate variability in the past, present and future.</p> <p>Following the same methodology, we apply the FEBE to regional scales: estimating model and forcing parameters to produce climate projections at 2.5<sup>o</sup>x2.5<sup>o</sup> resolutions. We compare the spatial patterns of climate sensitivity and projected warming between the FEBE and CMIP5/6 GCMs. </p> </div>

2020 ◽  
Author(s):  
Roman Procyk ◽  
Shaun Lovejoy ◽  
Raphael Hébert

Abstract. We produce climate projections through the 21st century using the fractional energy balance equation (FEBE) which is a generalization of the standard EBE. The FEBE can be derived either from Budyko–Sellers models or phenomenologically by applying the scaling symmetry to energy storage processes. It is easily implemented by changing the integer order of the storage (derivative) term in the EBE to a fractional value near 1/2. The FEBE has two shape parameters: a scaling exponent H and relaxation time τ; its amplitude parameter is the equilibrium climate sensitivity (ECS). Two additional parameters were needed for the forcing: an aerosol re-calibration factor α to account for the large aerosol uncertainty, and a volcanic intermittency correction exponent ν. A Bayesian framework based on historical temperatures and natural and anthropogenic forcing series was used for parameter estimation. Significantly, the error model was not ad hoc, but was predicted by the model itself: the internal variability response to white noise internal forcing. The 90 % Confidence Interval (CI) of the shape parameters were H = [0.33, 0.44] (median = 0.38), τ = [2.4, 7.0] (median = 4.7) years compared to the usual EBE H = 1, and literature values τ typically in the range 2–8 years. We found that aerosols were too strong by an average factor α = [0.2, 1.0] (median = 0.6) and the volcanic intermittency correction exponent was ν = [0.15, 0.41] (median = 0.28) compared to standard values α = ν = 1. The overpowered aerosols support a revision of the global modern (2005) aerosol forcing 90 % CI to a narrower range [−1.0, −0.2] W m−2 compared with the IPCC AR5 range [1.5, 4.5] K (median = 3.2 K). Similarly, we found the transient climate sensitivity (TCR) = [1.2, 1.8] K (median = 1.5 K) compared to the AR5 range TCR = [1.0, 2.5] K (median = 1.8 K). As commonly seen in other observational-based studies, the FEBE values are therefore somewhat lower but still consistent with those in IPCC AR5. Using these parameters we made projections to 2100 using both the Representative Carbon Pathways (RCP) and Shared Socioeconomic Pathways (SSP) scenarios and shown alongside the CMIP5/6 MME. The FEBE hindprojections (1880–2019) closely follow observations (notably during the hiatus, 1998–2015). Overall the FEBE were 10–15 % lower but due to their smaller uncertainties, their 90 % CIs lie completely within the GCM 90 % CIs. The FEBE thus complements and supports the GCMs.


2021 ◽  
Vol 12 (2) ◽  
pp. 489-511
Author(s):  
Shaun Lovejoy

Abstract. In Part 1, I considered the zero-dimensional heat equation, showing quite generally that conductive–radiative surface boundary conditions lead to half-ordered derivative relationships between surface heat fluxes and temperatures: the half-ordered energy balance equation (HEBE). The real Earth, even when averaged in time over the weather scales (up to ≈ 10 d), is highly heterogeneous. In this Part 2, the treatment is extended to the horizontal direction. I first consider a homogeneous Earth but with spatially varying forcing on both a plane and on the sphere: the new equations are compared with the canonical 1D Budyko–Sellers equations. Using Laplace and Fourier techniques, I derive the generalized HEBE (the GHEBE) based on half-ordered space–time operators. I analytically solve the homogeneous GHEBE and show how these operators can be given precise interpretations. I then consider the full inhomogeneous problem with horizontally varying diffusivities, thermal capacities, climate sensitivities, and forcings. For this I use Babenko's operator method, which generalizes Laplace and Fourier methods. By expanding the inhomogeneous space–time operator at both high and low frequencies, I derive 2D energy balance equations that can be used for macroweather forecasting, climate projections, and studying the approach to new (equilibrium) climate states when the forcings are all increased and held constant.


2021 ◽  
Author(s):  
Shaun Lovejoy

<p>The highly successful Budyko-Sellers energy balance models are based on the classical continuum mechanics heat equation in two spatial dimensions. When extended to the third dimension using the correct conductive-radiative surface boundary conditions, we show that surface temperature anomalies obey the (nonclassical) Half-order energy balance equation (HEBE, with exponent H = ½) implying heat is stored in the subsurface with long memory. </p><p> </p><p>Empirically, we find that both internal variability and the forced response to external variability are compatible with H ≈ 0.4.  Although already close to the HEBE and classical continuum mechanics, we argue that an even more realistic “effective media” macroweather model is a generalization: the fractional heat equation (FHE) for long-time (e.g. monthly scale anomalies).  This model retains standard diffusive and advective heat transport but generalize the (temporal) storage term.  A consequence of the FHE is that the surface temperature obeys the Fractional EBE (FEBE), generalizing the HEBE to 0< H ≤1.  We show how the resulting FEBE can be been used for monthly and seasonal forecasts as well as for multidecadal climate projections.  We argue that it can also be used for understanding and modelling past climates.</p>


2020 ◽  
Author(s):  
Roman Procyk ◽  
Shaun Lovejoy ◽  
Lenin Del Rio Amador

<p>The conventional energy balance equation (EBE) is a first order linear differential equation driven by solar, volcanic and anthropogenic forcings.  The differential term accounts for energy storage usually modelled as one or two “boxes”.  Each box obeys Newton’s law of cooling, so that when perturbed, the Earth’s temperature relaxes exponentially to a thermodynamic equilibrium.</p><p>However, the spatial scaling obeyed by the atmosphere and its numerical models implies that the energy storage process is a scaling, power law process, a consequence largely of turbulent, hierarchically organized oceans currents but also hierarchies of land ice, soil moisture and other processes whose rates depend on size.</p><p>Scaling storage leads to power law relaxation and can be modelled via a seemingly trivial change - from integer to fractional order derivatives - the Fractional EBE (FEBE): with temperature derivatives order 0 < H  < 1 rather than the EBE value H = 1.  Mathematically the FEBE is a past value problem, not an initial value problem.    Recent support for the FEBE comes from [Lovejoy, 2019a] who shows that the special H = 1/2 case (close to observations), the “Half-order EBE” (HEBE), can be analytically obtained from classical Budyko-Sellers energy balance models by improving the boundary conditions.</p><p>The FEBE simultaneously models the deterministic forced response to external (e.g. anthropogenic) forcing as well as the stochastic response to internal forcing (variability) [Lovejoy, 2019b].  We directly exploit both aspects to make projections based on historical data estimating the parameters using Bayesian inference.  Using global instrumental temperature series, alongside CMIP5 and CMIP6 standard forcings, the basic FEBE parameters are H ≈ 0.4 with a relaxation time ≈ 4 years.  </p><p>This observation-based model also produces projections for the coming century with forcings prescribed by the CMIP5 Representative Concentration Pathways scenarios and the CMIP6 Shared Socioeconomic Pathways.</p><p>We compare both generations of General Circulation Models (GCMs) outputs from CMIP5/6 alongside with the projections produced by the FEBE model which are entirely independent from GCMs, contributing to our understanding of forced climate variability in the past, present and future.  When comparing to CMIP5 projections, we find that the mean projections are about 10- 15% lower while the uncertainties are roughly half as large.  Our global temperature projections are therefore within the  CMIP5 90% confidence limits and thus give them strong, independent support.</p><p> </p><p><strong>References</strong></p><p>Lovejoy, S., The half-order energy balance equation, J. Geophys. Res. (Atmos.), (submitted, Nov. 2019), 2019a.</p><p>Lovejoy, S., Fractional Relaxation noises, motions and the stochastic fractional relxation equation Nonlinear Proc. in Geophys. Disc., https://doi.org/10.5194/npg-2019-39, 2019b.</p>


2020 ◽  
Author(s):  
Shaun Lovejoy ◽  
Lenin Del Rio Amador ◽  
Roman Procyk

<p>The conventional 1-D energy balance equation (EBE) has no vertical coordinate so that radiative imbalances between the earth and outer space are redirected in the horizontal in an ad hoc manner.  We retain the basic EBE but add a vertical coordinate so that the imbalances drive the system by imposing heat fluxes through the surface.   While this is theoretically correct, it leads to (apparently) difficult mixed boundary conditions.  However, using Babenko’s method, we directly obtain simple analytic equations for (2D) surface temperature anomalies T<sub>s</sub>(x,t): the Half-order Energy Balance Equation (HEBE) and the Generalized HEBE, (GHEBE) [Lovejoy, 2019a].  The HEBE anomaly equation only depends on the local climate sensitivities and relaxation times.  We analytically solve the HEBE and GHEBE for T<sub>s</sub>(x,t) and provide evidence that the HEBE applies at scales >≈10km.  We also calculate very long time diffusive transport dominated climate states as well as space-time statistics including the cross-correlation matrix needed for empirical orthogonal functions.</p><p>The HEBE is the special H = 1/2 case of the Fractional EBE (FEBE) [Lovejoy, 2019b], [Lovejoy, 2019c] and has a long (power law) memory up to its relaxation time t.  Several papers have empirically estimated H ≈ 0.5, as well as t ≈ 4 years, whereas the classical zero-dimensional EBE has H = 1 and t ≈ 4 years.   The former values permit accurate macroweather forecasts and low uncertainty climate projections; this suggests that the HEBE could apply to time scales as short as a month.  Future generalizations include albedo-temperature feedbacks and more realistic treatments of past and future climate states.</p><p><strong>References</strong></p><p> </p><p>Lovejoy, S., The half-order energy balance equation, J. Geophys. Res. (Atmos.), (submitted, Nov. 2019), 2019a.</p><p>Lovejoy, S., Weather, Macroweather and Climate: our random yet predictable atmosphere, 334 pp., Oxford U. Press, 2019b.</p><p>Lovejoy, S., Fractional Relaxation noises, motions and the stochastic fractional relxation equation Nonlinear Proc. in Geophys. Disc., https://doi.org/10.5194/npg-2019-39, 2019c.</p>


2020 ◽  
Author(s):  
Shaun Lovejoy

Abstract. In part I, we considered the zero-dimensional heat equation showing quite generally that conductive – radiative surface boundary conditions lead to half-ordered derivative relationships between surface heat fluxes and temperatures: the Half-ordered Energy balance Equation (HEBE). The real Earth – even when averaged in time over the weather scales (up to ≈ 10 days) – is highly heterogeneous, in this part II, we thus extend our treatment to the horizontal direction. We first consider a homogeneous Earth but with spatially varying forcing. Using Laplace and Fourier techniques, we derive the Generalized HEBE (the GHEBE) based on half-ordered space-time operators. We analytically solve the homogeneous GHEBE, and show how these operators can be given precise interpretations. We then consider the full inhomogeneous problem with horizontally varying diffusivities, thermal capacities, climate sensitivities and forcings. For this we use Babenko's operator method which generalizes Laplace and Fourier methods. By expanding the inhomogeneous space-time operator at both high and low frequencies, we derive 2-D energy balance equations that can be used for macroweather forecasting, climate projections and for studying the approach to new (thermodynamic equilibrium) climate states when the forcings are all increased and held constant.


Author(s):  
Shaun Lovejoy ◽  
Roman Procyk ◽  
Raphael Hébert ◽  
Lenin Del Rio Amador

Sign in / Sign up

Export Citation Format

Share Document