observational data analysis (Supplementary Material)

Similar documents
GAMINGRE 8/1/ of 7

Identifying the MJO Skeleton in Observational Data

particular regional weather extremes

Technical note on seasonal adjustment for M0

Annual Average NYMEX Strip Comparison 7/03/2017

Free and convectively coupled equatorial waves diagnosis using 3-D Normal Modes

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.

Time Series Analysis

ATSR SST Observations of the Tropical Pacific Compared with TOPEX/Poseidon Sea Level Anomaly

Jayalath Ekanayake Jonas Tappolet Harald Gall Abraham Bernstein. Time variance and defect prediction in software projects: additional figures

Interannual variation of MODIS NDVI in Lake Taihu and its relation to climate in submerged macrophyte region

Stratospheric planetary wave reflection and its influence on the troposphere

Average 175, , , , , , ,046 YTD Total 1,098,649 1,509,593 1,868,795 1,418, ,169 1,977,225 2,065,321

Average 175, , , , , , ,940 YTD Total 944,460 1,284,944 1,635,177 1,183, ,954 1,744,134 1,565,640

AdvAlg9.7LogarithmsToBasesOtherThan10.notebook. March 08, 2018

UNCLASSIFIED. Environment, Safety and Health (ESH) Report

In Centre, Online Classroom Live and Online Classroom Programme Prices

Record date Payment date PID element Non-PID element. 08 Sep Oct p p. 01 Dec Jan p 9.85p

WHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and Rainfall For Selected Arizona Cities

The Spectrum of Broadway: A SAS

Mr. XYZ. Stock Market Trading and Investment Astrology Report. Report Duration: 12 months. Type: Both Stocks and Option. Date: Apr 12, 2011

Technical note on seasonal adjustment for Capital goods imports

From El Nino to Atlantic Nino: pathways as seen in the QuikScat winds

South Eastern Australian Rainfall in relation to the Mean Meridional Circulation

Predictability of Sudden Stratospheric Warmings in sub-seasonal forecast models

Global Climates. Name Date

Modal view of atmospheric predictability

Variability and trend of the heat balance in the southeast Indian Ocean

ENGINE SERIAL NUMBERS

Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC

PICES XIV (2005) in Vladivostok

2003 Water Year Wrap-Up and Look Ahead

SYSTEM BRIEF DAILY SUMMARY

Momentum Budget Analysis of the Migra5ng Diurnal Tide in WACCM4: Seasonal Varia5on

Modulation of the Semiannual Oscillation Induced by Sudden Stratospheric Warming Events

SPECIMEN. Date Morning/Afternoon. A Level Geography H481/01 Physical systems Sample Question Paper. Time allowed: 1 hour 30 minutes PMT

Dynamics of the middle atmosphere at low, mid and high latitudes observed by the microwave wind radiometer WIRA

Eddy and Chlorophyll-a Structure in the Kuroshio Extension Detected from Altimeter and SeaWiFS

Jackson County 2018 Weather Data 67 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center

WindPRO version Dec 2005 Printed/Page 11/04/2008 1:54 PM / 1. Meteo data report, height: 66.0 Feet

2.1 Inductive Reasoning Ojectives: I CAN use patterns to make conjectures. I CAN disprove geometric conjectures using counterexamples.

Fig.3.1 Dispersion of an isolated source at 45N using propagating zonal harmonics. The wave speeds are derived from a multiyear 500 mb height daily

The Effects of Weather on Urban Trail Use: A National Study

Variability and trends in daily minimum and maximum temperatures and in diurnal temperature range in Lithuania, Latvia and Estonia

SYSTEM BRIEF DAILY SUMMARY

Jackson County 2013 Weather Data

Weak Constraints 4D-Var

SHADOW - Main Result. windpro CUMULTATIEVE EFFECTEN SLAGSCHADUW HERENTALS. EDF Luminus Markiesstraat Brussel

Long-term Water Quality Monitoring in Estero Bay

Algae and Dissolved Oxygen Dynamics of Landa Lake and the Upper Spring Run

Winter Season Resource Adequacy Analysis Status Report

GTR # VLTs GTR/VLT/Day %Δ:

Changes in spatial distribution of chub mackerel under climate change: the case study using Japanese purse seine fisheries data in the East China Sea

July 2017 LOGISTICAL HARMONY

Probabilistic predictions of monsoon rainfall with the ECMWF Monthly and Seasonal Forecast Systems

Meteo data report, height: 20.0 m. Name of meteo object: Ugashik

How are adding integers and subtracting integers related? Work with a partner. Use integer counters to find 4 2. Remove 2 positive counters.

Optimal Spectral Decomposition (OSD) for Ocean Data Analysis

Climatography of the United States No

In this activity, students will compare weather data from to determine if there is a warming trend in their community.

Salem Economic Outlook

El Niño, South American Monsoon, and Atlantic Niño links as detected by a. TOPEX/Jason Observations

Dates and Prices ICAEW - Manchester In Centre Programme Prices

Global Ocean Monitoring: Recent Evolution, Current Status, and Predictions

Calculations Equation of Time. EQUATION OF TIME = apparent solar time - mean solar time

FEB DASHBOARD FEB JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC

REPORT ON LABOUR FORECASTING FOR CONSTRUCTION

DEPARTMENT OF THE ARMY MILITARY SURFACE DEPLOYMENT AND DISTRIBUTION COMMAND (SDDC) 1 SOLDIER WAY SCOTT AFB, IL 62225

CWV Review London Weather Station Move

Habitat Suitability for Forage Fishes in Chesapeake Bay

Climate Variability in South Asia

DAILY QUESTIONS 28 TH JUNE 18 REASONING - CALENDAR

Winter Forecast. Allan Huffman RaleighWx

BRADSHAW'S RAILWAY GUIDE : accessible copies

Computing & Telecommunications Services

= observed volume on day l for bin j = base volume in jth bin, and = residual error, assumed independent with mean zero.

Quantification of energy losses caused by blade icing and the development of an Icing Loss Climatology

Computing & Telecommunications Services Monthly Report January CaTS Help Desk. Wright State University (937)

Climate Change and Arizona s Rangelands: Management Challenges and Opportunities

REGIONALIZATION OF PRECIPITATION ANOMALIES ASSOCIATED WITH EL NIÑO EVENTS IN SOUTHERN SOUTH AMERICA

Life Cycle of Convective Systems over Western Colombia

Hazard assessment based on radar-based rainfall nowcasts at European scale The HAREN project

Climatography of the United States No

Climatography of the United States No

ENSO Outlook by JMA. Hiroyuki Sugimoto. El Niño Monitoring and Prediction Group Climate Prediction Division Japan Meteorological Agency

OVERVIEW OF IMPROVED USE OF RS INDICATORS AT INAM. Domingos Mosquito Patricio

Orbit Design Marcelo Suárez. 6th Science Meeting; Seattle, WA, USA July 2010

Spatiotemporal variations of alpine climate, snow cover and phenology

Weather Products for Decision Support Tools Joe Sherry April 10, 2001

Suan Sunandha Rajabhat University

BAYESIAN PROCESSOR OF ENSEMBLE (BPE): PRIOR DISTRIBUTION FUNCTION

How to find Sun's GHA using TABLE How to find Sun's Declination using TABLE 4...4

Status SC23 Predictions SC24

The effect of ocean mixed layer depth on climate in slab ocean aquaplanet ABSTRACT

The MODIS Cloud Data Record

Where Was Mars At Your Birth?

Climatography of the United States No

Climatography of the United States No

Climatography of the United States No

Climatography of the United States No

Transcription:

1 Assessing the equatorial long-wave approximation: asymptotics and observational data analysis (Supplementary Material) 3 by H. Reed Ogrosky and Samuel N. Stechmann 4 1. Kelvin and MRG waves 5 6 The analysis of section 6 is repeated here for the φ and φ 1 components of equation (5.) from the main text corresponding to the Kelvin, MRG, and EIG waves. 7 a. Definition of wave variables 8 1) KELVIN WAVE 9 The φ component of (5.) is a single PDE governing a Kelvin wave, t r + x r =. (S1) 1 ) MRG AND EIG WAVES 11 The φ 1 component of (5.) is a system of two coupled PDEs governing an MRG and EIG wave, 1 t r 1 + x r 1 v =, t v + r 1 =. (Sa) (Sb) 13 Each of the variables in (S) can be expressed as a superposition of plane-wave ansatzes ˆr 1 (k,ω)e i(kx ωt), ˆv (k,ω)e i(kx ωt) ; (S3) 1

14 substituting (S3) into (S) gives i(k ω) 1 1 iω ˆr 1 ˆv =. (S4) 15 We are interested in finding the eigenmodes of (S4); its characteristic equation is ω kω 1 =. (S5) 16 17 There are two solutions, ω j, for j {MRG,EIG } to (S5). Each eigenvalue ω j for j {MRG,EIG } is associated with an eigenvector of the form ( ê j = i ) T, 1. (S6) k ω j 18 19 The resulting eigenvectors are shown in Fig. S1 after normalization. The degree to which the spatial structure of each of these two waves is seen in reanalysis data can be calculated by the projection technique described in section 6. Using this approach, the 1 Fourier coefficients of the MRG and EIG wave structures are defined as MRG(k) = ê MRG (k) ˆr 1 (k), ÊIG (k) = ê EIG (k) ˆv (k) ˆr 1 (k) ˆv (k), (S7) 3 respectively, where crosses denote the conjugate transpose. This spectral data may then be trans- formed back into physical space through an inverse Fourier transform. 4 b. Long-wave theory with wave variables 5 1) KELVIN WAVE 6 The Kelvin wave is unaffected by the long-wave approximation.

7 ) MRG AND EIG WAVES 8 Projecting (5.7), i.e. with δ 1, onto φ 1 results in two coupled PDEs, t r 1 + x r 1 v =, δ t v + r 1 = ; (S8a) (S8b) 9 In the limit of small δ, the system (S8) can be expressed as an eigenvalue problem (written here in terms of ˆv = δ ˆv ), 3 iδ(k ω) 1 1 iδω ˆr 1 ˆv =. (S9) 31 We are interested in finding the eigenmodes of (S9); its characteristic equation is δ ω δ kω 1 =. (S1) 3 33 34 There are again two solutions, ω j, for j {MRG,EIG } to (S1), but in the limit δ both of these roots are singular. Each eigenvalue ω j is associated with an eigenvector of the form ( ê j = i ) T, δ. (S11) k ω j In the long-wave limit, approximate eigenvalues may be found by expanding in powers of δ, ω MRG = δ 1 + O(1), ω EIG = δ 1 + O(1), (S1a) (S1b) 35 36 After a phase shift so that the r 1 component is positive and real, the normalized long-wave eigen- vectors are, to leading order in δ, given by MRG : ê MRG = ( 1, i ) T, EIG : ê EIG = ( 1, i ) T. (S13a) (S13b) 37 Note that the structure of these long-wave eigenvectors is independent of k. 3

38 c. Observational data analysis 39 1) KELVIN WAVE 4 41 4 43 A Hovmoller plot of the Kelvin wave structure isolated in reanalysis data is shown in Fig. S(a) for the one year period 1 July 9 through 3 June 1. Note the abundance of information propagating rapidly to the east, reminiscent of Kelvin waves. The power spectrum of the Kelvin wave structure is shown in the main text in Fig. 18(a). 44 ) MRG AND EIG WAVES 45 46 47 48 49 5 51 5 53 54 Fig. S(b-c) shows the non-long-wave MRG and EIG wave structures isolated in reanalysis data. The MRG wave exhibits the strongest activity between 18 and 9W, similar to the location of highest activity in the meridional winds v (see Fig. 4b in the main text). This can be anticipated in light of the MRG wave being comprised almost entirely of v at moderate to high wavenumbers, as in Fig. S1(a). Note also the large discrepancy between the MRG and EIG structures in amplitude: the EIG wave contains very little power at all frequencies and wavenumbers, while the MRG wave contains significant power at moderate wavenumbers. Both waves contain very little low-wavenumber information. These figures demonstrate that reanalysis data projects weakly onto the MRG and EIG waves over spatial and temporal scales where the long-wave approximation holds. 4

55 56 LIST OF FIGURES Fig. S1. Eigenvector components of the (a) MRG and (b) EIG wave structures........ 6 57 58 Fig. S. Hovmoller plot of (a) Kelvin, (b) MRG, and (c) EIG anomalies from a seasonal cycle. Time period shown is 1 July 9 through 3 June 1.............. 7 5

.5 1 (a) Components of ê MRG Re(r 1 ) Im(v ).5 1 (b) Components of ê EIG -.5 -.5-1 -1 1 k -1-1 1 k FIG. S1. Eigenvector components of the (a) MRG and (b) EIG wave structures. 6

E I G M RG K Jun May Apr Mar Feb Jan Dec Nov Oct Sep Aug 59 6 Jul 9E 18 9W (a).5.5 9E 18 9W.5 (b).5 (c).5 9E 18 9W.5 F IG. S. Hovmoller plot of (a) Kelvin, (b) MRG, and (c) EIG anomalies from a seasonal cycle. Time period shown is 1 July 9 through 3 June 1. 7