GISS Modele Progress and Plans
Total Page:16
File Type:pdf, Size:1020Kb
GISS ModelE Progress and Plans WGCM, Princeton, Nov 2016 Gavin Schmidt and team GISS Post-CMIP5 Progress Goddard Institute for Space Studies MJO variability and prediction skill Self-generated QBO Enhancements to forcings (irrigation, volcanic, solar) Better use of single forcing runs Greatly improved ocean/sea ice simulations Planned GISS CMIP6 Configurations Goddard Institute for Space Studies Multiple configs w/variations for DECK runs: 1. GISS-E2.1 (ready) Variations: OMA vs MATRIX; R vs H ocean; L40 vs L96/102 2. GISS-E3 (mid-2017) C90+L96/102, same oceans; self-generated QBO, MATRIX aerosols, M&G cloud microphysics, cold pool convection 3. GISS-E4 (2018?) C180+L96/102, GO2 (GISS Ocean 2) (cubed sphere/ALE vertical) 1JULY 2012 K I M E T A L . 4645 Newly resolved modes I: MJO 1JULY 2012 K I M E T A L . 4645 1JULY 2012GISS$E K I M E T A L .GISS$E2 GISS$E2.1))))4651 GISS$E2.1)))) FIG. 3. Space–time spectrum of the 158N–158S symmetric component of precipitation divided by the background spectrum for (a) GPCP, (b) AR4a, (c) AR5a, (d) C_AR5a, (e) AR5c, (f) AR5a_Ent1, (g) AR5a_Ent1_Re, and (h) AR5a_Ent2_Re. Superimposed are the dispersion curves of the odd-numbered meridional mode equatorial waves for the equivalent depths of 12, 25, and 50 m. data are from the International Best Track Archive the Kelvin, equatorial Rossby wave, and MJO bands for Climate Stewardship (IBTrACS) dataset (Knapp (the last being defined as wavenumbers 1–3, periods 30– et al. 2010). 60 days) that can be found in the observations (Fig. 3a). The most significant improvement that AR5a has com- 3. Simulations of the MJO using Model E2 pared to AR4a is its simulation of the Kelvin mode. The Kelvin mode in AR5a is similar to that in observations in a. Simulations of the MJO in AR4 and AR5 versions (Kim et al, 2012) both its amplitude and phase speed; the implied equiv- of Model E2 22 FIG. 9. Phase–longitude diagram of OLR [contour plotted every 3 W m , positive (green) and negativealent (purple)] depth and is surface about evap- 25 m. Compared to AR5a, AR4a oration (W m22)/ 925-hPa moisture convergenceFollowing (kg Wheeler kg21 s2 and1) for Kiladis (a), (c) observations,(1999), wavenumber– and (b), (d) AR5a_Ent1.has a much Phases weaker are from and the faster Kelvin mode, which is also MJO life cycle composite; values arefrequency averaged between diagrams 108S are and constructed108N. to determine the mostly confined to high frequencies (i.e., periods less FIG. 3. Space–time spectrum ofcapability the 158N–15 of8S the symmetric models component to simulate of convectivelyprecipitation divided cou- bythan the background 7 days). Despite spectrum these for improvements, AR5a still (a) GPCP, (b) AR4a, (c) AR5a, (d)pled C_AR5a, equatorial (e) AR5c, waves (f) AR5a_Ent1, and the MJO. (g) AR5a_Ent1_Re, Figure 3 shows and the (h) AR5a_Ent2_Re.lacks the MJO Superimposed mode. are the dispersion curves of the odd-numbered meridional mode equatorial waves for the equivalent depths of 12, 25, and 50 m. a2DphasespaceofthetwoleadingPCsfromtheCEOFsymmetric wavenumber–frequencytropospheric power temperature spectra [nor- increasesFigure (Fig. 3 also 12b). contains The rel- the symmetric components of analysis. In this 2D phase space,malized distance by estimated from the background or- ative humidity power, Wheeler (Fig. 12c) and andthe precipitable wavenumber–frequency water (Fig. spectra of equatorial pre- igin represents amplitude of theKiladis MJO. (1999)] The strong-MJO of equatorial12d) precipitation also increase. from These obser- systematiccipitation changes from thecaused different by versions of Model E2. The eventdata occurs are from during the March–April International 2000 Best in the Track simulation. Archive thethe decreasing Kelvin, equatorial entrainment Rossby rate wave, (from and AR5a_Ent1 MJO bands to for Climate Stewardship (IBTrACS)vations and dataset several (Knapp versions of(the Model last being E2. Our defined focus as is wavenumbersC_AR5a (Fig. 1–3, 3d) periods represents 30– a version of Model E2 that Hovmo¨ ller diagrams of total,on anomalous the signals (deviationsdistinct from theAR5a) background can be characterized spectrum in asuses enhanced higher stability horizontal in the resolution than that in AR5a by fromet al. the 2010). seasonal cycle) and 20–100-day filtered equa- tropics.60 days) This that result can be is consistent found in the with observations those of Kim (Fig. et 3a).al. torial (158S–158N) precipitation show the eastward pro- (2011b),The most who significant showed improvementthat models with that stronger AR5a has MJOs com- pagation3. Simulations of organized of the precipitation MJO using anomalies Model E2 with phase alsopared had to a AR4a cold bias is its in simulation the upper troposphereof the Kelvin relative mode. The to 21 Kelvin mode in AR5a is similar to that in observations in speeda. Simulations;5ms during of the this MJO period in AR4 (Fig. and 11). AR5 Daily versions restart both its amplitude and phase speed; the implied equiv- filesof are Model saved during E2 the period of this event and used to initialize the 30-day integrations of AR5a. Note that we alent depth is about 25 m. Compared to AR5a, AR4a useFollowing restart files Wheeler during February–May and Kiladis (1999), 2000 to wavenumber– encompass has a much weaker and faster Kelvin mode, which is also thefrequency whole strong-MJO diagrams are period. constructed to determine the mostly confined to high frequencies (i.e., periods less capabilityDuring the of the course models of the to simulate 30-day integration,convectively the cou- than 7 days). Despite these improvements, AR5a still AR5apled equatorialversion systematically waves and deviates the MJO. from Figure the AR5a_Ent1 3 shows the lacks the MJO mode. version.symmetric Figure wavenumber–frequency 12a shows the composite power deviations spectra of [nor- the Figure 3 also contains the symmetric components of troposphericmalized by temperature estimated background from the first power, day of simulation. Wheeler and It the wavenumber–frequency spectra of equatorial pre- indicatesKiladis that (1999)] the tropical of equatorial atmosphere precipitation becomes from stabilized obser- cipitation from the different versions of Model E2. The (warmervations and upper/colder several versions lower of troposphere) Model E2. Our gradually focus is C_AR5a (Fig. 3d) represents a version of Model E2 that untilon the day signals 30. The distinct warming from the aloft background is greater spectrum than the in uses higher horizontal resolution than that in AR5a by cooling below so that the mass-weighted average of FIG. 10. Schematic diagram of the reinitialization experiment. Cold pool parameterization: Formed'from'downdrafts,'used'to'restrict'occurrence Goddard Institute for Space Studies of'weakly'entraining'plumes GISS$E2 YOTC'20=day'MJO'rain hindcast'Hovmöller diagrams (0.63'20=day correlation'with'TRMM TMI'with'cold'pool'vs. 0.70'TMI=Radar correlation) GISS$E3 Self-generated stratospheric QBO Tropical zonal mean winds ERA40 GISS$E2.1(L102) Rind et al (2014) Ocean model improvements Reformulation of GM eddy parameterisation GISS-Vertical Mixing Scheme Inclusion of GM vertical dependence G. Danabasoglu et al. / Ocean Modelling 73 (2014) 76–107 83 Evaluation with stand-alone ocean CORE-I/II protocol ClimatologyGISS V2 GISS V2 Model Obs Fig. 6. Time-mean meridional heat transports for the Atlantic Ocean. The black lines denoted by L&Y09 represent implied time-mean transport calculated by Large and Yeager, 2009 with shading showing the implied transport range in individual years for the 1984–2006 period. Direct estimates with their uncertainty ranges from the RAPID Zonaldata (square; Johns mean et al., 2011) and SST from Bryden andgradients Imawaki, 2001 (triangle; B&I01) are also shown. AMOC Heat flux 26ºN in stark contrast with the other models and observationally-based ticipating models, Fig. 6 implies many model differences in details data. The latitudinalSouthern variations in MHT Ocean for MRI-A reflect its AMOC of surfaceCORE heat fluxes, resulting II:" primarilyDanabasoglu from differences in sim- et al (2014) structure. Such variations seem to be common in the MHT distribu- ulated SSTs. One example is the much larger heat gain in BERGEN tions obtained with some other data assimilation products as well between 10°N and 30°N in contrast with most of the other models (see Munoz et al., 2011). We believe that, as discussed in Msadek where much smaller heat gains or even losses are suggested. The et al., 2013, errors in representations of the NADW cell and, partic- oceanic heat gain evident in most models between 45°N and ularly, in the vertical structure of h (see Fig. 11), are largely respon- 55°N – as indicated by the positive MHT slopes – is associated with sible for the substantially lower MHTs in all model simulations the surface heat fluxes acting to damp the cold SST biases present compared to observational estimates even in simulations with real- in these models (see Fig. 8) due to the incorrect path of the North istic overturning strengths. Although much smaller in its contribu- Atlantic Current (NAC) (e.g., Danabasoglu et al., 2012). tion to MHT, errors in the gyre components can also explain some of As hinted at above, AMOC is the dominant contributor to the the differences (Msadek et al., 2013). We note that non-eddy- Atlantic Ocean MHT (Böning et al., 2001; Msadek et al., 2013). resolving horizontal resolutions of the present models can contrib- The relationship between AMOC and MHT is presented in Fig. 7, ute to low MHTs due to changes in the mean rather than the eddy considering the scatter plot of the maximum AMOC transport heat transport (Kirtman et al., 2012). against MHT at 26.5°N. Here and in subsequent scatter plots show- At equilibrium, there is negligible storage so the positive and ing AMOC strength at 26.5°N, we also include the RAPID data for negative MHT slopes with respect to latitude in Fig.