Sensitivity of the Wave Field to High Time-Space Resolution Winds during a Tropical Cyclone

: The impact of the high space-temporal variability of the wind ﬁeld during the moderate and intense storm stages of a tropical cyclone on the wave ﬁeld as computed by the numerical model WaveWatch III is investigated in this work. The realistic wind ﬁelds are generated by a high-resolution implementation of the HWRF model in the Gulf of Mexico and stored over 15 min intervals. The spatial structure of the wind ﬁeld computed by HWRF is highly variable in space and time, although its mean structure is very similar to that described for parametric hurricanes already speciﬁed in the previous studies. The resulting storm-generated wave ﬁelds have a persistent structure, with wave maxima present in the forward quadrants of the storm and in the rear quadrant II. This structure is determined by the strong winds and the extended fetch condition in quadrants I and II, as well as by the translation speed of the storm. When a shorter time interval is analyzed (e.g., a 3 h period, when the storm becomes a category 1 hurricane), the structure of the mean wind ﬁeld may differ greatly from the mean ﬁeld calculated with a sufﬁciently longer period; however, the spatial distribution of the wave ﬁeld around the hurricane tends to maintain its typical spatial structure. The use of wind ﬁelds with reduced time variability (e.g., with a 3 h moving average) does not change the structure of the mean wave ﬁeld, but reduces the mean wave height values by up to 10%.


Introduction
Tropical cyclones (TCs) are extreme low-pressure meteorological systems developed over the tropical oceans [1]. They are characterized by intense winds, irregular trajectories and fast translation speed [2]. Due to their vortex-like structure, the TCs present cyclonic winds which rapidly vary in space and time [3], thus modifying the intrinsic wind field and the eye of the tropical cyclone [4,5]. These events have a direct impact on the generation and evolution of the wave field. During the TCs' landfall, the wind and wave patterns become even more complex [4,5].
Several works have been made to study the behavior of the wind, the wave field response to a TC using buoys and altimetry measurements [6][7][8][9][10] and the wind field from parametric models and numerical simulations to analyze the features of the wave field [4,11,12], as well as parametric models of the wind field [11,13]. This information has been used to determine the average features of the wind and wave fields during tropical cyclone conditions. So, based on those results, it has been observed that the strongest winds are located to the right of the eye of the TCs, with respect to its translation movement; whereas, the weakest winds are found to the left of it [6,8,11,14]. This asymmetry emerges as a consequence of the TC curvature [4,[14][15][16]. Moreover, the spatial distribution of the mean significant wave height (Hs) shows large values to the right of the TC with respect to its propagation direction [4,6,9,10] and, when compared to the wind field, the spatial distribution of mean Hs also exhibits a higher asymmetry [10].
A common way to estimate the mean wind field for the analysis of waves during TCs, with a regular time and space distribution, is through the use of global numerical models. Nonetheless, these fields are not correctly solved by the global atmospheric circulation models, mainly because the wind field during a TC is a relatively small vortex [17]. In other studies, more realistic wind fields are built with higher time resolution by combining wind fields from both reanalysis and hurricane models as the HWIND [12]. The wind fields from HWIND are obtained from the NOAA/HRD [18] and are constructed with information obtained from buoys and ships, as well as from TC reconnaissance aircraft. The wind fields have a 6 km × 6 km resolution, thus covering an 8 • × 8 • region around the hurricane center, and also have a time resolution of 3 or 6 h. Another wind model used in the study of hurricanes is the Holland parametric model [19]. This model generates a simple vortex, and the wind field is then generated from a radial pressure profile representation of a TC to obtain a symmetrical wind field which, once the translation velocity is incorporated, approaches an asymmetric field characteristic of a tropical cyclone. By contrast, there are atmospheric models with physical parameterizations and specific schemes for hurricanes, which are initialized so that the initial state of the model's atmosphere contains a realistic cyclonic vortex with features closely resembling the observations [20]. The MM5 (Fifth Generation Mesoscale Model [21,22] is among these models. However, the lack of a refinement scheme for the vortex-following grid makes the high-resolution integration (lower than 5 km grid spacing) unrealistic for periods of a few hours [23]. Another commonly used atmospheric model is the WRF-NMM (Non-hydrostatic Mesoscale Model), by which it is possible to follow a vortex and to obtain a better simulation of the TC. Nonetheless, the corresponding vortex wind profile is generated from an analytic Rankine model [20]. Objectively, to simulate the behavior of the TCs with more exactitude, vortex formation schemes with realistic intensity and positioning are needed since this directly affects the trajectory, wind intensity and pressure errors [24]. On the other hand, to simulate wind waves during extreme events such as the tropical cyclones, it is recommended to use as an atmospheric forcing agent-a wind field with high resolution in time and space [25].
Since a high spatial resolution is required to adequately simulate the wind field during TCs, some novel nesting techniques of mobile grids have been incorporated into the models, thus allowing them to register the evolution of the wind field structure during tropical storms and hurricanes. The HWRF atmospheric model is a good example in which the user can obtain high-resolution wind fields throughout space and time around the TCs, besides from monitoring the cyclone trajectory. The high resolution is essential too if we wish to capture the small-scale structures observed in the surroundings of the cyclone's eyewall, where strong wind velocity gradients are present [17]. This, in turn, may have an important effect on the description of the wave field evolution [25].
Moon et al. [4] provided one of the first attempts to simulate the directional wave spectra under hurricane conditions by using high-resolution winds. In their work, wind data from the NOAA National Hurricane Center (NHC) were used to generate 30 min wind fields, which were used as a forcing to the spectral wave model WAVEWATCH III (WW3). Their results show that the high-resolution wind forcing, as well as the generated wave field, have a dependence on the radius of maximum wind and the cyclone's translation speed. When the translation speed is low, the wave direction is determined by the distance to the center of the storm; whereas, with high translation speed the dominant waves are rather determined by an effect of prolonged wave exposure to the forcing wind (extended fetch). Thus, the waves in the right quadrants of the cyclone are exposed to the wind for a longer period, remaining trapped within this region.
In a recent study, [26] a high spatial resolution model that includes an atmospherewave-ocean coupling has been used. They study the impact of the previously referred interaction in the modeling of the wind-wave field structure during hurricane Bejisa (2014), west of the Indian Ocean, and show the complexity of the wind field structure obtained from their numerical experiment. In the present study, the atmospheric model HWRF is used to obtain realistic wind velocity fields of the tropical cyclone Isaac (2012), with high spatial and time resolution. Based on the numerical results, the spatial structure of the average wind field and the corresponding time variability are analyzed. These winds are used as a forcing to the spectral wave model WW3, so the impact of the TCs' time variability in the evolution of the space-time structure of the wave field can be analyzed. This paper is organized as follows: In Section 2, the general description of the simulated TC is presented together with the main features of the implemented numerical models and the corresponding parameterizations. Within Section 3, model results validation is described through the analyses of the simulation results and the measurements from buoys. Finally, the conclusions are given in Section 5. The wind field for hurricane Isaac was obtained using the atmospheric model for hurricanes HWRF version 3.8. This model is a specialized version of the Weather Research and Forecasting (WRF) model, which is focused on ocean-atmosphere interactions [27]. This model solves the dynamic equations over an Arakawa E-grid using a hybrid sigma vertical coordinate system, also possessing the capability of using data assimilation and specific procedures for the relocation of the main vortex. The coverage and resolution of the grids used for the present implementation are shown in Table 1 and Figure 1. Further details on the features and formulations of the HWRF can be found in Gopalakrishnan et al. [28] and Biswas et al. [27]. Table 1. Characteristics of the grids used for the HWRF model configuration. The fixed grid is centered in the initial position of the tropical cyclone, whereas the mobile grids are centered in the instant position of the tropical cyclone eye.

Domain
Resolution Grid In this work, HWRF is set up to use the physical parameterizations indicated in Table 2. Although the HWRF model is able of coupling ocean-atmosphere systems, neither ocean coupling nor assimilation is considered here. As initial and boundary conditions for the domain D01, historical data from the NCEP FNL dataset (Research Data Archive: NCAR (https://rda.ucar.edu/datasets/ds083.3/dataaccess/ accessed on 22 September 2022) are used, with space and time resolution equal to 1.0 degree and 6-hourly, respectively. The D02 and D03 mobile grids are initialized using the corresponding downscale variables from the parent domain mean while vortex characteristics (i.e., location, size, intensity and structure), provided by NHC (TCvitals files), are incorporated during the run. The hurricane simulation was carried out from 26 August at 00 UTC, and the model run was set for a 120 h period. The results of the simulation were recorded with 15 min intervals. WaveWatch III ® (WW3; [33]) is a third generation model operated by the National Oceanic and Atmospheric Administration (NOAA), together with the National Centers for Environmental Predictions (NCEP). The model WW3 describes the evolution of the wave directional spectrum based on the numerical solution of the spectral action balance equation. Most of the results analyzed here refer to a TC developed over deep waters. Under these conditions, the evolution of the directional spectrum is mainly controlled by three processes: the energy transfer between the atmosphere and the waves (wind input), the energy dissipation by deep water wave breaking (whitecapping) and the resonant quadruplet wave-wave interactions. This is treated in our model configuration using the ST6 package (see Table 3). The governing equation is solved using a regular grid placed over the GoM between 18 • N and 31 • N, and between 79 • W and 98 • W ( Figure 2a). The time steps used in this implementation are shown in Table 3, and they are chosen taking into account the spatial and spectral discretization, so the CFL criteria is fulfilled. The bathymetry data were obtained from GEBCO (General Bathymetric Chart of the Oceans). In this implementation, data from boundaries are not taken into account because the predominant waves over the region during the analyzed period are those generated by the TC. In order to generate the initial conditions, the model is implemented over a 3-week period using wind information from the ERA5 reanalysis database. Then, starting from 26 August at 00 UTC, the model is forced using wind data from the HWRF output. The results of the WW3 are then recorded with a time interval of 15 min.

Characteristic Parameters of the TC
As mentioned above, the main objective of this paper is to determine and study the average features and time variability of the wind and wave fields during the occurrence of a TC numerical analysis. Hence, it is necessary to confirm that the results provided by the model are adequate by means of a qualitative and quantitative inspection. The validation is performed with information from different databases. From the NHC and HURDAT datasets, the following parameters of the TC were obtained: minimum pressure in the eye of the TC (P min ), intensity of the maximum sustained wind (V max ) and the radius of maximum wind (R max ). All these data are available at the 6 h sampling interval. From the operative buoys of the NDBC (National Data Buoy Center), the wind speed at 10 m, U 10 , and the significant wave height, Hs, were used.

Structure and Trajectory
Historically, V max and P min , as well as the position of the storm center are considered as key parameters to describe the TC and its evolution. These parameters are also used to evaluate the atmospheric model's ability to reproduce the TC's main features. Results from domain 3 (see Table 1) simulated with the HWRF model are compared to the trajectory ( Figure 2a, blue line) and characteristics reported by NHC (black line). It is possible to observe that both trajectories are similar. Nonetheless, the storm translation speed, V t simulated with the HWRF model is higher than the observations. The V t values computed by the HWRF and those reported by the NHC are shown in Figure 2d. A similar behavior on the trajectory and translation speed of hurricane Isaac is reported by Curcic et al. [37]. The V max simulated by the model (Figure 2b) shows high variability along the TC trajectory, particularly during the tropical storm stage, whereas those values from the NHC are typically smooth. During the first 18 h of simulation, the V max computed with the HWRF is slightly higher than reported by the NHC. According to the numerical results, the hurricane category is achieved during 28 August at 18:30 UTC (segmented line in Figure 2b), reaching maximum winds of 39.6 ms −1 . In contrast, according to the NHC report, the category of hurricane was achieved during 28 August, at 12:00 UTC presenting V max of 33.4 ms −1 with an increasing tendency. This tendency is well reproduced by the model. During the tropical storm stage, the P min values determined by the HWRF model are very similar to those reported by the NHC (Figure 2c). A clear decreasing tendency is observed, which induces an intensification of the storm and further transition from tropical storm to hurricane. Once this category is reached, the P min values reported by the NHC tend to decrease, as well as those from the HWRF. As for the V t values computed with the HWRF (Figure 2d, solid blue line), very high variability is observed in contrast to the NHC time series. Additionally, the TC simulated by the HWRF model travels about 2 ms −1 faster than those reported by NHC. However, it is possible to observe an interesting feature within the values given by the model: the decreasing tendency of the instantaneous values of V t is highly correlated to the reported information.

Fixed Observations
Recordings of U 10 and Hs from two NDBC fixed buoys (42001 and 42003) encountered in the surroundings of the TC trajectory were used for additional validations. In order to evaluate the performance of the model, three statistical parameters were used: The determination coefficient, where C xy indicates the covariance between the observed data, x and the model results, y; whereas, σ x and σ y indicate their corresponding standard deviation.
The root mean squared error (RMSE), and whereȳ is the average value obtained from the model results, andx is the average value from the observation.
Due to the difference between the observed and computed TC track, the position of the buoy relative to the storm eye may significantly vary. So, the comparison between wind and wave data from buoys 42001 and 42003 was made against the model results for the same relative positions as the one reported by the NHC. The buoy 42001 is located near the central GoM, at 89 • 39 25 W and 25 • 56 31 N (see Figure 2a). At some point, the shortest distance between this buoy and the hurricane's eye was about 200 km, over their left quadrants. In Figure 3a, we show the comparison between the wind intensity observed with the buoy (black line) and those values computed by HWRF (blue line). The model adequately reproduces the observed values until 28 August, then the model results tend to overestimate the values measured by the buoy. An objective comparison between the observations and the model (Table 4) yields an R 2 = 0.75, with a RMSE of 0.98 ms −1 and a bias of S = −0.22 ms −1 . The overestimation in the wind intensity is reflected in the larger Hs values computed by WW3 with respect to those measured by the buoy (see Figure 3b). In this case, a relatively high R 2 of 0.85 is computed with an RMSE of 0.5 m and an S = 0.32 m.  The buoy 42003 is located at 85 • 36 53 W and 25 • 55 31 N (Figure 2a). The shortest distance between this buoy and the hurricane's eye is approximately 50 km, over its left quadrants. Figure 3c shows the comparison between the buoy (black line) and the HWRF model (blue line) data for the wind intensity values. The closest approach between the buoy and the eye of the TC was around 16:00 UTC of 27 August, when the buoy reported a wind intensity minimum. During the following 4 h, the buoy shows increasing wind intensity, until reaching a maximum of 26 ms −1 , likely associated with the storm intensification. Then, U 10 values from the buoy tend to decrease to about 16 ms −1 ; whereas those from the HWRF model are clearly overestimated (up to 10 ms −1 ). As observed from the modeled U 10 time series, higher variability is depicted in comparison to the buoy observations, probably due to differences in the calculation of the eye size at those specific times. According to HWRF results, the intensification of the registered winds is noticeably higher than the observed, but the decreasing tendency is similar. The values of U 10 computed by HWRF at the vicinity of the storm eye indicate a model tendency to overestimate its magnitude, which is also consistent with the translation speed overestimation [38]. This, in turn, would explain the tendency of the WW3 to produce much higher Hs values at the relative position of the buoy. The objective comparison between the model results and observations is not as good as that of buoy 42001. As for the HWRF values, the computed determination coefficient is R 2 = 0.2, having a RMSE of 5.4 ms −1 with a S = 3.2 ms −1 . Finally, for the Hs computed with the WW3, R 2 = 0.47 with a RMSE of 1.7 m and a bias of 1.5 m. It is clear that the model does not reproduce the exact evolution of the TC, although it can be seen that the model results tend to reproduce the general pattern observed in the measurements.

Effects of the Time Resolution of the Wind Field on the Wave Field
The numerical results obtained from the HWRF model show a very high temporal variability, characterized by abrupt changes in the spatial structure of the wind field. According to Janssen [39], good quality of wind data is essential to improve the performance of the wave models and, when it comes to extreme events such as a TC, this fact becomes more relevant [40]. With poor time and space resolution, the wind models may underestimate the wind behavior, particularly at the storm center [41], thus inefficiently representing the wave field evolution. Tolman and Alves [25] also emphasize the relevance of a high time resolution in the characterization of the wind field. In order to analyze the effect of the temporal variability of the wind field on the wave spatial structure computed by WW3, two experiments are performed: one with the wind fields originally calculated by the HWRF model and the other with the same original wind fields but filtered using a 3 h running average. Both wind fields are fed to the wave model with a time resolution of 15 min. The wind data are grouped as a function of P min [10] and divided into two stages: moderate storm, when P mim > 980 mbar; and the intensification stage, when P min ≤ 980 mbar (Figure 2c). Prior to the calculation of average values, the instantaneous fields are normalized with their corresponding maximum value. The wind and wave fields are analyzed following the standardization procedure proposed by Young [14]. The procedure consists of normalizing the space dimension with the R max computed for each field and rotating it along the TC direction of propagation. This procedure allows the positive y-axis to be aligned with the direction of advancement, thus permitting a spatial discretization of the TC parameters into quadrants. The front quadrants of the TC are labeled as I (right) and IV (left), whereas the rear quadrants are identified as II (right) and III (left).

Average Wind and Hs Fields
In this section, the full-resolution wind fields computed by HWRF and the corresponding wave fields calculated by WW3 are analyzed. Figure 4 shows the mean wind field considering the full time period (18 h), moderate storm (9 h) and storm intensification stage (9 h) (Figure 4a,c,e, respectively). In general, the wind field presents a clearly-defined storm eye, which is typical from a well-developed storm. Generally, in the 18 h average (Figure 4a) and moderate storm stage (Figure 4c) an asymmetric behavior of the wind field is observed, with the quadrants I and II presenting the higher mean values. In addition, the wind maximum is present to the right of the eye, with strong winds distributed over I and II quadrants, between 1 and 2R max ; these strong winds are readily apparent in these quadrants and between 1 and 2R max during the 9 h moderate storm stage analyzed (Figure 4c). Large intensity values are obtained in quadrant IV during the full period 18 h average. It is possible to observe large mean wind values in regions further away from the storm center distance large than 1R max in quadrants I and II. This behavior is reported by Esquivel-Trava et al. [9] in quadrant II, and they associated it with the presence of concentric walls in the hurricane. The spatial pattern observed in Figure 4c could be associated with the spatial variability of the spatial structure of the hurricane simulated by HWRF, as well as to the possible presence of concentric walls, although a more detailed analysis is needed to further confirm this statement.
Furthermore, the left quadrant (III), as observed in the 18 h average and during the moderate storm stage (Figure 4c), is characterized by the presence of lower mean wind values. This asymmetry is caused by the cyclonic motion of the storm. Particularly, over the right quadrants, the translation speed of the TC is added to the wind speed near the storm eye vicinity, which allows wind intensification over this region [4,42]. This mean pattern is similar to previous observations from Young [14], Esquivel-Trava et al. [9] and Tamizi and Young [10] using in situ measurements, and it is also similar to simulations from Moon et al. [4], Hu and Chen [43] and Mora Escalante [44]. During the intensification stage (Figure 4e), the asymmetry of the wind field is reduced and moderate mean winds are present over the left quadrant (III and IV) at 1R max .
In Figure 4b,d,f, the Hs mean fields are presented as obtained from the corresponding averaged wind field (Figure 4a [11] and Tamizi and Young [10]. In all panels of Figure 4, the mean direction of the wind (blue arrows) and waves (white arrows) are included. Over the right quadrants (I and II), where the wind maxima are observed, wind and waves have nearly the same direction. The difference between them increases in the regions further away from the storm center. Over these quadrants, near the storm eye, the wind waves tend to be local, i.e., actively forced by winds [45]; although Holthuijsen et al. [8] mention that young swell traveling in the same direction as the storm may be present. Away from the regions of maximum winds, the difference between the mean direction of the wind and waves increases, mostly because the waves over these regions are predominantly non-local (swell). As for the left quadrant (III and IV), the waves propagate following a straight line to the right of the wind direction, and they form an angle near to 90 • . This quadrant is dominated by previously generated waves and has features different from the local wind [45]. Holthuijsen et al. [8] define this kind of pattern as crossing swell (quadrant IV) and opposing swell (quadrant III). In general, these features of the wind and waves' direction are in accordance with observations [8,9] and numerical model results [4]. It is worth noting that an important feature of the averaged wave field is the position of the minimum Hs, which lies approximately 0.5R max south (according to the storm track) of the eye location. This feature is similar to the Hs spatial structure shown by Moon et al. [4] in their Figure 8c, although its cause is not discussed. We can argue that this feature is produced by the translation speed of the storm and the response time of the wave field to the forcing; however, its precise cause needs to be investigated.

Mean Fields of Filtered Winds and the Corresponding Wind Waves
Tolman et al. [41] suggest using data with high spatial resolution in order to get better estimations of the wave field, and particularly to avoid wind underestimation near the storm core. Filtered data with a 3 h running mean were used in order to analyze the effect of low time resolution in U 10 and Hs, although the 15 min interval was maintained. In general, the filtered wind fields (Figure 5a,c,e) show a similar spatial pattern to that obtained when the original fields (Figure 4a,c,e) are used. The average values of the filtered wind fields are closer to one, which is the result of their lower spatial variability and slightly lower instantaneous maxima (about 3 ms −1 ) than those calculated using the original HWRF results. The mean Hs fields calculated with the filtered wind (Figure 5b,d,f) have a very similar spatial structure to that of the mean fields calculated with the original wind fields. However, the Hs values calculated with the latter wind data tend to be higher than those calculated with the filtered wind data. This difference is most obvious during the intensification stage of the storm when the Hs computed by the model reaches values greater than 9.5 m. Figure 6 shows the average spatial distribution of the differences between the 97.5 percentile values of the Hs calculated using the original winds minus those calculated using the filtered winds. The differences are reported in percentages with respect to the values calculated using the original winds. It is possible to observe that the large values of Hs computed using the original winds are up to 15% higher than those calculated with the filtered wind data. The largest differences are observed in the region around the storm eye wall, mainly in quadrants II and quadrants I and IV. As will be shown in Section 3.3, it is in this region that the greatest temporal wind variability is observed. This type of result is in agreement with the fact that high wind variability influences high Hs values [46].

Particular Case: Change from Tropical Storm to Hurricane
In order to assess the effect of wind variability on the wave field structure over a short period, a 3 h interval centered on the date when the storm reaches hurricane status (blue dashed line in Figure 2b) is analyzed. It is worth noting that this is a period of sudden increase in maximum wind intensity and translation speed of the TC (see Figure 2b,d).
In this case, the wave fields are calculated using both the original and filtered wind fields. Figure 7a,c, shows the average wind field; whereas, Figure 7b,d shows the average Hs fields. As in the previous examples, each field is normalized with its own maximum value. The average wind field computed with original data (Figure 7a) shows a very different structure with respect to the one computed over a larger period (see Figure 4a). This exhibits the high time-space variability of the wind field, only visible when the TC is reproduced with enough time and space resolution. It is possible to observe the presence of intense winds both on the right and left quadrants. In fact, the maximum winds are present in quadrants III and IV, at a distance of 1R max . This behavior differs much with respect to the typical structure obtained from the parametric hurricanes described in the literature. The mean Hs field (Figure 7b) shows a spatial structure similar to that computed with a larger set of data (see, e.g., Figure 4b). The higher Hs average values are seen in the right quadrants of the TC, at distances greater than R max . Large values are also observed in quadrants III and IV, which are associated with the wind maxima for the same regions (Figure 7a). By contrast, using the filtered data, U 10 and Hs average fields (Figure 7c,d, respectively) show the same patterns as those seen with the original wind data. Nonetheless, the U 10 mean fields have normalized values greater than those obtained with the unfiltered data (Figure 7a). This is due to the maximum wind values used to normalize the corresponding fields which are smaller. This behavior is not observed in the averaged Hs fields computed with the smoothed U 10 field, because a reduced time variability of the U 10 field will induce lower Hs values.
The comparison between average fields computed over a relatively long period against those computed over short periods indicate a high variability in the TC wind field. Although the structure of the U 10 field within the TC may show a similar behavior to the one described in Figure 7, the average field tends to be similar to that of Figure 4 if it is analyzed through a sufficiently long period. Nonetheless, the mean Hs field shows a rather persistent spatial structure, with maximum values located mainly in the front (I and IV) and second quadrants. This behavior is associated with an extended fetch condition over the right quadrants [4,11], which implies that the waves in this sector are forced by the wind over a longer period. In this case, the storm translation speed, V t , plays an important role in the generation and development of Hs, and also on the TC waves distribution [10]. For this condition to be present, V t must be comparable to the wave group velocity. This process allows the wind, by constantly blowing over the right side of the TC, to favor wave trapping inside the cyclone and to generate larger local waves [11].
Based on the direction of the spectral peak values computed by the wave model, it is possible to obtain the waves' average group velocity C g at the right quadrants of the TC (I and II), and to make a comparison with the V t computed using the HWRF output. These results are shown in Figure 8. On average, V t (7.2 ms −1 ) and C g (8.3 ms −1 ) are nearly equal. Particularly, during the period in which the TC reaches the hurricane category (Figure 7), large wave heights are observed covering-almost totally-quadrant I, regardless of the fact that intense winds are present in quadrants III and IV. The presence of high Hs values in quadrant II is mostly associated with the intensity of U 10 , as well as to the extended fetch condition; whereas, the high values in the front quadrants seem to be linked to the wind intensity and to the greatest C g values with respect to the mean V t .

Temporal Variability of the Wind and Significant Wave Height Fields
The temporal variability of the wind and significant wave height fields are determined by computing the standard deviation normalized by its respective average field. The normalization process is considered a measure of the data dispersion, which allows the determination of the variability between them, in relation to their mean value [47]. The normalized standard deviation values are expressed in percentage. Three important variability regions are present within the 18 h average wind field (Figure 9a), as well as during the moderate storm ( Figure 9c) and intensification (Figure 9e) stages. The first region is associated with low variability (less than 27%) and is located in the right quadrants (I and II), at a distance of about 1R max , thus belonging to the zone characterized by the maximum winds. During the moderate stage, the region with a low variability values is spatially enhanced to 3R max in quadrants I and IV, but the variability percentage reduces from 27% to values even lower than 7%. During the intensification stage (Figure 9e), the variability percentage rises up to 27%, thus being very close to the value obtained when a total of 18 h is considered (Figure 9a).
The second region with characteristics variability is located at quadrant III, away from the eye of the TC. Using the 18 h scheme, the variability oscillates between 21-48%; when the storm is moderate (Figure 9c), the variability reduces to 27% and, when the storm is intensified, it increases to 34%. Finally, the third variability region is around the eye of the TC. This area is defined by Weatherford and Gray [48] as the inner core of the storm, reaching to R max (∼111 km), and presenting the higher variability percentage.
As for the 18 h interval, the variability ranges from 41% to 109%; during the moderate stage, it is slightly reduced together with its spatial distribution and, when the storm is intensified, there is a considerable rise in the variability which reaches a maximum of 109%. At every stage of the storm, the minimum variability region is located in quadrants I and II, so this region is representative of the average features. Very low variability values at these regions is indicative of a nearly-constant wind direction. The third region characterized by high variability is also seen in the vicinity of the storm eye, which is likely to be generated by strong velocity gradients as well as rapid variations in the wind direction at the region adjacent to the cyclone's eyewall [40].
The U 10 variability field computed with the filtered data is presented in Figure 9b,d,f. As in the unfiltered case, three main regions are distinguished but with the expected variability reduction due to the filtering process. Particularly, during the 18 h interval ( Figure 9b) the minimum variability region located right of the TC eye is maintained below 14%; that value is notably reduced to almost zero during the moderate storm stage (Figure 10d), and extends up to 3R max over the four quadrants of the storm. During the intensification stage (Figure 9f), the variability is slightly higher than that of the moderate storm stage, with approximately 27%. The second variability region, located at quadrant II, has percentage values similar to those computed using the unfiltered data (27-41%). As for the vicinity of the TC eye, the variability percentage was reduced, ranging from 41% to 68%. Figure 10 shows the Hs time variability field using the original and filtered data for the different storm stages. As in the wind fields from Figure 9, a region of minimum variability is observed to the right of the storm center during the three stages. The variability level fluctuates from 11% during the 18 h interval to 16% during the moderate storm stage, then reduces to 5% during the intensification stage. In addition, during the moderate stage, this region of low variability expands over the four quadrants of the TC, and then it is spatially reduced during the intensification stage. Unlike the wind field, the Hs variability level does not show high percentage values around the storm eye. As seen in Figure 10e, the greatest variability is found beyond R max . Figure 10b,d,f, shows the Hs variability computed with wind filtered data. The results are similar to those obtained with the original wind data, specifically, those low-variability values in quadrants I and IV, as well as an expansion directed away from the center. Additionally, the level of Hs variability computed using the unfiltered wind data is somewhat maintained.
When analyzing the 3 h interval in which the storm becomes a category 1 hurricane, the wind variability (Figure 11a,c) is observed to be similar as in the distinct storm stages (Figure 9a,c,e). The high variability is located around the TC center and extends to R max . The lower variability percentage (14%) is still present in quadrant III. The time variability of Hs (Figure 11b,d) during this period shows values below 10% within 1R max , mainly over quadrant I and III. The highest variability (∼10%) is observed in quadrant II, between 2 and 3R max . In general, Hs presents slight variations which become more evident when reduced time intervals are considered.

Discussion
Our results show that the atmospheric model HWRF has the capability to reproduce the trajectory, V max , and P min of the TC Isaac in agreement with the NHC reports. However, the values of V t simulated with the HWRF (Figure 2d, solid blue line) show high variability in contrast with the NHC observations. The variations of V t with respect to the NHC observations could be attributed to the appropriate simulation of the internal parameters of the storm, particularly of V max , which has been reported to have a positive correlation with V t [49,50]. In addition, the good estimation of the internal parameters of the TC could be a result of using the HWRF physically configured for hurricanes [20]. The latter allows the generation of high time resolution wind data, unlike the HWIND wind fields used by NHC, which have a time resolution of 3 or 6 h.
The variability of U 10 (see Figure 3c) from the model time series is larger with respect to the buoy observations, probably due to differences in the calculation of the eye size at those specific times. According to the HWRF results, the further intensification of the registered wind is faster than the observed, but the decreasing tendency is similar. The values of U 10 (Figure 3c) computed by the HWRF at the vicinity of the storm eye indicate a model tendency to overestimate its magnitude, which is consistent with the translation speed overestimation [38]. By eliminating the 3 h variability of the original wind and Hs fields, a reduced spatial distribution of variability level was observed when compared to the original data, together with the reduction in the Hs values. In addition, when averaging over a short time lapse (i.e., 3 h, when the storm had a hurricane category as in Figure 7), the maximum winds are reproduced in quadrants III and IV, at a distance of 1 R max . This behavior is importantly different with respect to the typical structure obtained from the parametric hurricanes described in the literature. The mean Hs field shows a very persistent structure, with the maximum values located mainly in the front (I and IV) and second quadrants associated with an extended fetch condition. The comparison between the fields computed with a relatively long period, against those computed with short periods indicate a high variability in the TC wind field. Moreover, in Figure 10, the Hs presents slight variations which become more evident when reduced time intervals are considered. The observed variability in the wind field and Hs shows the importance of using high time resolution wind fields as forcing agents of the wind wave models, as was already suggested by Ponce and Ocampo-Torres [46], Tolman and Alves [25] as well as Collins et al. [51].

Conclusions
In this work, space-time variability of the wind field under tropical cyclone conditions was determined using numerical simulations of realistic wind fields with high space and time resolution. The atmospheric model HWRF was used to generate the wind fields and was configured to reproduce the evolution of hurricane Isaac during its passage through the GoM in the last few days of August 2012. Based on the realistic wind field obtained with the HWRF model, the wind wave field was computed using the third generation spectral model WaveWatch III. The main objective of this work is to determine the extent to which the high time and space variability of the wind field impacts on the evolution of the wind waves field during the development of these extreme events. The effect of the time resolution is analyzed by using the HWRF-generated data with an interval of 15 min (original data), as well as a 3 h filtered version of the same database (keeping ∆t = 15 min). Both datasets are used to simulate the wave field. Although reproducing the exact characteristics of hurricane Isaac is not intended here, a comparative analysis of the wind and wave fields is carried out using observations and the model-generated results in order to validate the simulated output. The numerical results indicate that the structure of the wind field has a high time and space variability, although it has an overall structure very similar to those reported for parametric hurricanes, as long as the number of averaged hours is sufficiently high. If a shorter time interval is analyzed (as in a 3 h study case, in which period the storm becomes a category 1 hurricane), the structure of the mean wind field may importantly differ from the mean field computed with a sufficiently longer period.
The wind waves response shows that, although the extended fetch is present in the TC, it becomes more important during the intensification stage of the storm, since it allows a considerable wave growth to the right of the TC when the wind maximum is on the left of the TC. For this condition to be developed, the V t and C g values must be of the same order; here, V t reached 7.2 ms −1 and C g was of 8.3 ms −1 . This phenomena causes the waves located to the right of the storm to be trapped for a prolonged period and distance, with the presence of previously generated swell and local wind waves. While to the left of the TC the waves are local and not trapped, wave growth is also inhibited because the wind blows for a shorter period. The implications of the obtained results are the following: • The wave field generated by the WaveWatch III model does not respond instantly to the variations in the wind field structure. More studies and specialized measurements are required to determine if the described behavior is either physically realistic, or is a deficiency of the wave model.

•
The results indicate that the structure of the wave field is strongly determined by the extended fetch process, which leads to the occurrence of high waves over quadrants I and II of the storm. The high wave values over the frontal quadrant IV are related to a translation speed of the TC smaller than the group velocity of the waves generated in quadrants I and II. • The wind field generated by parametric hurricanes successfully reproduces the wave field structure because the mean wind field has a tendency to produce a typical asymmetric spatial structure. Funding: This research was funded by Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE).