Application of a dynamical two-box surface-atmosphere model to the Mount Pinatubo cooling event
|
|
- Sarah McKenzie
- 6 years ago
- Views:
Transcription
1 v. 2, 1/21/09 Application of a dynamical two-box surface-atmosphere model to the Mount Pinatubo cooling event Robert S. Knox and Kevin J. LaTourette* Department of Physics and Astronomy University of Rochester, Rochester, NY Abstract We analyze the global temperature change due to the Mt. Pinatubo eruption using a simple two-layer model of the atmosphere and surface to obtain results consistent with satellite data. Through analytic and numerical analysis we find a principal characteristic response time of 5 to 8 months and a climate sensitivity of 0.17 to 0.20 C/(W/m^2), corresponding to a negative instantaneous feedback. Our solutions were fit to the data, reproducing the results of a one-box model, and providing somewhat more detailed information about the feedbacks related to surface layer temperature. The formalism for coupling of the surface layer to the thermocline is set up but not applied. 1. Introduction A leading goal of climatology is to determine the effect of various forcings on Earth s climate system. A forcing is defined as the effective change of net radiative input at the top of the atmosphere, and has been adapted to all forms of disturbance [1]. We have a particular interest in the eruption of Mount Pinatubo because the event provides a *Present address: Department of Applied Mathematics, University of Arizona, Tucson, AZ rsk@pas.rochester.edu, klatourette@math.arizona.edu
2 substantially greater disruption to the climate when compared with other climatic disturbances. With its climax occurring June 15, 1991, the eruption created the greatest aerosol cloud of the 20 th century, releasing sufficient SO 2 into the atmosphere to reduce the average global temperature briefly by approximately 0.5 C. This powerful climatic event had the potential to [exceed] the accumulated forcing due to all anthropogenic greenhouse gasses added to the atmosphere since the industrial revolution began. [2] Occurring over a period relatively free of other climate forcings including the ENSO (El Niño-Southern Oscillation) effects and solar fluctuations, the Pinatubo event is an excellent candidate for analysis. Lindzen [3,4] studied volcano temperature anomalies by setting up a one-box model of the climate system with coupling to the thermocline. (The thermocline is an ocean layer in which temperature drops from its surface value. A layer of varying depth, of the order of dozens of meters, called the mixing layer, lies between the actual surface and the thermocline. The mixing layer temperature is uniform, as a result of turbulence caused by surface winds.) The forcing used was appropriate to that due to the 19thcentury Krakatoa eruption, and means for approximating the exact solution were taken, resulting in a prediction of the correct order of magnitude of global cooling, but with a fairly poor fit of the result to the data. It has been argued [5] that the reason for the poor fit was that the characteristic relaxation time τ in the box was assumed to be of the order of several years, rather than several months as the data imply. The resulting formalism implied a climate sensitivity of λ = 0.4 C/(W/m 2 ) (to be discussed below, sec. IV). Douglass and Knox [5] applied the Lindzen theory, the forcing calculations of Hansen [6,7], and the aerosol optical density estimates of Ammann et al. [8] to the Pinatubo case, producing a rather good fit to the data with a concomitant smaller τ of 6.8 months. They neglected the coupling to the thermocline. The climate sensitivity [1] predicted was λ = 0.18 C/(W/m 2 ), implying negative feedback, in apparent disagreement with the prevailing thinking about climate response. As a result, two comments were
3 published [9,10] whose salient claim was that the neglect of coupling to the thermocline had caused the work of [5] to be invalid, because its characteristic decay time was too short. The authors responded [11,12] with a revised calculation that included a crude estimated coupling, finding τ = 5.8 months, λ = 0.18 C/(W/m 2 ). Being concerned about the approximations involved in the thermocline coupling, Douglass et al. [13] solved the surface-thermocline coupling problem exactly. This introduces one new important physical parameter, namely, the eddy diffusion constant κ. They found results very similar to those of [5,11,12], but the value of κ required to fit the data was quite small, m 2 /s, as compared with the usually assumed value of m 2 /s. A one-box model for the surface and atmosphere is a serious oversimplification, as there are obviously many physically different layers and the interactions between these layers cannot be accounted for. Splitting the box into a surface layer and atmospheric layer [14,15] allows us to account for the different fluxes that occur between the two. The one-box model also lacks a temperature structure in the atmosphere, preventing it from handling the distinct radiative flows at the top and bottom of the atmosphere layer. Differences in the radiative flow at the top and bottom are related to the temperature differences, and our particular two-box model includes this climatic feature [14] without significantly increasing the complexity of the model. Adding a third box to the model will allow for the inclusion of the deep ocean (thermocline) coupled with the surface layer and provides an exact mathematical solution. In this preliminary report we develop the formalism for the thermocline coupling but do not proceed to the necessary numerical computations. Thus, our two-box numerical results are to be compared with the results of the first Pinatubo calculation [5]. In Section II we describe our dynamical model for a generalized forcing. In Section III we apply the Pinatubo forcing, in Section IV climate sensitivity and feedbacks are discussed, and finally Section V is a summary.
4 II. The dynamical model We begin with a static model [14] for globally averaged quantities. This will enable us to establish values of some important parameters. Our model has an incoming solar radiation, S 0 = 342 W/m 2, averaged over the surface of Earth, and proposes the following energy balance equations between the two layers in the steady state (see [14] and Figure 1): and ( 1+ b)εf A εf S = AS 0 + S NR (1) bεf A + F S = BS 0 S NR Q, (2) representing the atmosphere and surface layers, respectively. F A and F S are the ideal Stefan-Boltzmann fluxes given respectively by σt A 4 and σts 4, and σ = W/m 2 /K 4. The long-wave imbalance parameter b is the ratio of the net downward flux at the lowest portion of the atmosphere to the net upward flux at the top of the layer. This acts as a proxy for the lapse rate (temperature gradient), which is negative with increasing altitude, leading to an imbalance of long-wave radiation. The temperature at the bottom of the atmosphere is b 1/4 T A. The long-wave emissivity ε also acts as absorptivity in the term εf S through Kirchhoff s Law. Parameters A and B are the fractions of incoming solar radiation (S 0 ) absorbed by the atmosphere and surface, respectively. Imbedded within A and B are corrections for the multiple reflections between the two layers, but A and B may also be understood as independent adjustable parameters under the condition that the overall system reflectivity (planetary albedo) is α = 1 (A+B) [14]. The nonradiative flux, S NR, is a result of turbulent exchanges of heat and latent heat between the surface and atmosphere, and represents a substantial fraction of the up-flow from the surface (approximately 100 W/m 2 ). We also wish to incorporate feedbacks, and will
5 handle them in detail below. Q represents the net heat flux from the surface to the deep ocean at the top of the thermocline. It is assumed zero in the steady state. All parameter values are chosen so that T S0 = 288K and α 0.30, both of which are accepted values. This forces the parameters into a limited possible range and ensures that their values reflect their natural states. Without these constraints, a generalization of the energy balance model would be unrealistic as it would not necessarily be even a crude representation of Earth s climate system. We now make this model time dependent by writing (see Appendix A): c A du A dt + K AA u A + K AS u S = Δ[A 0 (t)s 0 (t)] (3) and c S du S dt + K SA u A + K SS u S = Δ[B 0 (t)s 0 (t)] Q(t), (4) where u A = ΔT A and u S = ΔT S are anomalies, i. e., temperature deviations from steadystate values. c A and c S are, respectively, the effective area specific heats of the atmosphere and surface boxes. The relaxation and coupling factors K IJ, defined in Eqs. (A14), result from expanding Eqs. (1) and (2) in Taylor series, dropping all but terms linear in the temperature anomalies and the forcing variations ΔA 0 S 0 and ΔB 0 S 0. In our application of the model, subscript A refers to that part of the atmosphere above 2 km, containing about 76% of its mass, and S refers to the surface plus the lower 2 km of atmosphere. Thus c A = J/m 2 /K and c S is an effective specific heat that includes the lower atmosphere contribution along with that of the actual surface, or c S = J/m 3 /K + h m C S. Here h m is an effective mixed-layer depth [13] and C S = 4.1 x 10 6 J/m 3 /K, the volume specific heat of seawater. We have defined the A layer in this manner because the temperature data available refers to satellite measurement at roughly 2 km, and it can therefore be represented by b 1/4 T A. The addition of a third box
6 allows for the inclusion of the thermocline into our model (see Figure 1), and is accounted for using the method of Douglass et al. [13]. The flux at the interface of the thermocline and surface layer is Q = C V κ ( u T / x) x=0, where the thermocline anomaly u T varies with depth x, measured downward. Here κ is the effective eddy diffusion coefficient and the specific heat C V is the volume specific heat of sea water (the coupling at the thermocline interface, which is below the mixed layer at a depth of dozens of meters, does not involve the complication of land and atmosphere). Under Laplace Transformation with Laplace variable p, Eqs. (3,4) become and (c A p + K AA ) u A ( p) + K AS u S ( p) = L{Δ[A 0 (t)s 0 (t)]} (5) K SA u A ( p) + ( c S p + K SS + C v κ 1 2 p 1 2 ) u ( p) = L{Δ[B (t)s (t)]}, (6) S 0 0 where L{...} indicates LaPlace transform. In the transforms we have assumed that solutions of interest are zero at t = 0. Here, we will treat the case of constant solar flux, so the driving terms become S 0 A (p) and S 0 B (p). Note that the thermocline effect, represented by the term containing κ, appears in the surface equation. This term includes a factor (p/κ) 1/2 that originates in the factor ( u T (x, p)/ x) x=0 of Q. The thermocline anomaly u T itself has disappeared from the equations because its relevant value, u T (x = 0), is equal to u S. III. Application to Pinatubo A. Case of volcano forcing As in the earlier work [5], we take the quantitative effect of the Pinatubo event as a globally averaged reduction of solar intensity caused by an increased aerosol reflectivity, resulting in a forcing that is proportional to the increased aerosol optical density (AOD). This is represented in our model by
7 ΔF(t) = 0.439A(t/t V )exp( t/t V ) W/m 2, (7) where A has been theoretically derived [6,7], and t V is the time of peak AOD (7.6 months after t = 0, the beginning of the event). For historical reasons and without danger of confusion we retain A as the constant on the right side of (7), although it is also our symbol for the atmospheric fraction of solar absorption. For its value, we use the most recent calculated by Hansen and coauthors [7], A = 21. With the notation k V = 1/t V, the Laplace transform of this forcing is Δ F( p) = 0.439Ak V ( p + k V ) 2. (8) The forcing under consideration originates in the excess reflection due to aerosols and is expressed in the standard way in terms of flux change at the top of the atmosphere. This does not necessarily mean that it is the atmosphere layer that receives the full effect of the forcing. We assume that the effect will be felt in proportion to the absorption of solar energy by each of the layers, i. e., a fraction φ A is assigned to the atmosphere layer and φ S to the surface layer, with φ A + φ S = 1. The fractions will be in the ratio φ A / φ S = A/B, where A and B are the static model parameters (see eqs. 1 and 2 and Fig. 1). Thus the forcings to appear in eqs. (5,6) replacing the entire right hand sides are Δ A 0 ( p) = φ A Δ F( p) and Δ B 0 ( p) = φ S Δ F( p). compact form and With these forcings and with the notation k IJ = K IJ /c I, Eqs. (5) and (6) take the ( p + k AA ) u A ( p) + k AS u S ( p) = (φ A /c A )Δ F( p) (9) k SA u A ( p) + ( p + k SS + κ 1 2 p 1 2 ) u ( p) = (φ /c )Δ F( p). (10) S S S The additional abbreviation introduced. κ 1 2 = κ 1 2 (C V /c S ) with dimension s 1/2 has been
8 In the case κ = 0 the functions u A (t) and u B (t) are elementary linear combinations of exponential functions. Otherwise a numerical attack on the inverse transform must be mounted. For the applications in this report, confined to the case κ = 0, the exponential solutions have been obtained and are set out in Appendix B. B. Parameter determination Let us review the several parameters involved in fitting the volcano-induced temperature data with our model. The eigenvalues λ 1 and λ 2 (Eq. B2) are obtained from the k s of Eq. (B1), which come from the area specific heats and the kinetic coefficients given in Eqs. (A7). Because of the fairly large number of parameters, we must adopt a systematic approach. One subset of parameters, called static, is set by reference to steady-state properties. The static set consists of {S 0, b, ε, A, B, S NR }, which may be seen in Eqs. (1,2), Fig.1 and Table 1. These parameters determine F A and F S or equivalently T A and T S. As explained in [14,18], a reasonable static set is chosen by imposing several conditions: T A = 288 K, albedo in the range , S NR approximately 100 W/m 2, ε in the range , and a value of b that produces a lower-atmosphere temperature consistent with observed lapse rates in conjunction with the predicted T A. The solar constant is always taken to be 342 W/m 2. Since there are many more static parameters than associated observables, we must be prepared to test the sensitivity (or robustness ) of any subsequent results to the particular set of static parameter values. Here we do this by running the entire calculation for two distinct static sets (see Table 1). From Eqs. (A14-16) we see that of the static quantities only parameters {b, ε } and quantities {T A, T S } will be needed in the dynamic problem. Here T A and T S 3 3 determine q A = 4σT A and qs = 4σT S. Since TS is fixed at 288 K, we have q S = 5.42 W/m 2 K. The other three quantities have limited ranges within the criterion of reasonableness. For static set 1, we have centered our dynamical calculations on the values b = 1.65 and ε = Whenever these were changed, the others were adjusted to
9 maintain T S = 288 K and the satisfaction of other static criteria. These adjustment covered ranges of 1.55 b 1.80 and ε and did show that the results of the dynamical parameter fit were robust. Static set 2 differs considerably and still yields similar dynamical results. Returning to the fits of time-dependent data, we define the dynamical set of parameters as {c A, c S, [ f ]}, where [ f ] refers to the set of four feedbacks in Eqs. (A7). Because of our definition of the model and the atmosphere and surface layers (see Sec. II), c A is fixed and c S depends only on h m, an effective ocean mixing-layer depth. The set of dynamical parameters is thus further reduced to {h m, [ f ]}. It is clearly futile to vary all four feedbacks with a single data set under analysis, so we arbitrarily focus on f SS and f AS, namely, those induced by surface temperature changes. C. Data Fitting Our primary data source is the global monthly satellite Microwave Sounding Unit lower troposphere temperature (TLT) anomaly data set [16], modified with El Niño and solar irradiance cycle perturbations effectively removed (see the discussion in [5], section 2.2). All reported analysis is based on this modified data set, denoted TLTm. Furthermore, because of the noise obviously existing in the set at times beyond 8.4 t V, only the 61 points between t = 0 and 8.4 t V were used in the data fits (see any of Figs. 2-4). Previously, TLTm has been regarded as the temperature change at the undifferentiated surface layer, although measurement takes place in the lower levels of the atmosphere. The present model allows us to be a bit more precise about this, and so, as discussed above, we fit the quantity b 1/4 u A to the data. In what follows, this quantity will be referred to as u B.
10 Static parameter set 1 The fitting began with all feedbacks assumed to be zero. In this case, it being the only other variable parameter in our protocol, h was varied, starting with the arbitrary value 15 m, until the sum of the quantity (u TLTm) 2 was minimized. This least-squares method produced the value h m = 15 m and an R 2 = 0.71 (when the fit was u = u S ) and 0.69 (when the fit was to u = b 1/4 u A ). Values of h m as high as 30 m were considered. We note that the eigenvalues of the kinetic matrix at this stage were s 1 and s 1, corresponding to lifetime parameters of 0.056t V (13 days) and 1.17t V (8.9 months). From the solutions we can see that the larger eigenvalue (shorter lifetime) tends to contribute a negligible amount in the earlier part of the TLTm time course (see any of Figs. 2-4). To next refine the fit, the feedback f SS was introduced and treated in the same way, minimizing the sum of squares and maximizing R 2, holding h constant at its previously determined value. The result was a visibly better fit. This process was repeated with f AS only, then with both feedbacks, switching back and forth until a maximum R 2 was obtained. The resulting best fits, both with R 2 = 0.74, are shown in Figs (h m was again varied between alternation of feedbacks, with no changes in its required value.) It was not possible to assign a better fit to either u B (Fig. 3) or u S (Fig. 4), since both produced the same value of R 2. As explained earlier, we consider the u B fit to be the more realistic because of the nature of the data. Clearly, Fig. 4 shows that either of the solution sets is adequate to explain the data for the parameters that deliberately make u B fit. Value of the mixed layer depth. As we adjust the values of the mixed layer h m, both temperature minima (u S and u B ) occur later as h m increases, but only slightly. The two most notable changes induced by modifying h m are in the amplitude of the minimum temperature and the characteristic response time. By doubling h m to a value of 30 m, we see the amplitude decrease by 0.1 ± 0.05 C, and the climate response time doubles as h m
11 doubles. When a value significantly greater than 20 m is used for the mixed layer depth, however, the response time does not follow the trend of the data, growing to the order of years not months. The best value for h m was 17.0 ± 4.0 m. Varying the static parameters to check the robustness of the steady-state model made only slight changes in the value of R 2 (of order ). Finally, we note that all fits were done with surface forcing (φ S = 0.97, φ A = 0.03). There was virtually no effect of small variations of these fractions. Results are shown in Table 2. Of particular interest are the feedbacks, 0.31 and 0.62, to be discussed in the section IV. 2. Static parameter set 2. Parameter set 2 differs from set 1 principally in that a larger fraction of the solar absorption is allotted to the atmosphere. As a result, somewhat different values of the other parameters are required to produce the usual surface temperature and albedo, and the temperature of the atmosphere layer is increased. Within the philosophy of using elementary climate models, these parameter variations are acceptable if further results based on them are essentially unchanged. We find this to be the case here, in that nearly the same dynamical parameters are required to fit the Pinatubo data as well as does parameter set 1. The best fit with parameter set 2 was found to be that for the bottom temperature b 1/2 u A. The solutions fitted were virtually identical to those of parameter set 1 and are not shown here. The 16% larger value of h m, therefore the effective heat capacity of the surface box, probably reflects the need for slowing down the kinetics during the rising phase of the signal, since the increased involvement of the atmosphere brings in the effect of its smaller heat capacity at earlier times.
12 IV. Climate sensitivity and feedback The present analysis, to the extent that it is successful, determines the feedback parameters required in computing the parameter called sensitivity, defined as the steadystate temperature shift per unit step-function forcing, λ = ΔT S /ΔF. For purely solar forcing, the no-feedback sensitivity for our case is given by [17] λ S0 = ΔT S ΔF = ΔT S = (1 α)δs 0 T S0 1 4(1 α)s 0 1 γ, where γ = S NR [(1+ b) + ba]s 0. (13) The sensitivity in the presence of feedback is λ S = λ S0 /(1 f S eff), where f S eff is a straightforward but rather complicated mixture of the feedbacks (Appendix A, Eq. A18). The value of λ S0 is commonly considered to be 0.30 C/(W/m 2 ) [1] but has been corrected to 0.36 within the context of the present model [14,17]. The dynamics of the one-layer Pinatubo treatment [5] contained a single feedback parameter f that can be identified with our f eff. A negative effective feedback was indicated in that work, and, as seen above, we also require negative feedbacks. Any attempt to connect the value found in [5] to our model introduces many new undetermined parameters, as seen in Appendix A. Therefore, as an exploratory example, let us assume that the temperature dependence of non-radiative transfer, represented by q NR = ds NR /dt S, is a dominant feedback cause, to the exclusion of all others. It is not hard to imagine that the non-radiative flux will increase with surface temperature, and therefore the fitting-predicted negative feedbacks are consistent with our model. According to Eqs. C15 and C17, f AS = q NR /ε 0 q S and f SS = q NR /q S. These equations predict 2.0 and 1.43 W/m 2 /K, respectively, for the value of q NR. Comparison may be made with two estimates in the literature. The two-layer model work of Barker and Ross [19] predicts q NR = 0.33 W/m 2 /K. The thermodynamic treatment of Hartmann [20] predicts that the sensible heat part of q NR will be 12 W/m 2 /K. This wide disparity calls for much more careful analysis of the model and the meaning of S NR. For the non-zero feedbacks found here, we have (Eq. C20)
13 f S eff ( = 1+ b ) f SS bε f AS 1+ b bε. (14) Using the static parameters employed in obtaining the feedbacks, one predicts f S eff = 0.93 (static set 1) and 1.02 (static set 2). These are in remarkably good agreement with the value 1.0 ± 0.4 found in the earlier work [5]. Looking further into the feedback signs, it can also be argued that the term B 1 =( B/ T S ) 0, which is a positive-signed component of f SS, also has a negative value, as follows: as the surface temperature increases, cloud cover will also increase, which will raise the absorption A in the atmosphere and increase reflectivity, both of which will tend to decrease the absorption B of the surface layer. V. Summary/Conclusion Our solutions of the dynamical two-box model applied to the Mt. Pinatubo effect are entirely consistent with earlier work that modeled the system with one box. The parameter fits, while necessarily very crude because of the lack of multiple data sets and the number of parameters required, pass tests of reasonableness and provide more detailed information about surface-temperature-induced feedbacks. These feedbacks are introduced here by reference to temperature dependence of parameters in a two-box surface-atmosphere system. We believe that this treatment is new and can increase the applicability and value of the simple model.
14 Appendix A. Basis of the time dependent model Equations (1) and (2) are generalized to the time-dependent case as follows: de A dt = c A dδt A dt = ( 1+ b)εf A + εf S + AS 0 + S NR (A1) and de S dt = c S dδt S dt = bεf A F S + BS 0 S NR Q, (A2) where E A and E S are the energy content per unit area of the A and S boxes, respectively, and c A and c S are the area specific heats of the atmosphere and surface layers, respectively. Eqs. (1) and (2) of the text refer to the case in which E A and E S are constants. We have assumed, as in other EBM treatments (e. g., [15]), that in the region of temperatures concerned, temperature-independent specific heats exist such that de A = c A dt A and de S = c S dt S. Assuming that the temperatures make small departures from their steady-state values T A0 and T S0, such that T A = T A0 + ΔT A and T S = T S0 + ΔT S, we have, by Taylor expansion, assuming that ΔT A and ΔT S are small quantities: F A = F A0 + q A ΔT A and F S = F S0 + q S ΔT S, (A3,4) where F A0 and F S0 are the Stefan-Boltzmann fluxes σt 4 A0 and σ T 4 S0, and q A = 4σT 3 A0, q S 3 = 4σT S0. In addition to these explicitly temperature-dependent quantities, nearly all the parameters in Eqs. (A1,2) may well depend on T A or T S or both. The equations are therefore generally nonlinear. However, for small anomalies the Taylor expansion method may be extended to the parameters as well as the F s. We therefore write A(t) = A 0 + A 1 ΔT A (t), B(t) = B 0 + B 1 ΔT S (t), S NR (t) = S NR0 + q NR ΔT S (t), (A5,6,7)
15 where the new parameters A 1, B 1 and q NR are adjustable unknowns or may be estimated from consideration of Taylor expansion coefficients. For example, q NR = S NR T S, with other variables constrained as required by the context. In practice, we are interested in the effect of imposed variations of the solar flux and of the composition of the atmosphere. In the former case we have S(t) = S 0 + ΔS(t). (A8) In this paper we attribute the effect of Pinatubo s eruption to extrinsic changes in A 0 and B 0, ignoring the aerosol s effect on ε. In other applications, where changes of greenhouse gas concentrations are of interest, there would be some small effects on A and B, we assume these to be negligible at this modeling stage. For ε, we consider the possibility of both an externally imposed change Δε 0 (t) and a system-induced temperature dependence, so that ε(t) = ε 0 + Δε(t), (A9) where Δε(t) = Δε 0 (t) + ε 1A ΔT A (t) + ε 1S ΔT S (t). (A10) Finally, the imbalance parameter b must be considered. Following a suggestion by R. Henry [21], we assume that b depends primarily on ε and write b = b 0 + b 1 Δε, (A11) which, through Eq. (A10), gives b an implicit temperature dependence as well as a dependence on Δε 0 (t). b. Construction of the working equations. Upon substituting Eqs. (A3-11) into Eqs. (A1, 2) one obtains terms containing Δ factors to the zeroth, first, and second power. Only those with first power Δ factors are retained. Those with no Δ factors comprise the expression for the steady state of the system (because the left hand side is zero). Those with two or more Δ factors are neglected in order to be consistent with
16 the neglect of the higher terms in F A and F S, in keeping with the general philosophy of linearization. These assumptions must of course be revisited if the solutions appear to be driven out of the linear regime. In what follows, for economy and without likelihood of confusion, we drop the subscripts 0 from b 0, ε 0, A 0, B 0, F A0, F S0, and S 0. The working equations become c A dδt A (t) dt c S dδt S (t) dt where the K coefficients are defined by = K AA ΔT A (t) K AS ΔT S (t) + Δn A (t) (A12) = K SA ΔT A (t) K SS ΔT S (t) + Δn S (t) Q (A13) K AA = (1+ b)εq A (1 f AA ), K AS = εq S (1 f AS ), K SA = bεq A (1 f SA ), K SS = q S (1 f SS ), (A14a-d) and the driving terms are Δn A (t) = [F S (1+ b + b 1 ε)f A ]Δε 0 (t) + AΔS(t) + SΔA(t), Δn S (t) = (b + b 1 ε)f A Δε 0 (t) + BΔS(t) + SΔB(t). (A15) (A16) As pointed out in the text, the coupling term Q is proportional to ΔT S and is easily incorporated into the linearized equations. In the K s, all factors of the form (1 f) become unity when there is no parameter temperature dependence, i. e., when there is no feedback. A factor f XY may be thought of as the feedback on the energy balance of box X due to a change in temperature T Y. One may have the impression that by choosing arbitrary values of the f s nearly anything could be predicted with these equations. This is true, except that the f s are not in fact arbitrary. They are rather complicated combinations of the physical parameters introduced above:
17 f AA = [F S (1+ b + b 1 ε)f A ]ε 1A + SA 1 (1+ b)εq A, (A17a) f AS = [(1+ b + b 1 ε)f A F S ]ε 1S q NR εq S, (A17b) f SA = (b + b 1 ε)f A ε 1A + SA 1 bεq A, (A17c) f SS = (b + b 1 ε)f A ε 1S q NR + SB 1 q S, (A17d) These feedbacks are expected to lie in the range { < f < 1}, which excludes oscillatory behavior. In some contexts the individual feedback parameters are usefully combined into the following overall effective feedback, particularly in the case of the climate sensitivity. An expression such as Eq. (10) of the text can be shown to acquire an additional factor 1/(1 f S eff ), where f S eff = 1 1 (1+ b)( f AA + f SS f AA f SS ) bε 0 ( f AS + f SA f AS f SA ) 1+ b bε 0 1 f AA (1+ b)b + f SA ba (1+ b)b + ba, (A18)
18 Appendix B. Solutions The Laplace-transformed solutions, Eqs. (11,12), may be converted readily into linear combinations of exponential decays in the case κ = 0. The time constants of interest are the inverse of the eigenvalues of the kinetic matrix K AA /c A K SA /c S K AS /c A K SS /c S = k k AA AS k SA k SS, (B1) having the values λ 1 λ 2 = (k + k )/2 (k k AA SS AA SS )2 /4 + k AS k SA (k AA + k SS )/2 + (k AA k SS ) 2 /4 + k AS k. SA (B2) The resulting solutions, after quite a bit of tedious linear algebra, may be expressed as u A (t) = 0.439A φ A ( λ 2 λ 1 )t V c A Ω (t) k 2 AS φ S c S Ω (t) 1 (B3) and u S (t) = 0.439A φ S ( λ 2 λ 1 )t V c S Ω (t) k 3 SA φ A c A Ω (t) 1, (B4) where Ω 1 (t) = exp( λ t) 1 (k V λ 1 ) exp( λ t) 2 2 (k V λ 2 ) + 1+ (k λ )t V 2 1+ (k λ V 1 )t 2 (k V λ 2 ) 2 (k V λ 1 ) 2 Ω 2 (t) = k λ SS 1 (k V λ 1 ) exp( λ t) k λ SS (k V λ 2 ) exp( λ t) exp( k t), (B5) V + (k λ )[1+ (k λ )t] SS 2 V 2 (k λ )[1+ (k λ )t] SS 1 V 1 (k V λ 2 ) 2 (k V λ 1 ) 2 exp( k t), V (B6) and Ω 3 (t) = k λ AA 1 (k V λ 1 ) exp( λ t) k λ AA (k V λ 2 ) exp( λ t)
19 (k λ )[1+ (k λ )t] AA 2 V 2 (k λ )[1+ (k λ )t] AA 1 V 1 (k V λ 2 ) 2 (k V λ 1 ) 2 exp( k t). V (B7) In the case φ S = 1, φ A = 0, and k AS = k SA = 0 (reducing the problem effectively to one layer), solution (B4) reduces to the one obtained by Douglass and Knox [5], noting the additional correspondences BS 0 = A, k SS = 1/τ = 1/(c S λ). Including coupling to the ocean. Eqs. (9,10) may be solved directly for the Laplace-transformed anomalies: and u A ( p) = 0.439A t V u S ( p) = 0.439A t V 1 ( ) 2 p + k V 1 p + k V ( p + k SS + κ p 1 2 )(φ A / c A ) k AS (φ S / c S ) (B8) ( p + k AA )( p + k SS + κ p 1 2 ) k AS k SA ( p + k AA )(φ S /c S ) k SA (φ A /c A ). (B9) ( κ p 1 2 ) k AS k SA ( ) 2 ( p + k AA ) p + k SS + Preliminary work has shown that these transformed expressions may be inverted successfully in the case κ > 0 by using the same contour as that use in ref. [13] and that this two-box version reproduces the one-box ocean-delayed solution found there.
20 Appendix C. The meaning of λ and τ in a two-box model In a one-box model, the dynamical equation for the temperature anomaly is usually written in the form hc dδt (t) dt = ΔT (t) λ + ΔF(t), (C1) where hc is the effective heat capacity of the box and where λ can be identified with the asymptotic ( equilibrium ) temperature shift associated with a step-function forcing ΔF 0 : ΔT ( ) = λδf 0. (C1) λ remains a parameter under different forcings and may be used to fit the data. However, in a model any having more than one box, λ is a derived parameter, i. e., it must be defined and extracted from the theory. Since the one-box ΔT is used to describe the surface temperature, in the two-box case it seems most reasonable to compute a quantity λ S as the response ΔT S ( ) to a step-function forcing. According to Eqs. (A12,13), K AA ΔT A ( ) + K AS ΔT S ( ) = φ A ΔF 0, K SA ΔT A ( ) + K SS ΔT S ( ) = φ S ΔF 0, (C3) (C4) from which we find λ S = b(1 f SA )φ A + (1+ b)(1 f AA )φ S (1+ b)(1 f AA )(1 f SS ) bε(1 f AS )(1 f SA ) 1 q S. (C5) The one-box task of fitting data by adjusting λ is seen to have ballooned into selecting all the feedbacks. This is what one is doing, effectively, when adjusting λ in the one-box model, where there is only one feedback parameter. When there is no feedback, Eq. (C5) yields λ S0 = bφ A + (1+ b)φ S 1+ b bε 1 q S = ba + (1+ b)b (1 α)(1+ b bε) 1 q S. (C6) Eq. (C6) is algebraically equivalent to text Eq. (13).
21 There is no easy connection between the two-box lifetime parameters and the onebox τ = hcλ, since the two-box ΔT S depends on two relaxation times. As a representative value of τ we may use the inverse of the smaller eigenvalue, λ 1 of Eq.(B2), which will dominate the behavior of both ΔT A and ΔT S except at short times.
22 References 1. K. P. Shine, Y. Fouquart, V. Ramaswamy, S. Solomon, and J. Srinivasan, Radiative Forcing. in Climate Change 1994, edited by J. T. Houghton et al., Cambridge Univ. Press, Cambridge, 1995, pp J. Hansen, A. Lacis, R. Ruedy, and M. Sato (1992), Potential climate impact of Mount Pinatubo eruption, Geophys. Res. Lett. 19, R. S. Lindzen (1994), Climate dynamics and global change, Ann. Rev. Fluid Mech. 26, R. S. Lindzen (1995), Constraining possibilities versus signal detection in Martinson-D-G, Natural Climate Variability on Decade-to-Century Time Scales (National Academy Press, Washington, DC), pp D. H. Douglass and R. S. Knox, Climate forcing by the volcanic eruption of Mount Pinatubo, Geophys. Res. Lett. 32 (5), L05710, doi: /2004gl (2005). Revised version archived as 6. J. Hansen, M. Sato, and R. Ruedy (1997), Radiative forcing and climate response, J. Geophys. Res. 102, Hansen, J., and 27 others (2002), Climate forcings in Goddard Institute for Space Studies SI2000 simulations, J. Geophys. Res., 107(D18), 4347, doi: /2001jd C. M. Ammann, G. A. Meehl, and W. W. Washington (2003), A monthly and latitudinally varying volcanic forcing data set in simulations of 20th century climate, Geophys. Res. Lett., 30(12), 1657, doi: /2003gl T. M. L. Wigley, C. M. Ammann, and B. D. Santer (2005), Using the Mount Pinatubo volcanic eruption to determine climate sensitivity: Comments on Climate forcing by the volcanic eruption of Mount Pinatubo by David H.
23 Douglass and Robert S. Knox, Geophys. Res. Lett. 32, L20711, doi: /2005gl A. Robock (2005), Using the Mount Pinatubo volcanic eruption to determine climate sensitivity: Comments on Climate forcing by the volcanic eruption of Mount Pinatubo by David H. Douglass and Robert S. Knox, Geophys. Res. Lett. L20709, doi: /2005gl D. H. Douglass and R. S. Knox (2005), Reply to comment by T. M. Wigley et al. on Climate forcing by the volcanic eruption of Mount Pinatubo, Geophys. Res. Lett. 32, L20710, doi: /2005gl D. H. Douglass and R. S. Knox (2005), Reply to comment by A. Robock on Climate forcing by the volcanic eruption of Mount Pinatubo, Geophys. Res. Lett. 32, L20712, doi: /2005gl D. H. Douglass, R. S. Knox, B. D. Pearson, and A. Clark, Jr. (2006), Thermocline flux exchange during the Pinatubo event, Geophys. Res. Lett. 33, L19711, doi: /2006gl R. S. Knox (2004), Non-radiative energy flow in elementary climate models, Phys. Lett. A 329, K. M. Shell and R. C. J. Somerville (2005), A generalized energy balance climate model with parameterized dynamics and diabatic heating, J. Climate 18, J. R. Christy, R. W. Spencer, and W. D. Braswell (2000), MSU tropospheric temperatures: dataset construction and radiosonde comparisons, J. Atmos. Oceanic Tech., 17, Updates available at This sensitivity is the response to a direct forcing at the top of the atmosphere due to a variation in the solar flux. It differs slightly from the standard IPCC sensitivity to radiative forcing in two ways. First, radiative forcing is defined by a two-stage process discussed in [1] generally and in [14] as it relates to the
24 present model. For our purposes the difference in λ is negligible. Second, the factor 1/(1 γ ) is not normally included in standard treatments, probably because it is considered a feedback factor. Without this factor, one obtains the standard Stefan-Boltzmann result ΔT S /ΔF = 0.30 C/(w/m 2 ). With the factor included, 0.30 changes to R. S. Knox (1999), Physical aspects of the greenhouse effect and global warming, Amer. J. Phys. 67, J. R. Barker and M. H. Ross (1999), An introduction to global warming, Amer. J. Phys. 67, D. L. Hartmann (1994), Global Physical Climatology (Academic Press, San Diego et al.), pages 114 and R. L. Henry (2008), private communication.
25 Table 1. Static fitting parameters defining the steady-state problem, as discussed in text. See also Figure 1. Units Static set 1* Static set 2* T S K T A K α = 1 (A + B) b ε S NR W/m A (atm. absorp. fraction) B (sfc. absorp. fraction) S 0 W/m q S W/m 2 /K q A W/m 2 /K φ A φ S * Quantities in bold fixed; others chosen to satisfy the criteria of reasonability [14,18]. Derived quantities needed for the dynamical equations are shown in italics Table 2. Dynamical fitting parameters, as discussed in text. The derived quantities λ S and τ are defined and discussed in Appendix C. Units Static set 1 Static set 2 ƒ AA ƒ AS ƒ SS ƒ SA ƒ eff h m m λ S C/(W/m 2 ) τ months A (forcing parameter) 21 21
26 Figure 1. The two-layer box model and its principal fluxes, with coupling to a box representing the deeper ocean ( thermocline ). For full description of the parameters in the atmosphere and surfaces boxes, see text following Eqs. (1) and (2). Edited version of Fig. 1 of Knox [14].
27 Figure 2. Time dependent atmosphere and surface temperature anomalies induced by the Mt. Pinatubo forcing using parameter set 1 (forcing mainly felt at the surface layer). Circles are the TLTm data set. Here a best fit of the atmosphere layer solution (u B ) to the data (blue curve) has been found, with feedbacks constrained to zero and by varying the effective surface layer parameter h m. Forcing held in Surface layer Change in temperature TLTm Atmosphere Surface t' (t in units of tv)
28 Figure 3. Plot of the time dependent atmosphere and surface temperature anomalies induced by the Mt. Pinatubo forcing using parameter set 1 (forcing mainly felt at the surface layer). Circles are the TLTm data set. Here a best fit has been found with the atmosphere layer solution (u B ) to the data (blue curve). Results with parameter set 2, not shown, with forcing felt in both layers, are identical Change in temperature TLTm Atmosphere Surface t' (t in units of tv)
29 Figure 4. Plot of the time dependent atmosphere and surface temperature anomalies induced by the Mt. Pinatubo forcing using parameter set 1 (forcing mainly felt at the surface layer). Circles are the TLTm data set. Here a best fit has been found with the surface-layer solution (u S ) (pink curve) to the data Change in temperature TLTm Atmosphere Surface t' (t in units of tv)
Climate sensitivity of Earth to solar irradiance: update
Paper presented at 2004 Solar Radiation and Climate (SORCE) meeting on Decade Variability in the Sun and the Climate, Meredith, New Hampshire, October 27-29, 2004 Climate sensitivity of Earth to solar
More informationEquation for Global Warming
Equation for Global Warming Derivation and Application Contents 1. Amazing carbon dioxide How can a small change in carbon dioxide (CO 2 ) content make a critical difference to the actual global surface
More informationClimate Feedbacks from ERBE Data
Climate Feedbacks from ERBE Data Why Is Lindzen and Choi (2009) Criticized? Zhiyu Wang Department of Atmospheric Sciences University of Utah March 9, 2010 / Earth Climate System Outline 1 Introduction
More informationStefan-Boltzmann law for the Earth as a black body (or perfect radiator) gives:
2. Derivation of IPCC expression ΔF = 5.35 ln (C/C 0 ) 2.1 Derivation One The assumptions we will make allow us to represent the real atmosphere. This remarkably reasonable representation of the real atmosphere
More informationAn Introduction to Coupled Models of the Atmosphere Ocean System
An Introduction to Coupled Models of the Atmosphere Ocean System Jonathon S. Wright jswright@tsinghua.edu.cn Atmosphere Ocean Coupling 1. Important to climate on a wide range of time scales Diurnal to
More informationFundamentals of Atmospheric Radiation and its Parameterization
Source Materials Fundamentals of Atmospheric Radiation and its Parameterization The following notes draw extensively from Fundamentals of Atmospheric Physics by Murry Salby and Chapter 8 of Parameterization
More informationConsequences for Climate Feedback Interpretations
CO 2 Forcing Induces Semi-direct Effects with Consequences for Climate Feedback Interpretations Timothy Andrews and Piers M. Forster School of Earth and Environment, University of Leeds, Leeds, LS2 9JT,
More informationRadiation in climate models.
Lecture. Radiation in climate models. Objectives:. A hierarchy of the climate models.. Radiative and radiative-convective equilibrium.. Examples of simple energy balance models.. Radiation in the atmospheric
More informationClimate Change: some basic physical concepts and simple models. David Andrews
Climate Change: some basic physical concepts and simple models David Andrews 1 Some of you have used my textbook An Introduction to Atmospheric Physics (IAP) I am now preparing a 2 nd edition. The main
More informationArctic Climate Change. Glen Lesins Department of Physics and Atmospheric Science Dalhousie University Create Summer School, Alliston, July 2013
Arctic Climate Change Glen Lesins Department of Physics and Atmospheric Science Dalhousie University Create Summer School, Alliston, July 2013 When was this published? Observational Evidence for Arctic
More informationInterannual variability of top-ofatmosphere. CERES instruments
Interannual variability of top-ofatmosphere albedo observed by CERES instruments Seiji Kato NASA Langley Research Center Hampton, VA SORCE Science team meeting, Sedona, Arizona, Sep. 13-16, 2011 TOA irradiance
More informationSensitivity of climate forcing and response to dust optical properties in an idealized model
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jd007198, 2007 Sensitivity of climate forcing and response to dust optical properties in an idealized model Karen
More informationGlobal temperature record reaches one-third century
Dec. 16, 2011 Vol. 21, No. 7 For Additional Information: Dr. John Christy, (256) 961-7763 john.christy@nsstc.uah.edu Dr. Roy Spencer, (256) 961-7960 roy.spencer@nsstc.uah.edu Global temperature record
More informationLecture 7: The Monash Simple Climate
Climate of the Ocean Lecture 7: The Monash Simple Climate Model Dr. Claudia Frauen Leibniz Institute for Baltic Sea Research Warnemünde (IOW) claudia.frauen@io-warnemuende.de Outline: Motivation The GREB
More informationClimate Dynamics Simple Climate Models
Climate Dynamics Simple Climate Models John Shepherd School of Ocean & Earth Science Southampton Oceanography Centre 1) Basic facts and findings Overview : 4 Lectures The global energy balance Zero-dimensional
More informationA Review of Soden et al: Global Cooling After the Eruption of Mount Pinatubo: A Test of Climate Feedback by Water Vapor.
Suvi Flagan ESE/Ge 148a A Review of Soden et al: Global Cooling After the Eruption of Mount Pinatubo: A Test of Climate Feedback by Water Vapor. By: BJ Soden, RT Wetherald, GL Stenchikov, and A Robock.
More informationRecent Climate History - The Instrumental Era.
2002 Recent Climate History - The Instrumental Era. Figure 1. Reconstructed surface temperature record. Strong warming in the first and late part of the century. El Ninos and major volcanic eruptions are
More informationTHE EXOSPHERIC HEAT BUDGET
E&ES 359, 2008, p.1 THE EXOSPHERIC HEAT BUDGET What determines the temperature on earth? In this course we are interested in quantitative aspects of the fundamental processes that drive the earth machine.
More informationLet s make a simple climate model for Earth.
Let s make a simple climate model for Earth. What is the energy balance of the Earth? How is it controlled? ó How is it affected by humans? Energy balance (radiant energy) Greenhouse Effect (absorption
More informationTemperature response of Earth to the annual solar irradiance cycle
Physics Letters A 323 (2004) 315 322 www.elsevier.com/locate/pla Temperature response of Earth to the annual solar irradiance cycle David H. Douglass a, Eric G. Blackman a,b, Robert S. Knox a,b, a Department
More informationThe Planck Blackbody Equation and Atmospheric Radiative Transfer
The Planck Blackbody Equation and Atmospheric Radiative Transfer Roy Clark Ventura Photonics There appears to be a lot of confusion over the use of the terms blackbody absorption and equilibrium in the
More information1. Weather and climate.
Lecture 31. Introduction to climate and climate change. Part 1. Objectives: 1. Weather and climate. 2. Earth s radiation budget. 3. Clouds and radiation field. Readings: Turco: p. 320-349; Brimblecombe:
More informationGlobal Energy Balance Climate Model. Dr. Robert M. MacKay Clark College Physics & Meteorology
Global Energy Balance Climate Model Dr. Robert M. MacKay Clark College Physics & Meteorology rmackay@clark.edu (note: the value of 342 W/m 2 given in this figure is the solar constant divided by 4.0 (1368/4.0).
More informationLecture 3a: Surface Energy Balance
Lecture 3a: Surface Energy Balance Instructor: Prof. Johnny Luo http://www.sci.ccny.cuny.edu/~luo Total: 50 pts Absorption of IR radiation O 3 band ~ 9.6 µm Vibration-rotation interaction of CO 2 ~
More information2. Meridional atmospheric structure; heat and water transport. Recall that the most primitive equilibrium climate model can be written
2. Meridional atmospheric structure; heat and water transport The equator-to-pole temperature difference DT was stronger during the last glacial maximum, with polar temperatures down by at least twice
More informationEarth s Radiation Budget & Climate
Earth s Radiation Budget & Climate Professor Richard Allan University of Reading NERC Advanced Training Course Earth Observations for Weather & Climate Studies 5 9 September 2016 Quantify the main terms
More informationEarth: the Goldilocks Planet
Earth: the Goldilocks Planet Not too hot (460 C) Fig. 3-1 Not too cold (-55 C) Wave properties: Wavelength, velocity, and? Fig. 3-2 Reviewing units: Wavelength = distance (meters or nanometers, etc.) Velocity
More informationCHAPTER 4. THE HADLEY CIRCULATION 59 smaller than that in midlatitudes. This is illustrated in Fig. 4.2 which shows the departures from zonal symmetry
Chapter 4 THE HADLEY CIRCULATION The early work on the mean meridional circulation of the tropics was motivated by observations of the trade winds. Halley (1686) and Hadley (1735) concluded that the trade
More informationRadiation in the atmosphere
Radiation in the atmosphere Flux and intensity Blackbody radiation in a nutshell Solar constant Interaction of radiation with matter Absorption of solar radiation Scattering Radiative transfer Irradiance
More information2018 Science Olympiad: Badger Invitational Meteorology Exam. Team Name: Team Motto:
2018 Science Olympiad: Badger Invitational Meteorology Exam Team Name: Team Motto: This exam has 50 questions of various formats, plus 3 tie-breakers. Good luck! 1. On a globally-averaged basis, which
More informationThe Structure and Motion of the Atmosphere OCEA 101
The Structure and Motion of the Atmosphere OCEA 101 Why should you care? - the atmosphere is the primary driving force for the ocean circulation. - the atmosphere controls geographical variations in ocean
More informationIntroduction to Climate ~ Part I ~
2015/11/16 TCC Seminar JMA Introduction to Climate ~ Part I ~ Shuhei MAEDA (MRI/JMA) Climate Research Department Meteorological Research Institute (MRI/JMA) 1 Outline of the lecture 1. Climate System (
More informationSensitivity of Tropical Tropospheric Temperature to Sea Surface Temperature Forcing
Sensitivity of Tropical Tropospheric Temperature to Sea Surface Temperature Forcing Hui Su, J. David Neelin and Joyce E. Meyerson Introduction During El Niño, there are substantial tropospheric temperature
More informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More informationAbstract: The question of whether clouds are the cause of surface temperature
Cloud variations and the Earth s energy budget A.E. Dessler Dept. of Atmospheric Sciences Texas A&M University College Station, TX Abstract: The question of whether clouds are the cause of surface temperature
More informationTwentieth-Century Sea Surface Temperature Trends M.A. Cane, et al., Science 275, pp (1997) Jason P. Criscio GEOS Apr 2006
Twentieth-Century Sea Surface Temperature Trends M.A. Cane, et al., Science 275, pp. 957-960 (1997) Jason P. Criscio GEOS 513 12 Apr 2006 Questions 1. What is the proposed mechanism by which a uniform
More informationAPPLICATIONS WITH METEOROLOGICAL SATELLITES. W. Paul Menzel. Office of Research and Applications NOAA/NESDIS University of Wisconsin Madison, WI
APPLICATIONS WITH METEOROLOGICAL SATELLITES by W. Paul Menzel Office of Research and Applications NOAA/NESDIS University of Wisconsin Madison, WI July 2004 Unpublished Work Copyright Pending TABLE OF CONTENTS
More informationThe Energy Balance Model
1 The Energy Balance Model 2 D.S. Battisti 3 Dept. of Atmospheric Sciences, University of Washington, Seattle Generated using v.3.2 of the AMS LATEX template 1 ABSTRACT 5 ad 2 6 1. Zero-order climatological
More informationClimate 1: The Climate System
Climate 1: The Climate System Prof. Franco Prodi Institute of Atmospheric Sciences and Climate National Research Council Via P. Gobetti, 101 40129 BOLOGNA SIF, School of Energy, Varenna, July 2014 CLIMATE
More informationData and formulas at the end. Exam would be Weds. May 8, 2008
ATMS 321: Science of Climate Practice Mid Term Exam - Spring 2008 page 1 Atmospheric Sciences 321 Science of Climate Practice Mid-Term Examination: Would be Closed Book Data and formulas at the end. Exam
More informationExplaining Changes in Extremes and Decadal Climate Fluctuations
Explaining Changes in Extremes and Decadal Climate Fluctuations Gerald A. Meehl Julie Arblaster, Claudia Tebaldi, Aixue Hu, Ben Santer Explaining changes implies attributing those changes to some cause
More informationGlobal warming trend and multi-decadal climate
Global warming trend and multi-decadal climate variability Sergey Kravtsov * and Anastasios A. Tsonis * Correspondence and requests for materials should be addressed to SK (kravtsov@uwm.edu) 0 Climate
More informationATMO 551a Intro to Optical Depth Fall τ υ,z. dz = di υ. B[ v,t(z) ]e
Atmospheric Radiative Transfer We need to understand how energy is transferred via radiation within the atmosphere. We introduce the concept of optical depth. We will further show that the light moves
More informationData and formulas at the end. Real exam is Wednesday May 8, 2002
ATMS 31: Physical Climatology Practice Mid Term Exam - Spring 001 page 1 Atmospheric Sciences 31 Physical Climatology Practice Mid-Term Examination: Would be Closed Book Data and formulas at the end. Real
More informationSimple energy balance climate models
Chapter 2 Simple energy balance climate models Supplemental reading: 1 Budyko (1969) Held and Suarez (1974) Lindzen and Farrell (1977) North (1975) Sellers (1969) My initial purpose in beginning with an
More information1) The energy balance at the TOA is: 4 (1 α) = σt (1 0.3) = ( ) 4. (1 α) 4σ = ( S 0 = 255 T 1
EAS488/B8800 Climate & Climate Change Homework 2: Atmospheric Radiation and Climate, surface energy balance, and atmospheric general circulation Posted: 3/12/18; due: 3/26/18 Answer keys 1. (10 points)
More information3. Carbon Dioxide (CO 2 )
3. Carbon Dioxide (CO 2 ) Basic information on CO 2 with regard to environmental issues Carbon dioxide (CO 2 ) is a significant greenhouse gas that has strong absorption bands in the infrared region and
More informationLecture 3a: Surface Energy Balance
Lecture 3a: Surface Energy Balance Instructor: Prof. Johnny Luo http://www.sci.ccny.cuny.edu/~luo Surface Energy Balance 1. Factors affecting surface energy balance 2. Surface heat storage 3. Surface
More informationLecture 2 ENSO toy models
Lecture 2 ENSO toy models Eli Tziperman 2.3 A heuristic derivation of a delayed oscillator equation Let us consider first a heuristic derivation of an equation for the sea surface temperature in the East
More informationNOTES AND CORRESPONDENCE. On the Radiative and Dynamical Feedbacks over the Equatorial Pacific Cold Tongue
15 JULY 2003 NOTES AND CORRESPONDENCE 2425 NOTES AND CORRESPONDENCE On the Radiative and Dynamical Feedbacks over the Equatorial Pacific Cold Tongue DE-ZHENG SUN NOAA CIRES Climate Diagnostics Center,
More informationUnderstanding the Greenhouse Effect
EESC V2100 The Climate System spring 200 Understanding the Greenhouse Effect Yochanan Kushnir Lamont Doherty Earth Observatory of Columbia University Palisades, NY 1096, USA kushnir@ldeo.columbia.edu Equilibrium
More informationFOLLOW THE ENERGY! EARTH S DYNAMIC CLIMATE SYSTEM
Investigation 1B FOLLOW THE ENERGY! EARTH S DYNAMIC CLIMATE SYSTEM Driving Question How does energy enter, flow through, and exit Earth s climate system? Educational Outcomes To consider Earth s climate
More information( ) = 1005 J kg 1 K 1 ;
Problem Set 3 1. A parcel of water is added to the ocean surface that is denser (heavier) than any of the waters in the ocean. Suppose the parcel sinks to the ocean bottom; estimate the change in temperature
More informationSeeking a consistent view of energy and water flows through the climate system
Seeking a consistent view of energy and water flows through the climate system Robert Pincus University of Colorado and NOAA/Earth System Research Lab Atmospheric Energy Balance [Wm -2 ] 340.1±0.1 97-101
More informationEffect of Exclusion of Anomalous Tropical Stations on Temperature Trends from a 63-Station Radiosonde Network, and Comparison with Other Analyses
2288 JOURNAL OF CLIMATE VOLUME 16 Effect of Exclusion of Anomalous Tropical Stations on Temperature Trends from a 63-Station Radiosonde Network, and Comparison with Other Analyses JAMES K. ANGELL NOAA
More informationChapter 6: Modeling the Atmosphere-Ocean System
Chapter 6: Modeling the Atmosphere-Ocean System -So far in this class, we ve mostly discussed conceptual models models that qualitatively describe the system example: Daisyworld examined stable and unstable
More informationSome remarks on climate modeling
Some remarks on climate modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Selected overheads by Doug Nychka Outline Hierarchy of atmospheric modeling strategies
More informationWater Vapor Multiplier of Carbon Dioxide
Water Vapor Multiplier of Carbon Dioxide Harvey S. H. Lam October 10, 007 Abstract The temperature rise of the earth due to the direct greenhouse effects of atmospheric carbon dioxide can be amplified
More informationTeaching Energy Balance using Round Numbers: A Quantitative Approach to the Greenhouse Effect and Global Warming
Teaching Energy Balance using Round Numbers: A Quantitative Approach to the Greenhouse Effect and Global Warming Brian Blais Science and Technology Department Bryant College bblais@bryant.edu August 29,
More informationRadiative equilibrium Some thermodynamics review Radiative-convective equilibrium. Goal: Develop a 1D description of the [tropical] atmosphere
Radiative equilibrium Some thermodynamics review Radiative-convective equilibrium Goal: Develop a 1D description of the [tropical] atmosphere Vertical temperature profile Total atmospheric mass: ~5.15x10
More informationXV. Understanding recent climate variability
XV. Understanding recent climate variability review temperature from thermometers, satellites, glacier lengths and boreholes all show significant warming in the 2th C+ reconstruction of past temperatures
More informationExternally forced and internal variability in multi-decadal climate evolution
Externally forced and internal variability in multi-decadal climate evolution During the last 150 years, the increasing atmospheric concentration of anthropogenic greenhouse gases has been the main driver
More informationA New Basic 1-Dimension 1-Layer Model Obtains Excellent Agreement With the Observed Earth Temperature
A New Basic 1-Dimension 1-Layer Model Obtains Excellent Agreement With the Observed Earth Temperature Rainer Link and Horst-Joachim Lüdecke EIKE, European Institute for Climate and Energy, PO.Box 11011,
More informationSIMPLE CLIMATE MODELING. Ka Kit Tung
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES Volume 7, Number 3, May 2007 pp. 651 660 SIMPLE CLIMATE MODELING Ka Kit Tung Department of Applied Mathematics University
More informationClimate Variability Natural and Anthropogenic
Climate Variability Natural and Anthropogenic Jim Renwick NIWA Climate Research j.renwick@niwa.co.nz Climate equilibrium and climate forcings Natural forcings Anthropogenic forcings Feedbacks Natural variability
More informationGE510 Physical Principles of the Envt
GE510 Physical Principles of the Envt Radiative Forcing of Climate Radiative Forcing of Climate 1. What is radiative forcing of climate? Conceptually - A simple definition - Physical properties of molecules
More informationClouds in the Climate System: Why is this such a difficult problem, and where do we go from here?
Clouds in the Climate System: Why is this such a difficult problem, and where do we go from here? Joel Norris Scripps Institution of Oceanography CERES Science Team Meeting April 29, 2009 Collaborators
More informationGuest Blog for Climate Science: Roger Pielke Sr. Research Group News SPPI Reprint Series
Has the IPCC inflated the feedback factor? By Christopher Monckton Guest Blog for Climate Science: Roger Pielke Sr. Research Group News SPPI Reprint Series Has the IPCC inflated the feedback factor? In
More informationKey Feedbacks in the Climate System
Key Feedbacks in the Climate System With a Focus on Climate Sensitivity SOLAS Summer School 12 th of August 2009 Thomas Schneider von Deimling, Potsdam Institute for Climate Impact Research Why do Climate
More informationEarth s Energy Balance and the Atmosphere
Earth s Energy Balance and the Atmosphere Topics we ll cover: Atmospheric composition greenhouse gases Vertical structure and radiative balance pressure, temperature Global circulation and horizontal energy
More information2/18/2013 Estimating Climate Sensitivity From Past Climates Outline
Estimating Climate Sensitivity From Past Climates Outline Zero-dimensional model of climate system Climate sensitivity Climate feedbacks Forcings vs. feedbacks Paleocalibration vs. paleoclimate modeling
More informationCHAPTER 8. AEROSOLS 8.1 SOURCES AND SINKS OF AEROSOLS
1 CHAPTER 8 AEROSOLS Aerosols in the atmosphere have several important environmental effects They are a respiratory health hazard at the high concentrations found in urban environments They scatter and
More informationCLIMATE CHANGE: THE SUN S ROLE HUGH S 80 TH!
CLIMATE CHANGE: THE SUN S ROLE Gerald E. Marsh FOR HUGH S 80 TH! 1 BACKGROUND MATERIALS IPCC: Climate Change 2001: Working Group I: The Scientific Basis: http://www.grida.no/climate/ipcc_tar/wg1/index.htm
More informationLecture 10. Surface Energy Balance (Garratt )
Lecture 10. Surface Energy Balance (Garratt 5.1-5.2) The balance of energy at the earth s surface is inextricably linked to the overlying atmospheric boundary layer. In this lecture, we consider the energy
More informationEffects of Black Carbon on Temperature Lapse Rates
Effects of Black Carbon on Temperature Lapse Rates Joyce E. Penner 1 Minghuai Wang 1, Akshay Kumar 1, Leon Rotstayn 2, Ben Santer 1 University of Michigan, 2 CSIRO, 3 LLNL Thanks to Warren Washington and
More informationTorben Königk Rossby Centre/ SMHI
Fundamentals of Climate Modelling Torben Königk Rossby Centre/ SMHI Outline Introduction Why do we need models? Basic processes Radiation Atmospheric/Oceanic circulation Model basics Resolution Parameterizations
More informationReply to Lockwood and Fröhlich The persistent role of the Sun in climate forcing. Svensmark, H. and Friis-Christensen, E.
Reply to Lockwood and Fröhlich The persistent role of the Sun in climate forcing Svensmark, H. and Friis-Christensen, E. Danish National Space Center Scientific Report 3/2007 Reply to Lockwood and Fröhlich
More informationSolar radiation - the major source of energy for almost all environmental flows
Solar radiation - the major source of energy for almost all environmental flows Radiation = electromagnetic waves Different types of heat transfer: Heat conduction by molecular diffusion (no large-scale
More information2. Energy Balance. 1. All substances radiate unless their temperature is at absolute zero (0 K). Gases radiate at specific frequencies, while solids
I. Radiation 2. Energy Balance 1. All substances radiate unless their temperature is at absolute zero (0 K). Gases radiate at specific frequencies, while solids radiate at many Click frequencies, to edit
More informationSeparation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I. El Niño/La Niña Phenomenon
International Journal of Geosciences, 2011, 2, **-** Published Online November 2011 (http://www.scirp.org/journal/ijg) Separation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I.
More informationThe frequency response of temperature and precipitation in a climate model
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:0.029/20gl048623, 20 The frequency response of temperature and precipitation in a climate model Douglas G. MacMynowski, Ho Jeong Shin, 2 and Ken Caldeira 2 Received
More informationPoint-to-point response to reviewers comments
Point-to-point response to reviewers comments Reviewer #1 1) The authors analyze only one millennial reconstruction (Jones, 1998) with the argument that it is the only one available. This is incorrect.
More informationClimate Models & Climate Sensitivity: A Review
Climate Models & Climate Sensitivity: A Review Stroeve et al. 2007, BBC Paul Kushner Department of Physics, University of Toronto Recent Carbon Dioxide Emissions 2 2 0 0 0 0 7 6 x x Raupach et al. 2007
More informationLecture # 04 January 27, 2010, Wednesday Energy & Radiation
Lecture # 04 January 27, 2010, Wednesday Energy & Radiation Kinds of energy Energy transfer mechanisms Radiation: electromagnetic spectrum, properties & principles Solar constant Atmospheric influence
More informationThe Canadian Climate Model 's Epic Failure November 2016
The Canadian Climate Model 's Epic Failure November 2016 By: Ken Gregory The Canadian Centre for Climate Modeling and Analysis located at the University of Victoria in British Columbia submitted five runs
More informationEarth s Energy Budget: How Is the Temperature of Earth Controlled?
1 NAME Investigation 2 Earth s Energy Budget: How Is the Temperature of Earth Controlled? Introduction As you learned from the reading, the balance between incoming energy from the sun and outgoing energy
More informationName(s) Period Date. Earth s Energy Budget: How Is the Temperature of Earth Controlled?
Name(s) Period Date 1 Introduction Earth s Energy Budget: How Is the Temperature of Earth Controlled? As you learned from the reading, the balance between incoming energy from the sun and outgoing energy
More informationGlaciology HEAT BUDGET AND RADIATION
HEAT BUDGET AND RADIATION A Heat Budget 1 Black body radiation Definition. A perfect black body is defined as a body that absorbs all radiation that falls on it. The intensity of radiation emitted by a
More informationBoundary layer equilibrium [2005] over tropical oceans
Boundary layer equilibrium [2005] over tropical oceans Alan K. Betts [akbetts@aol.com] Based on: Betts, A.K., 1997: Trade Cumulus: Observations and Modeling. Chapter 4 (pp 99-126) in The Physics and Parameterization
More informationGEOL 437 Global Climate Change 2/1/18: Solar radiation and the annual cycle
GEOL 437 Global Climate Change 2/1/18: Solar radiation and the annual cycle Why are there seasons? How does the climate respond to the radiative annual cycle? How does the climate respond to changes in
More informationLecture 3: Global Energy Cycle
Lecture 3: Global Energy Cycle Planetary energy balance Greenhouse Effect Vertical energy balance Latitudinal energy balance Seasonal and diurnal cycles Solar Flux and Flux Density Solar Luminosity (L)
More informationShort-Term Climate Variability (Ch.15) Volcanos and Climate Other Causes of Holocene Climate Change
Short-Term Climate Variability (Ch.15) Volcanos and Climate Other Causes of Holocene Climate Change Volcanos and Climate We learned in Chapter 12 that the volanos play an important role in Earth s climate
More informationATMO 551a Fall The Carnot Cycle
What is a arnot ycle and Why do we care The arnot ycle arnot was a French engineer who was trying to understand how to extract usable mechanical work from a heat engine, that is an engine where a gas or
More informationLecture 3: Atmospheric Radiative Transfer and Climate
Lecture 3: Atmospheric Radiative Transfer and Climate Solar and infrared radiation selective absorption and emission Selective absorption and emission Cloud and radiation Radiative-convective equilibrium
More informationThe Sun Approaches Its 11 Year Minimum and Activity Cycle 24
The Sun Approaches Its 11 Year Minimum and Activity Cycle 24 Tom Woods, Laboratory for Atmospheric and Space Physics, University of Colorado, woods@lasp.colorado.edu Judith Lean, Naval Research Laboratory,
More informationClimate Sensitivity to Increasing Greenhouse Gases
Chapter 2 Climate Sensitivity to Increasing Greenhouse Gases James E. Hansen, Andrew A. Lacis, David H. Rind, and Gary L. Russell INTRODUCTION Climate changes occur on all time scales, as illustrated in
More informationLecture 2 Global and Zonal-mean Energy Balance
Lecture 2 Global and Zonal-mean Energy Balance A zero-dimensional view of the planet s energy balance RADIATIVE BALANCE Roughly 70% of the radiation received from the Sun at the top of Earth s atmosphere
More informationATMS 321: Sci. of Climate Final Examination Study Guide Page 1 of 4
ATMS 321: Sci. of Climate Final Examination Study Guide Page 1 of 4 Atmospheric Sciences 321: Final Examination Study Guide The final examination will consist of similar questions Science of Climate Multiple
More informationAn Introduction to Climate Modeling
An Introduction to Climate Modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Outline What is Climate & why do we care Hierarchy of atmospheric modeling strategies
More informationAtmospheric Thermodynamics
Atmospheric Thermodynamics R. Wordsworth February 12, 2015 1 Objectives Derive hydrostatic equation Derive dry and moist adiabats Understand how the theory relates to observed properties of real atmospheres
More information