Hurricanes are destructive natural phenomena. On average a half of dozen or so form each year over the warm tropical and subtropical waters of the North Atlantic and track westward for thousands of kilometers. Historically, hurricanes account for a majority of the costliest weather disasters in the United States. They rival earthquakes in destructive potential and loss of life. Despite technological advances in monitoring and prediction, hurricanes retain their potential to cause severe damage and numerous deaths (Arguez and Elsner 2001). During 1998, hurricane Mitch became a grim reminder that hurricanes can quickly kill thousands of people. In the United States, population and demographic shifts toward the coast are making the problem worse as development flourishes in areas prone to hurricane strikes: the warm subtropical shorelines and islands of the Atlantic Ocean and Gulf of Mexico.
Knowledge of their past occurrence, even if it is incomplete, provides clues about their future frequency and intensity that goes beyond what present numerical climate models are capable of. This understanding is important for land-use planning, emergency management, hazard mitigation, (re)insurance applications and, potentially, the long-term weather derivative market. Climatologists have been issuing seasonal hurricane activity forecasts for the North Atlantic using since 1984 (Gray, 1984; Elsner et al., 1996). These forecasts provide total basin estimates of annual counts using linear and generalized linear regression models. Forecasts are issued several months in advance of the season, which runs from June through November, and updated as the season approaches.
Some understanding has been achieved in solving the climate puzzle with regard to the question of where hurricanes are likely to go based on conditions a month or two in advance of the season. Lehmiller et al. (1997) outline the problem and show the potential for specific forecast models based on regional and large-scale climate factors. Elsner et al. (2000) demonstrate a link between the probability of a major hurricane strike on the East Coast and the strength of the North Atlantic oscillation (NAO). Physically it is reasoned that the strength of the NAO is an indication of the position and intensity of the subtropical high, which steers the hurricanes toward the coast. Taking advantage of a unique wind speed data set from hurricane landfalls in the United States, Jagger et al. (2001) develop a model for hurricane probabilities conditioned on climate anomalies, including El Niño and the NAO.
But current seasonal forecast models for the entire North Atlantic basin do not incorporate spatial and temporal information. Thus they fail to provide specific seasonal activity forecasts for different geographic regions of the hurricane basin, which includes the Gulf of Mexico and Caribbean Sea. Statistical models need to be created that combine spatial and temporal correlation in the data to generate regional forecasts. This paper introduces a class of space-time statistical models for count data that can be used for seasonal hurricane prediction. Ultimately the goal is a model that can predict the likely tracks of hurricanes for an entire season.
The paper is organized as follows. Section two describes the model grid and data. Model formulae and justification are presented in section three along with a description of the model predictors. Section four describes increasingly sophisticated models as a way to understand the full model. Final model selection procedures are given in section five. Some issues related to the estimation procedure and a single hindcast are given in section six. The paper ends with a summary in section seven.
To develop the model we divide the North Atlantic basin into a 6° by 6° latitude and longitude grid. Grid choice is a compromise between sample size and resolution. As noted later, model coefficients are estimated using a Monte Carlo procedure, and the procedure fails to converge for smaller grid boxes because between grid box correlations are too large.
Data are annual hurricane counts in the grid boxes. A tropical cyclone that records a position at hurricane intensity within the box is counted once. Hurricanes are tropical cyclones with maximum sustained winds reaching 65 kt or greater. A hurricane that loops around and reenters the box is counted as a single hurricane. Hurricane positions and intensities are obtained from the best-track records (Neumann et al., 1999), which are a compilation of the six-hourly information of all tropical cyclones back to 1886. Records are most reliable after 1943, particularly for the weaker, shorter lived cyclones. The period of record for the present study is 1900-93 (94 yr).
We remove grids having mostly land or historically low hurricane activity, leaving grid S with 40 cells as shown in Fig. . Total hurricane occurrences over the 94-yr period are shown for each region. Historically, hurricane activity is most pronounced over the Bahamas extending north and eastward toward Bermuda (see Elsner and Kara 1999). Count values in the grid boxes outside region S are used as boundary conditions for the model. This improves the estimates of model constants over the case where the boundary values must be assumed or estimated. For generating sample forecasts we set the boundary values equal to zero to approximate annual climatology.
Count values in the grid boxes obtained from the best-track data and used in model development are displayed as a series of space-time images in Fig. . Each image represents an 11 by 8 grid (6° by 6° latitude/longitude), covering the western half of the North Atlantic for a single season. Gray levels in the images indicate the number of hurricanes whose centers pass into the region during the year. Activity in the grids range from zero (light gray) to four (black). Low latitude hurricane activity dominated the region during 1916 and 1933 whereas high latitude activity was more pronounced during 1963 and 1969. Note that there are no hurricanes in the best-track data for the years 1907 and 1914. These annual maps of seasonal hurricane counts are used as the model response as detailed next.
Hurricane frequencies over time and space form a space-time counts
process. The dependence structure of this type of data can be
modeled by a conditional probability approach (Bartlett, 1968;
Whittle, 1963; Besag, 1974; Gilks et al., 1996). Besag (1974)
introduce conditionally specified auto-Poisson models for spatial
counts data, which link observation of a Poisson process at a
given location with those at its spatial neighborhoods. But the
auto-Poisson model proposed by Besag (1974) has restrictions on
the parameter space making it applicable only to spatial data in
which the interaction coefficients are non-positive. Here we
consider a class of space-time regression models for hurricane
activity based on the right-truncated Poisson distribution
|
We restrict the specification to a spatially invariant, nearest neighbor, first order autoregressive model, with predictors. The space-time neighborhood consists of five sites: north, south, east, west, and the previous season's activity at the given site as shown in Fig. . Note that using the TPSTAR model on this neighborhood requires three coupling parameters g0h,g0v,g1c for the east-west, the north-south and the temporal lag-one coupling, respectively.
Each grid box (t,ij) represents a 6° by 6° region
for one year. Let the response values Ht,ij, be the number of
distinct hurricanes passing into any portion of the grid (i,j) Î S during year t. The longitude of the location centers are
given by -107° + 6°·i, for i=1¼11 and the
latitude centers by 6°+6°·j, for j=1¼8.
We can describe the distribution model conditionally as:
|
|
We chose M=10 for the grid boxes, since it is larger than the maximum number of hurricanes observed in any grid box per year. The value of M is small enough to observe positive coupling, yet large enough to compare the results of this model with models based on the Poisson distribution.
The model makes use of three predictors types:
Predictors of type two represent the regional effects from predictors that are not measured at each location, such as sea surface temperature, location to land mass, and the variation in surface area from region to region. The parameters are given by aij. Without additional information the site offset gives an indication of average yearly activity in the associated grid box. As aij represents one parameter for every grid box, this can be a considerable number of parameters. Thus, we constrain aij to be the sum of two factors, longitude ai and latitude bj, where aij=ai+bj. This reduces the number of parameters associated with offsets to 14 in the present case.
Predictors of type three represent the climate covariates as identified by previous research studies (Jagger et al., 2001). Here we consider the influence of the covariate to be the same across grid boxes, and allow the value of the covariate to vary in intensity from year to year. The number of covariates is restricted by data availability. The climate covariates used in the TPSTAR model include:
To better examine hurricane activity as a count process over the basin, we present results from three different models, a time-series only model, a space-time model without instantaneous coupling parameters, and a full model with coupling parameters. The full model given in (3) represents a version of TPSTAR. The first two models are fit to a Poisson generalized linear model (GLM) with dispersion. To fit the full model we try both a maximum Poisson likelihood estimator (MPLE) and a Monte Carlo maximum likelihood estimator (MCMLE).
The MPLE extends the GLM model by adding neighborhood observations as covariates to the models. If the covariate matrix is of full rank, the MPLE using the canonical link function always produces a parameter estimate. If the model is shift invariant with finite range, the parameter estimates of the MPLE are asymptotically consistent with increasing domains (Winkler 1995). However, even in the shift invariant case, the MPLE is not necessarily efficient, and does not provide standard error estimates. We use the MCMLE method for consistent and efficient parameter estimates with consistent parameter covariance estimates as described in Geyer (1994) and used by Wu (1994). The MCMLE method was extended by Jagger (2000) to handle the autoregressive coupling in the TPSTAR model.
The MCMLE method has several problems. For one, it is computationally intensive as it uses a Markov Chain Monte Carlo method for estimating the log likelihood function. For another, the method fails to converge, unless the initial parameter estimates are close to the actual parameter values. Though the MPLE is biased in our case, the convergence problem is somewhat alleviated by using the MPLE for initial parameter estimates.
For the time-series model, we analyze annual hurricane counts over the entire region S. The time series for 95 years, Ht, t Î 1889 ¼1993 is shown in Fig. . Ht represents the number of separate hurricanes passing through any portion of the region in year t. Although not all predictors are significant in the time-series formulation, we compare the parameter values and their errors to values obtained after dividing the region into 6° by 6° grid boxes.
Table shows the parameter estimates from the time series model. The formulation is a Poisson GLM with dispersion. A Poisson GLM is used on U.S. hurricane activity in Elsner and Bossak (2001). We note that the model is somewhat under dispersed. The dispersion is estimated from the Pearson statistic åt=1n[(Ht-ft)/SE(ft)]2 where ft is the fitted value for Ht and SE(ft) is the standard error of ft. Dakar rainfall, Azores pressure, and the warm phase of the El Niño are significant at a = .05, but the autoregressive coefficient is not.
Table shows the parameter estimates for a Poisson GLM with dispersion after modifying the time-series model by dividing the basin into grids, and adding the offset factors for latitude and longitude. The associated t-values on the model coefficients indicate that every factor is significant at a = .05. The fact that the lag-one autoregressive coefficient is significant is unexpected, since it was insignificant in the time-series model. Variances for the global predictors are significantly smaller than the variance of the time series model. Thus, the t-values on the global predictors in the second model are about three times those in the first model.
The conclusions about the parameters and their standard errors are not entirely valid, as the GLM assumes no instantaneous coupling; that is, the conditional distribution given the past is independent for each site. This assumption does not hold because our data consist of hurricanes that pass between adjacent grid boxes creating correlations in annual counts between adjacent boxes. This suggests adding spatial structure to the model in the form of coupling parameters. Note that the GLM formulation does not require temporal independence, since the assumed MLE in the GLM estimator is the same as the actual MLE for an autoregressive time series. Thus, in the absence of instantaneous coupling, we can use the Poisson GLM treating past values of the response as covariates in the model.
Model estimates for the final model are derived using the results of a pseudo likelihood estimator as inputs to the MCML estimator. The MCML estimator is then run iteratively with 1000 samples at each stage. Table shows the MCML estimates after four iterations of the estimator. Parameter estimates converge with changes between the third and fourth iteration of less than 0.6 times the estimated standard error. The root mean square change is 0.20 s.
Results from the full TPSTAR model indicate all couplings are positive and significant. Note again that the lag-one coupling, which was not significant in the time-series model, is significant in the space-time model. This is a new finding that provides evidence for hurricane path persistence over successive years. Locations that were threatened by a hurricane one year are more likely to be threatened again in the next year. As with the second model, the climate predictors are significant.
Estimates from the MCML estimator appear to be reasonable. The model takes into consideration both spatial and temporal couplings. For example, the parameter estimates are the same sign, but smaller in the final model with instantaneous couplings as compared to the Poisson GLM without instantaneous couplings, or the time-series model. This makes sense, since a positive coupling causes the expected value of any statistic to be more sensitive to changes in the predictor value than would be expected in the absence of the coupling. The parameters and the estimated standard errors from the modified MCML estimator are smaller than those obtained with the first two models. This reduction in standard error might be real owing to the addition of instantaneous coupling and offsets, but it is more likely an artifact of the apparent increase in total hurricane counts from 435 in the time series model to 1676 in the TPSTAR model.
Table shows the estimated correlations between the coupling parameters and model intercept. The correlations between the temporal and spatial coupling parameters are small, whereas the correlation between the east-west and north-south coupling parameters are considerably larger. Surprisingly the intercept is not strongly correlated to the coupling parameters.
Backward elimination is applied to the TPSTAR model to arrive at a final model of the spatial/temporal variations in seasonal hurricane activity. The procedure makes use of the estimated changes in the Schwartz's Bayesian Information Criteria (SBC) and the Akaike Information Criteria (AIC). In fact, we need to use changes in criterion levels because the TPSTAR model generates rough estimates of the deviance or predicted values. Thus, common statistics for model selection cannot be calculated accurately. The deviance requires knowledge of the normalizing constant for the distribution, which cannot be estimated accurately in the presence of strong coupling.
We test significance with the Wald test, and use this to compare the difference in SBC or AIC between two models. Both AIC and SBC are used for model selection. Although no asymptotes for either of these statistics exist for our model, for a stationary Gaussian time series, minimizing the AIC gives the model with the smallest predictive error, whereas minimizing the SBC gives a consistent model (Brockwell 1991).
Assume we have two nested models with the second model generated from the first by adding an additional factor with q levels, for a total of p parameters. Let q=[q1, ¼, qq] be the parameter vector with q Î Rp and qq Î Rq for the added factor. Let Sq be the covariance matrix for the added factor, then the Wald test is H0: qq=0 and the statistic W = qq¢Sq-1qq, and asymptotically in n, W has a c2 distribution with q degrees of freedom. Since AIC=-2logl((q)) + 2p and the asymptotic distributions for W and 2logl(q)-2logl(q1) are equal, we can approximate the change in AIC for adding the factor with q levels as DAICF » -W+2·q. Since SBC=-2logl(q)) +p·log(n), we can calculate the change in SBC for adding a factor with q levels as DSBCF » -W+q·log(n). Now, if we remove a factor with q levels then the changes in the information are DAICB » W-2·q and DSBCB » W-q·log(n).
Let us consider four versions of the TPSTAR model with instantaneous coupling. The first version is without offsets, so it has three coupling parameters, the intercept, and five climate covariates giving a total of p1=9 parameters. Version two adds the longitude factor increasing the number of parameters to p2=p1+8=17. Version three adds the latitude factor to model one increasing the number of parameters to p3=p1+8=14. Version four adds both factors to the original model for a total of p4=p1+8+5=22 parameters. First, we test the estimated change in SBC excluding the latitude factor or longitude factor. If the change is negative for either factor we remove this factor from the model. Then we rerun the model and test the change in SBC and AIC by removing the other offset factor. Table shows these results, which indicate removing the latitude factor based on SBC. Results for the final selected model are shown in Table .
Removal of the latitude factor reduces the global covariate parameters to a small degree, while increasing g1c. The parameter variances also increase slightly. These results are expected since we are moving the variance explained by the latitude factor into the unexplained variance, which shows as an increase in parameter variance. Note that although removing the latitude factor reduces the SBC it increases the AIC, so we keep the full model when using it to make forecasts.
We experimented with larger spatial neighborhoods by adding the
current period's four diagonal neighbor sites using a single
parameter to the model. These sites are the NE, SE, SW and NW
neighbor regions, with parameter g0d (the diagonal
term). We have the same model described by Eq. 2, with
Eq. 3 changed to
|
We run the modified MCML estimator for several iterations, using the previous values of the parameter estimates from the full model in Table , and an initial MPLE estimate for g0d of 0.0 We estimate that g0d=0.090,sg0d=0.021 ,t=4.3 . Adding the diagonal parameter to the model significantly affects only the estimates for g0h and g0v reducing them to 0.319 and 0.209 respectively, while increasing their standard errors to 0.017 and 0.026, respectively. The correlation matrix of the coupling parameters and the intercept show that while the other terms have correlations less than 0.07 with each other g0d is significantly correlated. This term is not only significant, but improves the model. If this term is added both the estimated AIC and SBC are reduced by 17.2 and 10.9 respectively. Thus, future models should consider both larger spatial and temporal neighborhoods.
There are concerns about increasing the size of the neighborhood. If the new couplings are positive the estimator may fail to converge. Because we have more site values, the model will required additional boundary values and parameters. One may reduce the number of parameters by using combinations of various canonical parameters. For example, the parameter g0d is really the sum of two parameters one for the NW-SE diagonal and another for the NE-SW diagonal.
The model can be used to forecast hurricane activity when lagged values of the covariates are included. At each stage of the MCML estimator, we generate samples from the distribution of Xt conditioned on {Xs : s < t } using the observations { Xt-1,¼Xt-p }, the covariates at time t, zk(t), k=1¼P the model parameters and estimated boundary values. For our application we used the full model parameters given in Table , so as to have the smallest prediction error. Also, we estimate the boundary values using the mean for each cell.
As a single test case, we generate 103 sample forecasts of hurricane activity for 1994 in each 6° by 6° region of grid S using the 1994 values of the five global covariates. Hurricane activity during 1994, which was below the long-term average, is not included at any phase of model development. The spatial distribution and intensity of hindcast values are plotted in Fig. . Table shows the mean and standard deviation for each region along with a comparison to the actual 1994 values for that region. In this table, if at least one 1994 hurricane entered the grid box the actual value is one. Results indicate reasonable agreement with observations. In general, the hindcast indicates a greater probability in grid boxes that were actually hit.
Results are summarized in Fig. which shows boxplots of the predicted probability as a function of actual occurrence. Seven grid boxes were affected by hurricanes during 1994. The average hindcast probability over those sites is 0.66 with a median value of 0.59. This compares with an average and median probability of 0.43 and 0.41, respectively for the 33 sites not affected by hurricanes during 1994. Overall this single case study supports the contention that the TPSTAR model, or similar spatial count models, might be useful tools in predicting regional hurricane activity over the North Atlantic basin.
We introduce and apply a space-time count process model to North Atlantic hurricane activity. The model uses the best-track data consisting of historical hurricane positions and intensities together with climate variables to determine local space-time coefficients of a truncated Poisson process. The model, referred to as a truncated Poisson space-time autoregressive (or TPSTAR) model, is motivated by first examining a time-series model for the entire domain. Then a Poisson generalized linear model is considered that uses grids boxes within the domain and adds offset factors for latitude and longitude. A natural extension is then made that includes instantaneous local and autoregressive coupling between the grids. A final version of the model is found by backward selection of the predictors based on values of SBC and AIC. A single hindcast is performed on the 1994 hurricane season using a model having five nearest neighbors and statistically significant couplings. The parameters in the TPSTAR model are estimated using MPLE. The model showed promise as a potential forecast tool.
Several conclusions concerning the application of the TPSTAR model to seasonal North Atlantic hurricane activity are reached:
Acknowledgements. The study was partially funded
by the National Science Foundation (ATM-0086958) and the Risk
Prediction Initiative of the Bermuda Biological Station for
Research (RPI-99-001).
|
|
|
|
|
|
Value | Standard Errors | t-value | units | |
Intercept | 198.475 | 68.996 | 2.88 | |
g1c | 0.037 | 0.0194 | 1.88 | |
Cold | 0.193 | 0.108 | 1.78 | |
Warm | -0.302 | 0.126 | -2.38 | |
Dakar | 0.594 | 0.250 | 2.40 | years per meter |
Azores | -0.149 | 0.049 | -3.05 | per millibar |
Iceland | -0.046 | 0.027 | -1.71 | per millibar |
Value | Standard Errors | t-value | units | |
g1c | 0.094 | 0.028 | 3.34 | per hurricane |
Intercept | 244.823 | 34.009 | 7.19 | |
lat1 | 0.101 | 0.040 | 2.53 | |
lat2 | 0.097 | 0.020 | 4.92 | |
lat3 | 0.070 | 0.016 | 4.44 | |
lat4 | 0.035 | 0.013 | 2.62 | |
lat5 | -0.067 | 0.016 | -4.16 | |
long1 | 0.034 | 0.069 | 0.50 | |
long2 | 0.052 | 0.035 | 1.49 | |
long3 | 0.023 | 0.022 | 1.07 | |
long4 | 0.016 | 0.015 | 1.10 | |
long5 | 0.012 | 0.011 | 1.04 | |
long6 | -0.014 | 0.010 | -1.49 | |
long7 | -0.038 | 0.009 | -4.23 | |
long8 | -0.068 | 0.009 | -7.36 | |
Warm | -0.391 | 0.064 | -6.12 | |
Cold | 0.293 | 0.052 | 5.66 | |
Dakar | 1.260 | 0.120 | 10.47 | year per meter |
Azores | -0.193 | 0.029 | -8.12 | per millibar |
Iceland | -0.049 | 0.013 | -3.66 | per millibar |
| MCML Estimator | |||||||
|
| Standard Error |
| |||||
g0h | 0.389 | 0.358 | 0.014 | 24.90 | per hurricane | |||
g0v | 0.250 | 0.287 | 0.017 | 16.75 | per hurricane | |||
g1c | -0.002 | 0.079 | 0.020 | 4.03 | per hurricane | |||
intercept | 57.412 | 148.537 | 18.473 | 8.04 | ||||
lat1 | -0.072 | -0.010 | 0.038 | -0.27 | ||||
lat2 | 0.012 | 0.023 | 0.017 | 1.35 | ||||
lat3 | 0.039 | 0.022 | 0.013 | 1.74 | ||||
lat4 | 0.044 | 0.007 | 0.011 | 0.60 | ||||
lat5 | -0.012 | -0.047 | 0.013 | -3.62 | ||||
long1 | -0.197 | -0.116 | 0.076 | -1.53 | ||||
long2 | -0.056 | -0.003 | 0.035 | -0.09 | ||||
long3 | -0.047 | -0.009 | 0.020 | -0.45 | ||||
long4 | -0.042 | -0.009 | 0.013 | -0.68 | ||||
long5 | -0.024 | -0.013 | 0.010 | -1.32 | ||||
long6 | -0.024 | -0.032 | 0.009 | -3.73 | ||||
long7 | -0.018 | -0.039 | 0.008 | -4.77 | ||||
long8 | -0.037 | -0.055 | 0.008 | -6.79 | ||||
Warm | -0.124 | -0.230 | 0.038 | -6.02 | ||||
Cold | 0.072 | 0.192 | 0.027 | 7.15 | ||||
Dakar | 0.368 | 0.768 | 0.068 | 11.31 | year per meter | |||
Azores | -0.042 | -0.120 | 0.013 | -9.14 | per millibar | |||
Iceland | -0.015 | -0.028 | 0.007 | -3.96 | per millibar |
g0,h | g0,v | g1,c | Intercept | |
g0h | 1.000 | -0.763 | -0.075 | -0.022 |
g0v | -0.763 | 1.000 | 0.003 | -0.084 |
g1c | -0.075 | 0.003 | 1.000 | 0.103 |
Intercept | -0.022 | -0.084 | 0.103 | 1.000 |
Factor | Resulting | Number of | Information Changes | ||
Model | Removed | Model | levels (q) | DAICB | DSBCB |
4 | latitude | 2 | 5 | 10.9 | -20.3 |
4 | longitude | 3 | 8 | 95.2 | 45.3 |
2 | longitude | 1 | 8 | 92.5 | 42.6 |
Value | Standard Errors | t-value | units | |
g0h | 0.359 | 0.015 | 24.78 | per hurricane |
g0v | 0.289 | 0.017 | 16.64 | per hurricane |
g1c | 0.084 | 0.020 | 4.09 | per hurricane |
Intercept | 143.278 | 19.335 | 7.41 | |
long1 | -0.111 | 0.077 | -1.44 | |
long2 | -0.014 | 0.034 | -0.40 | |
long3 | -0.008 | 0.020 | -0.41 | |
long4 | -0.006 | 0.013 | -0.44 | |
long5 | -0.015 | 0.010 | -1.54 | |
long6 | -0.033 | 0.009 | -3.83 | |
long7 | -0.038 | 0.008 | -4.64 | |
long8 | -0.054 | 0.008 | -6.73 | |
Warm | -0.222 | 0.041 | -5.37 | |
Cold | 0.183 | 0.028 | 6.48 | |
Dakar | 0.740 | 0.072 | 10.29 | years per meter |
Azores | -0.115 | 0.014 | -8.40 | per millibar |
Iceland | -0.027 | 0.008 | -3.66 | per millibar |
Region | Sample | Actual | Region | Sample | Actual | ||||
Center | Mean | Std | Value | Center | Mean | Std | Value | ||
83 | °W 15 °N | 0.398 | 0.616 | 0 | 71 | °W 27 °N | 0.864 | 1.121 | 0 |
77 | °W 15 °N | 0.408 | 0.663 | 1 | 65 | °W 27 °N | 0.728 | 0.972 | 0 |
71 | °W 15 °N | 0.408 | 0.617 | 0 | 59 | °W 27 °N | 0.553 | 0.763 | 0 |
65 | °W 15 °N | 0.495 | 0.803 | 0 | 53 | °W 27 °N | 0.369 | 0.642 | 0 |
59 | °W 15 °N | 0.408 | 0.601 | 0 | 47 | °W 27 °N | 0.301 | 0.575 | 0 |
53 | °W 15 °N | 0.243 | 0.514 | 0 | 77 | °W 33 °N | 0.699 | 0.968 | 0 |
47 | °W 15 °N | 0.214 | 0.517 | 0 | 71 | °W 33 °N | 1.000 | 1.358 | 1 |
95 | °W 21 °N | 0.427 | 0.651 | 0 | 65 | °W 33 °N | 0.854 | 1.175 | 1 |
89 | °W 21 °N | 0.369 | 0.642 | 0 | 59 | °W 33 °N | 0.592 | 0.785 | 0 |
83 | °W 21 °N | 0.592 | 0.810 | 0 | 53 | °W 33 °N | 0.505 | 0.862 | 0 |
77 | °W 21 °N | 0.515 | 0.815 | 1 | 47 | °W 33 °N | 0.282 | 0.584 | 0 |
71 | °W 21 °N | 0.524 | 0.765 | 0 | 71 | °W 39 °N | 0.641 | 0.838 | 0 |
65 | °W 21 °N | 0.437 | 0.763 | 0 | 65 | °W 39 °N | 0.835 | 1.103 | 0 |
59 | °W 21 °N | 0.330 | 0.584 | 0 | 59 | °W 39 °N | 0.563 | 1.026 | 1 |
53 | °W 21 °N | 0.233 | 0.597 | 0 | 53 | °W 39 °N | 0.476 | 0.739 | 0 |
47 | °W 21 °N | 0.165 | 0.422 | 0 | 47 | °W 39 °N | 0.291 | 0.588 | 0 |
95 | °W 27 °N | 0.495 | 0.655 | 0 | 65 | °W 45 °N | 0.369 | 0.610 | 0 |
89 | °W 27 °N | 0.379 | 0.643 | 0 | 59 | °W 45 °N | 0.262 | 0.641 | 0 |
83 | °W 27 °N | 0.592 | 0.785 | 1 | 53 | °W 45 °N | 0.262 | 0.523 | 0 |
77 | °W 27 °N | 0.670 | 0.974 | 1 | 47 | °W 45 °N | 0.204 | 0.531 | 0 |