Investigation of HWIND Simulation Wind Speeds and a Methodology for Simulating Distributions of Extreme Winds in Hurricane Environments
by
Joseph B. Dannemiller
A Dissertation
In
Wind Science & Engineering
Submitted to the Graduate Faculty of Texas Tech University in Partial Fulfillment of the Requirements for the Degree of
DOCTOR OF PHILOSOPHY
Approved
Douglas A. Smith Co-Chair of Committee
Stephen M. Morse Co-Chair of Committee
John L. Schroeder
Kishor C. Mehta
Mark Sheridan Dean of the Graduate School
May 2019
Copyright 2018, Joseph Dannemiller
ii Texas Tech University, Joseph Dannemiller, May 2019
ACKNOWLEDGEMENTS My thanks to Dr. Douglas Smith and Dr. Stephen Morse for chairing my committee, and Drs. John Schroeder and Kishor Mehta for their help and assistance. Also, thank you to Dr. Audra Morse, Dr. Anna Young, and Dr. Rich Krupar for their guidance, friendship and support in matters pertaining to wind, research and life. I would also like to thank all of the people associated with the Wind Science & Engineering Research Center, and the National Wind Institute, for their help and assistance along the way. Having been an IGERT Fellow I would like to thank the National Science Foundation for the IGERT program for funding me during this process and allowing me the freedom to pursue this research. I would like to thank, most of all, my wife Sandra, my children Mark, Alexandra and James, my family including Mom, Dad, Katherine and Christopher, and all of my friends for supporting me because without their love and support I would never have found the fortitude to complete this journey. During this process I have met, interacted, and made friends with so many people that I am not mentioning here specifically, so to everyone that helped me, in any way, thank you.
iii Texas Tech University, Joseph Dannemiller, May 2019
TABLE OF CONTENTS ACKNOWLEDGEMENTS ...... iii ABSTRACT ...... xiii LIST OF TABLES ...... xiv LIST OF FIGURES ...... xvi DEFINITION OF VARIABLES ...... xxxii I. INTRODUCTION ...... 1 II. BACKGROUND ...... 3 Major Hurricane Damage ...... 3 Loads and Resistances ...... 6 Wind Loads ...... 9 Observational Systems ...... 17 Texas Tech University Hurricane Research Team ...... 19 Wind Speed Standardization ...... 28 Gust Factors ...... 33 Density of Surface Observation Deployments ...... 51 HWIND Hurricane Wind Field Model ...... 52 HWIND in SPA ...... 58 III. HWIND AND TTUHRT DATA PROCESSING ...... 60 Data Selection ...... 60 HWIND Data Processing ...... 60 Texas Tech University Hurricane Research Team Data Processing Methods ...... 71 Processed Texas Tech University Hurricane Research Team Data ...... 85 IV. TTUHRT VERSUS HWIND COMPARISON ...... 95 Differences Between TTUHRT and HWIND as a Function of Distance From Storm Center ...... 102 Comparison of TTUHRT and HWIND Gust Factors ...... 109
iv Texas Tech University, Joseph Dannemiller, May 2019
Comparison of TTUHRT and HWIND Mean Values ...... 113 Comparison of Recorded Peak Wind Speeds and Maximum Sustained Wind Speeds from HWIND Reconstruction ...... 115 V. ESTIMATING PARAMETERS FOR DISTRIBUTIONS OF EXTREME WIND SPEED IN 600s WINDOWS USING TTUHRT DATA ...... 118 Distribution of Extreme Winds Location and Scale Parameters ...... 118 Parent Wind Distribution Mean Parameters Versus Extreme Wind Distribution Location Parameters...... 122 Estimating Scale Parameters for Distributions of Extreme Winds ...... 130 Extreme Value Location Parameters Separated into Smoother than Open, Open and Rougher than Open Exposure Categories ...... 142 Extreme Value Scale Parameters Separated into Smooth, Open and Rough Exposure Categories ...... 167 Estimating Values of Extreme Wind Speed Distribution Location Parameters ...... 176 Estimating Values of Extreme Wind Speed Distribution of Scale Parameters ...... 180 Iterative Procedure to Generate Distributions of Extreme Winds ...... 186 VI. CASE STUDIES ...... 193
Case Study One – Smoother Than Open Exposure (z0<0.03m) ...... 193
Case Study Two –Open Exposure (0.3m≤z0≤0.07m) ...... 209 Case Study Three –Rougher Than Open Exposure (z0>0.7m) ...... 226 Case Study Results ...... 243 VII. FUTURE POSSIBILITIES FOR SIMULATING DISTRIBUTIONS OF EXTREME WIND SPEEDS ...... 245 VIII. CONCLUSIONS ...... 253 Future Work ...... 256 BIBLIOGRAPHY ...... 258
v Texas Tech University, Joseph Dannemiller, May 2019
APPENDICES ...... 263 A1 HWIND ASSIGNED SURFACE ROUGHNESS VALUES FOR 30° SECTORS AT ALL TTHURT DEPLOYMENTS ...... 263 A2 HWIND SIMULATION OUTPUT DATA FIGURES ...... 295 A3 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 1998 LANDFALL OF HURRICANE BONNIE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 349 A4 TEXAS TECH UNIVERSITY HURRICANE RESEARCH TEAM DEPLOYMENT SUMMARY STATISTICS ...... 455 A5 HWIND SIMULATION OUTPUT DATA TABLES ...... 480 A6 TTUHRT VERSUS HWIND COMPARISON PLOTS ...... 514 A7 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR ALL 600S WINDOWS ...... 539 A8 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING 3s MOVING AVERAGE DATA FOR ALL 600S WINDOWS ...... 550 A9 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING 60s MOVING AVERAGE DATA FOR ALL 600S WINDOWS ...... 561 A10 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR 600S WINDOWS WITH A MEAN ABOVE 15m/s ...... 572 A11 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING 3s MOVING AVERAGE DATA FOR 600s WINDOWS WITH A MEAN ABOVE 15m/s ...... 583 A12 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE PARAMETERS USING 60s MOVING AVERAGE DATA FOR 600s WINDOWS WITH A MEAN ABOVE 15m/s ...... 594 A13 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 LESS THAN 0.03m ...... 605
vi Texas Tech University, Joseph Dannemiller, May 2019
A14 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 BETWEEN 0.03m AND 0.07m ...... 616 A15 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 GREATER THAN 0.07m ...... 624 A16 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE LOCATION PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS SEPARATED INTO SMOOTH (z0<0.03M) OPEN (0.03M
vii Texas Tech University, Joseph Dannemiller, May 2019
A23 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE PARAMETERS USING 3s MOVING AVERAGE DATA FOR 600s WINDOWS WITH A MEAN ABOVE 15m/s ...... 728 A24 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE PARAMETERS USING 60s MOVING AVERAGE DATA FOR 600s WINDOWS WITH A MEAN ABOVE 15m/s ...... 743 A25 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE SCALE PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 LESS THAN 0.03m ...... 755 A26 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE SCALE PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 BETWEEN 0.03m AND 0.07m ...... 770 A27 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE SCALE PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS WITH A SURFACE ROUGHNESS VALUE z0 GREATER THAN 0.07m ...... 787 A28 CONDITIONAL DISTRIBUTIONS OF EXTREME VALUE PARAMETERS USING RAW DATA FOR ALL 600s WINDOWS SEPARATED INTO SMOOTH (z0<0.03M) OPEN (0.03M
viii Texas Tech University, Joseph Dannemiller, May 2019
A32 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #2 DURING THE 2002 LANDFALL OF HURRICANE LILI BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 947 A33 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2003 LANDFALL OF HURRICANE ISABEL BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1047 A34 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #2 DURING THE 2003 LANDFALL OF HURRICANE ISABEL BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1190 A35 SUMMARY STATISTICS RECORDED BY TTUHRT PMT BLACK DURING THE 2003 LANDFALL OF HURRICANE ISABEL BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1385 A36 SUMMARY STATISTICS RECORDED BY TTUHRT PMT CLEAR DURING THE 2003 LANDFALL OF HURRICANE ISABEL BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1519 A37 SUMMARY STATISTICS RECORDED BY TTUHRT PMT WHITE DURING THE 2003 LANDFALL OF HURRICANE ISABEL BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1660 A38 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2004 LANDFALL OF HURRICANE FRANCES BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 1801 A39 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #2 DURING THE 2004 LANDFALL OF HURRICANE FRANCES BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 2090 A40 SUMMARY STATISTICS RECORDED BY TTUHRT PMT BLACK DURING THE 2004 LANDFALL OF HURRICANE FRANCES BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 2567 A41 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2005 LANDFALL OF HURRICANE DENNIS BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 2769 A42 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2005 LANDFALL OF HURRICANE IVAN BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 2877
ix Texas Tech University, Joseph Dannemiller, May 2019
A43 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2005 LANDFALL OF HURRICANE KATRINA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 2930 A44 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #2 DURING THE 2005 LANDFALL OF HURRICANE KATRINA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3003 A45 SUMMARY STATISTICS RECORDED BY TTUHRT PMT BLACK DURING THE 2005 LANDFALL OF HURRICANE KATRINA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3190 A46 SUMMARY STATISTICS RECORDED BY TTUHRT PMT CLEAR DURING THE 2005 LANDFALL OF HURRICANE KATRINA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3311 A47 SUMMARY STATISTICS RECORDED BY TTUHRT PMT WHITE DURING THE 2005 LANDFALL OF HURRICANE KATRINA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3454 A48 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #1 DURING THE 2005 LANDFALL OF HURRICANE RITA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3711 A49 SUMMARY STATISTICS RECORDED BY TTUHRT WEMITE #2 DURING THE 2005 LANDFALL OF HURRICANE RITA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 3877 A50 SUMMARY STATISTICS RECORDED BY TTUHRT PMT BLACK DURING THE 2005 LANDFALL OF HURRICANE RITA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4026 A51 SUMMARY STATISTICS RECORDED BY TTUHRT PMT CLEAR DURING THE 2005 LANDFALL OF HURRICANE RITA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4046 A52 SUMMARY STATISTICS RECORDED BY TTUHRT PMT WHITE DURING THE 2005 LANDFALL OF HURRICANE RITA BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4260 A53 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 101A DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4332
x Texas Tech University, Joseph Dannemiller, May 2019
A54 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 102B DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4385 A55 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 103A DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4386 A56 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 104B DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4424 A57 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 105A DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4522 A58 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 106B DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4676 A59 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 107A DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4758 A60 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 108B DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4845 A61 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 109A DURING THE 2008 LANDFALL OF HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 4929 A62 SUMMARY STATISTICS RECORDED BY TTUHRT STICKNET 110A DURING THE 2008 LANDFALL OF
xi Texas Tech University, Joseph Dannemiller, May 2019
HURRICANE IKE BROKEN INTO SEQUENTIAL 600s WINDOWS ...... 5059
xii Texas Tech University, Joseph Dannemiller, May 2019
ABSTRACT The wind speeds simulated by the National Oceanic and Atmospheric Administration’s HWIND hurricane simulation model are compared to data gathered by Texas Tech University’s Hurricane Research Team. The errors associated with the HWIND reported wind speeds are quantified for data from 32 deployments, using 15 platforms, during the landfall of 10 hurricanes, from 1998 to 2005. The relationships between parent wind distribution parameters and extreme wind field parameters is explored and a methodology for simulating extreme wind speed distributions in hurricane environments is advanced. Utilizing the proposed methodology, hurricane wind field model data can be used to generate distributions of extreme winds for direct use in assessing structural performance during hurricane events.
xiii Texas Tech University, Joseph Dannemiller, May 2019
LIST OF TABLES
2.1 Gust factor computed in (Durst, 1970) ...... 42 2.2 Gust factors computed in (Krayer and Marshal, 1992) ...... 43 2.3 Gust factors, Cg3600,2, computed in (Vickery and Skerlj, 2005) ...... 45 2.4 Gust factors, Cg600,60, computed using (Vickery and Skerlj, 2005) data ...... 47 2.5 Gust factors, Cg600,2, computed using (Vickery and Skerlj, 2005) data ...... 48 2.6 Percent errors between simulated and observed peak wind speeds as reported by (Vickery et al., 2000) ...... 58 3.1 HWIND reconstruction data at the location of the TTUHRT WEMITE #1 during the 2005 landfall of Hurricane Katrina ...... 70 3.2 Breakdown of 600s windows from TTUHRT time histories ...... 82 3.3 Data for assessing the surface roughness of the upwind wind field for the TTUHRT StickNet 101A deployment during the 2008 landfall of Hurricane Ike ...... 84 4.1 Magnitude and percent differences between TTUHRT and HWIND maximum sustained wind speeds broken down by year...... 108 5.1 Frequency Table for μj60s, 600s vs σj60s, 600s histogram ...... 137 5.2 Expected Frequency Table for μj60s, 600s vs σj60s, 600s histogram ...... 139 5.3 Chi-Squared values for μj60s, 600s vs σj60s, 600s histogram ...... 140 5.4 GEV parameters for conditional location parameter distributions computed from TTUHRT Raw Data ...... 177 5.5 GEV parameters for conditional location parameter distributions computed from TTUHRT 60s MA Data...... 178 5.6 GEV parameters for conditional location parameter distributions computed from TTUHRT 3s MA Data ...... 179 5.7 GEV parameters for conditional scale parameter distributions computed from TTUHRT Raw data ...... 181 5.8 GEV parameters for conditional scale parameter distributions computed from TTUHRT 3s Data ...... 183
xiv Texas Tech University, Joseph Dannemiller, May 2019
5.9 GEV parameters for conditional scale parameter distributions computed from TTUHRT 60s MA Data...... 185 5.10 Mean and standard deviation of differences between the location and scale parameters recorded by TTUHRT raw data and the means of simulated distributions using only 600s windows with smoother than open exposure ...... 188 5.11 Mean and standard deviation of differences between the location and scale parameters recorded by TTUHRT raw data and the means of simulated distributions using only 600s windows with open exposure ...... 188 5.12 Mean and standard deviation of differences between the location and scale parameters recorded by TTUHRT raw data and the means of simulated distributions using only 600s windows with rougher than open exposure ...... 188 6.1 Assigned surface roughness values, surface roughness values upwind and downwind of the last roughness change, and distances to the last change in surface roughness at the TTUHRT WEMITE #1 deployment location during the 1998 landfall of Hurricane Bonnie in each 30 sector ...... 195 6.2 Assigned surface roughness values, surface roughness values upwind and downwind of the last roughness change, and distances to the last change in surface roughness at the TTUHRT PMT Black deployment location during the 2004 landfall of Hurricane Francis in each 30 sector ...... 211 6.3 Assigned surface roughness values, surface roughness values upwind and downwind of the last roughness change, and distances to the last change in surface roughness at the TTUHRT WEMITE #2 deployment location during the 2003 landfall of Hurricane Isabel in each 30 sector ...... 228
xv Texas Tech University, Joseph Dannemiller, May 2019
LIST OF FIGURES 2.1 Bolivar Peninsula before Hurricane Ike Landfall, September 9, 2008 (FEMA 2009) ...... 5 2.2 Bolivar Peninsula after Hurricane Ike Landfall, September 15, 2008 (FEMA 2009) ...... 5 2.3 Normal distribution probability density function (PDF) ...... 15 2.4 Relationship between the mean and mode for positively skewed distributions ...... 16 2.5 Comparison of normal and skewed distributions ...... 17 2.6 (a) Automated Surface Observation System (ASOS) located in Austin, TX (www.ncdc.noaa.gov, 2017), (b) Weather buoy stationed in the Gulf of Mexico (www.ndbc.noaa.gov, 2017), (c) Geostationary Operational Environmental Satellite R Series (GOES-R) (www.lockheedmartin.com, 2017), (d) National Oceanic and Atmospheric Administration (NOAA) WP-3D Orion weather research and hurricane intercept aircraft (www.esrl.noaa.gov, 2017), (e) National Center for Atmospheric Research (NCAR) dropsonde device in flight (www2.ucar.edu, 2017), (f) ship mounted weather stations ...... 18 2.7 WEMITE #1 near Vacherie, LA prior to the landfall of Hurricane Katrina, 2005 (TTUHRT, 2006) ...... 20 2.8 WEMITE #2 at Stennis International Airport following the landfall of Hurricane Katrina, 2005 (TTUHRT, 2006) ...... 21 2.9 PMT Clear at Slidell Municipal Airport prior to the landfall of Hurricane Katrina, 2005 (TTUHRT, 2006) ...... 22 2.10 StickNet 101, A probe, deployed prior to the landfall of Hurricane Sandy ...... 23 2.11 StickNet 102, B Probe, deployed for testing at Reese Technology Center, Lubbock, TX ...... 24 2.12 Wind speed time history, recorded at 10m by WEMITE #1 during Hurricane Bonnie (1998) ...... 26 2.13 Wind direction time history, recorded at 10m by WEMITE #1 during Hurricane Bonnie (1998) ...... 27 2.14 Hurricane Bonnie, WEMITE #1 deployment site ...... 29 2.15 Comparison of boundary layer velocities over rough and smooth terrain ...... 30
xvi Texas Tech University, Joseph Dannemiller, May 2019
2.16 The boundary layer transition past a change in roughness, rough to smooth ...... 31 2.17 Probability of exceedance relationship in relation to parent wind speed parameters ...... 36 2.18 Comparison of gust factor curves from ESDU (1983), Durst (1960) and Krayer and Marshal (1992), reproduced from (Vickery and Skerlj, 2005)...... 44 2.19 Cg600,2 gust factors computed using TTUHRT platform data where z0 is between 0.03m and 0.07m, reproduced from (Paulsen and Schroeder, 2005) ...... 48 2.20 Cg600,3 gust factors computed using TTUHRT and FCMP platform data, reproduced from (Giammanco et.al., 2012) ...... 49 2.21 Locations of TTUHRT, Florid Coastal Monitoring Project and Louisiana Monroe atmospheric measurement platforms during the 2006 landfall of Hurricane Katrina (Giammanco et al, 2006) ...... 51 2.22 Locations of Automated Surface Observation Stations (ASOS) and Automated Weather Observation Stations (AWOS) ...... 52 2.23 Observations aggregated into the HWIND wind field simulation of Hurricane Katrina at 6:00 UTC on August 29, 2005, NOAA ...... 54 2.24 NOAA HRD created figure showing the HWIND reconstruction for Hurricane Katrina at 06:00 UTC on August 29, 2005 ...... 56 3.1 NOAA HRD created figure showing the HWIND reconstruction for Hurricane Katrina at 00:00 UTC on August 29, 2005 ...... 62 3.2 Grid points employed in HWIND reconstruction ...... 63 3.3 First enlarged view of the grid points employed in HWIND reconstruction...... 64 3.4 Second enlarged view of the grid points employed in HWIND simulation ...... 64 3.5 Wind speeds in m/s from the NOAA HRD HWIND reconstruction of Hurricane Katrina on 00:00UTC August 29, 2005 ...... 65 3.6 Storm relative quadrants in a hurricane ...... 67
xvii Texas Tech University, Joseph Dannemiller, May 2019
3.7 Wind speeds in m/s from the NOAA HRD HWIND reconstruction of Hurricane Katrina on 12:00UTC August 26, 2005 ...... 68 3.8 Wind speed time history for the raw data (blue) and the 60s MA data (red) recorded from 09:15UTC to 09:25UTC by TTUHRT WEMITE #1 during the 1998 landfall of hurricane Bonnie ...... 73 3.9 Peak wind speed values for 60-sec segments ...... 76 3.10 Aerial image of the deployment site for TTUHRT StickNet 101A captured January 2008 (© Google) ...... 83 3.11 Parent and extreme wind distribution parameters for the raw data recorded by TTUHRT WEMITE #1 during the 2005 landfall of Hurricane Katrina ...... 86 3.12 Parent and extreme wind distribution parameters for the 3s MA data recorded by TTUHRT WEMITE #1 during the 2005 landfall of Hurricane Katrina ...... 88 3.13 Parent and extreme wind distribution parameters for the 60s MA data recorded by TTUHRT WEMITE #1 during the 2005 landfall of Hurricane Katrina ...... 89 3.14 Parent and extreme wind distribution parameters for the raw data recorded by TTUHRT PMT Black during the 2005 landfall of Hurricane Katrina ...... 91 3.15 Parent and extreme wind distribution parameters for the 3s MA data recorded by TTUHRT PMT Black during the 2005 landfall of Hurricane Katrina ...... 92 3.16 Parent and extreme wind distribution parameters for the 60s MA data recorded by TTUHRT PMT Black during the 2005 landfall of Hurricane Katrina ...... 93 4.1 TTUHRT and HWIND maximum sustained wind speeds, differences between TTUHRT and HWIND maximum sustained wind speeds, and the distances from platform to center of the storm, for the deployment of TTUHRT PMT Clear during the 2005 landfall of Hurricane Katrina ...... 97 4.2 TTUHRT and HWIND maximum sustained wind speeds, differences between TTUHRT and HWIND mean wind speeds, and the distances from platform to center of the storm, for the deployment of TTUHRT WEMITE #1 during the 1998 landfall of Hurricane Bonnie ...... 99
xviii Texas Tech University, Joseph Dannemiller, May 2019
4.3 TTUHRT and HWIND maximum sustained wind speeds, differences between TTUHRT and HWIND mean wind speeds, and the distances from platform to center of the storm, for the deployment of TTUHRT WEMITE #2 during the 2003 landfall of Hurricane Isabel ...... 100 4.4 TTUHRT and HWIND maximum sustained wind speeds, differences between TTUHRT and HWIND mean wind speeds, and the distances from platform to center of the storm, for the deployment of TTUHRT StickNet 108B during the 2008 landfall of Hurricane Ike ...... 101 4.5 (a) Differences (m/s) between observed TTUHRT xj,k t,Tw and HWIND U 60s, 600s data versus the radial distance from storm center to the location of the TTUHRT platform Rp, (b) %-differences between observed TTUHRT xj, k t, Tw xjraw, 600s and HWIND U 60s, 600s U 600s data versus the radial distance from storm center to the location of the TTUHRT platform, Rp ...... 103 4.6 (a) Distribution of differences between TTUHRT xj, k t, Tw and HWIND U 60s, 600s data, (b) Distribution of %- differences between TTUHRT xj,k t,Tw and HWIND U 60s, 600s data ...... 105 4.7 (a) Differences between TTUHRT xj,k t,Tw and HWIND U 60s, 600s data with respect to storm year, (b) Percent differences between TTUHRT xj,k t,Tw and HWIND U 60s, 600s data with respect to storm year...... 107 4.8 Comparison between TTUHRT gust factors (blue histogram) and the single gust factor employed by HWIND (red) using only TTUHRT gust factors with smoother than open exposure (z0<0.03m) upwind during the recording of 600s time histories ...... 110 4.9 Comparison between TTUHRT (blue histogram) gust factors and the single gust factor employed by HWIND (red line) using only TTUHRT gust factors with open exposure (0.03m≤z0≤0.07m) upwind during the recording of 600s time histories ...... 111 4.10 Comparison between TTUHRT (blue histogram) gust factors and the single gust factor employed by HWIND (red line) using only TTUHRT gust factors with rougher than open exposure (z0>0.07m) upwind during the recording of 600s time histories ...... 112
xix Texas Tech University, Joseph Dannemiller, May 2019
4.11 (a) Distribution of differences between TTUHRT xjraw, 600s and HWIND U 600s data, (b) Distribution of %-differences between TTUHRT xjraw, 600s and HWIND U 600s data ...... 114 4.12 (a) Distribution of differences between TTUHRT xj60s, 600s and HWIND U 60s, 600s , (b) Distribution of %-differences between TTUHRT xj60s, 600s and HWIND U 60s, 600s data ...... 116 5.1 Distribution of EV location parameters from 60s MA windows ...... 119 5.2 Distribution of EV location parameters from 3s MA windows ...... 120 5.3 Distribution of EV scale parameters from 60s MA windows ...... 121 5.4 Distribution of EV scale parameters from 3s MA windows ...... 122 5.5 Linear relationship between mean and location parameters from TTUHRT 60s MA data ...... 123 5.6 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using Equation 5.2 as a function of xj60s, 600s ...... 124 5.7 Percent of differences between the TTUHRT μj60s, 600s and estimates using Equation 5.2 as a function of xj60s, 600s ...... 125 5.8 Linear relationship to estimate μj60s, 600s parameters ...... 126 5.9 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using Equation 5.4 as a function of xjraw, 600s ...... 127 5.10 Percent of differences between the TTUHRT μj60s, 600s and estimates using Equation 5.4 as a function of xjraw, 600s ...... 127 5.11 Linear relationship to estimate μj3s, 600s parameters ...... 128 5.12 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using Equation 5.4 as a function of xjraw, 600s ...... 129 5.13 Percent of differences between the TTUHRT μj3s, 600s and estimates using Equation 5.5 as a function of xjraw, 600s ...... 129 5.14 Relationship between mean and scale parameters from 60s MA windows ...... 131
xx Texas Tech University, Joseph Dannemiller, May 2019
5.15 Relationship between mean and scale parameters from 3s MA windows ...... 131 5.16 Relationship between EV location and scale parameters from 60s MA windows ...... 132 5.17 Relationship between EV location and scale parameters from ...... 132 5.18 Two parameter histogram plotting the correlated frequencies of occurrence for μj60s, 600s and σj60s, 600s values, view 1 ...... 134 5.19 Two parameter histogram plotting the correlated frequencies of occurrence for μj60s, 600s and σj60s, 600s values, view 2 ...... 134 5.20 Two parameter histogram plotting the correlated frequencies of occurrence for μj60s, 600s and σj60s, 600s values, view 3 ...... 135 5.21 Two parameter histogram plotting the correlated frequencies of occurrence for μj60s, 600s and σj60s, 600s values, view 4 ...... 135 5.22 Two parameter histogram plotting the correlated frequencies of occurrence for μj60s, 600s and σj60s, 600s values, heat map...... 136 5.23 Relationship between parent mean and EV location using TTUHRT 60s MA data ...... 143 5.24 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using 5.8 for smoother than open exposure data as a function of xj60s, 600s...... 144 5.25 Percent differences between the TTUHRT μj60s, 600s and estimates using 5.8 for smoother than open exposure data as a function of xj60s, 600s ...... 145 5.26 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using 5.8 for open exposure data as a function of xj60s, 600s ...... 146 5.27 Percent differences between the TTUHRT μj60s, 600s and estimates using 5.8 for open exposure data as a function of xj60s, 600s ...... 147 5.28 Magnitudes of differences between the TTUHRT μj60s, 600s and estimates using 5.8 for rougher than open exposure data as a function of xj60s, 600s...... 148
xxi Texas Tech University, Joseph Dannemiller, May 2019
5.29 Percent differences between the TTUHRT μj60s, 600s and estimates using 5.8 for rougher than open exposure data as a function of xj60s, 600s ...... 149 5.30 Relationship between parent mean and EV location using TTUHRT 3s MA data ...... 150 5.31 Magnitudes of differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for smoother than open exposure data as a function of xj3s, 600s ...... 151 5.32 Percent differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for smoother than open exposure data as a function of xj3s, 600s ...... 152 5.33 Magnitudes of differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for open exposure data as a function of xj3s, 600s ...... 153 5.34 Percent differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for open exposure data as a function of xj3s, 600s ...... 154 5.35 Magnitudes of differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for rougher open exposure data as a function of xj3s, 600s ...... 155 5.36 Percent differences between the TTUHRT μj3s, 600s and estimates using Equation 5.9 for rougher than open exposure data as a function of xj3s, 600s ...... 156 5.37 Relationship between parent mean and EV location using TTUHRT raw data ...... 157 5.38 Magnitudes of differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for smoother than open exposure data as a function of xjraw, 600s ...... 158 5.39 Percent differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for smoother than open exposure data as a function of xjraw, 600s ...... 159 5.40 Magnitudes of differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for open exposure data as a function of xjraw, 600s ...... 160 5.41 Percent differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for open exposure data as a function of xjraw, 600s ...... 161
xxii Texas Tech University, Joseph Dannemiller, May 2019
5.42 Magnitudes of differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for rougher open exposure data as a function of xjraw, 600s ...... 162 5.43 Percent differences between the TTUHRT μjraw, 600s and estimates using Equation 5.10 for rougher than open exposure data as a function of xjraw, 600s ...... 163 5.44 Two parameter histogram plotting the correlated frequencies of occurrence for xjraw, 600s and μjraw, 600s values for multiple roughness regimes, view 1 ...... 164 5.45 Two parameter histogram plotting the correlated frequencies of occurrence for xjraw, 600s and μjraw, 600s values for multiple roughness regimes, view 2 ...... 165 5.46 Two parameter histogram plotting the correlated frequencies of occurrence for xjraw, 600s and μjraw, 600s values for multiple roughness regimes, view 3 ...... 165 5.47 Two parameter histogram plotting the correlated frequencies of occurrence for xjraw, 600s and μjraw, 600s values for multiple roughness regimes, view 4 ...... 166 5.48 Two parameter histogram plotting the correlated frequencies of occurrence for xjraw, 600s and μjraw, 600s values, heat map...... 166 5.49 Relationship between EV location and scale parameters using TTUHRT 60s MA data ...... 168 5.50 Relationship between EV location and scale parameters using TTUHRT 3s MA data ...... 169 5.51 Relationship between EV location and scale parameters using TTUHRT raw data ...... 169 5.52 Two parameter histogram plotting the correlated frequencies of occurrence for μjraw, 600s and σjraw, 600s values for multiple roughness regimes, view 1 ...... 170 5.53 Two parameter histogram plotting the correlated frequencies of occurrence for μjraw, 600s and σjraw, 600s values for multiple roughness regimes, view 2 ...... 171 5.54 Two parameter histogram plotting the correlated frequencies of occurrence for μjraw, 600s and σjraw, 600s values for multiple roughness regimes, view 3 ...... 171
xxiii Texas Tech University, Joseph Dannemiller, May 2019
5.55 Two parameter histogram plotting the correlated frequencies of occurrence for μjraw, 600s and σjraw, 600s values for multiple roughness regimes, view 4 ...... 172 5.56 Two parameter histogram plotting the correlated frequencies of occurrence for μjraw, 600s and σjraw, 600s values, heat map...... 172 5.57 Histograms for smooth, open and rough exposure σjraw, 600s values ...... 173 5.58 Coefficients of variation for the GEV distributions of σjraw, 600s values in all μjraw, 600s bins ...... 174 5.59 Differences in magnitude between the peak 3s wind speed recorded by the TTUHRT platforms and peak 3s wind speeds computed using the mean observed by the TTUHRT platforms and a gust factor equal to 1.44 ...... 189 5.60 Differences in percent between the peak 3s wind speed recorded by the TTUHRT platforms and peak 3s wind speeds computed using the mean observed by the TTUHRT platforms and a gust factor equal to 1.44 ...... 190 5.61 Differences in magnitude between the peak 60s wind speed recorded by the TTUHRT platforms and peak 60s wind speeds computed using the mean observed by the TTUHRT platforms and a gust factor equal to 1.18 ...... 191 5.62 Differences in percent between the peak 60s wind speed recorded by the TTUHRT platforms and peak 60s wind speeds computed using the mean observed by the TTUHRT platforms and a gust factor equal to 1.18 ...... 191 6.1 Deployment location for the TTUHRT WEMITE #1 platform during the 1998 landfall of Hurricane Bonnie ...... 194 6.2 Wind speed time history recorded by the TTUHRT WEMITE #1 platform during the 1998 landfall of Hurricane Bonnie ...... 195 6.3 Simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 197 6.4 Simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 197 6.5 Differences in magnitude between the location parameter recorded by TTUHRT WEMITE #1 and the values in the
xxiv Texas Tech University, Joseph Dannemiller, May 2019
simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 198 6.6 Differences in percent between the location parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 199 6.7 Differences in magnitude between the scale parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 200 6.8 Differences in magnitude between the scale parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 201 6.9 Distribution of mean differences in magnitude between the location parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window...... 202 6.10 Distribution of mean differences in percent between the location parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window...... 203 6.11 Distribution of mean differences in magnitude between the scale parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 204 6.12 Distribution of mean differences in percent between the scale parameter recorded by TTUHRT WEMITE #1 and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Bonnie raw data for a single 600s window ...... 205
xxv Texas Tech University, Joseph Dannemiller, May 2019
6.13 Differences in magnitude between the peak 3s wind speed recorded by the TTUHRT WEMITE #1 during the 1998 landfall of Hurricane Bonnie and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Bonnie time history ...... 207 6.14 Differences in percent between the peak 3s wind speed recorded by the TTUHRT WEMITE #1 during the 1998 landfall of Hurricane Bonnie and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Bonnie time history ...... 207 6.15 Differences in magnitude between the peak 60s wind speed recorded by the TTUHRT WEMITE #1 during the 1998 landfall of Hurricane Bonnie and peak 60s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.18 using all 600s windows in the Hurricane Bonnie time history ...... 208 6.16 Differences in percent between the peak 60s wind speed recorded by the TTUHRT WEMITE #1 during the 1998 landfall of Hurricane Bonnie and peak 60s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.18 using all 600s windows in the Hurricane Bonnie time history ...... 209 6.17 Deployment location for the TTUHRT PMT Black platform during the 2004 landfall of Hurricane Francis ...... 210 6.18 Simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 212 6.19 Simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 213 6.20 Differences in magnitude between the location parameter recorded by TTUHRT PMT Black and the values in the simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 214 6.21 Differences in percent between the location parameter recorded by TTUHRT PMT Black and the values in the simulated distribution of location parameters for the smooth
xxvi Texas Tech University, Joseph Dannemiller, May 2019
roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 215 6.22 Differences in magnitude between the scale parameter recorded by TTUHRT PMT Black and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 216 6.23 Differences in magnitude between the scale parameter recorded by TTUHRT PMT Black and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 217 6.24 Distribution of mean differences in magnitude between the location parameter recorded by TTUHRT PMT Black and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 218 6.25 Distribution of mean differences in percent between the location parameter recorded by TTUHRT PMT Black and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 219 6.26 Distribution of mean differences in magnitude between the scale parameter recorded by TTUHRT PMT Black and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 220 6.27 Distribution of mean differences in percent between the scale parameter recorded by TTUHRT PMT Black and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Francis raw data for a single 600s window ...... 221 6.28 Differences in magnitude between the peak 3s wind speed recorded by the TTUHRT PMT Black during the 2004 landfall of Hurricane Francis and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Francis time history ...... 223 6.29 Differences in percent between the peak 3s wind speed recorded by the TTUHRT PMT Black during the 2004
xxvii Texas Tech University, Joseph Dannemiller, May 2019
landfall of Hurricane Francis and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Francis time history ...... 223 6.30 Differences in magnitude between the peak 60s wind speed recorded by the TTUHRT PMT Black during the 2004 landfall of Hurricane Francis and peak 60s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.18 using all 600s windows in the Hurricane Francis time history ...... 224 6.31 Differences in percent between the peak 60s wind speed recorded by the TTUHRT PMT Black during the 2004 landfall of Hurricane Francis and peak 60s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.18 using all 600s windows in the Hurricane Francis time history ...... 225 6.32 Deployment location for the TTUHRT WEMITE #2 platform during the 2003 landfall of Hurricane Isabel ...... 227 6.33 Simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 229 6.34 Simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 230 6.35 Differences in magnitude between the location parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 231 6.36 Differences in percent between the location parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distribution of location parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 232 6.37 Differences in magnitude between the scale parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 233
xxviii Texas Tech University, Joseph Dannemiller, May 2019
6.38 Differences in percent between the scale parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distribution of scale parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 234 6.39 Distribution of mean differences in magnitude between the location parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 235 6.40 Distribution of mean differences in percent between the location parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distributions of location parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 236 6.41 Distribution of mean differences in magnitude between the scale parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 237 6.42 Distribution of mean differences in percent between the scale parameter recorded by TTUHRT WEMITE #2 and the values in the simulated distributions of scale parameters for the smooth roughness case study using TTUHRT Hurricane Isabel raw data ...... 238 6.43 Differences in magnitude between the peak 3s wind speed recorded by the TTUHRT WEMITE #2 during the 2003 landfall of Hurricane Isabel and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Bonnie time history ...... 240 6.44 Differences in percent between the peak 3s wind speed recorded by the TTUHRT WEMITE #2 during the 2003 landfall of Hurricane Isabel and peak 3s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.44 using all 600s windows in the Hurricane Bonnie time history ...... 241 6.45 Differences in magnitude between the peak 60s wind speed recorded by the TTUHRT WEMITE #2 during the 2003 landfall of Hurricane Isabel and peak 60s wind speeds computed using the mean observed by the TTUHRT platform
xxix Texas Tech University, Joseph Dannemiller, May 2019
and a gust factor equal to 1.18 using all 600s windows in the Hurricane Bonnie time history ...... 242 6.46 Differences in percent between the peak 60s wind speed recorded by the TTUHRT WEMITE #2 during the 2003 landfall of Hurricane Isabel and peak 60s wind speeds computed using the mean observed by the TTUHRT platform and a gust factor equal to 1.18 using all 600s windows in the Hurricane Bonnie time history ...... 243 7.1 Relationship between parent mean and turbulence intensity using TTUHRT data from all 600s windows...... 245 7.2 Relationship between parent mean and EV TI using TTUHRT data from all 600s windows, broken into roughness regimes (smooth, open and rough) ...... 246 7.3 Relationship between parent mean and EV TI using only TTUHRT data from 600s windows with a mean wind speed greater than 15m/s, broken into roughness regimes (smooth, open and rough) ...... 247 7.4 Relationship between parent mean and EV COV using TTUHRT data from all 600s windows, broken into roughness regimes (smooth, open and rough) ...... 248 7.5 Relationship between parent mean and EV COV using only TTUHRT data from 600s windows with a mean wind speed greater than 15m/s, broken into roughness regimes (smooth, open and rough) ...... 249 7.6 Relationship between TI and EV TI using TTUHRT data from all 600s windows, broken into roughness regimes (smooth, open and rough) ...... 250 7.7 Relationship between TI and EV COV using TTUHRT data from all 600s windows, broken into roughness regimes (smooth, open and rough) ...... 250 7.8 Relationship between TI and EV TI using only TTUHRT data from 600s windows with a mean wind speed greater than 15m/s, broken into roughness regimes (smooth, open and rough) ...... 251 7.9 Relationship between TI and EV COV using only TTUHRT data from 600s windows with a mean wind speed greater than 15m/s, broken into roughness regimes (smooth, open and rough) ...... 251
xxx Texas Tech University, Joseph Dannemiller, May 2019
xxxi Texas Tech University, Joseph Dannemiller, May 2019
DEFINITION OF VARIABLES
𝐶 – Expected gust factor 𝑓𝑟𝑒𝑞𝑢 – Sampling frequency for TTUHRT instrument platform, Hz 𝑔 𝑡 – Peak factor for computing gust factors
𝑔 , 𝑡, T – Segment gust factor of duration t seconds in segment k, from T second duration window 𝑗 of a time history 𝑖 – Index for time history data, from 1 to 𝑛
𝑗 – Index for windows in time history records, from 1 to 𝑤
𝑘 – Index for segmenting windows, from 1 to 𝑠 𝑙 – Largest lateral, longitudinal or diagonal distance for a structural component used in computing wind gust averaging times 𝑛 – Total number of data points per time history record
𝑛 – Total number of data points in a 𝑇 second segment
𝑛 – Total number of data points in a 𝑇 second window
𝑠 – Total number of segments in a time history window 𝑠𝑑 – Standard deviation 𝑠𝑑 𝑡, 𝑇 – Standard deviation of t second duration gusts wind speeds in T second time history
𝑠𝑑 𝑡, T – Standard deviation of t second gust wind speeds in the T second window 𝑗 of a time history
𝑠𝑑 𝑟𝑎𝑤, 600𝑠 – Standard deviation of wind speeds from window 𝑗 in a TTUHRT raw data time history, where the gust duration depends on the frequency of capture, and the record duration is equal to 600s
𝑠𝑑 3𝑠, 600𝑠 – Standard deviation of wind speeds from window 𝑗 in a TTUHRT 3s MA data time history, where the gust duration equals 3s, and the record duration is equal to 600s
𝑠𝑑 60𝑠, 600𝑠 – Standard deviation of wind speeds from window 𝑗 in a TTUHRT 60s MA data time history, where the gust duration equals 60s, and the record duration is equal to 600s 𝑡 – Time duration of gust wind speeds, second
𝑤 – Total number of complete 𝑇 second windows in a time history record
xxxii Texas Tech University, Joseph Dannemiller, May 2019
𝑥 𝑡 – Single t second gust wind speed value 𝑥̅ – Distribution mean value 𝑥̅ 𝑡, 𝑇 – Mean of t second duration gusts wind speeds in a T second time history
𝑥 ̅ 𝑡, T – Mean wind speed of t second duration gusts in the T second window 𝑗 of a time history
𝑥 ̅ 𝑟𝑎𝑤, 600𝑠 – Mean wind speed from window 𝑗 in a TTUHRT raw data time history, where the gust duration depends on the frequency of capture, and the record duration is equal to 600s
𝑥 ̅ 3𝑠, 600𝑠 – Mean wind speed from window 𝑗 in a TTUHRT 3s MA data time history, where the gust duration equals 3s, and the record duration is equal to 600s
𝑥 ̅ 60𝑠, 600𝑠 – Mean wind speed from window 𝑗 in a TTUHRT 60s MA data time history, where the gust duration equals 60s, and the record duration is equal to 600s
𝑥̅ 𝑧 𝑎𝑛𝑑̅ 𝑥 𝑧 – mean wind speeds used when standardizing time history data 𝑥 𝑡, 𝑇 – Peak t second gust wind speed in a T second time history
𝑥 𝑡, T – Peak t second gust wind speed in the T second window 𝑗 of a time history
𝑥 , 𝑡, T – Peak t second gust wind speed, in the T second segment,k, from the j-th T second window of a time history
⏞𝑥 , 𝑡, 𝑇 – Estimate of the peak t second gust wind speed in the T second window 𝑗 of a time history computed using 𝑥 ̅ 𝑡, T multiplied by the HWIND gust factor equal to 1.18
𝑧 – Surface roughness, meters
𝑧 and 𝑧 – Heights used when standardizing time history data, meters
𝑧 and 𝑧 – Surface roughness values used when standardizing time history data, meters
𝑧 , – Surface roughness upwind of a roughness change, meters
𝑧 , – Surface roughness downwind of a roughness change, meters
𝑧 – Surface roughness where data is recorded, in a region containing a transitional boundary layer above the site, meters
xxxiii Texas Tech University, Joseph Dannemiller, May 2019
𝑧 – Height data is recorded by an observational platform, meters 𝐴 – Area over which wind induced pressures act 𝐶 – Product of coefficients accounting for phenomena affecting the wind field when computing wind induced pressures
𝐷 – Difference between x raw, 600s and U 600s , reported as both [m/s] and [%] 𝐺 𝑡,𝑇 – Gust factor for t second gust in a T second duration record 𝐺 60𝑠, 600𝑠 – Gust factor used in HWIND to compute the peak 60𝑠 of wind in a 600𝑠 window
𝐺 𝑡, T – Gust factor of t second duration gusts in the T second window 𝑗 of a time history
𝐻 – Height of an atmospheric observational system, used in HWIND wind speed standardization, meters
𝑆𝐷𝐸𝑉 𝑡, T – Standard deviation of extreme t second duration gusts wind speeds in a T second time history
𝑀𝐸𝑉 𝑡, T – Mean of extreme t second duration gusts wind speeds in a T second time history 𝑂 – Outcome of structural performance assessment, pass/fail or computed probability 𝑄 – Load acting on a component/system used in structural performance assessment 𝑅 – Resistance of a component/system used in structural performance assessment
𝑅 – Distance from observational system to the center of a hurricane
𝑅 – Radius of maximum winds, distance from the location of the maximum wind in a storm relative quadrant to the center of the hurricane
𝑇 – Time duration of time history record, seconds
𝑇 – Time duration of segment, seconds
𝑇 – Time duration of window, seconds
𝑇𝐻 – TTUHRT measurement station time history data where gust duration time depends on the frequency at which data was captured
xxxiv Texas Tech University, Joseph Dannemiller, May 2019
𝑇𝐻 – TTUHRT measurement station 3s moving average time history data
𝑇𝐻 – TTUHRT measurement station 60s moving average time history data
𝑇𝐼 𝑡, T – Turbulence intensity ratio of t second duration gusts in the T second window 𝑗 of a time history 𝑈 600𝑠 – HWIND computed 600s mean wind speed at 10-m in open exposure 𝑈 60𝑠, 600𝑠 – HWIND computed peak 60s wind speed in a 600s window, at 10-m over open exposure 𝑉 – Wind speed used to compute wind induced pressures 𝑉 – Mean wind speed used in the Time, Velocity, Length model for computing wind gust averaging times 𝑋 – Distance downwind from a change in roughness to where data is recorded, meters
𝑎 – Estimate of extreme value Type-I location parameter 𝜇 – Distribution location parameter 𝜇 𝑡, 𝑇 – Extreme Value Type-I location parameter for t second duration gusts wind speeds in a T second time history
𝜇 𝑡, T – Extreme Value Type-I location parameter for t second duration gusts wind speeds in a t second duration time history, unadjusted back to the windows length
𝑢 𝑡, T –Location parameter of t second duration gusts in the T second window 𝑗 of a time history, adjusted to the window length
𝜇 𝑟𝑎𝑤, 600𝑠 – Extreme Value Type-I location parameter from window 𝑗 in a TTUHRT raw data time history, where the gust duration depends on the frequency of capture, and the record duration is equal to 600s
𝜇 3𝑠, 600𝑠 – Extreme Value Type-I location parameter from window 𝑗 in a TTUHRT raw data time history, where the gust duration depends on the frequency of capture, and the record duration is equal to 600s
𝜇 60𝑠, 600𝑠 – Extreme Value Type-I location parameter from window 𝑗 in a TTUHRT 60s MA data time history, where the
xxxv Texas Tech University, Joseph Dannemiller, May 2019
gust duration equals 60s, and the record duration is equal to 600s 𝜌 – Air density
𝛿 – Maximum height of a well established boundary layer after an upwind change in roughness
𝛿 – Maximum height of the transition boundary layer after an upwind change in roughness 𝜎 – Distribution scale parameter
𝜎 𝑡, T – Extreme Value Type-I scale parameter for t second duration gusts wind speeds in a T second time history, unadjusted back to the windows length
𝜎 𝑡, T – Extreme Value Type-I scale parameter for t second duration gusts wind speeds in a T second time history, adjusted back to the windows length
𝜎 𝑟𝑎𝑤, 600𝑠 – Extreme Value Type-I scale parameter from window 𝑗 in a TTUHRT raw data time history, where the gust duration depends on the frequency of capture, and the record duration is equal to 600s
𝜎 3𝑠, 600𝑠 – Extreme Value Type-I scale parameter from window 𝑗 in a TTUHRT 3s MA data time history, where the gust duration equals 3s, and the record duration is equal to 600s
𝜎 60𝑠, 600𝑠 – Extreme Value Type-I scale parameter from window 𝑗 in a TTUHRT 60s MA data time history, where the gust duration equals 60s, and the record duration is equal to 600s
𝜎 – Log-Logistic standard deviation parameter 𝛽 – Distribution shape parameter 𝛤 – Gamma function
xxxvi Texas Tech University, Joseph Dannemiller, May 2019
CHAPTER 1
INTRODUCTION
Hurricanes are responsible for billions of dollars in damages and place at risk over 6 million homes in the United States each year. In some cases, the devastation following hurricane landfall is so severe little-to-nothing of the built environment remains in entire regions. In cases like these, on site forensic investigation is no longer viable and alternative methods, like modeling, must be employed to assess the cause, sequence and probabilities associated with structural failures. Such an analysis can be done using deterministic or stochastic methods, but stochastic methods offer the advantage of accounting for the variability of random variables affecting both the loads and resistances factoring into the outcome of a failure assessment. Of course, such an analysis can only be employed if the true variability of each random variable can be accounted for. This requires understanding the true nature of each random variable and, in the case of hurricane wind loads, modeling the distribution of extreme winds becomes paramount. To correctly model the distribution of extreme winds at a site where failure assessment is sought, high resolution, high fidelity, high precision wind records are ideal as they provide data applicable directly to the site of interest. Unfortunately, while many groups do gather, high resolution, high fidelity, high precision wind data, the network of observation stations is sparse leaving a majority of the geography affect by hurricane winds without suitable data useful in structural performance analysis (SPA). As an alternative, hurricane wind field models can be used to provide wind data for any location in the region affected by a hurricane. Hurricane wind models are most useful if two things can be done: (1) the values produced by the model can be verified using high resolution, high fidelity, high precision wind data, and (2) a methodology is advanced to take hurricane wind field modeled data and generate distributions of extreme winds suitable for stochastic SPA.
1 Texas Tech University, Joseph Dannemiller, May 2019
The objectives of this work are to investigate the hurricane wind field model developed by the National Oceanic and Atmospheric Administration Hurricane Research Division’s HWIND. To accomplish this investigation, HWIND reconstructions are compared against wind data gathered by Texas Tech University’s Hurricane Research Team (TTUHRT) collected for 32 total deployments, with 15 different platforms during 10 separate hurricane landfall events. The same TTUHRT data is then used to form relationships between parent and extreme wind fields. The established relationships and hurricane model wind field data can then be used together to simulate distributions of extreme winds for use in SPA to assess the cause, sequence and probabilities associated with structural failures leading to total structural collapse.
2 Texas Tech University, Joseph Dannemiller, May 2019
CHAPTER 2
BACKGROUND
Major Hurricane Damage Of the 314 million people living in the United States in 2013, 85 million, or 27%, live in hurricane prone coastal counties between Maine and Texas (NOAA, 2013). In these areas, around 6.6 million homes are at risk of being damaged or destroyed during a hurricane’s landfall (Botts et al, 2015). Of the 6.6 million, 2.8 million homes are on the Gulf coast, Texas to Florida, with the remaining 3.8 million on the Atlantic coast, between Florida and Maine (Botts et al., 2015). Major hurricanes cause widespread devastation, disrupt communities and result in injury and death to many people in the affected regions. Examples of costly, major hurricanes to hit the Gulf coast include the Galveston Hurricane of 1900, Hurricane Camille in 1969, Hurricane Katrina in 2005 and Hurricane Ike in 2008. The most recent major hurricane to devastate the Atlantic coast is the 2012 landfall of Hurricane Sandy and in decreasing order, Hurricanes Katrina, Sandy, and Ike are the first, second and fifth costliest natural disasters in United States history. Together, these three hurricanes have resulted in a combined total of $110.4 billion in insured losses (FEMA, 2006), (FEMA, 2009), (FEMA, 2013). Between the Galveston Hurricane of 1900 and Hurricane Katrina in 2000, 10,000 people lost their lives making these two storms two of the deadliest natural disasters in United States history. Economically, the cost to rebuild all 6.6 million wood frame homes at risk on the Gulf and Atlantic coast is estimated by (Botts et al., 2015) to be over $1.5 trillion. The homes used by Botts et al. (2015) for this estimate are only those at risk of experiencing damage from hurricane surge, and does not include homes at higher elevations that are still vulnerable to damages caused by hurricane winds. It is unlikely the full $1.5 trillion estimate will ever be realized as a hurricane making landfall in North Carolina is unlikely to also affect homes on the Gulf coast of Texas. However, the economics associated with the structures at risk (the 6.6 million wood
3 Texas Tech University, Joseph Dannemiller, May 2019 frame houses plus engineered structures and infrastructure) makes assessing how and why (wind or surge) structures fail during hurricane landfall an important issue for policy makers, home and business owners, insurers, researchers, tax payers, and specifically any member of the public with at least one asset insured against such losses. Structural damages due to hurricane winds span from minor shingle and siding loss to being completely wiped off their foundations. Situations where structures are completely wiped off their foundations are often referred to as “slab-claims”, since the only remnant of the structure post event is the concrete foundation slab itself. For perspicuity, the term slab-claim will henceforth refer to a situation where only a foundation remains. In slab-claim cases it can be difficult to determine what structural component failed first, the sequence of structural failures, or make a conclusive assessment regarding whether hurricane winds or storm surge caused any one specific failure, where both hazards are present. In fortuitous cases, similarly constructed structures near a slab-claim survive providing examples of damage patterns that can be extrapolated to the analysis of a slab-claim. However, in the case of a slab-claim, where no similar structures in proximity survive the same event, observational analysis methods are no longer conclusive. Other methods, such as simulation, must be utilized to determine the most likely first component failure, the sequence of component failures, and whether wind or surge led to any, or all, of the component failures. Bolivar Peninsula, Texas is an example of a region devastated by hurricane winds and/or surge, leaving an entire region of slab-claims after the 2008 landfall of Hurricane Ike. Before and after aerial photographs are shown in Figure 2.1 and Figure 2.2 to illustrate the magnitude of the destruction.
4 Texas Tech University, Joseph Dannemiller, May 2019
Figure 2.1 Bolivar Peninsula before Hurricane Ike Landfall, September 9, 2008 (FEMA 2009)
Figure 2.2 Bolivar Peninsula after Hurricane Ike Landfall, September 15, 2008 (FEMA 2009)
The arrows in the before and after photographs in Figure 2.1 and Figure 2.2 identify the locations of two wood frame structures for reference. Most of the other
5 Texas Tech University, Joseph Dannemiller, May 2019 structures between the referenced structures and the coastline are gone. These sites represent slab-claims and onsite forensic investigation is no longer conclusive for determining the hazard leading to failure, the progression of failures, or what ultimately led to total structural collapse. Attributing individual component failures to wind or surge and determining the order of component failures becomes a complex issue.
Loads and Resistances Assessing failures of structural components requires relating loads and strengths with deterministic, or stochastic, analysis. Both load and resistance values are random phenomena, each a function of other random variables, making assessments of failure, or non-failure, a random variable. For failure to occur, the loads acting on a component must exceed the strength of the component. In the opposite case, where the strength of a component exceeds the load, the component does not fail. Deterministic analysis is the process of using a single value for each random variable contributing to the production of a single outcome. The single values used in deterministic analysis are computed leveraging knowledge of the underlying distributions governing each random phenomenon. Stochastic analysis involves generating distributions for each random variable and then repeatedly sampling to compute a distribution of outcomes. With such a distribution, the probabilities associated with individual component failures can be assessed. Deterministic analysis is preferable where representative cases are sought to illustrate a common, or likely, result for a random outcome. In deterministic analysis, it is common for strength values to be underestimates of the mean strength of a sample of similar components, while the values of load are commonly overestimates of the mean of the load acting on a component resulting in a “worst case” condition i.e., a weak structure meets a strong load (Salmon et al., 2008). Structures survive this “worst-case” i.e., the structure can
6 Texas Tech University, Joseph Dannemiller, May 2019 still perform its intended function, when actual loads are lower than the overestimate, or when actual strengths are stronger than the underestimate. The alternative to deterministic analysis is stochastic analysis where distributions of random variables are simulated, and then sampled to compute instances of loads and strengths related together to compute probabilities of failure. The distributions of all random variables associated with the strengths and loads must be known to properly assess whether failure has occurred. The process of independently assessing whether failure has occurred is similar to deterministic analysis, except that in stochastic analysis the process is repeated many times using values sampled from each of the underlying distributions. The number of times this process is repeated, or repetitions, is governed by how many realizations are necessary to converge on a statistically significant solution, or one with an acceptable level of confidence. Once the independent assessments have been completed, the probability that failure, under the utilized loading condition, can be computed. An in-depth discussion of stochastic analysis specifically as it pertains to structures, is given in (Nowak and Collins, 2000). As stated above, the process of assessing whether failure has occurred in a structural system requires understanding and relating two phenomena, loads and strengths. The strength of a system, henceforth referred to as resistance, is governed by the properties of the materials used, member geometries and the type of structural system. The loads acting on the system are a function of the environment, and during hurricane landfall the two of principle concern are loads due to storm surge and hurricane winds. Storm surge is left for another discussion and only hurricane winds are discussed herein. Where loads and resistances are known, they can be related to determine if a failure has occurred using Equation 2.1.
O R Q 2.1
7 Texas Tech University, Joseph Dannemiller, May 2019
Where O denotes the outcome, R denotes the resistance of the system, and Q denotes the load acting on the system. If the outcome is less than zero, the load exceeds the system’s ability to resist and failure occurs. Instead of referring to outcomes less than or equal to zero as failures, some literature refers to these scenarios as limit states. Failure of a limit state is the point where a structural component/system stops performing its intended function. Limit states fall into two main categories, ultimate limit states (ULS) and serviceability limit states (SLS). When ULS are reached, the properties of a system physically change. In most cases, a structural system can no longer resist the loading environment and the structure is no longer sound. In other words, the strength of a system can be reduced to a fraction of its original capacity, or all the way to zero. Serviceability limit states refer to vibration, deflection and localized deformations. When serviceability limit states have been surpassed, a structure can still be deemed sound, but may not be fit due to human perception. As an example, excessive vibrations and deflections can make structures undesirable to some human occupants. Occupant perceptions can lead to illness, or to some occupants losing faith in a structure as some equate deflections with failures. When a structure is wiped off its foundation, ULS have been exceeded and SLS are no longer a governing concern. The first ULS reached (first failure), as well as the sequence in which ULS occurred (damage sequence), is difficult to impossible to identify. Where damage sequences are difficult to ascertain, the systems comprising an entire structure can be analyzed, one at a time, to assess the conditional probabilities associated with all possible damage sequences. A scientist that understands the engineering and physics applicable to a specific assessment can reduce the number of possible damage sequences investigated when assessing causal failure. This can make the assessment more manageable in scope, but even analyses that are less than holistic can be long and cumbersome. The process of assessing the conditional probabilities, for all possible ULS damage sequences, for a single structure, is henceforth referred to as Structural Performance Analysis (SPA).
8 Texas Tech University, Joseph Dannemiller, May 2019
Wind Loads As hurricane wind loads are the primary concern herein, it is important to understand how winds impart loads on structures. Wind induced loads come from dynamic pressures exerted on structural systems, acting over the area of a system, or on a localized sub-system. An example of a sub-system representing the “weak link” in a structural system is a fastener used to attach metal roofing to a main structural element. The wind loads acting on the fastener are highly localized and are higher in magnitude per unit area than the wind loads acting on the panel, per unit area, that fastener attaches to the structure. The differentiation in the area used to compute wind pressures is important because wind loads do not act fully correlated, or can be lessened due to pressure coefficients applicable over differentiated areas. The area over which wind loads are of concern for any system can be computed using the Time, Velocity and Length (TVL) model from (Lawson, 1980). The TVL method is used to determine the size of wind gust that acts fully correlated over the area of a component. The size of gust is measured in time, t, and is computed using Equation 2.2 below.
tV 4.5l 2.2
Where the longest length between any two points bounding the area for a component is denoted as l, and the mean wind velocity is denoted as V . The wind gust of duration t is then used to compute the wind loads acting on a component based on its dimensions. As the mean wind speed for finite intervals changes, using the magnitude of the mean wind speed to determine the governing gust duration size is impractical in application. However, considerable work has been completed in developing structural codes which establish standard gust durations for: (1) entire structural systems, and (2) for smaller components of structures. For larger structural systems, (Cook, 1985b) uses a gust duration of 60s. In the United States, the gust durations established by ASCE 7-10 Design of Loads for Buildings and Other Structures (ASCE 7, 10) are 3s as 3s gust durations are recorded by the National Weather Service’s observational
9 Texas Tech University, Joseph Dannemiller, May 2019 systems (ASCE &, 10). The work herein will utilize raw wind data, as well as 60s and 3s gust wind data. Computing wind loads can vary based mostly on two things, how wind loads are assessed and how the phenomena affecting the wind field are considered. In terms of assessment, (Cook, 1985a) breaks the “problem” of assessing wind loads into categories that couple variations of three pieces: the wind climate, the atmospheric boundary layer, and a structure. These three pieces are separated because each represents a major factor governing when, and how, peak loads act on a structure. The wind climate refers to wind events that occur on the order of several days to years. Events occurring at this frequency that produce the highest wind speeds are thunderstorms, downbursts, hurricanes, tornados, and several other unique weather phenomena with winds higher than those produced by the daily diurnal cycle. The boundary layer is the second piece of how (Cook, 1985a) breaks apart the “problem” of assessing wind loads. Boundary layer considerations account for how the wind field is modeled near the surface of the earth. The shear caused by interactions with the surface of the earth itself, plus the man-made structures on the surface, slow wind speeds close to the surface. The last piece in assessing wind loads is the effects of a structure itself. Structural characteristics such as the height, geometry and material makeup affect what percentage of the dynamic pressures produced by the wind field translate into a structure, versus the percentage that breaks around a structure like water around a rock. Combining these three pieces together, a peak wind speed can be determined based on: (1) a type of wind event, and (2) our understanding of how such a wind field could be modeled near the surface. Once a wind speed has been selected, wind loads can be computed, for any assessment method (static, quasi-static, quasi- steady, pseudo-steady) described in (Cook, 1985a), using the same generic equation, Equation 2.3.
1 Q ∗ρ∗ V ∗C∗A 2.3 2
10 Texas Tech University, Joseph Dannemiller, May 2019
Where the wind load, denoted as Q, is computed as a function of the air density, denoted as ρ, the square of a wind velocity, denoted as V, a set of coefficients accounting for the random variables affecting the wind velocity, denoted as C, and the area over which the wind pressure acts, denoted as A. Some phenomena included in computing values of C include localized pressure coefficients, upwind terrain roughness, terrain elevation changes, turbulence intensity, gust ratios, structural height, etc. Structures in the United States are governed by (ASCE 7,10). ASCE 7- 10, which uses both quasi-static theory, where fluctuations in loads are directly correlated to fluctuations in boundary layer winds; and the quasi-steady theory, where wind loads are a function of both boundary layer induced fluctuations in the wind field, and building induced fluctuations (Cook, 1985). Simulating wind loads in SPA requires knowledge of ρ, C, A, but most importantly V. The magnitude of the wind velocity is much larger than the other random variables, and then it is squared. This makes accurate estimates of wind speeds vital when computing wind loads. In stochastic SPA, it is important for the distribution of wind speeds to correctly model the highest wind speeds occurring in finite periods. The distribution that best models the highest wind speeds is an extreme value (EV) type-I distribution. But, before the EV distribution can be discussed further, the nature of the distribution of parent winds must be explored further The parent wind field contains every wind speed measurement in a time history record. The statistical distribution that best fits the parent wind field, in the hurricane environment, is the Weibull distribution (Hennessey, 1978), (Justus, 1978), (Simiu and Scanlan, 1996) and (Seguro and Lambert, 2000). The Weibull distribution is defined by location, scale and shapes parameters providing for translation of the distribution itself, and an adaptability in the shape of the distribution as the scale and shape parameters change. Some research has fit parent wind data with a Rayleigh distribution (Conradsen and Nielsen, 1984), where the Weibull shape parameter equals two, making this a specialized fit using a Weibull distribution by constraining the
11 Texas Tech University, Joseph Dannemiller, May 2019 shape parameter to a single value. The Probability Density Function (PDF) for the three parameter Weibull distribution is expressed in Equation 2.4.
β 𝑥 μ 𝑓 𝑥|μ,σ,β ∗ ∗𝑒 2.4 σ σ
Where 𝑓 𝑥|μ, σ, β denotes the probability, f(x), of a value x occurring given values of location, scale, and shape parameters, denoted as μ, σ, and β, respectively. The mean, denoted as 𝑥̅, and variance, denoted as 𝑣 , of the Weibull distribution are computed using numerical software and fit Equations 2.5 and 2.6.
1 𝑥̅ μ σ∗Γ 1 2.5 β