Aug 22, 2018

Harmonic Feedback

Plenty of harmony in the pattern analysis over the past month. Interesting to note the uptick in correlation between RRWT projected and observed 500mb hPa anomaly. Attached is the 46-50d outlook verification. More verification maps here; http://www.consonantchaos.com/a-all.html



If there are any questions, comments, or suggestions on the material presented please let me know. Thank you for reading.

Jul 19, 2018

RRWT Methods Listing

Abstract
The impact of the closely related nature of planetary and synoptic waves is the motivation for this Solution. This Solution attempts to identify periodic upper atmospheric modes that subscribe to repeated weather occurrences.

Introduction
The recurring behavior of observed upper air patterns is reflected as future weather based on the empirical data via statistical modeling. In this pattern correlation phase space, each point represents a possible instantaneous state of the system. A solution of the governing equations is represented by a point traveling along a trajectory in the phase space. A single point in phase space determines the entire future trajectory, providing a composite of potential future weather and climate outcomes. [1]

Background
Studies have attempted to identify recurrent modes that contribute to the alternation between active and inactive periods of subseasonal regimes. Generally, two components of variability were found to be dominant on an intraseasonal time scale: the “short-term” component, with a period of less than 1 month, and the “long-term” component, with a period of 1–2 months but apparently less than a season. Both can be referred to as the intraseasonal oscillations (ISOs) since they are clearly separated from the high-frequency synoptic-scale (less than 10 days) features and the slowly varying seasonal cycle as well. [2]

Methods
Recurrent modes can be well observed from geopotential height (500 hPa) data. Here the method uses observations of daily 00Z 500mb from 2012 to present with a spatial resolution of 1° latitude from 45N at 10° longitude increments from 0E around the globe. Satellite observations of daily 500 hPa data provided by Navy Operational Global Atmospheric Prediction System (NOGAPS; available online at erdNogaps1D500mb). The dataset has a spatial resolution of 1° latitude 1° longitude with global coverage. [3]

Daily frequencies are obtained by correlating “current” patterns with “previous” patterns at lengths of 20 to 110 days for seven ~50° longitudinal wide regions. These correlation data are then reconstructed with respect to a daily top 10 for any calendar year. These daily top 10 correlations then construct a daily running average for the given calendar year (GAxA).

With the above criteria a global daily top 10 is constructed from each of the seven regions. These top 10 are then constructed into a daily running average for the given calendar year and a 15 day mode from present (M15D). Along with these frequencies, a daily average of the regional running averages is also constructed (RAxD). To overcome an evident intraseasonal jet position bias the method also seeks the top correlated frequency interval from the previous 720 days based from the present RAxD frequency (BeOP).

These frequencies (GAxA, M15D, RAxD, BeOP) are then used to seek analogs of observed wave trains and extend the encompassing weather from the present day. Deriving analog dates is accomplished by adding the forecast range to the present day then subtracting the interval rate. A three and four week forecast generated on 3/26 would provide a forecast for 4/10-23. Using an interval rate of 39, for example, would provide analog dates of 3/2-15.

Results
A time table displaying daily top 10 pattern correlation shows consistent features within the results. High correlated patterns centered on the 30-60 day period with top correlation greater than a 0.70 coefficient. Some interesting features are found with the time table as well. The highest coefficients form harmonious modulation of a seasonal standing wave. In addition, the 30-60 day interval displays moderate flux in correlation power with biannual dips centered on the equinox. During this time a lower frequency, interval greater than 60 days, over powers the 30-60 day interval. Also observed are periods when higher frequency, lower than a 30 day interval, are the only active intervals. [Figure 1]

Conclusions
Planetary friction and blocking will cause temporal and spatial difference in sensible weather returns. These challenges continue to advocate trial and error experiments within the framework. One theorized action incorporated a high correlated synoptic wave of 10-20 days to aid in sudden wave change. This idea proved too delicate in long range projection as it relied heavily on persistence. Another theorized action extended beyond the intraseasonal timescale and sought correlation of the average interval rate within the past two annual cycles. This idea proved robust in long range projection as it maintained a long-term longwave with seasonal sensibility.

Improvements to the framework correlation domain are desired, such as greater data points, finer grid resolution, and able mathematics. Utilization of the recurrent nature of the atmosphere can provide insight to future waveguides. The novelty of this approach is the ability to anticipate long range patterns based on recently observed patterns.

References
1. Xubin Zeng, Roger A. Pielke and R. Eykholt. Chaos Theory and Its Application to the Atmosphere. Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado. 1993;74:631-44. [Bulletin of the American Meteorological Society]

2. Bin Guan and Johnny C. L. Chan. Nonstationarity of the Intraseasonal Oscillations Associated with the Western North Pacific Summer Monsoon. Laboratory for Atmospheric Research, Department of Physics and Materials Science, City University of Hong Kong, Hong Kong, China. 2006;19:622-28. [Bulletin of the American Meteorological Society]

3. NOGAPS 500mb data. Prior 2012 500mb data were gathered from NOAA/ESRL Radiosonde Database available online at https://ruc.noaa.gov/raobs/.

Figure 1

Jun 10, 2018

Bering Sea Rule Verification

For the Bering Sea Rule (Renken et al.) fans out there, I am slowly automating verification results. Not sure what metrics to correlate or at what time frequency, but it is a start. Attached is an example of a 5-day H5 anomaly for 6/2-6.



A popular way to visualize the Bering Sea Rule is to utilize the maps @crankywxguy generates (http://www.stormhamster.com/bsr/bsr.php) which showcases daily surface and 500MB. I may just run with those metrics for the time being and concentrate on the "pattern" scoring 5-day anomalies.

For those interested in the verification method; the verification data are generated from ESRL NetCDF files utilizing the scorr function within GrADS. A snippet of the GrADS script is below. Beware, my GrADS scripting is amatuer at best. Enjoy!

function script(args)
myVarA=subwrd(args,1)
myVarB=subwrd(args,2)
myVarC=subwrd(args,3)
myVarD=subwrd(args,4)

'reinit'
rc = gsfallow("on")
nrows=2
ncols=1
args=nrows%' '%ncols
'sdfopen Z:\VMSharedFolder\Projects\Verification\esrl-ncdeflect\'myVarC
'set mproj latlon'
panels(args)

p=1
_vpg.1
'set parea 1 9.8 1 7.5'
'set t 1'
'set grads off'
'set gxout shaded'
'set rgb 16 0 0 255'
'set rgb 17 55 55 255'
'set rgb 18 110 110 255'
'set rgb 19 165 165 255'
'set rgb 20 220 220 255'
'set rgb 21 255 220 220'
'set rgb 22 255 165 165'
'set rgb 23 255 110 110'
'set rgb 24 255 55 55'
'set rgb 25 255 0 0'
'set rgb 26 255 250 250'
'set rgb 27 125 125 125'
'set clevs -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 16 17 18 19 20 27 21 22 23 24 25'
'set lon 140 210'
'set lat 39 65'
'set frame off'
'set mpdraw on'
'set gxout shaded'
'set csmooth on'
'd hgt1'
'd scorr(hgt1,hgt2,lon=0,lon=75,lat=57.5,lat=90)'
myVar=result
'run cbarc.gs'
'draw string 1.6 .35 Bering Sea Rule Intialized: 'myVarB
'draw string 5.8 .35 Bering Sea Rule Forecast Dates: 'myVarA
'draw string 2.8 4.12 500mb GPH (m) Anomaly Spatial Correlation (Bering Sea to CONUS Teleconnection)'
'draw title 'myVar

p=2
_vpg.2
'set parea 1 9.8 1 7.5'
'set t 1'
'set grads on'
'set gxout shaded'
'set rgb 16 0 0 255'
'set rgb 17 55 55 255'
'set rgb 18 110 110 255'
'set rgb 19 165 165 255'
'set rgb 20 220 220 255'
'set rgb 21 255 220 220'
'set rgb 22 255 165 165'
'set rgb 23 255 110 110'
'set rgb 24 255 55 55'
'set rgb 25 255 0 0'
'set rgb 26 255 250 250'
'set rgb 27 125 125 125'
'set clevs -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 16 17 18 19 20 27 21 22 23 24 25'
'set lon 232.5 302.5'
'set lat 24.5 50'
'set frame off'
'set mpdraw on'
'set gxout shaded'
'set csmooth on'
'd hgt2'
'draw string 3.9 .35 CONUS Observed Dates: 'myVarA
'draw string 0.05 .27 @OSNW3'

'set vpage off'
'printim Z:\VMSharedFolder\Projects\Verification\result-images\'myVarD
'set gxout print'
'set prnopts %6.2f'
write('bsr-hgt.txt', myVarB' 'myVar, append)
;


Also, there is some pre-GrADS NetCDF magic to line up the two domains for use with the scorr function. If there are any questions, comments, or suggestions on the material presented please let me know. Thank you for reading.