Exploring Spatial Patterns of Virginia Tornadoes Using Kernel Density and Space-Time Cube Analysis (1960–2019)

: This study evaluates the spatial-temporal patterns in Virginia tornadoes using the National Weather Service Storm Prediction Center’s Severe Weather GIS (SVRGIS) database. In addition to descriptive statistics, the analysis employs Kernel Density Estimation for spatial pattern analysis and space-time cubes to visualize the spatiotemporal frequency of tornadoes and potential trends. Most of the 726 tornadoes between 1960–2019 occurred in Eastern Virginia, along the Piedmont and Coastal Plain. Consistent with other literature, both the number of tornadoes and the tornado days have increased in Virginia. While 80% of the tornadoes occurred during the warm season, tornadoes did occur during each month including two deadly tornadoes in January and February. Over the 60-year period, a total of 28 people were killed in the Commonwealth. Most tornado activity took place in the afternoon and early evening hours drawing attention to the temporal variability of risk and vulnerability. Spatial analysis results identify signiﬁcant, non-random clusters of tornado activity and increasing temporal frequency. While this study improves weather-related literacy and addresses a need in the Commonwealth, more research is necessary to further evaluate the synoptic and mesoscale mechanisms of Virginia tornadoes.


Introduction
Aimed to enhance the awareness of severe weather and climate hazards for the Commonwealth of Virginia, this study builds on similar regional tornado assessments [1]. Our general objective is to provide a spatial-temporal analysis of tornadoes across Virginia, with the specific objective of detecting and describing spatial patterns of tornadoes in Virginia, particularly to explore clusters or significant hot spots. Secondarily, we intend to explore available data for any trends, including spatial anomalies of tornado severity or emergent patterns over time. Our research question posits that, like intensity of tornado activity in the upper (inner) Coastal Plain and Piedmont of North Carolina [2], we may find a significant clustering of tornadoes in Virginia in the lee of the Blue Ridge and along the Fall Line where the Coastal Plain and Piedmont adjoin, or potentially on the outer coastal plain associated with passing tropical storms or severe air mass thunderstorms. Meteorological conditions generally supporting this include the availability of thermal troughing along the Piedmont front range, urban agglomerations supporting thermal forcing, and onshore flow from sea breezes and pre-frontal warm air advection and unstable tropical airmasses along the Coastal Plain. Combined with frontal uplift and leeward side vorticity associated with the Blue Ridge Mountains, the eastern half of the state should support elevated tornado activity with strong storms and squall lines [3,4].
While the influence of climate change on extreme events is well documented [10], the multitude of factors associated with severe weather indicates no clear trend as to an increasing frequency or intensity of tornadoes [11][12][13]. Recently, Elsner et al. (2019) noted an upward trend in the power of tornadoes due to changes in convective storm environments [14]. Other researchers have sought to document an east to southeasterly shift of the traditionally perceived Tornado Alley, sparking a debate over the existence and popular interest in the "Dixie Alley," in the Southeast of the U.S. [15,16]. Future changes in risk and exposure may increase the disaster potential [17], thus a comprehensive assessment of tornado climatology is necessary to inform stakeholders in the Commonwealth of Virginia.
A suite of spatial analysis techniques is available to explore and quantify the spatial pattern and timing (including annual and seasonal occurrence as well as climatological trends). Among the available techniques, this project applies linear spatial pattern analysis to quantify the spatial distribution of documented tornadoes, seeking to discern regional hotspots and potential clusters. Tornadoes and clusters of outbreaks can be mapped as points or choropleths at small scales to infer and communicate broad risk patterns between regions or states [18]. Cartographic visualization and Kernel Density Estimates (KDE) are appropriate techniques for exploring such spatial patterns and subsequent analytical inferences using symbology and visual techniques (e.g., line density plots, 3D, or surface analyses). KDE methods have a long history in spatial analysis and frequent application to point-based, discrete events such as crimes [19] or road accidents [20]. With respect to severe weather, Brown et al. (2016) explored the seasonal as well as spatial pattern of tornadoes by comparing three cities and aiming to control for the observational bias of observer density in urban areas [21]. Coleman and Dixon (2014) mapped the number of tornadoes, using cartographic visualization to explore patterns, derived isopleths of tornado counts, and calculated the density and an average annual path length using KDE techniques [22]. They further plotted different severity and path-lengths of tornado tracks to ascertain potential damage impacts and area affected. KDE techniques, however, are not without their faults, and sensitive parameters such as the spatial search distance (bandwidth selection) are advisable to use iteratively after visual inspection [23]. Further, "Space-Time Cubes" (STCs) have a history in their conceptualization dating to Hagestrand (1970) yet also newly accessible computational tools in Esri ArcGIS Pro [24,25]. STCs are applicable to discerning and evaluating spatial patterns intrinsically over time, or spatialtemporal "trajectories" of phenomena. STCs have shown promise for analytical potential where domain experts (e.g., tornado forecast meteorologists, climatologists, and emergency managers) can be empowered with deeper understanding of the spatial-temporal patterns relevant to tornado forecasting, climate-sensitivity, or disaster risk reduction. Research on STCs, however, has shown limited success where STCs have not provided a usercentered perspective and analytical capabilities have not been functionally demonstrated, and careful use of visual variables and depth cues has been suggested to improve their functionality [26].
With most of the Virginia population situated along and east of I-95, understanding where and when tornado activity has taken place may help inform hazard researchers interested in mitigating risk ( Figure 1). Currently, no comprehensive assessment of tornadoes for the Commonwealth of Virginia exists. This study also strengthens the weather-related information for the Commonwealth of Virginia, filling a knowledge gap within the Virginia Climate Office [27,28]. The development of such assessments helps build synergies across local and regional stakeholders beyond the Weather Forecast Office (WFO). Virginia's Department of Emergency Management conducts an annual Tornado Exercise in March, aligning with the National Weather Service's annual Severe Weather Awareness Week [29].
Having updated and spatially refined information on the pattern and timing of tornado hazards could thus inform the implementation of these efforts, such as providing forecasters and science communicators graphical educational materials, guiding counties or areas for media campaigns and outreach, or also providing preparatory background to forecasters and emergency responders working in areas predisposed to tornado events.
Weather GIS (SVRGIS) database (https://www.spc.noaa.gov/gis/svrgis/, accessed on 1 January 2021). The SVRGIS database [30] contains severe weather event information in tabular and GIS formats and includes verified tornado information on the track path, intensity, and associated impacts (e.g., categorized estimated property loss). In addition to enabling mapping and spatial analyses, these data enabled descriptive statistics to be compiled for tornadoes in Virginia . Other studies noted database deficiencies and tornado reporting bias [31][32][33]. The 60-year time period provides a long-term record for analyzing trends. The SVRGIS tornado data include unique tornado identifier, date of occurrence (local time, day, month, and year), U.S. state, EF-scale magnitude, track length, tornado width, and begin-and end-points (latitude and longitude).

Datasets
Tornado information was downloaded from the National Oceanic and Atmospheric Administration's (NOAA) National Weather Service Storm Prediction Center Severe Weather GIS (SVRGIS) database (https://www.spc.noaa.gov/gis/svrgis/, accessed on 1 January 2021). The SVRGIS database [30] contains severe weather event information in tabular and GIS formats and includes verified tornado information on the track path, intensity, and associated impacts (e.g., categorized estimated property loss). In addition to enabling mapping and spatial analyses, these data enabled descriptive statistics to be compiled for tornadoes in Virginia . Other studies noted database deficiencies and tornado reporting bias [31][32][33]. The 60-year time period provides a long-term record for analyzing trends. The SVRGIS tornado data include unique tornado identifier, date of occurrence (local time, day, month, and year), U.S. state, EF-scale magnitude, track length, tornado width, and begin-and end-points (latitude and longitude).

Tornado Track Mapping
The initial exploratory phase of the spatial analysis focused on mapping the extent of tornado tracks extracted from SVRGIS for the period of study. Tornado track data were extracted by a spatial selection of the full conterminous U.S. database, including a buffer of approximately 100 km and any tracks that intersected the boundaries of the Commonwealth of Virginia (including starting or ending within the state boundary). Including tornado tracks within a buffer distance ensured that no tracks having any extent within Virginia were excluded or would erroneously affect the density and proximity analysis. Including tornadoes proximate but not within the state is also justifiable to mitigate against potential neighborhood statistical calculations around the periphery of the state. As a step of quality assurance and exploratory analysis, we first visually analyzed all the retrieved tornado checks, observing for potential gross errors or spurious digitizing errors that may affect subsequent spatial analyses. One tornado track was found to be anomalous and oriented as a perfectly south-north arc; the tornado event was cross-checked and determined that a geographic coordinate was likely erroneously entered into SVRGIS for this unique track. The begin-and end-points suggested a long track and gross error near Middlesex County and between towns of Deltaville and Whitestone on the Middle Peninsula of the Chesapeake Bay. In addition, this single long track could not be verified from any media or WFO reported tornado events. Likely an erroneous input of the lift point coordinate, this datapoint was removed out of concern for the error propagation of the track in subsequent spatial analyses.
Inspection of the general pattern of tornado tracks suggested the use of linear kernel density techniques (to capture the line density across space, e.g., track length per unit area) and replication of the map to see if the pattern and density would hold for differing levels of tornado intensity with the Enhanced Fujita (EF) scale (the Enhanced F scale adopted by the National Weather Service in 2007) [34,35]. Visualizations of tornado frequency and density were created using kernel density analysis. Point kernel density estimation produces a continuous density surface based on the tornado tracks, highlighting areas of Virginia where tornado characteristics such as strength and length are more pronounced [36]. Using Virginia tornado track lines segmented by magnitude, we produced a 100-m cell resolution density surface that illustrated differences in frequency of higher-magnitude tornadoes. It was revealed that Northern Virginia near Washington, D.C., Richmond (central Virginia), and the Hampton Roads region (Southeastern Virginia) have experienced a higher frequency of strong tornadoes than other parts of the Commonwealth. We also examined the distribution of Virginia's longest tracking tornadoes. Using the Virginia tornado track line length, another 100-m cell resolution density surface that illustrated differences in the tornado travel distance was developed.

Spatial Statistics and Kernel Density Analysis
Our approach to visually assess observed tornadoes included both cartographic visualizations using a composite map of all recorded events and an index of tornado density interpolated by kernel line density estimation. To expand beyond choroplethic mapping and Modifiable Areal Unit Problem (MAUP) effects [37], we sought a spatial representation of the overall pattern and density (tornado track lengths per unit area) and opted to implement kernel line density interpolation. The selection of the search distance for KDE followed iterative visual analysis of results and computation of a nominal search distance [25]. The KDE output grid resolution was also evaluated after the methods of Scott (1992) and rule of thumb for KDE to evaluate the standard distance of measured observation [38]. Our method used a KDE approach with polylines, similar to a linear KDE analysis of road density presented by Tang et al. (2015) [39]. After iterative computation with KDE bandwidths between 1 km to 10 km, we ultimately implemented a 4 km bandwidth very close to the default 1/30th ratio of the maximum divided by the minimal length-width ratio of the study area (Commonwealth of Virginia) and similarly implemented by Cai et al. (2013) [40]. This selection of bandwidth provided for visual discrimination within counties and along potential patterning to capture topographic influences, such as the Blue Ridge Mountains, Ridge and Valley Province, and the Fall Zone running north-south. A coarser bandwidth may overly smooth the density of tornado patterns along these features, while a finer scale visually obscured a subregional trend.

Space-Time Cube
As a temporally variable pattern, a space-time analysis or pattern persistence metrics were deemed worth exploring. For this geographic extent, having evidently patchy areas of concentrated tornado activity and relatively sparse areas (even empty), we opted to explore the application of new Space-Time cubes.
Our space-time cube analysis workflow using the Esri's ArcGIS Pro software required the disaggregation of tornado track line data into point data evenly spaced at 5-km intervals along each original track line. Each resulting point contained key tornado attributes including magnitude, direction, and a unique date/time stamp. All data points were then binned into decadal time step intervals between 1960 and 2020. A net CDF format data cube was aggregated using the time-binned tornado points to produce hexbin polygons, each covering 250 square kilometers. The rationale for hexbinning (i.e., hexagonal tessellation) is to primarily optimize equidistant measurements between centroids of areal bins [41]. Square or rectangular tessellations also tend to break-up or distort curvilinear patterns owing to acute angles. There is also an inherent benefit to using hexagons to reduce the spatial sampling bias of phenomena, since they have a low perimeter-area ratio. Since our study included the relatively large area of Virginia (east-west~690 km), hexagons would also show less distortion owing to the curvature of the earth as compared to a fishnet or grid tessellation [41].
Development of the netCDF cube allowed us to visualize trends in tornado activity over time as a 2-dimensional hexagonal polygon surface wherein temporal trends for each polygon were exposed using a "trend description" attribute. Our approach sought to discover and describe spatial-temporal patterns of tornado activity, a research challenge appropriate to applying advances in geocomputation and visualization [42]. To facilitate analysis temporal trends in the spatial pattern of tornadoes, we opted to enhance the cartographic analysis of the decadal small multiples using space-time cube computation. This approach emphasized the discrete, event-based nature of tornadoes and referenced time-scale changes by decade. The cubic spatial-temporal compartmentalization would also provide for temporal assessments, such as changes in rates or patterns of tornado activity by decade across the study area. Rather than computing Markov Chains for discrete state transitions (e.g., classifying tornado activity into ranked levels of intensity), we instead chose to use the space-time cube inferentially to explore and assess potential gross changes in activity and statistical significance [43].

Statistical Descriptive Analysis
Of the 69 FEMA disaster declarations in Virginia, five are associated with tornadoes. One of these was linked to a tropical cyclone, Gaston in 2004 [44]. Over the course of 338 tornado days, a total of 727 tornadoes occurred over the time period (Figures 2 and 3). Most tornadoes occurred in the Piedmont and Coastal Plain region of Virginia.  Across the 60 years, 28 people were killed across 15 individual tornado days ( Table  1). It is worth noting that the most deadly tornado outbreak in Virginia's history occurred outside of the time period of this study, on 30 September 1959. On this day, 12 people died. With an estimated $20 million damage, the costliest tornado occurred on 15 April 2018 near Lynchburg, Virginia. In the last 5 years (2015-2019), total damages exceeded $65 million. Other studies highlight the potential causality and economic impacts of tornadoes [46,47]. For example, Fricker and Elsner (2020) draw attention to how community socioeconomic and demographic characteristics increase vulnerability to unusually devasting tornadoes [48]. Mere geographic location is not the only contributing factor to vulnerability.
Only two systems resulted in EF 4 tornadoes: the historic 6 August 1993 outbreak that transversed North Carolina and Southeast Virginia, killing 4 and injuring over 250, and the 24 September 2001 event which was part of a severe weather outbreak in the Washington, DC-Baltimore metro area causing over $100 million in damage. In all, nine days with more than 10 tornadoes occurred (Table 1). Associated with a well-forecast When comparing the earlier  to the more recent time periods (1990-2019), average tornadoes (tornado days) per year increased from 5.9 (4.4) to 18.3 (6.9) per year ( Figure 3). As shown in Figure 3, an upward trend is observed over the full study period. We may infer that decadal changes in tornadoes reflect increasing storminess from tropical systems affecting the coast or inland "backdoor" storms, that land use change and suburban sprawl could indirectly increase the observation and reporting of tornadoes, and/or that technological advances in weather radar across coastal, Piedmont, and mountain areas could impart the increased tornado detection and verification [32,33,45]. The precise atmospheric mechanics behind the changes are beyond the scope of this research and draw attention to the need for continued investigation. Across the 60 years, 28 people were killed across 15 individual tornado days (Table 1). It is worth noting that the most deadly tornado outbreak in Virginia's history occurred outside of the time period of this study, on 30 September 1959. On this day, 12 people died. With an estimated $20 million damage, the costliest tornado occurred on 15 April 2018 near Lynchburg, Virginia. In the last 5 years (2015-2019), total damages exceeded $65 million. Other studies highlight the potential causality and economic impacts of tornadoes [46,47]. For example, Fricker and Elsner (2020) draw attention to how community socioeconomic and demographic characteristics increase vulnerability to unusually devasting tornadoes [48]. Mere geographic location is not the only contributing factor to vulnerability.
Only two systems resulted in EF 4 tornadoes: the historic 6 August 1993 outbreak that transversed North Carolina and Southeast Virginia, killing 4 and injuring over 250, and the 24 September 2001 event which was part of a severe weather outbreak in the Washington, DC-Baltimore metro area causing over $100 million in damage. In all, nine days with more than 10 tornadoes occurred (Table 1). Associated with a well-forecast extratropical system, two of the longest tornadoes were part of a major tornado outbreak in the United States, resulting in over 175 confirmed tornadoes (14-16 April 2011). In terms of seasonal cycle, a bimodal distribution was found with spring (April) and autumn (September) peaks ( Figure 4). While 80% of the tornadoes occurred during the warm season (April to September), tornadoes did occur during each month including two deadly tornados in  January and February. Like other studies, most tornado activity took place in the afternoon and early evening hours ( Figure 5). extratropical system, two of the longest tornadoes were part of a major tornado outbreak in the United States, resulting in over 175 confirmed tornadoes (14-16 April 2011). In terms of seasonal cycle, a bimodal distribution was found with spring (April) and autumn (September) peaks ( Figure 4). While 80% of the tornadoes occurred during the warm season (April to September), tornadoes did occur during each month including two deadly tornados in January and February. Like other studies, most tornado activity took place in the afternoon and early evening hours ( Figure 5).    extratropical system, two of the longest tornadoes were part of a major tornado outbreak in the United States, resulting in over 175 confirmed tornadoes (14-16 April 2011). In terms of seasonal cycle, a bimodal distribution was found with spring (April) and autumn (September) peaks ( Figure 4). While 80% of the tornadoes occurred during the warm season (April to September), tornadoes did occur during each month including two deadly tornados in January and February. Like other studies, most tornado activity took place in the afternoon and early evening hours ( Figure 5).

Patterns, Trends and Spatial Statistics
A trend of longer tracking tornadoes was evidenced for areas immediately south of Richmond (central VA), Washington, DC, USA (Northern VA), and the Middle Peninsula region (Eastern VA) ( Figure 6). Three principal clusters are visually evident in the KDE density map surface of Figure 6, two large clusters are visible along the Fall Line in Northern Virginia and southeast of Petersburg near the Richmond Metro area, while a third is distinguishable on the outer Coastal Plain from Suffolk to Chesapeake. These derived patterns are inferred as significant clustering, supported by a statistically significant nearest neighbor ratio (N r = 0.5996, z-score -18.434, p < 0.000). The Piedmont/Fall Line clusters are somewhat expected, owing to leeward side troughing and frontogenesis where subtropical airmass along the coastal plain and Chesapeake Bay increase the baroclinicity and moisture advection, especially in the case of approaching cold fronts trailing Midwest and Atlantic storm track shortwave lows. The cluster evident in Southeast Virginia, the rural areas of Hampton Roads metropolitan area, is suggestive of tropical storm influences from passing, recurving tropical storms and hurricanes. Results also depict minor clusters in Southwest Virginia in the Southern Shenandoah Valley and the inner Piedmont (also where favorable lee-side troughing exists). These areas show a prevalence for short-  track tornadoes suggestive of somewhat terrain-controlled airmass thunderstorms in addition to frontal systems. We also found that developing high-resolution (100 m) density surfaces, as suggested by Deng et al. (2016), provided for the inference of clusters of tornado hazard [49]. The multi-modality of the clusters (Piedmont and Coastal) suggests a disjunct process, providing our hypothesis generation as to frontal, tropical, and combined tornadogenic events. However, quantifying the event conditions and origins requires a more careful consideration of the temporal circumstances as well as the synoptic and dynamic meteorological forcing at the scale of days or hours. Nonetheless, to further explore and develop empirical awareness of the tornado patterns, we turn to the results of the space-time cube analysis. hazard [49]. The multi-modality of the clusters (Piedmont and Coastal) suggests a disjunct process, providing our hypothesis generation as to frontal, tropical, and combined tornadogenic events. However, quantifying the event conditions and origins requires a more careful consideration of the temporal circumstances as well as the synoptic and dynamic meteorological forcing at the scale of days or hours. Nonetheless, to further explore and develop empirical awareness of the tornado patterns, we turn to the results of the spacetime cube analysis.

Space-Time Cube Results
Space time cube analysis techniques which offer an opportune visual way to study the relationship between events with time-and space-related variables [41,50]. Trend attribution includes either an up or down trend description with an associated trend confidence (99%, 95%, 90%). A trend description of "No Significant Trend" is also possible. These trend surface hexagonal visualizations confirmed that areas in Northern Virginia, the Richmond metropolitan area, and the Hampton Roads area all exhibited upward trends with a high degree of confidence. We found that using space-time cube methods in concert with kernel density analyses of tornado path locations, strength, and length, allowed for the delineation of a large area of Virginia that exhibits high tornado activity (Figure 7). Temporal trends, summarized in hexagon bins in Figure 7a, depict two general zones of marked upward frequency trends in the Northern and central Piedmont-Fall Zone transition (overlapping these clusters observed in Figure 6) and a moderate upward trend surrounding the inner Coastal Plain (shown as an arc of orange hexes along the inner coastal plain from the Middle Peninsula south to Williamsburg and Suffolk.) While these same areas are qualitatively considered exurban land use growth areas since the 1980s, they also are on the inner range of some of the larger passing tropical cyclones of recent decades. For example, Hurricane Isabel (2003) and Hurricane Irene (2011) produced several tornadoes in that area. The STC methodology produced an estimate of statistical significance of non-random temporal trends for each hexagon, with the highlighted colors having increasing confidence through the 1951-2018 time-period. The zonal west-east pattern, when accounting for the cumulative spatial patterns of KDE and STC analysis, lead us to suggest that Virginia has a modest to minor area of elevated tornado activity (highlighted in Figure 7b.) This zone is qualitatively described as having both increased prevalence of tornadoes spatially (KDE results) and an increasing trend in their abundance temporally (STC results).

Space-Time Cube Results
Space time cube analysis techniques which offer an opportune visual way to study the relationship between events with time-and space-related variables [41,50]. Trend attribution includes either an up or down trend description with an associated trend confidence (99%, 95%, 90%). A trend description of "No Significant Trend" is also possible. These trend surface hexagonal visualizations confirmed that areas in Northern Virginia, the Richmond metropolitan area, and the Hampton Roads area all exhibited upward trends with a high degree of confidence. We found that using space-time cube methods in concert with kernel density analyses of tornado path locations, strength, and length, allowed for the delineation of a large area of Virginia that exhibits high tornado activity (Figure 7). Temporal trends, summarized in hexagon bins in Figure 7a, depict two general zones of marked upward frequency trends in the Northern and central Piedmont-Fall Zone transition (overlapping these clusters observed in Figure 6) and a moderate upward trend surrounding the inner Coastal Plain (shown as an arc of orange hexes along the inner coastal plain from the Middle Peninsula south to Williamsburg and Suffolk.) While these same areas are qualitatively considered exurban land use growth areas since the 1980s, they also are on the inner range of some of the larger passing tropical cyclones of recent decades. For example, Hurricane Isabel (2003) and Hurricane Irene (2011) produced several tornadoes in that area. The STC methodology produced an estimate of statistical significance of non-random temporal trends for each hexagon, with the highlighted colors having increasing confidence through the 1951-2018 time-period. The zonal west-east pattern, when accounting for the cumulative spatial patterns of KDE and STC analysis, lead us to suggest that Virginia has a modest to minor area of elevated tornado activity (highlighted in Figure 7b.) This zone is qualitatively described as having both increased prevalence of tornadoes spatially (KDE results) and an increasing trend in their abundance temporally (STC results).

Discussion
This study provided the first comprehensive, spatial-temporal assessment of Virginia tornadoes . With most tornadoes occurring in Eastern Virginia, the geographic landscape offers insight as to the spatial patterns. Virginia experiences tropical storm activity which may serve as a catalyst for severe weather outbreaks. Leeward side dynamics associated with the Blue Ridge and Appalachian Mountains enhances vorticity. Onshore flow and warm air advection along the Coastal Plain provides unstable physical environments conducive for tornadic activity. Long-term and lagged effects of soil moisture deficit or drought with interspersed abundance in river valleys, reservoirs or topographic depressions may affect the seasonal tornado activity [51]. Davis et al. (1997) identified five synoptic patterns associated with Virginia tornadoes: thermodynamic/northwest flow air masses, a mixed thermodynamic/dynamic category, a polar front jet situation, cool season midlatitude cyclone events, and Great Plains type tornadic air masses [10,52]. Results indicate Eastern Virginia as a high-risk zone for tornadic activity, and a heightened focus for severe weather preparedness and warning coordination for the National Weather Service (Wakefield and Sterling weather forecast offices).

Discussion
This study provided the first comprehensive, spatial-temporal assessment of Virginia tornadoes (1960-2019). With most tornadoes occurring in Eastern Virginia, the geographic landscape offers insight as to the spatial patterns. Virginia experiences tropical storm activity which may serve as a catalyst for severe weather outbreaks. Leeward side dynamics associated with the Blue Ridge and Appalachian Mountains enhances vorticity. Onshore flow and warm air advection along the Coastal Plain provides unstable physical environments conducive for tornadic activity. Long-term and lagged effects of soil moisture deficit or drought with interspersed abundance in river valleys, reservoirs or topographic depressions may affect the seasonal tornado activity [51]. Davis et al. (1997) identified five synoptic patterns associated with Virginia tornadoes: thermodynamic/northwest flow air masses, a mixed thermodynamic/dynamic category, a polar front jet situation, cool season midlatitude cyclone events, and Great Plains type tornadic air masses [10,52]. Results indicate Eastern Virginia as a high-risk zone for tornadic activity, and a heightened focus for severe weather preparedness and warning coordination for the National Weather Service (Wakefield and Sterling weather forecast offices).
Studies indicate shifts in the frequency of tornadoes, particularly multi-tornado days, over the past few decades [8,14,53]. Gensini and Brooks (2018) found an upward trend a) b)  [54]. While an increasing trend in Virginia tornadoes was found, this study did not explicitly evaluate the causal mechanisms. Improved reporting, technological monitoring, and dataset bias may also contribute to the observed shifts [31][32][33]. While cold-season tornadoes seem uncommon, Childs et al. (2018) found an increasing trend associated with these storms [55]. The study found conditions during La Niña episodes, characterized by warm, moist conditions across the southeast and an anomalous trough across the Western U.S., to be more favorable for cold-season severe weather development. The study also noted a relationship with the Arctic Oscillation; future analysis may consider the role of large-scale global circulation and how teleconnection patterns influence Virginia tornadoes [56]. Other studies note database deficiencies and tornado reporting bias [32,33,45,56,57]. Untrained witnesses [45] and improved detection methods [58] certainly present limitations, particularly in the early time period of this study.
Broomwell et al. (2020) notes the knowledge gaps in public knowledge about hazardous weather risks, particularly overnight and wintertime storms [59]. Using geovisualization as a tool, this study increases severe weather salience in the Commonwealth of Virginia. This analysis fills a knowledge gap with respect to the spatial patterns, changes, and impacts of tornadoes in the Commonwealth of Virginia. Efforts such as NWS-Skywarn spotter training may benefit from the results of this research and the content may be used in K-16 educational networks.

Conclusions
This study was the first statewide evaluation of the Virginia tornadoes. Most tornadoes occurred (1) during the afternoon and evening hours, (2) during the warm season, and (3) along the eastern flank of the Piedmont and Coastal Plain. The physiographic features certainly play a role in the development of severe weather, but the precise impact of these environments is beyond the scope of this study and is an area of continued research [60,61]. Future research may consider the intersection of tornado tracks with social factors such as building material or socio-economic status [48,62,63]. While the relative number and severity of tornadoes may not rival that of the U.S. Midwest and Great Plains, the observational data explored are suggestive of a statistically significant concentration (clusters) of tornado activity as well as indications of temporal increases. Outlining the areas of higher tornado activity may facilitate further research on the meteorological mechanisms and potential for forecasting and communication improvements. The Commonwealth of Virginia and NWS host an annual Tornado Drill to exercise risk preparedness, communications, and awareness. Drawing attention to the evolving spatial and temporal patterns of tornado hazards and their causal weather and potential risks can improve the overall hazard management. Such information can help build a more detailed database for Virginia-related hazards and to inform various stakeholders, including emergency managers and the general public, as to the spatial-temporal variability of tornadoes.