Climatological analysis of tornado report counts using a hierarchical Bayesian spatiotemporal model
|
|
- Cecily Poole
- 5 years ago
- Views:
Transcription
1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. D4, 9005, doi:10.109/00jd00806, 003 Climatological analysis of tornado report counts using a hierarchical Bayesian spatiotemporal model Christopher K. Wikle Department of Statistics, University of Missouri, Columbia, Missouri, USA Christopher J. Anderson Department of Agronomy, Iowa State University, Ames, Iowa, USA Received 8 July 00; revised 18 March 003; accepted 18 April 003; published 31 October 003. [1] Recent evidence suggests that tornado report counts for monthly or longer periods may be correlated to climate indices in space and time. However, the climatological analysis of tornado reports is complicated by reporting errors, by the non-gaussian nature of count data and rare events, and by the presence of spatial and temporal correlation. Typically, such factors have not been formally included in the underlying statistical model used for inferential decisions. We present an analysis in which a statistical model is constructed so that the aforementioned characteristics of tornado reports are explicitly accommodated in the model structure. In particular, a hierarchical Bayesian framework is considered, in which the various complicated structures are considered in a series of conditional models, formally linked by basic probability rules. This formalism allows one to evaluate characteristics of the spatiotemporal variability of an underlying (unobservable) tornado count process given the noisy observations (tornado reports). In addition, one can consider the effects of exogenous climate processes on this tornado count process. We find that an index of El Niño activity is significantly associated with tornado reports over the continental U.S. and that there is substantial regional variability in this relationship. In addition, we find evidence of temporal trends in tornado counts, with spatial variation in the magnitude and sign of the trend. INDEX TERMS: 160 Global Change: Climate dynamics (3309); 310 Mathematical Geophysics: Modeling; 3399 Meteorology and Atmospheric Dynamics: General or miscellaneous; KEYWORDS: Bayesian, hierarchical, Poisson, uncertainty, spatiotemporal Citation: Wikle, C. K., and C. J. Anderson, Climatological analysis of tornado report counts using a hierarchical Bayesian spatiotemporal model, J. Geophys. Res., 108(D4), 9005, doi:10.109/00jd00806, Introduction [] Interest in tornado climatology stems from the destructive potential of tornadoes. Clearly, the loss of human life caused by tornadoes is reason enough to obtain greater understanding. In addition, the estimated cost of damage caused by tornadoes is close to one billion per year in 1999 U.S. dollars [Brooks and Doswell, 001] (see also R. A. Pielke Jr., Extreme Weather Sourcebook 001 Edition, even though most tornadoes occur in rural areas. It is important to gain a fundamental understanding of the climatological behavior of such phenomena in order to develop mitigation strategies for interests ranging from realistic assessment of tornado variability under different climate regimes to seasonal prediction on a regional scale. Because physical understanding of climate influences on the meteorological phenomena at the tornado scale is Copyright 003 by the American Geophysical Union /03/00JD00806$09.00 limited, it is useful to approach this problem from a statistical perspective Summary of Recent Literature [3] Climatological analysis of tornado counts is complicated by the quality of tornado reports archived by the National Weather Service (NWS). As summarized by Doswell and Burgess [1988], most reports come from relatively untrained witnesses, with only a few events prior to the late 1970s subject to on-site inspection by trained specialists. Thus there are two primary limitations of such a database: reliance on human observation and a variable rating procedure. That is, the extent to which this database represents the true occurrence of tornadoes rests entirely on the willingness of humans to look for and report tornadoes and on the ability of humans to assign accurate assessments of tornado strength. It is very likely that some tornadoes are not witnessed, especially those that occur in rural areas and/or at night, while others do not cause damage. Thus the analysis discussed herein constitutes an examination of tornado reports rather than tornado occurrence. STS 1-1
2 STS 1 - WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL [4] An overview of statistical characteristics of the severe weather report archive maintained by the NWS was provided by Brooks [1998]. He found that severe thunderstorm report intensities have increased dramatically since Tornado reports were found to exhibit similar changes over time, although regional differences were noted. For example, Florida and the front range counties of Colorado showed different reporting trends than other regions of the U.S. [5] An approach for estimating the daily climatological probability (i.e., average frequency over a long time period) of F F5 tornadoes (where the F stands for the Fujita scale, a rating system for tornado damage, in which F0 represents minimal damage and F5 represents complete destruction ) is proposed by Concannon et al. [000]. Their analysis was based on a more comprehensive database of F F5 tornadoes compiled by Grazulis [1993] for the period They found that (1) the maximum number of tornado days per century occurs in an L-shaped pattern, from southeastern Mississippi, to southern Oklahoma, to southwestern Iowa; () fluctuations as large as 5% of the climatological frequency have occurred, but a national trend is not evident, although it was noted that slightly more tornado days occurred in Oklahoma and Kansas during the early rather than late portion of the record; (3) the annual peak frequency exhibits less interannual variability at points in the central rather than eastern U.S.; and (4) peak frequency occurs in March and April in the southeast U.S. and progresses from Texas to Canada from April to July. Thus their results support the presupposition that climate fluctuations impact tornado report frequency. [6] Bruening et al. [00] apply the analysis procedure of Concannon et al. [000] to F0 F5 tornado reports in order to examine year-to-year variations of tornado counts. They documented a clearly increasing trend in national tornado reports of 14 tornadoes per year, based on simple linear regression results. These authors also identified significant interannual and intraannual variability in tornado reports. [7] Monfredo [1999] examined the association between tornado reports in the midsouth and southern plains of the U.S. and the Southern Oscillation index (SOI) using reports from the NWS log. Specifically, he used dew point temperature and precipitation to define two regions that had similar meteorological characteristics. The regions were labeled (1) southern low plains (SLP), which included Missouri, Arkansas, the portion of Kansas east of the 100th meridian, all of Oklahoma except the panhandle, and north central Texas, and () midsouth (MSO), which contained Kentucky, Tennessee, and the northern halves of both Alabama and Mississippi. He considered only F F5 tornadoes from February to July (his tornado season ) of each year. SOI was accumulated over the canonical El Niño period of March February, preceding the tornado season. Spearman s rank correlation test showed significantly lower counts of F F5 tornado reports during El Niño and significantly higher counts during La Niña in both the SLP and MSO regions. The time lag between canonical El Niño and tornado seasons somewhat obscures the interpretation of this result, but Monfredo s [1998] previous work suggests a physical rationale for such findings. [8] Marzban and Schaefer [001] provided the most comprehensive analysis to date relating tornado counts and El Niño. They noted that recent studies presented at conferences on this subject showed conflicting results. A concern was that in each of these studies, El Niño was based on indices that differed in duration, underlying variables, timing with respect to the tornado season, and the indicator of anomalous behavior. In order to mitigate ambiguity of this type, they averaged equatorial Pacific sea surface temperatures (SST) in four zones and examined the monthly contemporaneous correlation between these SST indices with both the number of tornadoes and the number of tornado days that occur nationally and within certain regions of the U.S. They considered data from the NWS log for the years , and they found increasing trends in national counts of all (F0 F5) tornado reports related to all SST regions and a decreasing trend in national counts of F F5 tornado reports. They speculated that the increasing trend in F0 F5 tornado counts was due to increasing population and improvements in detecting tornadoes and that the decreasing trend of F F5 counts may have been related to a procedural change in the manner by which F-scale ratings were assigned. [9] Marzban and Schaefer [001] used Kendall s t nonparametric measure of association to examine the relationship between tornado activity and the various SST indices. For tornado activity they considered separately total report counts as well as number of days with tornadoes. They considered two intensity levels, F0 F5 ( all tornadoes ) as well as F F5 ( strong and violent ). They found relatively weak correlations, in general, after removing seasonal effects. However, they did find statistically significant correlations between SST and tornado activity, with the level of significance varying with location over the U.S. and spatial region corresponding to the SST index. In particular, the most significant correlation ( p value = ) was between SST in the eastern equatorial Pacific and the number of tornadic days in the northeast U.S. Furthermore, the SST in the central equatorial Pacific was found to be significantly correlated with strong and violent tornado counts (and tornado days) in the northeast U.S., citing results from Hoerling et al. [1997] as physical evidence in support of this result. In particular, timeaveraged upper-level wind fields show that the jet stream over the northeastern U.S. is stronger when the equatorial Pacific SSTs are cooler than normal. The central and southeast U.S. showed less significant associations in general. They concluded that tornado report frequency and equatorial Pacific SST are weakly negatively correlated. 1.. Study Motivation [10] The review of recent literature presented in section 1.1 suggests that statistical analyses of tornado climatology should account for: (1) observational uncertainty present in tornado report counts; () spatial variability in the longterm frequency of tornado counts; (3) temporal trends that may vary with space; and (4) exogenous climatological factors, primarily tropical Pacific SST, that are associated with tornado reports, with the potential for such effects to vary spatially. [11] Indeed, the climatological analysis of tornado reports is not only complicated by reporting errors but by the non-
3 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL STS 1-3 Gaussian nature of count data and rare events and by the presence of spatial and temporal correlation. The analyses to date have estimated sample statistics, such as sample mean frequency or sample correlation between tornado counts and climate indices. Hypothesis tests have often been applied to these estimates of sample statistics. In most cases, either an underlying assumption is made that the observations are independent, or the observations are grouped into relatively large geographical regions, potentially smearing the signal. Although some studies have accounted for some of the complicating factors mentioned above (e.g., preprocessing the data to remove spatial and temporal correlation), all of these factors have not been formally included in the underlying statistical model used for inferential analysis. [1] We take an alternative approach, in which a statistical model is constructed so that the aforementioned characteristics of tornado reports are explicitly accommodated in the model structure. In particular, a hierarchical Bayesian framework is considered, in which the various complicated structures are considered in a series of conditional models, formally linked by basic probability rules. This formalism allows one to evaluate characteristics of the spatiotemporal variability of the underlying (unobservable) tornado count process given the noisy observations. Such a model can easily incorporate spatially varying climatological effects (e.g., overall spatial mean and trend) as well as spatially varying effects of temporally varying exogenous variables (e.g., El Niño) on the underlying tornado count process.. Data [13] Tornado reports submitted to the NWS from 1950 to 1995 are archived at the National Oceanic and Atmospheric Administration Storm Prediction Center. We have extracted all reports from 1953 to 1995 for our analysis,since tornado reports prior to 1953 were very infrequent. This was also the first year of issuance of severe weather watches by the NWS. A short list of information contained in each report includes the time and location of the initial and final observation of the tornado, estimated tornado path width, estimated tornado intensity, estimated dollar amount of damage produced by the tornado, and other forms of severe weather observed in the immediate vicinity of the tornado. The accuracy of estimated observations, such as tornado intensity or cost of damage, are open to debate since they rely upon interpretation by human observers with varying levels of expertise in assessing these characteristics [Doswell and Burgess, 1988]. In this study we examine the initial position of each tornado report, which is considered the most reliable entry. [14] We overlaid a grid (400 km 1700 km) of 50-km boxes on the continental U.S., with the upper-left grid box centered at 45.5 N, 105 W and with 54 grid boxes in the east-west direction and 34 in the north-south direction. Because tornado reports were very infrequent within the Rocky Mountains and along the West Coast, we considered only tornado report counts east of 105 W. For every year we tallied the number of initial observations of an F F5 tornado that occurred within each box (hereinafter referred to as tornado report counts), resulting in a time series of length 43 years for each of the 1836 grid boxes. Figure 1 shows the total counts observed in each grid box over the 43 years. Note that there are regions with no observations (e.g., oceanic areas). Furthermore, the NWS database does not include tornado report counts over Canada. In general, there are a significant number of zero counts over the domain (i.e., 90% of the observations are zero). [15] In this study we considered the spatiotemporal distribution of the tornado report counts and the relationship of these counts to an index of the El Niño/Southern Oscillation (ENSO) phenomenon, as given by the Niño3.4 index [e.g., Trenberth, 1997]. These data were obtained as monthly time series from the National Centers for Environmental Prediction Climate Prediction Center web site. Yearly average (mean) values of the index were used in the study so as to correspond to the yearly tornado report counts. 3. Methodology [16] Tornado counts are distinctly non-gaussian (i.e., counts are discrete), have significant sampling and observational variability, contain many zeros, and exhibit relatively complicated spatiotemporal dependence. These issues make traditional regression methods intractable. In contrast, Bayesian hierarchical models attempt to decompose the observed data into a series of conditional models, all linked together formally through basic probability relationships. In this way, one can account for the observational process error, the true process variability (e.g., the variation in true counts, which can never be observed), and parameter uncertainty (e.g., the strength of correlation among parameters). [17] The basic hierarchical model is then written in three primary stages: (1) data model ([datajprocess, parameters]); () process model ([processjparameters]); and (3) parameter model ([parameters]), where the brackets refer to the probability distribution and the vertical line indicates a conditional distribution. Thus the fundamental idea is to approach the complex statistical modeling problem by breaking it into pieces [e.g., Berliner, 1996; Wikle et al., 1998]. [18] The first stage is concerned with the observational process, or data model, which specifies the distribution of the data (e.g., observational counts of tornadoes) given the process of interest (e.g., the true tornado count over some geographical area) and parameters that describe the data model. The second stage then describes the process, conditional on other parameters. For example, this might be a regression model, relating tornado occurrence to climate indices, with parameters describing the strength of the association. Finally, the last stage accounts for the uncertainty in the parameters by assigning them distributions. For example, the effect of a climate index could be different at different geographical locations. Thus we might allow the regression parameters to be spatially varying and assign them a distribution that includes a spatial correlation. [19] Ultimately, our interest is with the distribution of the process and parameters, updated by the data. This so-called posterior distribution is obtained via Bayes s theorem [e.g., Gelman et al., 1995]. Specifically, the distribution [process, parametersjdata] is proportional to the product of the data model, process model, and parameter model from the three stages listed above: ½process; parametersjdata Š / ½datajprocess; parametersš ½processjparametersŠ½parametersŠ: ð1þ
4 STS 1-4 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL Figure 1. (top) Log counts of F F5 tornadoes summed over the years (bottom) Time series of F F5 tornado counts for each year summed over the spatial domain shown in Figure 1 (top). This formula serves as the basis for estimation in Bayesian hierarchical models. Although simple in principle, there are several critical issues. First, the specification of the component distributions in equation (1) can be challenging. However, this is the case with any statistical model, and the specification of each stage requires the modeler to address issues about uncertainty directly (rather than implicitly, as in many traditional statistical models). Much of the uncertainty quantification is subjective (e.g., choice of distributions, parameterizations, etc.), but the quantification of such subjective judgment is, in fact, the strength of the Bayesian approach in that it provides a coherent probabilistic framework in which to incorporate scientific reasoning, judgment, and experience explicitly in the model (see Carlin and Louis [000] for an overview of issues related to Bayesian modeling). [0] The second critical issue concerns the calculation of the posterior distribution of equation (1). In principle, this is simple, but in practice, it can be very difficult to obtain the normalizing constant for the posterior distribution due to the complexity and high dimensionality of the model. However, although we cannot usually compute the posterior directly, we can simulate from it using Markov Chain Monte Carlo (MCMC) computational procedures. Such procedures for Bayesian computation have become common in recent years (e.g., see Gilks et al. [1996], Robert
5 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL STS 1-5 and Casella [1999], and Carlin and Louis [000] for overviews) Hierarchical Model for Tornado Counts [1] Poisson regression is often desired when the response outcomes are counts, especially in cases where large counts are relatively rare. In addition, one wishes to relate the expected mean response (i.e., the Poisson intensities) to a set of covariates. The general class of models for such a process are known as generalized linear models (see McCullagh and Nelder [1989] for an overview and Yan et al. [00] for an illustrative meteorological example). This procedure can be extended to accommodate random effects as well as fixed effects, in which case the class of models is known as generalized linear mixed models (e.g., see McCulloch and Searle [001] for an overview). In recent years, investigators have considered complicated random effects such as those induced by spatial random fields. Although there have been some classical implementations of such models [e.g., Heagerty and Lele, 1998], typically, their added complexity has led researchers to consider hierarchical Bayesian implementations [e.g., Diggle et al., 1998; Royle et al., 00; Best et al., 00; Christensen and Waagepetersen, 00; Crainiceanu et al., 003; Wikle, 00]. In the case of tornado count data we would like to relate the mean of the underlying Poisson count process to a spatially varying temporal trend as well as a spatially varying ENSO effect, where the regression coefficients are themselves spatial random processes. However, examination of the data suggests that there are too many zeros to justify a Poisson model for counts. Rather, we consider a socalled zero-inflated Poisson model, in which there is a probability p that the observation is from a Poisson random variable and a probability 1-p that the only possible observation is 0 [e.g., Lambert, 199; Hall, 000; Nair et al., 001]. The hierarchical model is described in detail in sections Data Model [] Let Y(s i ; t) be the tornado counts in the geographical region (i.e., grid box) indexed by s i, i =1,..., n at times t = 1,..., T (representing years 1953,..., 1995). Given that there are too many zero counts to justify a Poisson model, we utilize a zero-inflated Poisson data model [e.g., Lambert, 199] given by Yðs i ; t 8 < POIðlðs i ; tþþ with probability ps ð i ; tþ; Þjlðs i ; tþ; ps ð i ; tþ : 0 with probability 1 ps ð i ; tþ ðþ for all i and t. That is, conditional on the Poisson intensity (or mean) l(s i ; t), which is a spatiotemporal process, and on a spatiotemporal probability process p(s i ; t), the tornado counts are either 0 or follow a Poisson distribution. Thus it is a mixture distribution of a point mass at 0 and a Poisson distribution. The observations are conditionally independent. This does not suggest that the counts at two locations/ times (e.g., Y(s i ; t) and Y(s j ; t)) are independent but rather that their spatiotemporal dependence is due to the underlying intensity process, l, and probability process, p, which are modeled at the next stage. We say that the Y process is dependent marginally but is independent conditionally Process Models [3] The process models describe the underlying spatiotemporal processes, l(s i ; t) and p(s i ; t). First, we consider the l process, as it is the process of primary interest meteorologically. [4] In the present case we believe that the underlying spatiotemporal variability of the l process is partially explained by a linear temporal trend and by an ENSO effect. Furthermore, there are spatiotemporal random effects that account for the additional variability. That is, logðlðs i ; tþþ ¼ b 1 ðs i Þx 1t þ b ðs i Þx t þ hðs i ; tþþðs i ; tþ; ð3þ where H t (h(s 1 ; t),..., h(s n ; t)) 0 is a spatially correlated (but temporally uncorrelated) random process with mean D h, correlation matrix R(q h ), and variance s h, e(s i ; t) iid N(0,s e ), and where b k (s i ) is a regression parameter at location s i for the kth covariate, x kt, k = 1,. That is, the log of the Poisson intensity is modeled by a linear regression with normal errors such that the regression coefficients are spatially varying so that each site has its own impact from the given covariate. In this analysis, covariates include a year index (x 1t = {1,..., 43}) and the ENSO climate index (x t ), although, as described in section 3., these covariates are centered and standardized in the implementation. Thus b 1 represents the temporal trend, and b represents the ENSO effect. The important point here is that these coefficients are spatially varying as well. The error process represented by the term h(s i ; t) accounts for spatiotemporal random effects representing unobserved or unknown climate effects that are spatially correlated, whereas e(s i ; t) represents such random effects that are not spatially correlated (e.g., see Best et al. [00] for similar formulation in the spatial context). That is, these processes account for variability in the spatiotemporal Poisson intensity process not directly modeled by the aforementioned covariates. For example, these noise processes can partially account for geographical differences caused by sampling biases (e.g., demographic biases) and/or by unspecified climatological and meteorological factors. [5] The zero-inflated Poisson probability process p(s i ; t) is traditionally modeled on the logit scale [Lambert, 199; Hall, 000; Nair et al., 001]. In our case, logitðps ð i ; tþþ ¼ nðs i Þþgðs i ; tþ; ð4þ where logit(z) exp(z)/(1 + exp(z)) is the logit function, n(s i ) is a spatial random effect (specified at the next level of the hierarchy) and g(s i ; t) iid N(0, s g ) is an uncorrelated spatiotemporal error process. In this case we know less about the probability process, p(s i ; t), a priori than we do about the l(s i ; t) process. However, we believe that the underlying process that forces counts to deviate from a true Poisson process is likely to be spatial rather than temporal in nature. For example, areas over water and over regions for which we have no tornado reports (e.g., Canada) obviously have too many zeros. We will account for this in the mean structure for the spatially varying v process, as shown in section Parameter Models for the Poisson Mean Process, L [6] At this stage of the hierarchical model the parameters are given distributions. For example, it is our prior belief,
6 STS 1-6 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL based on the published literature (see the literature review in section 1.1), that the bs in the process model vary spatially and that they should exhibit spatial coherence. That is, b k jq k ; s k N 0; s k R kðq k Þ ; for k ¼ 1; ð5þ where B k is an n 1 vector of rastorized and spatially indexed coefficients corresponding to the kth covariate, s k is the variance of this parameter process, and R k (q k )isa spatial correlation matrix that itself depends on the parameter q k. Again, the important point here is that our prior belief is that there may be spatial structure in the regression coefficients. That is, if the regression coefficient in one grid box is high, we expect that the adjacent grid boxes will be more likely to show high coefficients as well. In this way, estimates of any given b k (s i ) can borrow strength from their neighbors, so the effective number of parameters actually estimated for b k (s i ), i =1,..., n is much less than n. However, the model must be flexible enough that if the data suggest that this is not the case, it can be reflected in the posterior analysis. Hence the spatial dependence of each b is controlled by the parameter q k, itself a random parameter with a specified distribution (see equation (15)). [7] Similarly, as mentioned in section 3.1., the h process is spatially correlated with prior distribution H t jd h ; s h ; q h N D h ; s h R q h : ð6þ [8] In practice, the implementation of the hierarchical model in the case with parameter model equations (5) and (6) can be very difficult because of the high dimensionality of the spatial processes (recall that n = 1836 in our case, implying, for example, that equation (5) is a 1836-dimensional multivariate normal distribution and R k (q k ) is Alternatively, as shown by Wikle [00], Bayesian computation of Poisson hierarchical spatial models is facilitated by considering spectral representations of the spatial processes. Specifically, we let B k = A k, and H t = X t, where is an n n matrix of spectral basis functions. Thus we can rewrite equation (3) as logðl t Þ ¼ A 1 x 1t þ A x t þ x t þ t ; ð7þ where, x 1t, and x t are known and equations (5) and (6) become A k jq k ; s k N 0; s k D ð q kþ ; for k ¼ 1; ð8þ x t jd; q x ; s x N D; s x D q x ; for all t; ð9þ where D(q k ) 0 R(q k ) and D(q h ) 0 R(q x ). [9] In principle, any orthogonal basis set can be chosen for. For example, possible choices for basis functions in these settings include empirical orthogonal functions [e.g., Obled and Creutin, 1986], wavelets [e.g., Wikle et al., 001], and Fourier basis functions [e.g., Wikle, 00]. The Fourier and wavelet modes are advantageous computationally when data are gridded (as in the present case) due to the computational efficiency of discrete transform algorithms (note that if is a matrix of Fourier basis functions, then A is just the inverse discrete Fourier transform operation on the vector A and 0 y is just the discrete Fourier transform (DFT) of y). Furthermore, if the spatial process is stationary, then use of the Fourier basis set implies that the matrix D(q k ) is asymptotically diagonal [e.g., Shumway and Stoffer, 000, section T3.1]. [30] In the present case we assume the correlation structure of these processes is isotropic and stationary and that they follow an exponential model (r(h) = exp( qh), where h is the Euclidean distance between sites) but allow the dependence parameter q to be random (modeled at the next stage of the hierarchy). Thus, conditional on the dependence parameter, one can obtain an analytical expression for the matrix D(q k ), as shown by Wikle [00]. That is, for the exponential correlation function corresponding to a spatial process in two dimensions the corresponding spectral density function at frequency w is f (w) / 1/(q + w ) [see e.g., Stein, 1999, p. 49]. Thus if are Fourier basis functions, then the prior structure for D(q) is assumed to be diagonal (this is true asymptotically and is thought to be a reasonable prior assumption here) and is defined completely by the spatial dependence parameter q. That is, this matrix is diagonal, with the diagonal elements corresponding to the spectral densities given by f (w). In practice, one must make sure that the rastorization (that takes the two-dimensional spatial field into a vector) of the two-dimensional spectral density corresponds with the particular DFT and inverse DFT algorithms used in the implementation. In addition, one must account for edge effects (e.g., by zero padding) as is typical for DFT applications. Finally, we note that the class of covariance functions that can be considered in this way is much broader than that presented here. For instance, Wikle [00] discusses the implementation of the Matern class of covariance functions, in which the inherent smoothness (i.e., differentiability) of the correlation function is parameterized as well as the spatial dependence. [31] Typically, spectral transformations, such as the Fourier transform described above, can be truncated without substantial loss of information. In our case we do not truncate as we have little prior information to suggest which coefficients are needed, especially since we are not sure how much spatial dependence should be in the spatial random fields. In addition, the DFT algorithm is very efficient and requires the full-dimensional spatial field for implementation. We note, however, as shown below, that many of the a coefficients are very close to zero in the posterior mean. [3] One might question the choice of a stationary isotropic spatial model when the spatial process covers such a broad geographical region as is the case here. This may be a strong concern for the l process, but we have little prior knowledge to suggest that the b processes or the h process are nonstationary. To account for spatial nonstationarity, one could use wavelet basis functions to model the spatial processes, as in the work of Wikle et al. [001] and Nychka et al. [00]. This will be explored elsewhere. [33] The D process has prior distribution given by h i Djs d;1 ; s d; N 0; diag s d;1 s d; 10 n 1 ; ð10þ
7 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL STS 1-7 Table 1. Prior and Posterior Summary of Univariate Model Parameters Prior Prior Parameter Mean Variance Prior Parameters where diag[z 0 ] refers to a diagonal matrix with elements of the vector z on the main diagonal and 1 n 1 is an n 1 dimensional vector of ones. The first component of the d-process has variance s d,1, and the other n-1 elements have homogeneous variance given by s d,. The motivation for this is our prior experience that the first (mean) component of spectral random fields is often substantially larger than the remaining elements. We specify inverse gamma (IG) prior distributions for all of the variance components: s e IG ð q e; r e Þ; ð11þ s x IG q x; r x ; ð1þ s k IG ð q k; r k Þ; k ¼ 1; ð13þ s d;j IG q d;j; r d;j ; j ¼ 1; ð14þ where the q and r hyperparameters are specified as given in Table 1. [34] Finally, we specify prior distributions for the spatial dependence parameters: q k Uu ð 1 ; u Þ; k ¼ 1; ð15þ q x Uu ð 1 ; u Þ; ð16þ where U(u 1,u ) refers to a uniform distribution between u 1 and u and where the hyperparameters u 1, u are specified as given in Table 1. Small values of q correspond to strong spatial dependence, and large values correspond to weak spatial dependence. Also, we use a standardized distance unit for the exponential model, in which the distance between adjacent grid boxes (50 km) is considered one distance unit Parameter Models for Zero-Inflated Probability Process, p [35] The spatial random process associated with the logit of p(s i ; t) is modeled as nðs i Þja 1 ; a ; s n Na ð 1z 1 ðs i Þ Posterior Mean Posterior Variance s e q e =.001, r e = s x q x =.05, r x = s q 1 =.000, r 1 = s q =.000, r = q u 1 = 0.1, u = q u 1 = 0.1, u = q x u 1 = 0.1, u = s d, q d,1 = 11.0, r d,1 = s d, q d, =.0, r d, = s g 0. q g =.00, r g = s n q e =.005, r e = a ~a 1 ¼ 0:1; ~s a ¼ a ~a ¼ 0:01; ~s a ¼ þ a z ðs i Þ; s n ; ð17þ independently for all i and where z 1, z represent dummy variables related to data-rich and data-poor regions. Specifically, z 1 (s i ) takes the value of 1 if the grid box is in a data-rich region and zero otherwise; z (s i ) takes the value of zero for data-sparse grid boxes and 1 otherwise. A grid box is assumed to be data-sparse if it is in a boundary region (e.g., water, Canada, western part of domain) for which we have no tornado observations over the 43 years of data coverage. We let the a parameters have prior distributions a j N ~a j ; ~s a ð18þ for j = 1,, where ~a j ; ~s a are hyperparameters specified as given in Table 1. Similarly, s n IG ð q n; r n Þ; ð19þ where the hyperparameters q n, r n are specified as given in Table 1. In addition, the error variance associated with the logit( p(s i ; t)) spatiotemporal error process, g(s i ; t), is given a prior s g IG(q g, r g ), where the hyperparameters are specified as given in Table Model Implementation [36] Given the hierarchical representation presented in sections , one can then evaluate the posterior distribution of all of the processes and parameters given the observed counts: L1 ;...; L T ; p 1 ;...; p T ; A 1 ; A ; X 1 ;...; X T ; N; a 1 ; a ; s ; s x ; s d;1 ; s d; ; s g ; s n ; s 1; s ; q 1 ; q ; q x jy Y T Y nt / ½y t ðs i Þjlðs i ; tþ; ps ð i ; tþš: t ¼ 1 i ¼ 1 lðs i ; tþja 1 ; A ; X t ; s h i psi ð ; tþjnðs i Þ; s g Y T h i h i X t ¼ 1 tjd; q x ; s s h i x x q x Djs d;1 ; s d; h ih i s d;1 s d; Nja 1 ; a ; s n ½ a1 ; a Š s n ; Y A k ¼ 1 kjq k ; s k ½ qk Š s s k ; ð0þ where Y represents all observed counts and where L t, and p t are n 1 vectorizations of the n spatial locations at time t for the Poisson means and zero-inflated probabilities, respectively. One cannot evaluate this posterior distribution analytically and must resort to numerical simulation methods. In particular, MCMC methods can be used to generate samples from the posterior distribution. See Wikle et al. [1998], Berliner et al. [000], and Wikle et al. [001] for examples of MCMC applied to atmospheric science problems. For MCMC implementations of Poisson spatial models, see Diggle et al. [1998], Best et al. [00], Christensen and Waagepetersen [00], and Crainiceanu et al. [003]. For MCMC implementation of spatiotemporal Poisson models, see Wikle [001] and Wikle [003]. In the present paper we utilize a special case of MCMC known as the Gibbs sampler. This is described in more detail in
8 STS 1-8 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL Figure. Histogram of posterior samples for the b parameter variances (left) s 1 and (right) s. Appendix A. For general discussions of the Gibbs sampler and MCMC, see Robert and Casella [1999], Carlin and Louis [000], and Gilks et al. [1996]. [37] There are a few issues regarding implementation that should be mentioned in greater detail. First, as discussed by Gilks and Roberts [1996] and Crainiceanu et al. [003], the covariates (x 1t, x t ) were centered (and standardized) to facilitate convergence of the Markov chain. In addition, it is well known that hierarchical implementation of random effects models is facilitated by so-called hierarchical centering, wherein one models the mean of one stage as the mean of the random effect at the next stage [e.g., Gilks and Roberts, 1996]. Crainiceanu et al. [003] demonstrate that this is also critical for facilitating convergence in hierarchical Poisson regression models with random effects. Thus the spatial mean D is parameterized in the x process rather than in the model for log(l). Similarly, for the logit p process, the mean was deferred to the following stage and thus was modeled in the n process. [38] Specification of the hyperparameter values given in Table 1 was subjective, but these values suggest vague prior distributions so that the data have more influence. Furthermore, sensitivity analyses were performed, and in all cases the posterior estimates were not overly sensitive to the specified hyperparameter values. [39] Samples for the b and h processes were obtained by simply taking the inverse Fourier transform of the a and x processes, respectively. Scientific interpretation is easier to conceptualize in physical space rather than in spectral space. [40] The MCMC simulation of the posterior distribution was run for 00,000 iterations after a 50,000 iteration burn in. Convergence was evaluated using basic univariate measures and different starting values (e.g., see Robert and Casella [1999] for an overview of approaches). Note that the slowest mixing of the Markov chain occurred with the s and s x parameters (although this was not problematic after the 50,000 iteration burn in). This is not particularly surprising as there are some potential identifiability problems with these variance components due to the lack of posterior spatial dependence in the h (i.e., x) process (see section 4.1). However, the sum of the two variance components is in fact quite well behaved. [41] Posterior standard errors for spatial processes were calculated by the batch mean approach [e.g., Roberts, 1996]. For scalar parameters the effective degrees of freedom approach was used to obtain posterior standard errors [e.g., Roberts, 1996]. 4. Results 4.1. Model Results [4] Posterior means and variances are shown, along with the priors for the model s scalar parameters, in Table 1. First, we note that the values of the spatial dependence parameters (q k ) for the b process (i.e., a process in the model) are quite small, which indicates substantial spatial dependence (recall that our distances are standardized so that the 50-km distance between adjacent grid boxes is a distance of 1). By contrast, the spatial dependence parameter (q x ) for the random effects process (h, modeled spectrally as x) is quite large. This value suggests that there is virtually no spatial dependence in this random effects process and that a simple independent noise model would probably suffice. Note, however, that there was enough dependence to identify the associated variance components, with posterior means s = 0.06 and s x = 0.3, yet the mixing was rather slow, as mentioned in section 3.. The variance components for the b processes were relatively small by comparison to the sum of these noise components (e.g., s 1 = 0.06 and s =.03). However, we did not expect that the temporal trend and ENSO components would be the dominant components of variability in the tornado counts but rather that they would be important components. On the basis of the variance of these distributions and the posterior distributions shown in Figure, these variance components do represent a significant portion of the tornado count variability. The significance of the b estimates is discussed further in section 4.. [43] Although the zero-inflated probability process is considered more as a nuisance process in this analysis, it is instructive to examine the associated parameters. First, we note from the comparison of the posterior means of s g and s n that the spatial process N is not particularly important. However, examination of the spatial map of the N process
9 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL STS 1-9 Figure 3. Posterior mean of the n spatial process. Note that values >0 imply zero-inflated probabilities >0.5 and that values <0 imply zero-inflated probabilities <0.5. shown in Figure 3 suggests that the probability that the count is Poisson tends to be >0.5 in regions of largest counts and <0.5 in the Appalachian corridor and edge regions (for which there are few nonzero observations). Note that a value of n(s i ) > 0 suggests that there is a >0.5 probability that the process follows a Poisson distribution and that a value of n(s i ) < 0 suggests that the probability is <0.5 that the process is Poisson. Thus the model is performing as expected in that it is increasing the probability of zeros in areas where observations are limited or the process cannot be nonzero. However, we note that the differences in probabilities implied by this n process are not very large (e.g., the spatial variation in p implied by N ranges between 0.49 and 0.51). More sophisticated models could be developed for the n process, although such models are not considered here. We do note that the trend and ENSO covariates were considered in the logit of the p process in a preliminary analysis. These results did not show a strong relationship and thus were not considered further. 4.. Results Related to Tornado Climatology [44] The parameters of most meteorological importance related to tornado climatology are B 1, B, and H t. Recall that these are obtained directly from samples of A 1, A, and X t through the inverse DFT, respectively. For simplicity of interpretation we focus our results on the spatial domain rather than on the spectral. It is interesting to note, however, that the majority of the spectral coefficients for the a Figure 4. Posterior mean of the spatial mean parameters D h.
10 STS 1-10 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL Figure 5. (top) Posterior mean and (bottom) posterior standard deviation for B 1 (linear temporal trend) parameters. processes are near zero (i.e., 80% of the posterior mean for a 1 coefficients are between 0.0 and 0.0; 88% of the posterior means for a are within this interval). This is not surprising as the Fourier transform does a good job of dimension reduction when the spatial correlation is relatively large, which is the case for a 1 and a, as discussed in section 4.1. In contrast, only 4% of the x posterior mean coefficients are between 0.0 and 0.0, reflecting the fact that there is very little spatial correlation in the h random component, as discussed in section 4.1. That is, the x process is not much different than a white noise process (flat spectrum). [45] In general, posterior means of D h (i.e., the inverse Fourier transform of D) shown in Figure 4 reflect the spatial pattern within the raw counts. It shows relative minima in Missouri, in the western plains states, and in the Appalachian corridor. In addition, there are relative maxima in northeast Texas, Oklahoma, eastern Kansas, Iowa, Indiana, Alabama, Mississippi, and, perhaps surprisingly, western Connecticut. The spatial pattern of this process exhibits much more detail than the results of Concannon et al. [000] that were produced from the so-called NWSsmoothed log. The similarities between raw counts and other smoothed plots lend credibility to the model results.
11 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL STS 1-11 Figure 6. (top) Posterior mean and (bottom) posterior standard deviation for the B (Niño 3.4 index) parameters. [46] The maps of posterior means and posterior standard deviations for the b parameters are given in Figures 5 6. If one considers the posterior distribution of the bs, the spatial regions corresponding to large magnitudes in the posterior mean in Figures 5 and 6 have 95% credible regions that are unlikely to cover zero, as shown in Figure 7. Thus one can say that these regions have less than a 5% probability of showing no response and thus are significant in a Bayesian (credible) sense. In addition, comparison with the total report counts over all years shown in Figure 1 shows that, as one would expect, the largest standard deviations correspond to areas of little or no data. This can be seen more clearly in Figure 8, which shows a scatterplot of the posterior standard deviations for b (s i ) versus the total raw tornado counts for each grid box summed over all years. [47] It has been consistently reported that national counts of F F5 tornadoes have either remained nearly constant or slightly decreased over the period analyzed herein [Marzban and Schaefer, 001]. The spatial field for posterior means of b 1 (the linear trend coefficients) shown in Figure 5 exhibits substantial spatial variability. Relatively large positive values of b 1 (indicating positive trends over time) are found near urban areas, such as along the eastern seaboard and
12 STS 1-1 WIKLE AND ANDERSON: CLIMATOLOGICAL ANALYSIS USING A BAYESIAN MODEL note that one has to be careful in interpreting the maps near the boundaries due to edge effects and limited data. This result confirms findings from Marzban and Schaefer [001] and Monfredo [1999] and establishes that the spatial pattern of this association is very coherent, yet regionally variable. [49] Smaller regions with small positive and negative values throughout the central U.S. are inconsistent with results from Monfredo [1999]. Monfredo [1999] has reported significant seasonal lag correlation for F F5 tornado reports that were accumulated over Missouri, Arkansas, eastern portions of Kansas and Oklahoma, and northern Texas. A similar inconsistency with results from Marzban and Schaefer [001] is found in the northeast U.S. Our approach differs from both studies in that, by using annual tornado counts, we have simultaneously examined effects of both seasonal lag and contemporaneous dependence. Further research is needed to determine how our results are influenced separately by lag and contemporaneous dependence at less than annual timescales in order to fully explain these contradictions. [50] Positive association in the western plains states is a unique result. Previous studies have not considered the relationship between equatorial Pacific SST and tornado reports in this region. However, positive seasonal lag correlation between equatorial Pacific SST and summertime precipitation has been identified in this region [Ropelewski and Halpert, 1996; Ting and Wang, 1997]. Since nearly all summertime precipitation in this region is supplied by thunderstorms, increased precipitation implies an increase in either thunderstorm frequency or precipitation efficiency of thunderstorms. If it is shown that thunderstorm frequency increases, it might be reasonable to find an increase of F F5 tornado frequency as well since most F F5 tornadoes are spawned from thunderstorms. Figure 7. Spatial regions corresponding to Figures 5 and 6, with 95% credible intervals that do not cover zero, indicated by the filled regions. (top) b 1 (linear temporal trend). (bottom) b (Niño 3.4 index). near Denver, Colorado, whereas a large region of negative values (indicative of negative temporal trends) are located in the central U.S. and in the Florida panhandle region, where the spatial mean is relatively large. The large-scale spatial structure in this map suggests a southwest to northeast orientation to the positive and negative estimates. These results do not necessarily contradict previous results. It is possible that a large increase of tornado reports in regions where tornadoes have been and still are infrequent may be offset by decreases over the central U.S., where tornado frequency is much higher. However, it is clear that a national trend should not be considered representative of local trend. [48] Posterior means of b (Figure 6) indicate that dependence does exist between equatorial Pacific SST anomalies (Niño 3.4) and tornado report frequency in the U.S. An expansive area of negative values in the southeast U.S. suggests a reduced (increased) risk of F F5 tornadoes during El Niño (La Niña) in this region, while the opposite relationship is suggested for the western plains states. We 5. Discussion [51] The statistical model applied herein to tornado report counts is an exploratory tool. Confidence in the validity of our results leads us to pursue further our analysis in the form of physical diagnoses, with the goal of determining causal links between equatorial SST and tornado reports. Much has been learned in the past 5 years about mesoscale weather conditions that support tornadic thunderstorms. The question we are led to ask is, How might slowly varying climate forces create situations in which such mesoscale conditions repeatedly form or fail to form? This will require analyses of mean and transient flow during El Niño and La Niña events. Already, it is known that the wintertime jet stream over the northeastern U.S. is sometimes stronger during La Niña than during El Niño [Hoerling et al., 1997], which, as noted by Marzban and Schaefer [001], could provide favorable conditions for tornadic thunderstorms should the anomaly persist into spring and summer months. It is important to follow up this finding with more detailed analysis. That is, it is necessary to determine what role, if any, transient synoptic-scale systems play in producing environments supportive of tornadic thunderstorms. In addition, the immediate environment of tornadic thunderstorms should be examined to determine if the structure of vertical wind shear and thermal stability resemble that of tornadoes elsewhere in the U.S.
Weather and Climate Summary and Forecast March 2018 Report
Weather and Climate Summary and Forecast March 2018 Report Gregory V. Jones Linfield College March 7, 2018 Summary: The ridge pattern that brought drier and warmer conditions from December through most
More informationNational Wildland Significant Fire Potential Outlook
National Wildland Significant Fire Potential Outlook National Interagency Fire Center Predictive Services Issued: September, 2007 Wildland Fire Outlook September through December 2007 Significant fire
More informationJP4.17 TORNADO AND SEVERE WEATHER CLIMATOLOGY AND PREDICTABILITY BY ENSO PHASE IN THE NORTH CENTRAL U.S.: A COMPOSITING STUDY
JP4.17 TORNADO AND SEVERE WEATHER CLIMATOLOGY AND PREDICTABILITY BY ENSO PHASE IN THE NORTH CENTRAL U.S.: A COMPOSITING STUDY Barbara E. Mayes *1 NOAA/NWS WFO Quad Cities IA/IL, Davenport, Iowa Craig Cogil
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
More informationSOUTHERN CLIMATE MONITOR
SOUTHERN CLIMATE MONITOR MARCH 2011 VOLUME 1, ISSUE 3 IN THIS ISSUE: Page 2 to 4 Severe Thunderstorm Climatology in the SCIPP Region Page 4 Drought Update Page 5 Southern U.S. Precipitation Summary for
More informationForced and internal variability of tropical cyclone track density in the western North Pacific
Forced and internal variability of tropical cyclone track density in the western North Pacific Wei Mei 1 Shang-Ping Xie 1, Ming Zhao 2 & Yuqing Wang 3 Climate Variability and Change and Paleoclimate Working
More informationImperfect Data in an Uncertain World
Imperfect Data in an Uncertain World James B. Elsner Department of Geography, Florida State University Tallahassee, Florida Corresponding author address: Dept. of Geography, Florida State University Tallahassee,
More informationALLEY CATS: HOW EL NIÑO INFLUENCES TORNADO ALLEY AND THE THREAT OF CATASTROPHES
ALLEY CATS: HOW EL NIÑO INFLUENCES TORNADO ALLEY AND THE THREAT OF CATASTROPHES Kevin Van Leer - Sr. Product Manager, Model Product Management National Tornado Summit - Tuesday, March 1 st, 2016 1 EL NIÑO/SOUTHERN
More informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
More informationPRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Event Response
PRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK June 2014 - RMS Event Response 2014 SEASON OUTLOOK The 2013 North Atlantic hurricane season saw the fewest hurricanes in the Atlantic Basin
More information2013 ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Cat Response
2013 ATLANTIC HURRICANE SEASON OUTLOOK June 2013 - RMS Cat Response Season Outlook At the start of the 2013 Atlantic hurricane season, which officially runs from June 1 to November 30, seasonal forecasts
More informationStatistics Introduction to Probability
Statistics 110 - Introduction to Probability Mark E. Irwin Department of Statistics Harvard University Summer Term Monday, June 26, 2006 - Wednesday, August 16, 2006 Copyright 2006 by Mark E. Irwin Personnel
More informationEl Niño Update Impacts on Florida
Current Issues in Emergency Management (CIEM) Sessions 1 &2 October 12 th and 26 th, 2006 Florida Division of Emergency Management Tallahassee, Florida El Niño Update Impacts on Florida Bart Hagemeyer
More informationThe New Normal or Was It?
The New Normal or Was It? by Chuck Coffey The recent drought has caused many to reflect upon the past and wonder what is in store for the future. Just a couple of years ago, few agricultural producers
More informationSpatio-temporal precipitation modeling based on time-varying regressions
Spatio-temporal precipitation modeling based on time-varying regressions Oleg Makhnin Department of Mathematics New Mexico Tech Socorro, NM 87801 January 19, 2007 1 Abstract: A time-varying regression
More informationThe role of teleconnections in extreme (high and low) precipitation events: The case of the Mediterranean region
European Geosciences Union General Assembly 2013 Vienna, Austria, 7 12 April 2013 Session HS7.5/NP8.4: Hydroclimatic Stochastics The role of teleconnections in extreme (high and low) events: The case of
More informationUnited States Streamflow Probabilities based on Forecasted La Niña, Winter-Spring 2000
United States Streamflow Probabilities based on Forecasted La Niña, Winter-Spring 2000 contributed by Michael D. Dettinger 1, Daniel R. Cayan 1, and Kelly T. Redmond 2 1 U.S. Geological Survey, Scripps
More informationWeather and Climate Summary and Forecast February 2018 Report
Weather and Climate Summary and Forecast February 2018 Report Gregory V. Jones Linfield College February 5, 2018 Summary: For the majority of the month of January the persistent ridge of high pressure
More informationBayesian data analysis in practice: Three simple examples
Bayesian data analysis in practice: Three simple examples Martin P. Tingley Introduction These notes cover three examples I presented at Climatea on 5 October 0. Matlab code is available by request to
More informationDiagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS)
Diagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS) Christopher L. Castro and Roger A. Pielke, Sr. Department of
More informationWeather and Climate Summary and Forecast Winter
Weather and Climate Summary and Forecast Winter 2016-17 Gregory V. Jones Southern Oregon University December 5, 2016 Well after an October that felt more like November, we just finished a November that
More informationWeather and Climate Summary and Forecast Winter
Weather and Climate Summary and Forecast Winter 2016-17 Gregory V. Jones Southern Oregon University February 7, 2017 What a difference from last year at this time. Temperatures in January and February
More informationSummer 2018 Southern Company Temperature/Precipitation Forecast
Scott A. Yuknis High impact weather forecasts, climate assessment and prediction. 14 Boatwright s Loop Plymouth, MA 02360 Phone/Fax 508.927.4610 Cell: 508.813.3499 ClimateImpact@comcast.net Climate Impact
More informationClimate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014
Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014 Short Term Drought Map: Short-term (
More informationWeather and Climate Summary and Forecast Fall/Winter 2016
Weather and Climate Summary and Forecast Fall/Winter 2016 Gregory V. Jones Southern Oregon University November 5, 2016 After a year where we were seemingly off by a month in terms of temperatures (March
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season Forecast : June 11, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met Service)
More informationGlobal Temperature Is Continuing to Rise: A Primer on Climate Baseline Instability. G. Bothun and S. Ostrander Dept of Physics, University of Oregon
Global Temperature Is Continuing to Rise: A Primer on Climate Baseline Instability G. Bothun and S. Ostrander Dept of Physics, University of Oregon The issue of whether or not humans are inducing significant
More informationWeather and Climate Summary and Forecast April 2018 Report
Weather and Climate Summary and Forecast April 2018 Report Gregory V. Jones Linfield College April 4, 2018 Summary: A near Miracle March played out bringing cooler and wetter conditions to the majority
More information1. TEMPORAL CHANGES IN HEAVY RAINFALL FREQUENCIES IN ILLINOIS
to users of heavy rainstorm climatology in the design and operation of water control structures. A summary and conclusions pertaining to various phases of the present study are included in Section 8. Point
More informationFigure 1. Time series of Western Sahel precipitation index and Accumulated Cyclone Energy (ACE).
2B.6 THE NON-STATIONARY CORRELATION BETWEEN SAHEL PRECIPITATION INDICES AND ATLANTIC HURRICANE ACTIVITY Andreas H. Fink 1 and Jon M. Schrage 2 1 Institute for Geophysics 2 Department of and Meteorology
More informationExtreme Weather Events and Climate Change
Extreme Weather Events and Climate Change Robert M Rabin NOAA/National Severe Storms Lab Norman, Oklahoma Most material presented is from: Climate Change 2013: The Physical Science Basis. Intergovernmental
More informationTHE SOUTHERN OSCILLATION AND ITS EFFECT ON TORNADIC ACTIVITY IN THE UNITED STATES. Supported by NSF Grant No. ATM R.A. Pielke Sr.
THE SOUTHERN OSCILLATION AND ITS EFFECT ON TORNADIC ACTIVITY IN THE UNITED STATES Supported by NSF Grant No. ATM-9017849 R.A. Pielke Sr., PI THE SOUTHERN OSCILLATION AND ITS EFFECT ON TORNADIC ACTIVITY
More informationImpacts of the April 2013 Mean trough over central North America
Impacts of the April 2013 Mean trough over central North America By Richard H. Grumm National Weather Service State College, PA Abstract: The mean 500 hpa flow over North America featured a trough over
More informationClimate Forecast Applications Network (CFAN)
Forecast of 2018 Atlantic Hurricane Activity April 5, 2018 Summary CFAN s inaugural April seasonal forecast for Atlantic tropical cyclone activity is based on systematic interactions among ENSO, stratospheric
More informationWeather and Climate Summary and Forecast March 2019 Report
Weather and Climate Summary and Forecast March 2019 Report Gregory V. Jones Linfield College March 2, 2019 Summary: Dramatic flip from a mild winter to a top five coldest February on record in many locations
More informationTornado Hazard Risk Analysis: A Report for Rutherford County Emergency Management Agency
Tornado Hazard Risk Analysis: A Report for Rutherford County Emergency Management Agency by Middle Tennessee State University Faculty Lisa Bloomer, Curtis Church, James Henry, Ahmad Khansari, Tom Nolan,
More informationSemiblind Source Separation of Climate Data Detects El Niño as the Component with the Highest Interannual Variability
Semiblind Source Separation of Climate Data Detects El Niño as the Component with the Highest Interannual Variability Alexander Ilin Neural Networks Research Centre Helsinki University of Technology P.O.
More information2011 National Seasonal Assessment Workshop for the Eastern, Southern, & Southwest Geographic Areas
2011 National Seasonal Assessment Workshop for the Eastern, Southern, & Southwest Geographic Areas On January 11-13, 2011, wildland fire, weather, and climate met virtually for the ninth annual National
More informationEl Niño / Southern Oscillation
El Niño / Southern Oscillation Student Packet 2 Use contents of this packet as you feel appropriate. You are free to copy and use any of the material in this lesson plan. Packet Contents Introduction on
More informationWeather and Climate Summary and Forecast January 2018 Report
Weather and Climate Summary and Forecast January 2018 Report Gregory V. Jones Linfield College January 5, 2018 Summary: A persistent ridge of high pressure over the west in December produced strong inversions
More informationCorrection to Spatial and temporal distributions of U.S. winds and wind power at 80 m derived from measurements
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2004jd005099, 2004 Correction to Spatial and temporal distributions of U.S. winds and wind power at 80 m derived from measurements Cristina L. Archer
More informationAN ANALYSIS OF THE TORNADO COOL SEASON
AN ANALYSIS OF THE 27-28 TORNADO COOL SEASON Madison Burnett National Weather Center Research Experience for Undergraduates Norman, OK University of Missouri Columbia, MO Greg Carbin Storm Prediction Center
More informationClimate Variability and El Niño
Climate Variability and El Niño David F. Zierden Florida State Climatologist Center for Ocean Atmospheric Prediction Studies The Florida State University UF IFAS Extenstion IST January 17, 2017 The El
More informationAn ENSO-Neutral Winter
An ENSO-Neutral Winter This issue of the Blue Water Outlook newsletter is devoted towards my thoughts on the long range outlook for winter. You will see that I take a comprehensive approach to this outlook
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season Forecast Issued: May 10, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met
More informationThe Correlation Between Fall and Winter Temperature Anomalies in the Midwest during ENSO Events
The Correlation Between Fall and Winter Temperature Anomalies in the Midwest during ENSO Events Morgan E. Brown Iowa State University Mentors: Ray Arritt and Alan Czarnetzki ABSTRACT According to the National
More informationHighlight: Support for a dry climate increasing.
Scott A. Yuknis High impact weather forecasts, climate assessment and prediction. 14 Boatwright s Loop Plymouth, MA 02360 Phone/Fax 508.927.4610 Cell: 508.813.3499 ClimateImpact@comcast.net Climate Impact
More informationMonitoring and Prediction of Climate Extremes
Monitoring and Prediction of Climate Extremes Stephen Baxter Meteorologist, Climate Prediction Center NOAA/NWS/NCEP Deicing and Stormwater Management Conference ACI-NA/A4A Arlington, VA May 19, 2017 What
More informationWeather and Climate Summary and Forecast August 2018 Report
Weather and Climate Summary and Forecast August 2018 Report Gregory V. Jones Linfield College August 5, 2018 Summary: July 2018 will likely go down as one of the top five warmest July s on record for many
More informationModelling trends in the ocean wave climate for dimensioning of ships
Modelling trends in the ocean wave climate for dimensioning of ships STK1100 lecture, University of Oslo Erik Vanem Motivation and background 2 Ocean waves and maritime safety Ships and other marine structures
More informationCHAPTER 11 THUNDERSTORMS AND TORNADOES MULTIPLE CHOICE QUESTIONS
CHAPTER 11 THUNDERSTORMS AND TORNADOES MULTIPLE CHOICE QUESTIONS 1. A thunderstorm is considered to be a weather system. a. synoptic-scale b. micro-scale c. meso-scale 2. By convention, the mature stage
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 informationStatistics Introduction to Probability
Statistics 110 - Introduction to Probability Mark E. Irwin Department of Statistics Harvard University Summer Term Monday, June 28, 2004 - Wednesday, August 18, 2004 Personnel Instructor: Mark Irwin Office:
More informationWeather and Climate Summary and Forecast Summer 2017
Weather and Climate Summary and Forecast Summer 2017 Gregory V. Jones Southern Oregon University August 4, 2017 July largely held true to forecast, although it ended with the start of one of the most extreme
More informationNOTES AND CORRESPONDENCE. El Niño Southern Oscillation and North Atlantic Oscillation Control of Climate in Puerto Rico
2713 NOTES AND CORRESPONDENCE El Niño Southern Oscillation and North Atlantic Oscillation Control of Climate in Puerto Rico BJÖRN A. MALMGREN Department of Earth Sciences, University of Göteborg, Goteborg,
More informationP3.6 THE INFLUENCE OF PNA AND NAO PATTERNS ON TEMPERATURE ANOMALIES IN THE MIDWEST DURING FOUR RECENT El NINO EVENTS: A STATISTICAL STUDY
P3.6 THE INFLUENCE OF PNA AND NAO PATTERNS ON TEMPERATURE ANOMALIES IN THE MIDWEST DURING FOUR RECENT El NINO EVENTS: A STATISTICAL STUDY Dayton Vincent 2, Sam Lashley 1, Sam O Connor 2, Michael Skipper
More informationREGIONAL VARIABILITY OF CAPE AND DEEP SHEAR FROM THE NCEP/NCAR REANALYSIS ABSTRACT
REGIONAL VARIABILITY OF CAPE AND DEEP SHEAR FROM THE NCEP/NCAR REANALYSIS VITTORIO A. GENSINI National Weather Center REU Program, Norman, Oklahoma Northern Illinois University, DeKalb, Illinois ABSTRACT
More informationInterannual Variability of the South Atlantic High and rainfall in Southeastern South America during summer months
Interannual Variability of the South Atlantic High and rainfall in Southeastern South America during summer months Inés Camilloni 1, 2, Moira Doyle 1 and Vicente Barros 1, 3 1 Dto. Ciencias de la Atmósfera
More informationUsing modern time series analysis techniques to predict ENSO events from the SOI time series
Nonlinear Processes in Geophysics () 9: 4 45 Nonlinear Processes in Geophysics c European Geophysical Society Using modern time series analysis techniques to predict ENSO events from the SOI time series
More informationIntroduction of climate monitoring and analysis products for one-month forecast
Introduction of climate monitoring and analysis products for one-month forecast TCC Training Seminar on One-month Forecast on 13 November 2018 10:30 11:00 1 Typical flow of making one-month forecast Observed
More informationThe U. S. Winter Outlook
The 2018-2019 U. S. Winter Outlook Michael Halpert Deputy Director Climate Prediction Center Mike.Halpert@noaa.gov http://www.cpc.ncep.noaa.gov Outline About the Seasonal Outlook Review of 2017-18 U. S.
More informationThai Meteorological Department, Ministry of Digital Economy and Society
Thai Meteorological Department, Ministry of Digital Economy and Society Three-month Climate Outlook For November 2017 January 2018 Issued on 31 October 2017 -----------------------------------------------------------------------------------------------------------------------------
More informationSeasonal Climate Watch July to November 2018
Seasonal Climate Watch July to November 2018 Date issued: Jun 25, 2018 1. Overview The El Niño-Southern Oscillation (ENSO) is now in a neutral phase and is expected to rise towards an El Niño phase through
More informationCHAPTER 1: INTRODUCTION
CHAPTER 1: INTRODUCTION There is now unequivocal evidence from direct observations of a warming of the climate system (IPCC, 2007). Despite remaining uncertainties, it is now clear that the upward trend
More informationNorth Pacific Climate Overview N. Bond (UW/JISAO), J. Overland (NOAA/PMEL) Contact: Last updated: September 2008
North Pacific Climate Overview N. Bond (UW/JISAO), J. Overland (NOAA/PMEL) Contact: Nicholas.Bond@noaa.gov Last updated: September 2008 Summary. The North Pacific atmosphere-ocean system from fall 2007
More informationDid we see the 2011 summer heat wave coming?
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl051383, 2012 Did we see the 2011 summer heat wave coming? Lifeng Luo 1 and Yan Zhang 2 Received 16 February 2012; revised 15 March 2012; accepted
More informationWeather and Climate Summary and Forecast January 2019 Report
Weather and Climate Summary and Forecast January 2019 Report Gregory V. Jones Linfield College January 4, 2019 Summary: December was mild and dry over much of the west, while the east was much warmer than
More informationModels for models. Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research
Models for models Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Outline Statistical models and tools Spatial fields (Wavelets) Climate regimes (Regression and clustering)
More informationEl Niño Seasonal Weather Impacts from the OLR Event Perspective
Science and Technology Infusion Climate Bulletin NOAA s National Weather Service 41 st NOAA Annual Climate Diagnostics and Prediction Workshop Orono, ME, 3-6 October 2016 2015-16 El Niño Seasonal Weather
More informationRainfall variability over the Indochina peninsula during the Boreal Winter, Part I: Preliminary data analysis
Rainfall variability over the Indochina peninsula during the Boreal Winter, Part I: Preliminary data analysis Sirapong Sooktawee*, sirapong@deqp.go.th; Atsamon Limsakul, atsamon@deqp.go.th, Environmental
More informationAnuMS 2018 Atlantic Hurricane Season Forecast
AnuMS 2018 Atlantic Hurricane Season : August 12, 2018 by Dale C. S. Destin (follow @anumetservice) Director (Ag), Antigua and Barbuda Meteorological Service (ABMS) The *AnuMS (Antigua Met Service) is
More information2015: A YEAR IN REVIEW F.S. ANSLOW
2015: A YEAR IN REVIEW F.S. ANSLOW 1 INTRODUCTION Recently, three of the major centres for global climate monitoring determined with high confidence that 2015 was the warmest year on record, globally.
More informationWeather and Climate Summary and Forecast Summer 2016
Weather and Climate Summary and Forecast Summer 2016 Gregory V. Jones Southern Oregon University June 6, 2016 May 2016 continued the warm trend for portions of the west, while providing some relief for
More informationINVISIBLE WATER COSTS
Every Drop Every Counts... Drop Counts... INVISIBLE WATER COSTS Corn - 108.1 gallons per pound How much water it takes to produce... Apple - 18.5 gallons to grow Beef - 1,581 gallons per pound Oats - 122.7
More informationEnhancing Weather Information with Probability Forecasts. An Information Statement of the American Meteorological Society
Enhancing Weather Information with Probability Forecasts An Information Statement of the American Meteorological Society (Adopted by AMS Council on 12 May 2008) Bull. Amer. Meteor. Soc., 89 Summary This
More informationlecture 11 El Niño/Southern Oscillation (ENSO) Part II
lecture 11 El Niño/Southern Oscillation (ENSO) Part II SYSTEM MEMORY: OCEANIC WAVE PROPAGATION ASYMMETRY BETWEEN THE ATMOSPHERE AND OCEAN The atmosphere and ocean are not symmetrical in their responses
More information2015 Summer Forecast
2015 Summer Forecast Inside This Forecast Summer 2014 Recap 2 Implications of El Niño and the PDO on the Upcoming Summer Summary of Global Indices & Oscillations Analog Year: The Midwest a battleground
More informationClimate informed flood frequency analysis and prediction in Montana using hierarchical Bayesian modeling
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L05404, doi:10.1029/2007gl032220, 2008 Climate informed flood frequency analysis and prediction in Montana using hierarchical Bayesian
More informationCLIMATE EXTREMES AND GLOBAL WARMING: A STATISTICIAN S PERSPECTIVE
CLIMATE EXTREMES AND GLOBAL WARMING: A STATISTICIAN S PERSPECTIVE Richard L. Smith Department of Statistics and Operations Research University of North Carolina, Chapel Hill rls@email.unc.edu Statistics
More informationExtreme Rainfall in the Southeast U.S.
Extreme Rainfall in the Southeast U.S. David F. Zierden Florida State Climatologist Center for Ocean Atmospheric Prediction Studies The Florida State University March 7, 2016 Causes of Extreme Rainfall
More informationAnalysis Links Pacific Decadal Variability to Drought and Streamflow in United States
Page 1 of 8 Vol. 80, No. 51, December 21, 1999 Analysis Links Pacific Decadal Variability to Drought and Streamflow in United States Sumant Nigam, Mathew Barlow, and Ernesto H. Berbery For more information,
More informationPaul T. Schickedanz, Prameela V. Reddy, and Stanley A. Changnon
Illinois State Water Survey at the University of Illinois Urbana, Illinois SPATIAL AND TEMPORAL RELATIONSHIPS IN CROP-HAIL LOSS DATA by Paul T. Schickedanz, Prameela V. Reddy, and Stanley A. Changnon Final
More informationProceedings, International Snow Science Workshop, Banff, 2014
RELATING AVALANCHES TO LARGE-SCALE OCEAN ATMOSPHERIC OSCILLATIONS Scott Thumlert 1, Sascha Bellaire 2, Bruce Jamieson 1 1 Dept. of Civil Engineering, University of Calgary, AB, Canada 2 Institute of Meteorology
More informationSeasonal Climate Watch February to June 2018
Seasonal Climate Watch February to June 2018 Date issued: Jan 26, 2018 1. Overview The El Niño-Southern Oscillation (ENSO) is expected to remain in a weak La Niña phase through to early autumn (Feb-Mar-Apr).
More informationWilliam M. Frank* and George S. Young The Pennsylvania State University, University Park, PA. 2. Data
1C.1 THE 80 CYCLONES MYTH William M. Frank* and George S. Young The Pennsylvania State University, University Park, PA 1. Introduction myth: A traditional story accepted as history and/or fact. For the
More informationENSO Cycle: Recent Evolution, Current Status and Predictions. Update prepared by Climate Prediction Center / NCEP July 26, 2004
ENSO Cycle: Recent Evolution, Current Status and Predictions Update prepared by Climate Prediction Center / NCEP July 26, 2004 Outline Overview Recent Evolution and Current Conditions Oceanic NiZo Index
More informationThe North Atlantic Oscillation: Climatic Significance and Environmental Impact
1 The North Atlantic Oscillation: Climatic Significance and Environmental Impact James W. Hurrell National Center for Atmospheric Research Climate and Global Dynamics Division, Climate Analysis Section
More informationWhat Is the Weather Like in Different Regions of the United States?
Learning Set 1 What Is Weather, and How Is It Measured and Described? 1.3 Explore What Is the Weather Like in Different Regions of the United States? trends: patterns or tendencies you can see over a broad
More informationColorado CoCoRaHS. Colorado CoCoRaHS. Because Every Drop Counts! October 2014 Volume 2, Issue 10
U.S. PRECIPITATION (% OF AVERAGE) LOOKING BACK AT SEPTEMBER 2014 September was dry for much of the nation with a few exceptions; one of those being the central Rockies. In Colorado, it was wetter-than-average
More informationNational Wildland Significant Fire Potential Outlook
National Wildland Significant Fire Potential Outlook National Interagency Fire Center Predictive Services Issued: April 1, 2008 Next Issue: May 1, 2008 Wildland Fire Outlook April 2008 through July 2008
More informationEarly May Cut-off low and Mid-Atlantic rains
Abstract: Early May Cut-off low and Mid-Atlantic rains By Richard H. Grumm National Weather Service State College, PA A deep 500 hpa cutoff developed in the southern Plains on 3 May 2013. It produced a
More informationSEASONAL ENVIRONMENTAL CONDITIONS RELATED TO HURRICANE ACTIVITY IN THE NORTHEAST PACIFIC BASIN
SEASONAL ENVIRONMENTAL CONDITIONS RELATED TO HURRICANE ACTIVITY IN THE NORTHEAST PACIFIC BASIN Jennifer M. Collins Department of Geography and Geosciences Bloomsburg University Bloomsburg, PA 17815 jcollins@bloomu.edu
More informationExtremes Seminar: Tornadoes
Dec. 01, 2014 Outline Introduction 1 Introduction 2 3 4 Introduction 101: What is a tornado? According to the Glossary of Meteorology (AMS 2000), a tornado is a violently rotating column of air, pendant
More informationNorth Pacific Climate Overview N. Bond (UW/JISAO), J. Overland (NOAA/PMEL) Contact: Last updated: August 2009
North Pacific Climate Overview N. Bond (UW/JISAO), J. Overland (NOAA/PMEL) Contact: Nicholas.Bond@noaa.gov Last updated: August 2009 Summary. The North Pacific atmosphere-ocean system from fall 2008 through
More information37. Regional Climate Change in the Southern United States: The Implications for Wildfire Occurrence
1998 Springer-Verlag New York, Inc. The Productivity & Sustainability of Southern Forest Ecosystems in a Changing Environment by Mickler & Fox 37. Regional Climate Change in the Southern United States:
More informationSHORT COMMUNICATION EXPLORING THE RELATIONSHIP BETWEEN THE NORTH ATLANTIC OSCILLATION AND RAINFALL PATTERNS IN BARBADOS
INTERNATIONAL JOURNAL OF CLIMATOLOGY Int. J. Climatol. 6: 89 87 (6) Published online in Wiley InterScience (www.interscience.wiley.com). DOI:./joc. SHORT COMMUNICATION EXPLORING THE RELATIONSHIP BETWEEN
More informationP3.1 Development of MOS Thunderstorm and Severe Thunderstorm Forecast Equations with Multiple Data Sources
P3.1 Development of MOS Thunderstorm and Severe Thunderstorm Forecast Equations with Multiple Data Sources Kathryn K. Hughes * Meteorological Development Laboratory Office of Science and Technology National
More informationLecture 14 Bayesian Models for Spatio-Temporal Data
Lecture 14 Bayesian Models for Spatio-Temporal Data Dennis Sun Stats 253 August 13, 2014 Outline of Lecture 1 Recap of Bayesian Models 2 Empirical Bayes 3 Case 1: Long-Lead Forecasting of Sea Surface Temperatures
More informationWeather and Climate Summary and Forecast Summer 2016
Weather and Climate Summary and Forecast Summer 2016 Gregory V. Jones Southern Oregon University August 5, 2016 The month of July saw the continued flip, flop of the western US temperatures. It was warm
More informationA COMPREHENSIVE 5-YEAR SEVERE STORM ENVIRONMENT CLIMATOLOGY FOR THE CONTINENTAL UNITED STATES 3. RESULTS
16A.4 A COMPREHENSIVE 5-YEAR SEVERE STORM ENVIRONMENT CLIMATOLOGY FOR THE CONTINENTAL UNITED STATES Russell S. Schneider 1 and Andrew R. Dean 1,2 1 DOC/NOAA/NWS/NCEP Storm Prediction Center 2 OU-NOAA Cooperative
More information