Hurricane Florence Makes Landfall in the Southeast USA: Sensitive Dependence on Initial Conditions, Parameterizations, and Integrated Enstrophy

The overall purpose of this paper is to post-evaluate the predictability of Hurricane Florence using the Advanced Research Weather Research Forecast (WRF) (ARW) version of a mesoscale model. This was performed over the period from 0000 UTC 13 September 2018 through 0000 UTC 18 September 2018. The WRF ARW core resolution used here was the 27-km grid spacing chosen to in order to balance finer resolution against in house processing time and storage. The large-scale analysis showed that a change in the Northern Hemisphere flow regime, especially the flow in the western part of the Northern Hemisphere may have contributed partly to the reduced forward speed of the tropical cyclone. In order to measure the predictability of a system, we will use different convective and boundary layer schemes initialized from the same conditions. The results demonstrated that the sign of the local IRE tendency was similar to that of the Northern Hemisphere Integrated Enstrophy. The results also showed that when the boundary layer, convective, and cloud microphysical schemes of the model were varied, the areal coverage of heavy precipitation of Florence was under-forecast by approximately 10% or more, and the heaviest amounts were under-forecast by an average of about 20%.


Introduction
Hurricane Florence made landfall in the southeast United States during [14][15][16][17] September 2018 resulting in more than 15 deaths and causing over one billion dollars in damage [1]. The storm was associated with at least one report of more than 75 cm of rain over the four-day period, see the Community Collaborative Rain Hail and Snow network (CoCoRaHS-http://www.cocorahs.com/). Hurricane Florence began as a tropical depression over the eastern tropical Atlantic on 31 August 2018 [2]. By 11 September, Florence was classified as a major hurricane (category four) on the Saffir-Simpson Scale with maximum sustained winds of more than 70 m•s −1 (140 mph). However, the storm made landfall as a category one storm on 14 September after it was anticipated to make landfall as a major hurricane.
While the track of Florence was well-forecast by operational numerical models, the rapid loss in intensity as the storm approached the coast as well as the amount of rain produced by the storm was not well forecast. Hurricane track prediction has improved substantially in recent years, but models often fail to capture the intensity of very strong tropical cyclones [3]. They [3] demonstrated that the predictability of a rapidly developing tropical cyclone (Hurricane Patricia-October 2015) was associated with dynamics consistent with violent turbulence, or rough dependence on the initial conditions (RDIC) [4] [5]. Such conditions indicate that, in some cases, predictability of certain phenomena at some time scales is nearly impossible given our current understanding of fluid dynamics [6].
The quantity integrated enstrophy (IE) was first proposed by [7] and they demonstrated that in a barotropic atmosphere, IE is proportional to the system stability and predictability. Then, [8] showed that if a large region of the NH was used to calculate IE, the regional estimate and its change with time was similar to the hemispheric value. Recent studies [3] [9] [10] [11] used the quantity Integrated Regional Enstrophy (IRE) to examine predictability within mid-latitude events such as synoptic-scale cyclones and blocking events. The study of [9] examined the development of a cold-season mid-latitude cyclone over the southeast United States using two different models and varying the different boundary layer physics and convective precipitation parameterizations. They found that differences in forecasts using parameterization schemes within a model were smaller than those using two different models.
The work of [10] used IE to study the performance of the Global Ensemble Forecast System (GEFS) in simulating the dynamic behavior of two blocking events. They found that as expected, individual ensemble members performed better than the control or the ensemble mean, but that the IE was forecast reasonably. Then [11] demonstrated that the timing of onset and the location for blocking was forecast well by the GEFS ensemble mean fields. However, the block intensity was not forecasted well by the GEFS model. Additionally, the GEFS model could not sustain the blocking event for as long as they were ob-served to exist. The model had difficulty with the timing of block termination as well. Finally, [11] examined the relationship between block intensity and IRE during block onset. They found that the rate of precursor upstream rapid cyclogenesis was related to a peak in IRE near block onset followed by a peak in block intensity.
This study will investigate the ability of the in-house version of the Advanced Research Weather Research Forecast (WRF-ARW) model to represent the behavior of the large-scale conditions associated with the NH flow during the landfall of Hurricane Florence as well as the area and amount of the heaviest precipitation associated with the storm. This study will examine the dynamics of the large-scale NH flow at the time of landfall for Hurricane Florence in section three using teleconnection indexes and IRE as the diagnostics. It is postulated that changes in the hemispheric large-scale flow may have been partly responsible for the weakening and change in the forward motion of the storm. Additionally, the sensitivity of the IE diagnostic and precipitation character of Hurricane Florence to the WRF-ARW model parameterizations of convective and cloud microphysics schemes, as well as boundary layer schemes will be compared in section four. Section five will present the main conclusions.

Data
In this study, observational data were obtained from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analyses. These data are archived at the NCAR research facilities in Boulder, CO [12] [13], which can be used to provide large-scale meteorological data at various resolutions from 1˚ × 1˚ to 2.5˚ × 2.5˚ latitude-longitude grids. The characteristics of Hurricane Florence were obtained from [2]. Finally, teleconnection index information was obtained from the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center (CPC) teleconnections website (https://www.cpc.ncep.noaa.gov/products/precip/CWlink/). The daily Arctic Oscillation (AO), North Atlantic Oscillation (NAO), and Pacific North American (PNA) teleconnection indexes were examined here and are described more thoroughly in [14] or [15] and references therein. In order to produce an ensemble, the convective schemes and boundary layer schemes were varied here and these are described below.

Methodology
The work of [7] postulated that if the atmosphere is barotropic, the sum of the positive Lyapunov exponent(s) in the atmosphere can be expressed as the area of integrated enstrophy (IE). Then, [8] demonstrated that the IE summed over a large area of the mid-latitude Northern Hemisphere over time will correspond strongly to the hemispheric value, which they called IRE or Integrated Regional Enstrophy. This is expressed as; where λ i is the ith Lyapunov exponent that is greater than zero in a dynamic system, ζ is the vorticity, or the curl of the wind vector, and the quantity squared is called enstrophy, which is the dissipation tendency of a fluid. A constant proportional to unity is implied by the very high correlation [3] [6]. Also, [8] demonstrated the utility of IRE in identifying the onset and termination of atmospheric blocking, while [10] showed that this quantity could be used to study blocking dynamics. Additionally, [10] and [17] used this quantity to study flow regime transition. In (1), the vorticity is calculated via the geostrophic relationship using second-order finite differencing and IRE is calculated within the region defined above.

Methodology
In this study, a pseudo ensemble was produced in house by running the WRF-ARW model using several different convective, cloud microphysics, and boundary layer schemes from the same initial state. The initial state was chosen was 0000 UTC 13 September 2018 as the tropical cyclone was weakening and had begun to slow its forward motion. This was also the first time that the storm was in the research domain chosen for heavy precipitation studies in the southeast USA by the University of Missouri synoptic scale dynamics group. For the control simulation schemes WRF Single Moment 6-class was used for cloud microphysics, the Tiedke scheme for cumulus parameterization, the RRTMG for long and shortwave radiation parameterization, the Mellor-Yamada-Janjic (YMJ) for boundary layer parameterization, the Unified Noah for clay surface parameterization, and Eta Similarity for surface physics parameterization. Changes were made to one parameterization at a time to identify changes to the simulations and were analyzed as in the following text below. The Grell-Freitas convective (GFC) scheme features a scaling and aerosol-aware approach to convective parameterization and features a tri-modal ap-proach to handling convection, which matches the types of convection present in tropical systems [18]. This scheme is utilized seldom in recent tropical cyclone sensitivity studies in spite of being used heavily to study other tropical phenomena [19]. Previous studies have shown that cumulus physics schemes are more likely to affect intensity than storm track [20] [21]. However, given the unique track of Florence, the track produced by the GFC scheme was of interest as well. Because of the relatively new improvements made to this scheme, which includes the handling of mid-level convection and adjustments to rainfall estimation [22], it would be useful to know how well it predicts intensity and track adjustments as well as the predictability accuracy in comparison to other physics schemes.
The Betts Miller convective (BMC) scheme is a convective adjustment scheme, which includes both deep and shallow convection [23]. For deep convection, this scheme uses quasi-equilibrium thermodynamic profiles as a reference whose shape is based on observations [23]. The work of [24] made appropriate adjustments to the BMC (now BMJC) scheme to add the cloud efficiency to reduce the excessive precipitation [23].
The BMJC scheme strengths are its ability to operate well in moist environments without a capping inversion. It includes complicated features but does not require much computing resources making it ideal for use in house. The BMJC scheme demonstrated lower bias in precipitation over land and over the Gulf of Mexico region when combined with WRF single moment microphysics class5 (WSM5) microphysics option. The BMJC also provided better hurricane track forecasts with the Kessler (KESS) [25] [26] or the microphysics [27], compared with other schemes such as the Kain and Fritsch convective (KFC), no cumulus parameterization (NCP), and Grell-Devenyi (GD) [28] [29]. The BMJC scheme, however, produces more precipitation than Kain-Fritsch convective (KFC) scheme in the simulations of Hurricane Ivan [30]. The BMJC and KFC schemes both gave a clear representation of the intense structure of Tropical Cyclone Nargis, but the KFC scheme produced a realistic intensity of the storm [31]. This scheme generally deals with convection better than other cumulus parameterization schemes, thus, this scheme was chosen for this study.
The KFC parameterization scheme was chosen for use here in spite of the fact is not currently integrated into most tropical cyclone numerical modelling studies. However, it has been used in the past for studies since it does handle convection reasonably well both across the tropics and the mid-latitudes [28] [32]. A key reason why the KFC was chosen was because the mass-flux scheme that is embedded within this scheme handles convection in such a way that it accommodates both entrainment and detrainment for both updrafts and downdrafts, respectively [33]. This is particularly important for tropical cyclone simulations since spiral banding features within hurricanes have corresponding regions of updrafts and downdrafts associated with the inner and outer parts of the warm-core low. In addition, the KFC handles convective updraft generation such that these updrafts remove higher theta-e values from the lower tropos-phere and transport it aloft, a key component of energy transfer within the eyewall and inner core of a tropical cyclone [34]. The parameterization of moist downdrafts that is included in the KFC scheme, especially for stronger downdrafts, reduces compensatory subsidence and consequently generates more realistic convective warming and drying tendencies at the lower part of the convective boundary layer. Furthermore, the KFC has been designed for models with 20 to 50 km spatial resolution, making it a viable candidate for use here [32]. Additionally, the scheme also allows for the development of an effective secondary circulation that often generates a numerically modelled environment, which is more conducive for generating more intense warm-core low-pressure systems. This would be helpful for modelling features within an intense tropical cyclone such as (but certainly not limited to) the eye, eye-wall, a moat when and if one was present, spiral convective bands, and other features.
Numerical weather prediction models are known to have difficulty predicting precipitation. One of the reasons for this is the lack of a realistic depiction for clouds in these models. It is for this reason that integrating schemes that address the issue of cloud and precipitation forecasting is crucial. One such scheme would be the Ferrier scheme that has been included in the Eta scheme. The Ferrier (FER) scheme is ideal for use when looking at moisture transport and forecasting precipitation. This scheme has a vertical resolution of 38 layers and a horizontal resolution at 40 km. The FER scheme addresses the issue of forecasting clouds and precipitation by using separate schemes for clouds and precipitation within the code. For the cloud scheme, the primary feature is the calculation of cloud water and cloud ice, since the cloud water/ice mixing ratio is a key variable in the cloud formation. By including cloud ice as well as the horizontal advection of clouds, a large improvement on cloud precipitation and relative humidity is achieved. After calculating the cloud water/ice mixing ratio, the FER is used for the diagnostic calculation of precipitation. Few schemes or models have a specific approach for calculating precipitation directly, rather they get precipitation outcomes as a result of calculating a different parameter, this is one of the reasons the FER scheme is so useful. The precipitation parameterization accounts for two important features. The first being the gradual melting of snow in warm cloud regions allowing for snow and rain to coexist, and the second being it takes into account evaporation from precipitation.
There are many reasons for using the FER scheme to analyze the moisture transport associated with hurricane Florence. The first reason is that this scheme addresses the issue of excessive condensation of moisture over oceans. In order to account for this, the critical value of relative humidity is set as a constant. Six different precipitation processes are accounted for also in the scheme. These processes are observed in clouds to produce a more accurate precipitation forecast. They include, auto conversion of cloud water to rain, collection of cloud droplets by the falling rain drops, auto conversion of ice particles to snow, melting of snow below the freezing level and the evaporation of precipitation below cloud bases. The FER scheme also explicitly defines cloud levels, with the low cloud layer being from the surface to 640 hPa, the mid-level from 640 -350 hPa and the high level from 350 hPa to the top of the model. The cloud scheme has improved forecasts for precipitation at all thresholds and is most apparent at levels below 850 hPa and above 550 hPa (e.g. [28]).
The WDM6 scheme is a double moment, bulk microphysical scheme that predicts the mixing ratio of water vapor, cloud droplets, cloud ice, snow, rain, and graupel along with the concentrations for rain and cloud water (e.g. [35]).
Single moment schemes predict only the mixing ratios. As such, this scheme is ideal for predicting the number concentration of cloud and rainwater. The double-moment approach allows for an evolving mean diameter, leading to a more efficient simulation of microphysical processes. Additionally, this scheme includes prognostic equations that allow for investigating aerosol effects on clouds and precipitation process. Simulations of the WDM6 can be run at different resolutions. Coarse resolution simulations will not resolve supersaturation in one grid box, thus, higher resolution simulations should be employed if this variable is desired. A bulk microphysical parameterization is computationally efficient compared to a bin-resolving method. Single moment bulk schemes predict only the mixing ratios of hydrometeors using a distribution function for the hydrometeor size. A double-moment scheme has the advantage of, not only predicting the mixing ratio, but also the number concentration of hydrometeors. A double moment approach allows more size flexibility and the diameter of the hydrometeor can evolve improving the single moment approach. Double moment schemes require more computational time but have been shown to simulate precipitating convective clouds more accurately. Double-moment schemes rely on the accuracy of several microphysical processes. The WDM6 scheme predicts the concentration of cloud and rainwater with the CCN. The six water species types that the WDM6 predicts mixing ratios for are water vapor, cloud droplets, cloud ice, snow, rain, and graupel.
The WDM6 scheme was chosen because of its computational efficiency and its ability to simulate precipitation convective clouds compared to a single moment scheme. In the control run of the experiment the CONUS physics suite was used, since the focus of the experiment is the land-falling period of the hurricane. The CONUS physics suite uses the WRF Single moment 6 class Scheme (WSM6). In an attempt to improve the rainfall prediction, the WDM6 scheme was tested as a comparison to the control. Features of the WDM6 scheme include the ability to be used in warm rain processes, and the prediction of CCN, and the number concentrations of cloud and rain. Studies indicate that the rain number concentration is important in determining the precipitation rate and storm morphology because of the regulation of other microphysical processes particularly the evaporation rate. The study of [36] found that increased rain-drop concentrations in double moment schemes result in more evaporation in the stratiform region of the storm and lead to lower rainfall rates. In general [36] found that light precipitation activity is reduced compared to the WSM6 while moderate and heavy precipitation increases. However, a high-CCN condition can suppress conversion from cloud droplet to raindrops and decrease precipitation.
When forecasting weather phenomena, it is important to look at the physical processes within the planetary boundary layer (PBL) with the most important process being buoyancy. Buoyancy describes the thermodynamic and kinematic profiles that models incorporate in their forecasts. These processes continue to be poorly represented in the PBL, so forecasters use many different schemes to help combat this problem. One of these schemes is the Yonsei University (YSU) Scheme. This is an extension of the previous Medium Range Forecast (MRF) model PBL scheme (e.g., [37]) and uses the MM5 similarity theory. This scheme focuses on buoyancy transport and is a powerful scheme that many modelers can use to help better forecast. According to [38], the YSU is the first and most widely used PBL scheme.
The YSU scheme was selected primarily because of its ability to accurately portray buoyancy transport and boundary layer depth, as compared to other PBL schemes. Then, [39] has shown that the YSU scheme accurately depicts the depth of the PBL. The boundary layer height is calculated by using the bulk Richardson number method. A major advantage is that the YSU scheme represents entrainment at the top of the PBL explicitly. Previous studies, particularly with Hurricane Isabel, have shown radial winds to be consistent with observed data and the slope of the eyewall to be well predicted [37]. Furthermore, YSU tends to have a less cold and moist bias than other schemes. Since it is a non-local closure scheme, it is able to incorporate all interactions between the atmospheric quantities within all the layers of the PBL, meaning that vertical mixing from the entire depth of the PBL is considered.
The University of Washington (TKE) is a planetary boundary layer scheme specifically developed for use in a subtropical marine environment [40]. This model is based on a 1.5-order turbulent closure scheme and the [41] formulation. This boundary layer model was revised by [42] to include a new parameterization for shallow convection in cloud-topped boundary layers. The model was revised again by [43] to improve the treatment of moist turbulence in cumulus updrafts. This scheme was chosen since Hurricane Florence was located over the ocean at the time of the initial conditions.

Large-Scale Analysis
Hurricane Florence began as a tropical wave (Figure 1  September around 22˚N and 45˚W. During this time, the 500 hPa NH flow became less zonal (Figure 2), and this is supported by a decrease in the AO Index (Table 1). Also, both the NAO and PNA indexes show that a change in flow regime was occurring during this time. The 500 hPa height anomalies ( Figure 2) showed an amplified wavelike pattern across the Pacific, North America, and Atlantic at the start of September and by 7 September the flow pattern was weaker and less organized across the region. Examining the observed IRE and its time derivative (DIRE [10]) for September 2018 (Figure 3) suggests that IRE peaked at or around 3 September 2018 signaling a change in the NH flow regime adding further support for the analysis here. Additionally, the 500 hPa maps suggest that the subtropical Atlantic was dominated by anti-cyclonic flow which would steer Florence westward.
During 6 September, Florence weakened to a tropical storm and began to take a more westerly track (Figure 1). Then, the storm intensified after 9 September culminating in a maximum intensity of a Category 4 storm with wind speeds at 70 m•s −1 on 12 September. During this time, the tropical cyclone turned to move northwestward again toward the Carolinas in the USA. As the storm made landfall on 14 September, the forward speed slowed considerably and the storm had      Table 2 shows the change in IRE over the WRF domain from the initial time 1200 UTC 13 September 2018 for a 120 h period every 24-h. A comparison to Figure 3 shows that the IRE within the WRF domain was forecast well when comparing to the IE over the entire NH. In general, the NH values were decreasing until 17 September and increased for the last 24 h period. The 24-h model control forecasts were of the same sign for all but the last 24 h period. There was relatively little change (small increase) in NH IE and control model domain IRE from 48 h to 72 h, while the ensemble members forecast a modest decrease. However, in summary, the control model and the ensemble members forecast IRE were generally consistent with each other and overall followed the change in NH IE. This latter finding is consistent with the results of Athar and Lupo (2010) described above. This result implied that local predictability was increasing as Florence slowed and moved inland. Examining individual members (Table 2), the 24-h IRE tendencies in the WRF-ARW runs using the YSU (0.04) and WDM6 (0.08) boundary layer schemes were closest to the control run. These two schemes produced root-mean-square (RMS) values (0.04 and 0.08, respectively). The YSU RMS was the only run that was better than the ensemble mean. Additionally, the FER microphysics (0.14) and GFC convective schemes (0.18) produced RMS values that were superior to the remaining members.

Predictability during Landfall
Examining the pattern and rainfall amounts from each scheme also was performed. Figure 5 shows the WRF control simulation results for the heaviest 120-h precipitation (>300 mm) totals. Much of North Carolina and part of the South Carolina border regions were modeled to receive more than 100 mm of  precipitation (about 114,000 km 2 ) with isolated values up to 650 mm along the border region (Table 3). The observed precipitation totals ( Figure 5(b)) show that the WRF control model under-forecast the heaviest precipitation total amounts, but the general pattern of the heaviest amounts (areal coverage) was over-forecast in that the highest totals were found further inland in North Carolina and along the North Carolina-South Carolina state border. The heaviest observed precipitation amounts exceeded 750 mm with the greatest amount at 912.6 mm, while the areal coverage of amounts greater than 100 mm was approximately 149,000 km 2 ( Table 3). The areal coverage of precipitation was estimated by fitting a box to the precipitation area on a scaled map. Examining the individual ensemble member forecast ( Figure 6) shows that all the members made forecasts similar to the control in terms of the areal coverage of the heaviest precipitation. Thus, many of the areal coverage forecasts of heavy precipitation were under-forecast by more than 10% (Table 3) Table 3) compared most closely with the control (Figure 5(a)), while the FER and TKE schemes ( Figure 6(b), Figure 6(e), Table 3) was closest to the observed ( Figure 5(b)). The former schemes under-forecast the heaviest precipitation amounts, and the heaviest amounts in Figure 6(a), Figure 6(d) were near the coast as in the observed. These schemes also predicted local precipitation maxima in western North Carolina associated with the Appalachian Mountains as in the observed. This might be because the BJMC microphysics scheme was designed to handle precipitation associated with tropical systems and simulates tropical cyclone tracks well [25] [26] [27]. Also, the KFC scheme does simulate convection very well whether in the tropics or mid-latitudes [28] [31] [32]. Additionally, both schemes produced slightly higher maximum precipitation amounts than the control run and closer to the observed near the coast (Table 3) and the mountains (over-forecast). The TKE scheme performed very well in comparison to the observed precipitation since it was designed for marine boundary layers.
The heaviest observed precipitation was close to 1000 mm and these areas were on the coast especially near the border of the Carolinas. The schemes that produced similar amounts of precipitation were the GFC and KFC convective schemes ( Figure 6(c), Figure 6(d)) (Table 3), respectively, and the YSU planetary boundary layer scheme (Figure 6(g)). Each of these produced more than along the coast. Each of these schemes produces an elongated maximum that extends well into North Carolina, which was not observed. It is not surprising that these two convective schemes produce better precipitation amounts near the coast as these schemes include convection that is similar to observed tropical convection [18] [33].
In summary, the WRF-ARW modeled IRE captured the NH wide decrease of IE over the five-day model forecast. The WRF-ARW model domain covered the eastern two-thirds of North America and surrounding environment. The size of this domain is about one quarter that used by [8] and [10] for their studies of blocking anticyclones (40˚ latitude by 60˚ longitude) using IRE. These studies used the smaller domains as they showed that local IRE matched reasonably with the NH IE values. This seems to the case here as well, and may indicate the WRF-ARW is capturing the local NH flow regime transition to a more zonal state as shown in section three. This transition was likely responsible at least partly for the slowing of Florence as it came ashore. The WRF-ARW under-forecast the areal coverage of the heaviest precipitation for most variations of the model tested here, and generally under-forecast the maximum amount precipitation found near the coast and over-forecast the heavy precipitation found inland. Many of the model variations also captured the topologically enhanced maxima observed in western North Carolina. The model forecast using the FER and TKE schemes produced the best area coverage for the heaviest precipitation. Two of the three model runs producing the best maximum precipitation values (YSU and GFC in Table 3) also produced two of the four best IRE tendency forecasts. Thus, for the landfall of Florence, there is some correspondence between those WRF-ARW runs producing the best IRE tendencies (DIRE) and the area coverage or maximum precipitation amounts. However, more study would be required in order to determine if this correspondence is true more generally

Summary and Conclusions
A multi-scale study of the Hurricane Florence landfall during mid-September 2018 was carried out here using the WRF-ARW mesoscale model. This was a major storm, which produced record rainfalls along the Carolina Coast and inland over the southeast USA causing severe economic damage and casualties. The large-scale study was performed using the NCEP re-analyses. The mesoscale study was done using model output from the WRF-ARW over a 120-h period from 1200 UTC 13 to 1200 UTC 18 September 2018 in which the convective, cloud microphysics, or boundary layer schemes were varied. This period covered the hurricane landfall and stalling over the southeast USA and subsequent flooding.
The key findings were: • Hurricane Florence began as a classical eastern Atlantic storm and moved toward the northwest around 3 September 2018. The course of the storm became westward after a change in the NH flow regime from zonal to meridional; • The storm slowed its forward movement as it came ashore on 14 September 2018 when the NH flow regime was transitioning back to a zonal pattern characteristic of early September The NH IRE was decreasing during the latter stage of the meridional flow regime, but increased as the zonal regime was established; • The WRF-ARW model control run and ensemble mean captured well the observed NH IRE tendency. However, only the YSU model run matched well the mean IRE tendency and root-mean-square values; • The model generally under-forecasted the areal coverage of the heaviest precipitation by 10%. The FER and YSU cloud microphysics and boundary layer schemes best captured this area coverage; • The same model under-forecast the heaviest amounts of precipitation near the coast, and the GFC and KFC convective schemes and the YSU boundary layer) were closest to the coastal maximum; • Lastly, four of the best model forecasts of IRE tendency also produced either one of the four best areal coverage or coastal precipitation maxima during the forecast period.