hurricaneexposure
packageThis package allows users to explore and map data of county-level exposures to Atlantic-basin tropical storms between 1988 and 2018 for a number of storm hazards (e.g., wind, rain, flood events, distance from the storm track) for counties in the eastern half of the United States. Users can map exposures by county for a selected storm and can also identify all storms to which selected counties were exposed, based on user-specified thresholds (for example, the package allows the user to identify all storms that brought wind of 34 knots or higher to Miami-Dade County in Florida).
If you use this package and the data in the associated data package
(hurricaneexposuredata
) for research, please cite both
packages. In particular, be sure to include the version of the packages
that you used, as this will make your research project more
reproducible, and the data will likely be updated as we get further
years of data and find improved ways to measure tropical storm exposure.
The two packages can be cited as:
Anderson B, Yan M, Ferreri J, Crosson W, Al-Hamdan M, Schumacher A and Eddelbuettel D (2020). hurricaneexposure: Explore and Map County-Level Hurricane Exposure in the United States. R package version 0.1.1, <URL: https://github.com/geanders/hurricaneexposure>.
Anderson B, Schumacher A, Crosson W, Al-Hamdan M, Yan M, Ferreri J, Chen Z, Quiring S and Guikema S (2020). hurricaneexposuredata: Data Characterizing Exposure to Hurricanes in United States Counties. R package version 0.1.0, <URL: https://github.com/geanders/hurricaneexposuredata>.
To generate BibTex entries for the packages, you can use the
citation
function in R.
This package depends on data in a data package
(hurricaneexposuredata
) that is available through a
drat
repository on GitHub. To use the
hurricaneexposure
package, you will need to install
hurricaneexposuredata
on your computer. You can do that by
adding the drat
archive to the list of repositories your
system will query when adding and updating R packages. Once you do this,
you can install the hurricaneexposuredata
using the
install.packages
function (and later update it using the
update.packages
function):
You will want to have the latest version of the
hurricaneexposuredata
package. If you have previously
installed hurricaneexposuredata
, you may need to re-run the
above code if you update the hurricaneexposure
package, to
update to the latest data.
The hurricaneexposuredata
data package includes data
that characterizes county-level exposure to tropical storms in counties
in the eastern half of the United States between 1988 and 2015 (for some
hazards, exposure data is only included for a subset of these years).
Tropical storms that did not pass within at least 250 km of at least one
US county were excluded from these datasets. The following datasets are
included with the hurricaneexposuredata
data package:
county_centers
: Location of United States county
centers of populationhurr_tracks
: Storm tracks for Atlantic-basin storms,
1988-2018closest_dist
: Closest distances between counties and a
storm track, for Atlantic-basin storms, 1988-2018rain
: Rainfall for US counties during Atlantic basin
tropical storms, 1988-2011; daily rainfall is given from five days
before to three days after the storm’s closest approach to the
countystorm_winds
: Modeled county wind speeds for
Atlantic-basin storms, 1988-2018storm_events
: Listings from the NOAA Storm Events
database that occurred near in time and location to tropical storms,
1988-2018. This database changed the types of events it reported in
1996, which should be considered when using the data.ext_tracks_wind
: Estimated county wind speeds for
Atlantic-basin storms, 1988-2018, based on the wind radii listed in the
HURDAT2Once you’ve installed and loaded hurricaneexposuredata
,
you can load the included data using the data
function. For
example, you can access the data on hurricane tracks
(hurr_tracks
) using:
## # A tibble: 6 × 6
## storm_id usa_atcf_id date latitude longitude wind
## <chr> <chr> <chr> <dbl> <dbl> <dbl>
## 1 Alberto-1988 AL011988 198808051800 32 -77.5 20
## 2 Alberto-1988 AL011988 198808060000 32.8 -76.2 20
## 3 Alberto-1988 AL011988 198808060600 34 -75.2 20
## 4 Alberto-1988 AL011988 198808061200 35.2 -74.6 25
## 5 Alberto-1988 AL011988 198808061800 37 -73.5 25
## 6 Alberto-1988 AL011988 198808070000 38.7 -72.4 25
For each dataset included in hurricaneexposuredata
, you
can see the helpfiles for the data for more information (e.g.,
?hurr_tracks
to read documentation on the hurricane tracks
data shown above). This data documentation includes both explanations of
how the dataset was created, definitions of the variables included in
each dataset, and the sources of data included.
The following table shows the storms covered by the
hurricaneexposuredata
data package, as well as which hazard
metrics are available for each year. All storms passed within 250 km of
at least one U.S. county.
Year | Storms | Distance | Wind | Rain | Flood | Tornado |
---|---|---|---|---|---|---|
1988 | Alberto, Beryl, Chris, Florence, Gilbert, Keith, AL13, AL14, AL17 | x | x | x | x | |
1989 | Allison, Chantal, Hugo, Jerry | x | x | x | x | |
1990 | AL01, Bertha, Marco | x | x | x | x | |
1991 | Ana, Bob, Fabian, AL12 | x | x | x | x | |
1992 | AL02, Andrew, Danielle, Earl | x | x | x | x | |
1993 | AL01, Arlene, Emily | x | x | x | x | |
1994 | Alberto, AL02, Beryl, Gordon | x | x | x | x | |
1995 | Allison, Dean, Erin, Gabrielle, Jerry, Opal | x | x | x | x | |
1996 | Arthur, Bertha, Edouard, Fran, Josephine | x | x | x | x | x |
1997 | AL01, Ana, Danny | x | x | x | x | x |
1998 | Bonnie, Charley, Earl, Frances, Georges, Hermine, Mitch | x | x | x | x | x |
1999 | Bret, Dennis, AL07, Floyd, Harvey, Irene | x | x | x | x | x |
2000 | AL04, Beryl, AL09, Gordon, Helene, Leslie | x | x | x | x | x |
2001 | Allison, Barry, Gabrielle, Karen, Michelle | x | x | x | x | x |
2002 | Arthur, Bertha, Cristobal, Edouard, Fay, Gustav, Hanna, Isidore, Kyle, Lili | x | x | x | x | x |
2003 | Bill, Claudette, AL07, Erika, Grace, Henri, Isabel | x | x | x | x | x |
2004 | Alex, Bonnie, Charley, Frances, Gaston, Hermine, Ivan, Jeanne, Matthew | x | x | x | x | x |
2005 | Arlene, Cindy, Dennis, Emily, Katrina, Ophelia, Rita, Tammy, Twenty-Two, Wilma | x | x | x | x | x |
2006 | Alberto, Beryl, Chris, Ernesto | x | x | x | x | x |
2007 | Andrea, Barry, Erin, Gabrielle, Humberto, Ten, Noel | x | x | x | x | x |
2008 | Cristobal, Dolly, Edouard, Fay, Gustav, Hanna, Ike, Kyle, Paloma | x | x | x | x | x |
2009 | One, Claudette, Ida | x | x | x | x | x |
2010 | Alex, Two, Bonnie, Five, Earl, Hermine, Nicole, Paula | x | x | x | x | x |
2011 | Bret, Don, Emily, Irene, Lee | x | x | x | x | x |
2012 | Alberto, Beryl, Debby, Isaac, Sandy | x | x | x | x | |
2013 | Andrea, Dorian, Karen | x | x | x | x | |
2014 | Arthur | x | x | x | x | |
2015 | Ana, Bill, Claudette | x | x | x | x | |
2016 | Bonnie, Colin, Eight, Hermine, Julia, Matthew | x | x | x | x | |
2017 | Cindy, Emily, Harvey, Irma, Jose, Nate, Philippe | x | x | x | x | |
2018 | Alberto, Chris, Florence, Gordon, Michael | x | x | x | x |
In functions throughout the hurricaneexposure
package,
storms are identified based on their name and year (e.g., Hurricane
Floyd in 1999 is identified as “Floyd-1999”). Note that a few storms
satisfied the criteria to be included in the data but were unnamed
(e.g., “AL13” in 1988). These storms are identified in functions in the
hurricaneexposure
package based on the identifying name
listed in the above table, which is based on the first four characters
of the storm’s US ATCF code for the storm. For some research projects,
however, you may want to consider excluding these unnamed storms from
the analysis.
The hurricaneexposure
package has functions to interact
with the data in the hurricaneexposuredata
package. First,
there are functions that can be used to create some different maps of
hurricane exposures based on distance to the storm track and
rainfall.
For several of the storm hazard measurements (currently rain, wind,
and distance), you can use the map_counties
function in
hurricaneexposure
to plot a map for a specific storm
showing the continuous value of that metric for all eastern U.S.
counties. For example, you can map county-level rain during Hurricane
Floyd (1999) using the following call:
## Warning: `tbl_df()` was deprecated in dplyr 1.0.0.
## ℹ Please use `tibble::as_tibble()` instead.
## ℹ The deprecated feature was likely used in the hurricaneexposure package.
## Please report the issue at
## <https://github.com/geanders/hurricaneexposure/issues>.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning: `mutate_()` was deprecated in dplyr 0.7.0.
## ℹ Please use `mutate()` instead.
## ℹ See vignette('programming') for more help
## ℹ The deprecated feature was likely used in the hurricaneexposure package.
## Please report the issue at
## <https://github.com/geanders/hurricaneexposure/issues>.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
By default, this map shows the cumulative rain in each county for two
days before to one day after the date that the storm passed closest to
the county (see the help file for the closest_dist
dataset
in hurricaneexposuredata
to find out more about how this
date of closest approach was determined for each county-storm
combination). The period used to calculate cumulative rainfall for the
map can be customized using the days_included
option in the
map_counties
function. The two following examples show the
difference in rain maps for Tropical Storm Allison (2001) when using
rain only from the day before and the day the storm was closest to each
county (days_included = -1:0
) versus from five days before
to three days after the storm’s closest approach
(days_included = -5:3
):
map_counties(storm = "Allison-2001", metric = "rainfall", days_included = -1:0) +
ggplot2::ggtitle("Rain during Allison (2001) for day before and day of closest approach")
map_counties(storm = "Allison-2001", metric = "rainfall", days_included = -5:3) +
ggplot2::ggtitle("Rain during Allison (2001) for five days before to three days\nafter the day of closest approach")
You can also use the map_counties
function to plot the
maximum wind during the storm for each county. For this, you use the
argument metric = "wind"
. For example, you can plot
county-specific storm winds during Hurricane Katrina (2005) with the
call:
These county winds are determined based on a model of wind speeds,
using as input data from the hurricane tracks data. See the
documentation for the storm_winds
dataset
(?storm_winds
) to find out more about this data; further
details on the wind model are available in the vignettes for the stormwindmodel
package.
You can map a few other variables of wind with this function. For
example, you can map the duration of winds of 20 m / s or more by using
the argument wind_var = "sust_dur"
:
You can map estimated gust winds, rather than sustained winds (the
default), using the argument wind_var = "vmax_gust"
.
Further, we have included a second source of estimated winds in the data
available in hurricaneexposuredata
. These wind estimates
are based on the wind radii in the HURDAT2 dataset (you can find out
more about this data in the helpfile for the
ext_tracks_wind
dataset). These data provide estimates of
which counties were exposed during a storm to sustained winds in four
categories: 0–34 knots; 34–50 knots; 50–64 knots; and 64 knots or
higher. Therefore, these data provide a categorical rather than
continuous estimate of county wind speeds. However, they may be
preferable to the model wind speeds for some storms, especially storms
in extra-tropical transition. Further, you may find it interesting to
use these data in a sensitivity analysis, to compare if impact study
results are sensitive to whether these wind data or the modeled wind
data are used for exposure classification. You can create a wind map
based on this data using the wind_source = "ext_tracks"
argument in the map_counties
function. For example, to map
estimated winds during Hurricane Katrina based on the HURDAT2 wind
radii, you can run:
Finally, you can also use the map_counties
function to
plot the closest distance between the storm and each county. For this,
you use the argument metric = "distance"
. For example, the
following call plots county-level distances from the track of Hurricane
Sandy (2012):
See the documentation for the closest_dist
dataset
(?closest_dist
) for more details on how distances from the
storm track were calculated for each county.
You can also map binary county-level exposures, where each county is classified as “exposed” based on whether it meets a threshold for one of the hazard-based metrics. These maps can be created for the metrics shown in the continuous exposure maps above (wind, rain, and distance from the storm track). In addition, binary classification maps can also be created based on storm event listings for flood events and tornado events from the NOAA Storm Events database.
For example, you can map a binary variable of distance-based exposure
using map_distance_exposure
, with “exposure” defined as
that the county was within the number of kilometers given in
dist_limit
of the storm track. For example, you can map the
counties that were within 75 km of Hurricane Sandy’s storm track with
the call:
Similarly, you can map binary county exposure based on rain using
map_rain_exposure
. In the case of rain, you must define
exposure thresholds for both rainfall and distance. This is because the
rain data is observed data, so neglecting to limit “exposed” counties to
those close to the storm track may pick up counties far from the storm
that experienced a lot of rain over the same period. For example, you
can use the following call to map rain exposure to Tropical Storm
Allison (2001), where a county is defined as exposed to the storm if
cumulative rainfall from five days before to three days after the
storm’s closest approach was 175 mm or more and the county was within
500 km of the storm’s track:
As with the continuous maps of rain exposure, the number of days
included to calculate cumulative rain can be adjusted with the
days_included
option in map_rain_exposure
(see
the help file for map_rain_exposure
for more details).
You can use the map_wind_exposure
function to map
counties exposed to a certain winds during a storm. For example, to
identify all counties exposed to sustained winds of 20 m / s or higher
during Hurricane Katrina, you can run:
If you would like to base a wind threshold on knots rather than m /
s, you can use the convert_wind_speed
function from the
weathermetrics
package to do that. For example, to map
counties in which winds were 34 knots or higher during the storm, you
can run:
library(weathermetrics)
map_wind_exposure(storm = "Katrina-2005",
wind_limit = convert_wind_speed(34, "knots", "mps"))
Finally, you can map which counties were exposed to specific types of
events, as listed in the NOAA Storm Events Database, using the
map_event_exposure
function. For example, you can map
counties for which a flood event was listed for Hurricane Floyd (1999)
with the call:
Similarly, you can map which counties were exposed to tornado events during Hurricane Ivan (2004) with:
When using this function, note that the types of events reported in the NOAA Storm Events Database changed in 1996, so for many event types there will be no listings before 1996. Tornado listings do extend earlier, but flood listings do not. Further, when using this data, the user should be aware of and take into account the limitations of the database, especially for comparing storms across years. Reporting standards likely have changed over time; for example, it is possible that earlier years have more false negatives than more recent years.
Event listings from the NOAA Storm Events database were linked to
storms based on the start date of the event being within a five-day
window of the date the storm was closest to the county and the storm
coming within 500 km of the county. The noaastormevents
package was used to create this dataset. For more information on the
data, see the help file for the storm_events
dataset
(?storm_events
) or the vignette for the
noaastormevents
package.
If you would like to add a storm’s track to one of the maps of
county-level exposure, you can do so using the map_tracks
function. This function can be used either to map a storm track on a
blank U.S. map or to add a track to an existing map, if the map is a
ggplot
object (as are all maps generated in the previous
section).
First, the map_tracks
function can be used to create a
map of the hurricane tracks for one or more storms. For example, to plot
the tracks of Hurricane Floyd in 1999, you can run:
There are some different options you can use for the tracks’ appearance. For example, if you wanted to plot the tracks of several storms, and plot each point when the track locations were measured (typically every six hours), use some transparency so you can see all the tracks, and show the tracks in blue, you can run:
map_tracks(storms = c("Andrew-1992", "Katrina-2005", "Rita-2005"),
alpha = 0.5, plot_points = TRUE, color = "blue")
As another example, to map all tracks for storms in 2018, you can run:
library(dplyr)
library(tidyr)
storms_2018 <- hurr_tracks %>%
select(storm_id) %>%
distinct() %>%
mutate(year = stringr::str_extract(storm_id, "-[0-9].+")) %>%
filter(year == "-2018")
map_tracks(storms = storms_2018$storm_id)
The map_tracks
function is used in all the mapping
functions described in the previous section, to add the storm’s track to
the exposure maps. However, if you would like to customize the
appearance of the storm’s track on the map, you can do so by plotting
the exposure map without the storm track
(add_track = FALSE
) and then adding the track yourself with
the map_tracks
function. For example, to change the color
of the storm track in a map of flood event exposure during Hurricane
Floyd (1999) and to add points showing the available observations for
the storm in the Best Tracks hurricane tracking data, you can run:
The hurricaneexposure
package also has several functions
that can input a list of counties and output a list of all of the storms
to which each county was exposed over a certain period, based on one of
the hazard metrics. These exposure listings include the storm dates,
making them easy to integrate into time series datasets of health or
other impact data.
The county_rain
function takes a list of county FIPS
codes, bounds on the starting and ending years of the analysis, and
thresholds to define rain-based exposure (cumulative rainfall and
distance from the storm’s track, as explained for the binary rain maps
shown in an earlier section) and creates a list of all storms that met
these thresholds for the counties. For example, to get a dataset of all
the storms to which Orleans Parish (FIPS 22071), and Newport News,
Virginia (FIPS 51700), were exposed between 1995 and 2005, where
“exposed” means that the storm passed within 100 kilometers of the
county center and the rainfall over a three-day window of the date of
closest approach was 100 millimeters or more, you can run:
county_rain(counties = c("22071", "51700"), start_year = 1995, end_year = 2005,
rain_limit = 100, dist_limit = 100, days_included = c(-1, 0, 1))
## Key: <storm_id, fips>
## storm_id fips closest_date storm_dist tot_precip local_time
## <char> <char> <char> <num> <num> <char>
## 1: Bill-2003 22071 2003-06-30 47.74045 141.1 2003-06-30 18:15
## 2: Charley-2004 51700 2004-08-14 55.58752 136.2 2004-08-14 19:45
## 3: Cindy-2005 22071 2005-07-06 26.62602 113.2 2005-07-06 02:15
## 4: Floyd-1999 51700 1999-09-16 49.43163 207.5 1999-09-16 11:15
## 5: Isidore-2002 22071 2002-09-26 14.27431 249.0 2002-09-26 06:45
## 6: Katrina-2005 22071 2005-08-29 44.07046 196.2 2005-08-29 09:00
## 7: Matthew-2004 22071 2004-10-10 74.88275 123.2 2004-10-10 10:00
## closest_time_utc
## <char>
## 1: 2003-06-30 23:15
## 2: 2004-08-14 23:45
## 3: 2005-07-06 07:15
## 4: 1999-09-16 15:15
## 5: 2002-09-26 11:45
## 6: 2005-08-29 14:00
## 7: 2004-10-10 15:00
This function draws on the same dataset (rain
) from the
hurricaneexposuredata
package as the functions for mapping
rain exposure shown above. For this and other exposures, see the table
earlier in this vignette to see which years of data are available (for
rain, data is currently available for storms through 2011).
In addition to giving the names and closest dates of each storm for
each county (closest_date
– note, this is given based on
local time for the county; see documentation for the countytimezones
package to find out more about how local time was calculated from
UTC for each county), this function also gives you the distance between
the county and the storm’s track at the time when the storm was closest
to the county’s population weighted center (storm_dist
, in
kilometers) and the total precipitation over the included days
(tot_precip
). The returned dataframe also gives the time of
closest approach in both local time (local_time
) and UTC
(closest_time_utc
), based on 15-minute intervals along the
storm’s track.
Similar functions are available in the package to create listings of county exposures to storms based on wind, distance to the storm track, and events listings from the NOAA Storm Events Database (flood and tornado events). For example, to get a listing of all storms between 1988 and 2015 for which Miami-Dade county (FIPS: ) has experienced sustained winds of 34 knots or more (17.5 m / s), you can run:
## storm_id fips vmax_sust vmax_gust sust_dur gust_dur usa_atcf_id
## <char> <char> <num> <num> <num> <num> <char>
## 1: Andrew-1992 12086 56.60869 84.34695 615 960 AL041992
## 2: Mitch-1998 12086 18.79781 28.00873 0 285 AL131998
## 3: Harvey-1999 12086 23.01816 34.29706 75 285 AL101999
## 4: Irene-1999 12086 31.88419 47.50744 615 1290 AL131999
## 5: Frances-2004 12086 20.77793 30.95911 285 1830 AL062004
## 6: Jeanne-2004 12086 20.53852 30.60239 90 1065 AL112004
## 7: Katrina-2005 12086 33.14848 49.39124 525 1740 AL122005
## 8: Wilma-2005 12086 33.32450 49.65351 315 600 AL252005
## 9: Ernesto-2006 12086 18.02811 26.86188 0 765 AL062006
## 10: Fay-2008 12086 18.85232 28.08996 0 1305 AL062008
## closest_time_utc storm_dist local_time closest_date
## <char> <num> <char> <char>
## 1: 1992-08-24 09:00 30.48032 1992-08-24 05:00 1992-08-24
## 2: 1998-11-05 12:45 130.19217 1998-11-05 07:45 1998-11-05
## 3: 1999-09-21 20:45 48.90460 1999-09-21 16:45 1999-09-21
## 4: 1999-10-15 23:15 43.07201 1999-10-15 19:15 1999-10-15
## 5: 2004-09-05 02:30 156.03624 2004-09-04 22:30 2004-09-04
## 6: 2004-09-26 03:15 156.99903 2004-09-25 23:15 2004-09-25
## 7: 2005-08-26 00:30 12.33865 2005-08-25 20:30 2005-08-25
## 8: 2005-10-24 12:45 76.55817 2005-10-24 08:45 2005-10-24
## 9: 2006-08-30 09:15 69.47981 2006-08-30 05:15 2006-08-30
## 10: 2008-08-19 10:15 127.08899 2008-08-19 06:15 2008-08-19
The returned dataframe includes the estimated maximum values during
the storm for sustained winds (vmax_sust
) and gust winds
(vmax_gust
), both in m / s, as well as the duration for
which sustained and gust winds were over 20 m / s in the county
(sust_dur
and gust_dur
, respectively). It also
includes some of the variables on the date and time of closest approach
and the distance at the storm’s closest approach.
The county_wind
function allows some further
specifications when identifying the storms to which a county was exposed
based on wind estimates. First, it allows the use of alternative wind
variables. While the default is to use the maximum sustained wind speed
experienced in the county during the storm, it is also possible to base
the exposure metric on the maximum gust winds
(wind_var = "vmax_gust"
) or on the duration of either
sustained or gust winds at or above 20 m / s
(wind_var = "sust_dur"
and
wind_var = "gust_dur"
, respectively). For example, to
determine which storms brought Orleans Parish sustained winds of 34
knots or higher (17.5 m / s or higher) for an hour or more
(wind_limit = 60
; for wind durations, this wind limit
should be specified in minutes), you can run:
county_wind(counties = "12086", start_year = 1988, end_year = 2015,
wind_var = "sust_dur", wind_limit = 60)
## storm_id fips vmax_sust vmax_gust sust_dur gust_dur usa_atcf_id
## <char> <char> <num> <num> <num> <num> <char>
## 1: Andrew-1992 12086 56.60869 84.34695 615 960 AL041992
## 2: Harvey-1999 12086 23.01816 34.29706 75 285 AL101999
## 3: Irene-1999 12086 31.88419 47.50744 615 1290 AL131999
## 4: Frances-2004 12086 20.77793 30.95911 285 1830 AL062004
## 5: Jeanne-2004 12086 20.53852 30.60239 90 1065 AL112004
## 6: Katrina-2005 12086 33.14848 49.39124 525 1740 AL122005
## 7: Wilma-2005 12086 33.32450 49.65351 315 600 AL252005
## closest_time_utc storm_dist local_time closest_date
## <char> <num> <char> <char>
## 1: 1992-08-24 09:00 30.48032 1992-08-24 05:00 1992-08-24
## 2: 1999-09-21 20:45 48.90460 1999-09-21 16:45 1999-09-21
## 3: 1999-10-15 23:15 43.07201 1999-10-15 19:15 1999-10-15
## 4: 2004-09-05 02:30 156.03624 2004-09-04 22:30 2004-09-04
## 5: 2004-09-26 03:15 156.99903 2004-09-25 23:15 2004-09-25
## 6: 2005-08-26 00:30 12.33865 2005-08-25 20:30 2005-08-25
## 7: 2005-10-24 12:45 76.55817 2005-10-24 08:45 2005-10-24
Further, the county_wind
function allows you to pull
wind estimates from an alternate source. By default, this function uses
wind estimates for each county from a wind model (see the vignettes for
the stormwindmodel
package for much more detail on this
modeling process). However, you can also pull estimates based on the
wind radii from HURDAT2 by using the option
wind_source = "ext_tracks"
. See the help file for the
ext_tracks_wind
dataset that comes with
hurricaneexposuredata
for more details on this data source
and how the county-specific exposures were determined from the data
source for use in the county_wind
function. It is important
to note that this data results in wind estimates with breaks at 34
knots, 50 knots, and 64 knots, rather than continuous estimates. This
means that a sustained wind estimate from this data source of 17.4896 m
/ s (34 knots) is estimating that the county had maximum wind speeds of
34 knots or higher during the storm, but not as high as 50 knots.
Further, the dur_sust
value when using this data source is
based on number of minutes with winds at or above 34 knots (rather than
the 20 m / s value used for durations for the modeled wind data). Gust
durations can not be determined using the HURDAT2 wind radii.
The HURDAT2 wind radii are more tied to observations during a
specific storm than the modeled data. In certain cases, especially
storms in extratropical transition, this dataset might be preferable to
the modeled wind speeds, even though it provides less continuous
estimates. This source of wind data can be specified using the
wind_source
option in the county_wind
function. For example, to generate a list of all storms in Orleans
Parish with maximum sustained winds of 34 knots or more based on the
HURDAT2 wind radii, you can run:
county_wind(counties = "12086", start_year = 1988, end_year = 2015,
wind_var = "vmax_sust", wind_limit = 17.4, wind_source = "ext_tracks")
## storm_id fips vmax_sust vmax_gust sust_dur usa_atcf_id
## <char> <char> <num> <num> <num> <char>
## 1: Bonnie-2010 12086 17.4896 26.05950 180 AL032010
## 2: Ernesto-2006 12086 17.4896 26.05950 855 AL062006
## 3: Fay-2008 12086 17.4896 26.05950 1230 AL062008
## 4: Frances-2004 12086 17.4896 26.05950 1455 AL062004
## 5: Isaac-2012 12086 17.4896 26.05950 720 AL092012
## 6: Jeanne-2004 12086 17.4896 26.05950 570 AL112004
## 7: Katrina-2005 12086 32.9216 49.05318 720 AL122005
## 8: Wilma-2005 12086 32.9216 49.05318 795 AL252005
## closest_time_utc storm_dist local_time closest_date
## <char> <num> <char> <char>
## 1: 2010-07-23 15:30 27.25331 2010-07-23 11:30 2010-07-23
## 2: 2006-08-30 09:15 69.47981 2006-08-30 05:15 2006-08-30
## 3: 2008-08-19 10:15 127.08899 2008-08-19 06:15 2008-08-19
## 4: 2004-09-05 02:30 156.03624 2004-09-04 22:30 2004-09-04
## 5: 2012-08-26 15:15 248.34174 2012-08-26 11:15 2012-08-26
## 6: 2004-09-26 03:15 156.99903 2004-09-25 23:15 2004-09-25
## 7: 2005-08-26 00:30 12.33865 2005-08-25 20:30 2005-08-25
## 8: 2005-10-24 12:45 76.55817 2005-10-24 08:45 2005-10-24
The county_distance
function can similarly be used to
generate a listing of all storms that came within a certain distance of
specified counties. For example, to get a list of all storms that came
within 50 km of Orleans Parish (FIPS 22071) between 1988 and 2015, you
can run:
## storm_id fips closest_date storm_dist local_time
## <char> <char> <char> <num> <char>
## 1: Beryl-1988 22071 1988-08-09 3.379783 1988-08-09 04:45
## 2: Florence-1988 22071 1988-09-10 4.217016 1988-09-10 03:00
## 3: Allison-2001 22071 2001-06-11 23.289781 2001-06-11 02:30
## 4: Bertha-2002 22071 2002-08-05 16.190337 2002-08-05 03:45
## 5: Isidore-2002 22071 2002-09-26 14.274314 2002-09-26 06:45
## 6: Bill-2003 22071 2003-06-30 47.740448 2003-06-30 18:15
## 7: Cindy-2005 22071 2005-07-06 26.626022 2005-07-06 02:15
## 8: Katrina-2005 22071 2005-08-29 44.070459 2005-08-29 09:00
## 9: Bonnie-2010 22071 2010-07-25 40.631881 2010-07-25 05:30
## closest_time_utc
## <char>
## 1: 1988-08-09 09:45
## 2: 1988-09-10 08:00
## 3: 2001-06-11 07:30
## 4: 2002-08-05 08:45
## 5: 2002-09-26 11:45
## 6: 2003-06-30 23:15
## 7: 2005-07-06 07:15
## 8: 2005-08-29 14:00
## 9: 2010-07-25 10:30
Finally, you can use the county_events
function to get a
listing of all storms for which a county had a certain type of NOAA
Storm Events listing. For example, to get a list of all storms for which
Norfolk, Virginia, (FIPS: 51710) had a flood event listing, you can
run:
## fips storm_id usa_atcf_id closest_time_utc storm_dist
## 1 51710 Arthur-1996 AL011996 1996-06-20 08:00 164.49133
## 2 51710 AL01-1997 AL011997 1997-06-01 18:30 246.23488
## 3 51710 Dennis-1999 AL051999 1999-09-05 05:00 200.44453
## 4 51710 Floyd-1999 AL081999 1999-09-16 15:00 20.91885
## 5 51710 Irene-1999 AL131999 1999-10-18 07:30 234.49601
## 6 51710 Kyle-2002 AL122002 2002-10-12 08:15 149.04476
## 7 51710 Alex-2004 AL012004 2004-08-03 18:45 199.38999
## 8 51710 Alberto-2006 AL012006 2006-06-14 20:15 105.73478
## 9 51710 Ernesto-2006 AL062006 2006-09-01 23:45 69.60457
## 10 51710 Irene-2011 AL092011 2011-08-28 00:00 54.49080
## 11 51710 Sandy-2012 AL182012 2012-10-29 20:30 289.06525
## 12 51710 Claudette-2015 AL032015 2015-07-12 00:00 182.77412
## local_time closest_date
## 1 1996-06-20 04:00 1996-06-20
## 2 1997-06-01 14:30 1997-06-01
## 3 1999-09-05 01:00 1999-09-05
## 4 1999-09-16 11:00 1999-09-16
## 5 1999-10-18 03:30 1999-10-18
## 6 2002-10-12 04:15 2002-10-12
## 7 2004-08-03 14:45 2004-08-03
## 8 2006-06-14 16:15 2006-06-14
## 9 2006-09-01 19:45 2006-09-01
## 10 2011-08-27 20:00 2011-08-27
## 11 2012-10-29 16:30 2012-10-29
## 12 2015-07-11 20:00 2015-07-11
Some datasets covering storm impacts might have communities composed
of multiple counties rather than county-specific listings. To get a
dataframe listing the relevant storms for multi-county communities, you
can use the multi_county_*
family of functions in a similar
way to the county_*
family of functions. For example, to
get rain exposure listings in the case where one of the communities (New
York, NY) is comprised of multiple counties, you can run:
communities <- data.frame(community_name = c(rep("ny", 6), "no", "new"),
fips = c("36005", "36047", "36061",
"36085", "36081", "36119",
"22071", "51700"))
multi_county_rain(communities = communities, start_year = 1995, end_year = 2005,
rain_limit = 100, dist_limit = 100)
## # A tibble: 10 × 9
## # Groups: community_name [3]
## community_name storm_id closest_date local_time closest_time_utc mean_dist
## <chr> <chr> <chr> <chr> <chr> <dbl>
## 1 new Charley-20… 2004-08-14 2004-08-1… 2004-08-14 23:45 55.6
## 2 new Danny-1997 1997-07-24 1997-07-2… 1997-07-24 19:00 61.9
## 3 new Floyd-1999 1999-09-16 1999-09-1… 1999-09-16 15:15 49.4
## 4 no Allison-20… 2001-06-11 2001-06-1… 2001-06-11 07:30 23.3
## 5 no Bill-2003 2003-06-30 2003-06-3… 2003-06-30 23:15 47.7
## 6 no Cindy-2005 2005-07-06 2005-07-0… 2005-07-06 07:15 26.6
## 7 no Isidore-20… 2002-09-26 2002-09-2… 2002-09-26 11:45 14.3
## 8 no Katrina-20… 2005-08-29 2005-08-2… 2005-08-29 14:00 44.1
## 9 no Matthew-20… 2004-10-10 2004-10-1… 2004-10-10 15:00 74.9
## 10 ny Floyd-1999 1999-09-16 1999-09-1… 1999-09-17 00:15 41.2
## # ℹ 3 more variables: mean_rain <dbl>, max_rain <dbl>, min_dist <dbl>
The output from this function includes columns for the average
closest distance for any of the counties in the community
(mean_dist
), the average precipitation for all the counties
(mean_precip
), the highest precipitation for any of the
counties (max_rain
), and the smallest distance between the
storm track and any of the county population-weighted centers
(min_dist
).
Similar functions exist for wind exposure
(multi_county_wind
) and distance exposure
(multi_county_distance
). In some cases, these multi-county
functions may have less functionality than their county_*
counterparts. For example, the multi_county_wind
function
currently only allows use of the sustained wind speeds for the exposure
metric, while the county_wind
function allows the use of
sustained winds, gust winds, or durations of either.