Bibliography - Morris A Bender

1. Bender, Morris A., I Ginis, Robert E Tuleya, B Thomas, and Tim Marchok, December 2007: The operational GFDL coupled hurricane–ocean prediction system and a summary of its performance. Monthly Weather Review, 135(12), doi:10.1175/2007MWR2032.1.
   [ Abstract ]

   The past decade has been marked by significant advancements in numerical weather prediction of hurricanes, which have greatly contributed to the steady decline in forecast track error. Since its operational implementation by the U.S. National Weather Service (NWS) in 1995, the best-track model performer has been NOAA’s regional hurricane model developed at the Geophysical Fluid Dynamics Laboratory (GFDL). The purpose of this paper is to summarize the major upgrades to the GFDL hurricane forecast system since 1998. These include coupling the atmospheric component with the Princeton Ocean Model, which became operational in 2001, major physics upgrades implemented in 2003 and 2006, and increases in both the vertical resolution in 2003 and the horizontal resolution in 2002 and 2005. The paper will also report on the GFDL model performance for both track and intensity, focusing particularly on the 2003 through 2006 hurricane seasons. During this period, the GFDL track errors were the lowest of all the dynamical model guidance available to the NWS Tropical Prediction Center in both the Atlantic and eastern Pacific basins. It will also be shown that the GFDL model has exhibited a steady reduction in its intensity errors during the past 5 yr, and can now provide skillful intensity forecasts. Tests of 153 forecasts from the 2004 and 2005 Atlantic hurricane seasons and 75 forecasts from the 2005 eastern Pacific season have demonstrated a positive impact on both track and intensity prediction in the 2006 GFDL model upgrade, through introduction of a cloud microphysics package and an improved air–sea momentum flux parameterization. In addition, the large positive intensity bias in sheared environments observed in previous versions of the model is significantly reduced. This led to the significant improvement in the model’s reliability and skill for forecasting intensity that occurred in 2006.

2. Bender, Morris A., and I Ginis, 2000: Real-case simulations of hurricane-ocean interaction using a high-resolution coupled model: effects on hurricane intensity. Monthly Weather Review, 128(4), 917-946.
   [ Abstract PDF ]

   In order to investigate the effect of tropical cyclone-ocean interaction on the intensity of observed hurricanes, the GFDL movable triply nested mesh hurricane model was coupled with a high-resolution version of the Princeton Ocean Model. The ocean model had 1/6º uniform resolution, which matched the horizontal resolution of the hurricane model in its innermost grid. Experiments were run with and without inclusion of the coupling for two cases of Hurricane Opal (1995) and one case of Hurricane Gilbert (1988) in the Gulf of Mexico and two cases each of Hurricanes Felix (1995) and Fran (1996) in the western Atlantic. The results confirmed the conclusions sugggested by the earlier idealized studies that the cooling of the sea surface induced by the tropical cyclone will have a significant impact on the intensity of observed storms, particularly for slow moving storms where the SST decrease is greater. In each of the seven forecasts, the ocean coupling led to substantial improvements in the prediction of storm intensity measured by the storm's minimum sea level pressure. Without the effect of coupling the GFDL model incorrectly forecasted 25-hPa deepening of Gilbert as it moved across the Gulf of Mexico. With the coupling included, the model storm deepened only 10 hPa, which was much closer to the observed amount of 4 hPa. Similarly, during the period that Opal moved very slowly in the southern Gulf of Mexico, the coupled model produced a large SST decrease northwest of the Yucatan and slow deepening consistent with the observations. The uncoupled model using the initial NCEP SST's predicted rapid deepening of 58 hPa during the same period. Improved intensity prediction was achieved both for Hurricanes Felix and Fran in the western Atlantic. For the case of Hurricane Fran, the coarse resolution of the NCEP-SST analysis could not resolve Hurricane Edouard's wake, which was produced when Edouard moved in nearly an identical path to Fran four days earlier. As a result, the operational GFDL forecast using the operational SST's and without coupling incorrectly forecasted 40-hPa deepening while Fran remained at nearly constant intensity as it crossed the wake. When the coupled model was run with Edouard's cold wake generated by imposing hurricane wind forcing during the ocean initialization, the intensity prediction was significantly improved. The model also correctly predicted the rapid deepening that occurred as Fran began to move away from the cold wake. These results suggest the importance of an accurate initial SST analysis as well as the inclusion of the ocean coupling, for accurate hurricane intensity prediction with a dynamical model. Recently, the GFDL hurricane-ocean coupled model used in these case studies was run on 163 forecasts during the 1995-98 seasons. Improved intensity forecasts were again achieved with the mean absolute error in the forecast of central pressure reduced by about 26% compared to the operational GFDL model. During the 1998 season, when the system was run in near-real time, the coupled model improved the intensity forecasts for all storms with central pressure higher than 940 hPa although the most significant improvement (~60%) occurred in the intensity range of 960-970 hPa. These much larger sample sets confirmed the conclusion from the case studies, that the hurricane-ocean interaction is an important physical mechanism in the intensity of observed tropical cyclones.

3. Wu, C-C, Morris A Bender, and Yoshio Kurihara, 2000: Typhoon forecast with the GFDL Hurricane Model: Forecast skill and comparison of predictions using AVN and NOGAPS global analyses. Journal of the Meteorological Society of Japan, 78(6), 777-788.
   [ Abstract PDF ]

   A hurricane model developed at GFDL, NOAA, was combined with each of AVN and NOGAPS global analyses to construct typhoon prediction systems GFDS and GFDN, respectively. The GFDS system performed 125 (178) forecast experiments for 16 (24) storms in the western North Pacific basin during 1995 (1996). It exhibited considerable skill in the forecast of tropical cyclone tracks. The average forecast position errors as 12, 24, 36, 48 and 72 h in 1995 (1996) were 95 (108), 146 (178), 193 (227), 249 (280), and 465 (480) km. The improvement with GFDS in the typhoon position forecast over CLIPER was roughly 30%. The reduction of position errors in both average and standard deviations indicates superior forecast accuracy and consistency of GFDS, although there existed systematic northward bias in the forecast motion at low latitudes. On the other hand, intensity forecast was not satisfactory, showing a tendency to overpredict weak storms and underpredict strong storms, similar to the tendency in the Atlantic Two sets of forecasts performed in the 1996 season, the one by GFDS and the other by GFDN, were compared with each other. Forecast skills of the storm position with the two systems were comparable. However, the two forecast positions tended to be systematically biased toward different directions. As a result, when the two forecasts were averaged, the mean error was 10% smaller than that of each forecast. Also, overall improvement in track forecast was obtained in supplemental experiments in which individual forecasts were corrected for systematic biases. Though systematic biases is not steady, there may be ways to utilize it for improvement of tropical cyclone forecasts.

4. Kurihara, Yoshio, Robert E Tuleya, and Morris A Bender, 1998: Application and improvement of the GFDL Hurricane Prediction System In Research Activities in Atmospheric and Oceanic Modelling, WMO/TD No. 865, Geneva, Switzerland, World Meteorological Organization, 5.31.
5. Kurihara, Yoshio, Robert E Tuleya, and Morris A Bender, 1998: The GFDL Hurricane Prediction System and its performance in the 1995 hurricane season. Monthly Weather Review, 126(5), 1306-1322.
   [ Abstract PDF ]

   The Geophysical Fluid Dynamics Laboratory (GFDL) Hurricane Prediction System was adopted by the U.S. National Weather Service as an operational hurricane prediction model in the 1995 hurricane season. The framework of the prediction model is described with emphasis on its unique features. The model uses a multiply nested movable mesh system to depict the interior structure of tropical cyclones. For cumulus parameterization, a soft moist convective adjustment scheme is used. The model initial condition is defined through a method of vortex replacement. It involves generation of a realistic hurricane vortex by a scheme of controlled spinup. Time integration of the model is carried out by a two-step iterative method that has a characteristic of frequency-selective damping. The outline of the prediction system is presented and the system performance in the 1995 hurricane season is briefly summarized. Both in the Atlantic and the eastern Pacific, the average track forecast errors are substantially reduced by the GFDL model, compared with forecasts by other models, particularly for the forecast periods beyond 36 h. Forecasts of Hurricane Luis and Hurricane Marilyn were especially skillful. A forecast bias is noticed in cases of Hurricane Opal and other storms in the Gulf of Mexico. The importance of accurate initial conditions, in both the environmental flow and the storm structure, is argued.

6. Bender, Morris A., 1997: The effect of relative flow on the asymmetric structure in the interior of hurricanes. Journal of the Atmospheric Sciences, 54(6), 703-724.
   [ Abstract PDF ]

   Asymmetric structure of tropical cyclones simulated by the Geophysical Fluid Dynamics Laboratory high-resolution triply nested movable-mesh hurricane model was analyzed. Emphasis was placed on the quasi-steady component of the asymmetric structure in the region of the eyewall. It was found that the asymmetry was primarily caused by the relative wind, that is, the flow entering and leaving the storm region relative to the moving storm. A set of idealized numerical experiments was first performed both with a constant and a variable Coriolis parameter (f) and the addition of basic flows that were either constant or sheared with height. Analysis was then made for one case of Hurricane Gilbert (1988) to demonstrate that the quasi-steady asymmetric structure analyzed in the idealized studies could be identified in this real data case. Vorticity analysis in the variable f experiment indicated that quasi-steady asymmetries resulted in the eyewall region through the effect of vorticity advection due to differences between the beta gyre flow in the lower free atmosphere and the storm motion. This was roughly matched with a persistent area of divergence and vorticity compression in the lower free atmosphere ahead of the storm and enhanced convergence and vorticity stretching to the rear. An asymmetric structure in the upward motion and accumulated precipitation, when averaged over a sufficiently long period of time, exhibited a corresponding maximum in the eyewall's rear quadrant. With the addition of an easterly basic flow, a pronounced change in the asymmetry of the time-averaged boundary layer convergence resulted, with maximum convergence located ahead of the storm. However, the asymmetries in the average vertical motion in the middle troposphere and accumulated precipitation were more affected by the convergence field in the lower free atmosphere produced by the relative flow there. The relative flow depended on both the basic and beta gyre flow. With the addition of an easterly vertical shear to the easterly basic flow, the storm moved faster than the lower-level winds, and strong relative wind was from the front to the rear in the lower free atmosphere and from the opposite direction in the outflow layer aloft. As a result, the upward motion was significantly increased in the front of the storm and reduced in the rear, and the precipitation maximum shifted to the left front quadrant Overall, analysis results suggest that the flow relative to the storm motion is an important factor contributing to the formation of quasi-steady asymmetries in the convergence and vertical motion fields, as well as in the mean precipitation pattern of tropical cyclones.

7. Bender, Morris A., C-C Wu, M A Rennick, and Yoshio Kurihara, 1997: Comparison of the GFDL Hurricane Model prediction in the Western Pacific using the NOGAPS and AVN Global Analysis In 22nd Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 615-616.
8. Kurihara, Yoshio, Morris A Bender, and Robert E Tuleya, 1997: For hurricane intensity forecast: Formulation of a new initialization method for the GFDL Hurricane Prediction Model In 22nd Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 543-544.
9. Kurihara, Yoshio, Robert E Tuleya, and Morris A Bender, 1997: Improvement of the GFDL Hurricane Prediction System In CAS/JSC Working Group on Numerical Experimentation - Research Activities in Atmospheric & Oceanic Modelling, Report No. 25, WMO/TD-No. 792, Geneva, Switzerland, World Meteorological Organization, 5.22.
10. Bender, Morris A., Robert E Tuleya, Yoshio Kurihara, and S J Lord, 1996: Results of the operational GFDL hurricane model at NCEP In 11th Conference on Numerical Weather Prediction, Boston, MA, American Meteorological Society, 393-395.
11. Ginis, I, and Morris A Bender, 1996: Coupled tropical cyclone-ocean modeling In Research Activities in Atmospheric and Oceanic Modelling, CAS/JSC Working Group on Numerical Experimentation, Report No. 23 WMO/TD No. 734, World Meteorological Organization, 9.10-9.11.
12. Kurihara, Yoshio, Robert E Tuleya, and Morris A Bender, 1996: Simulation studies of tropical cyclones In Research Activities in Atmospheric and Oceanic Modelling, CAS/JSC Working Group on Numerical Experimentation, Report No. 23 WMO/TD No. 734, World Meteorological Organization, 5.17-5.18.
13. Tuleya, Robert E., Morris A Bender, and Yoshio Kurihara, 1996: Prediction of hurricane landfall using the GFDL model In 11th Conference on Numerical Weather Prediction, Boston, MA, American Meteorological Society, 407-408.
14. Bender, Morris A., 1995: Numerical study of the asymmetric structure in the interior of tropical cyclones In 21st Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 600-602.
15. Kurihara, Yoshio, Morris A Bender, and Robert E Tuleya, 1995: Performance evaluation of the GFDL Hurricane Prediction System in the 1994 hurricane season In 21st Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 41-43.
16. Kurihara, Yoshio, Morris A Bender, Robert E Tuleya, and R Ross, 1995: Improvements in the GFDL Hurricane Prediction System. Monthly Weather Review, 123(9), 2791-2801.
    [ Abstract PDF ]

    The hurricane model initialization scheme developed at GFDL was modified to improve the representation of the environmental fields in the initial condition. The filter domain defining the extent of the tropical cyclone in the global analysis is determined from the distribution of the low-level disturbance winds. The shape of the domain is generally not circular in order to minimize the removal of important nonhurricane features near the storm region. An optimum interpolation technique is used to determine the environmental fields within the filter domain. Outside of the domain, the environmental fields are identical to the original global analysis. The generation process of the realistic and mode-compatible vortex has also undergone some minor modifications so that reasonable vortices are produced for various data conditions. The upgraded hurricane prediction system was tested for a number of cases and compared against the previous version and yielded an overall improvement in the forecasts of storm track. The system was run in an automated semioperational mode during the 1993 hurricane season for 36 cases in the Atlantic and 36 cases in the eastern Pacific basin. It demonstrated satisfactory skill in the storm track forecasts in many cases, including the abrupt recurvature of Hurricane Emily in the Atlantic and the landfall of Hurricane Lidia onto the Pacific coast of Mexico.

17. Bender, Morris A., I Ginis, and Yoshio Kurihara, 1993: Numerical simulations of hurricane-ocean interaction with a high-resolution coupled model In 20th Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 563-566.
    [ PDF ]
18. Bender, Morris A., I Ginis, and Yoshio Kurihara, 1993: Numerical simulations of tropical cyclone-ocean interaction with a high-resolution coupled model. Journal of Geophysical Research, 98(D12), 23,245-23,263.
    [ Abstract PDF ]

    The tropical cyclone-ocean interaction was investigated using a high-resolution tropical cyclone ocean coupled model. The model design consisted of the NOAA Geophysical Fluid Dynamics Laboratory tropical cyclone prediction model which was coupled with a multilayer primitive equation ocean model. Coupling between the hurricane and the ocean models was carried out by passing into the ocean model the wind stress, heat, and moisture fluxes computed in the hurricane model. The new sea surface temperature (SST) calculated by the ocean model was then used in the tropical cyclone model. A set of idealized numerical experiments were performed in which a tropical cyclone vortex was embedded in both easterly and westerly basic flows of 2.5, 5, and 7.5 m s^-1 with a fourth experiment run with no basic flow specified initially. The profile of the trangential wind for Hurricane Gloria at 1200 UTC 22, September 1985 was used as the initial condition of the tropical cyclone for each of the experiments. The model ocean was initially horizontally homogenous and quiescent. To clarify the impact of the ocean response to the hurricane's behavior, analogous experiments were also carried out with the SST kept constant (control cases). The experiments indicated that the cooling of the sea surface induced by the tropical cyclone resulted in a significant impact on the ultimate storm intensity due to the reduction of total heat flux directed into the tropical cyclone above the regions of decreased SST. The sea surface cooling produced by the tropical cyclones was found to be larger when the storms moved slower. In the experiments run without an initial basic flow,the maximum SST anomaly was about -5.6°C with a resulting difference in the minimum sea level pressure and maximum surface winds of 16.4 hPa and -7 m s^-1, respectively. In contrast, in the experiments run with the 7.5 ms^-1 basic flow, the maximum SST anomalies ranged from about 2.6° to 3.0°C with a difference in the minimum sea level pressure and maximum surface winds of about 7.3 hPa and -2.7 m s^-1. The tropical cyclone-ocean coupling significantly influenced the storm track only for the case with no basic flow and the 2.5 m s^-1 easterly flow. In these cases the storm with the ocean interaction turned more to the north and east (no basic flow) or the north (2.5 m s^-1 easterly flow) of the experiments with constant SST. In the first case, the storm by 72 hours was located over 70 km to the east-southeast of the control case. A possible explanation for this track deviation is related to a systematic weakening of the mean tangential flow at all radii of the storm due to the interaction with the ocean and resulting alteration of the beta drift.

19. Bender, Morris A., R Ross, Robert E Tuleya, and Yoshio Kurihara, 1993: Improvements in tropical cyclone track and intensity forecasts using the GFDL initialization system. Monthly Weather Review, 121(7), 2046-2061.
    [ Abstract PDF ]

    The initialization scheme designed at GFDL to specify a more realistic initial storm structure of tropical cyclones was tested on four real data cases using the GFDL high-resolution multiply nested movable mesh hurricane model. Three of the test cases involved Hurricane Gloria (1985) in the Atlantic basin; the fourth involved Hurricane Gilbert (1988) in the Gulf of Mexico. The initialization scheme produced an initial vortex that was well adapted to the forecast model and was much more realiztic in size and intensity than the storm structure obtained from the NMC T80 global analysis. As a result, the erratic storm motion seen in previous intergrations of the GFDL model has been nearly eliminated with dramatic improvements in track forecasts during the first 48 h of the prediction. Using the new scheme, the average 24-h and 48-h forecast error for the four test cases was 58 and 94 km, respectively, compared with 143 and 191 km for the noninitialized forecasts starting from the global analysis. The average National Hurricane Center operational forecast error at 24 and 48 h was 118 and 212 km for the same four cases. After 48 h the difference in the average track error became small between the integrations starting from the global analysis and the forecasts starting from the fields obtained by the initialization scheme With accurate specification of the initial vortex structure, changes in the storm intensity were also well predicted in these cases. The model correctly forecasted the rapid intensification of Hurricane Gloria just after the system was first upgraded to a hurricane. The model storm intensification also ceased at approximately the same time as observed, with gradual weakening as the storm moved north and approached the east coast of the United States. In the forecast of Hurricane Gilbert, the model storm initially weakened as it moved over the Yucatan Peninsula and underwent only moderate reintensification after moving over the Gulf of Mexico, in good agreement with observations Finally, in the case where the track of Hurricane Gloria was well forecast, the distribution of the maximum low-level winds was accurately predicted as the storm moved up the east coast of the United States. During this period the model successfully reproduced many observed features such as large asymmetries in the wind field, with strongest winds occurring well east of the storm center, and a sharp decrease of the wind speed at the coast. Although asymmetry in the wind distribution was reproduced to a first order in the forecast starting with the global analysis, the agreement with observations was much better with the specified vortex, primarily due to a more realistic radius of maximum wind and storm intensity.

20. Kurihara, Yoshio, Morris A Bender, and R Ross, 1993: An initialization scheme of hurricane models by vortex specification. Monthly Weather Review, 121(7), 2030-2045.
    [ Abstract PDF ]

    A scheme is presented to improve the representation of a tropical cyclone in the initial condition of a high-resolution hurricane model. In the proposed method, a crudely resolved tropical cyclone in the large-scale analysis is replaced by a vortex that is properly specified for use in the prediction model. Appropriate filters are used to remove the vortex from the large-scale analysis so that a smooth environmental field remains. The new specified bogus vortex takes the form of a deviation from the environmental field so that it can be easily merged with the latter field at the correct position. The specified vortex consists of both axisymmetric and asymmetric components. The symmetric component is generated by the time integration of an axisymmetric version of the hurricane prediction model. This ensures dynamical and thermodynamical consistency in the vortex structure, including the moisture field, and also compatibility of the vortex with the resolution and physics of the hurricane model. In the course of the integration of the axisymmetric model, the tangential wind component is gradually forced to a target wind profile determined from observational information and empirical knowledge. This makes the symmetric vortex a good approximation to the corresponding real tropical cyclone. The symmetric flow thus produced is used to generate an asymmetric wind field by the time integration of a simplified barotropic vorticity equation, including the beta effect. The asymmetric wind field, which can make a significant contribution to the vortex motion, is then added to the symmetric flow. After merging the specified vortex with the environmental flow, the mass field is diagnosed from the divergence equation with an appropriately controlled time tendency. The wind field remains unchanged at this step of initialization. Since the vortex specified by the proposed method is well adapted to the hurricane prediction model, problems of initial adjustment and false spinup of the model vortex, a long-standing difficulty in the dynamical prediction of tropical cyclones, are alleviated. It is anticipated that the improvement of the initial conditions can reduce the error in hurricane track forecasting and extend the feasibility of tropical cyclone forecasting to intensify change.

21. Kurihara, Yoshio, Robert E Tuleya, Morris A Bender, and R Ross, 1993: Advanced modeling of tropical cyclones In Tropical Cyclone Disasters, Proceedings of ICSU/WMO International Symposium, October 12-16, 1992, Beijing, China, Peking University Press, 190-201.
    [ Abstract ]

    Advanced tropical cyclone models of sufficiently fine resolution are capable of representing important internal structure of the vortex. In the model, a vortex should interact with ocean and land in a realistic manner. Interaction with the ocean can significantly moderate the storm intensity. Inclusion of the heat budget of the soil layer retards the storm intensity over land. How to improve the treatment of deep convection is an issue which is wide open for future study. Specification of a realistic, yet model-adapted vortex in the initial condition of the model is essential for improvement of tropical cyclone track and intensity prediction.

22. Kurihara, Yoshio, Morris A Bender, Robert E Tuleya, and R Ross, 1993: Hurricane forecasting with the GFDL automated prediction system In 20th Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 323-326.
23. Bender, Morris A., R Ross, Yoshio Kurihara, and Robert E Tuleya, 1991: Improvements in tropical cyclone track and intensity forecasts using a bogus vortex In 19th Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 324-325.
    [ PDF ]
24. Kurihara, Yoshio, R Ross, and Morris A Bender, 1991: Toward improvement of the dynamical prediction of tropical cyclones: A hurricane model initialization scheme In 19th Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Geophysical Union, 326-327.
25. Kurihara, Yoshio, Morris A Bender, Robert E Tuleya, and R Ross, 1990: Prediction experiments of Hurricane Gloria (1985) using a multiply nested movable mesh model. Monthly Weather Review, 118(10), 2185-2198.
    [ Abstract PDF ]

    The prediction capability of the GFDL triply nested, movable mesh model, with finest grid resolution of 1/6 degree, was investigated using several case studies of Hurricane Gloria (1985) during the period that the storm approached and moved up the east coast of the United States. The initial conditions for these experiments were interpolated from an NMC T80 global analysis at 0000 UTC 25 September and 1200 UTC 22 September. The integrations starting from 0000 UTC 25 September were run 72 h, while those starting on 1200 UTC 22 September were run 132 h. The lateral boundary conditions were obtained from either an integration of the NMC T80 forecast model or the T80 global analysis, or were fixed to the initial value. The model's predicted track of Gloria for each integration was compared against the best track determined by the National Hurricane Center (NHC). For the case starting from 0000 UTC 25 September using a forecasted boundary condition, the model successfully forecasted significant acceleration of the storm's movement after 48 h. The 72 h forecast error was about 191 km, compared to 480 km for the official track forecast made by the NHC. To examine the model's skill in simulating the storm structure, distributions of the low level maximum wind and total storm rainfall during passage of the model storm are shown and compared with observed values. The model successfully reproduced many observed features such as the occurrence of strong winds well east of the storm center, with an abrupt decrease of the wind field along the coastline. When the storm track was accurately forecasted, the total storm rainfall amounts agreed well with the observed values. In both the model integration and observations, a significant structural change took place as the storm accelerated toward the north with little significant precipitation occurring south of the storm center and heavy precipitation spreading well north of the storm. It appears that the gross features of the structure of the storm's outer region resulted from the interaction of the vortex with its environment. Sensitivity of the model forecast to the lateral boundary condition and the horizontal resolution was also investigated. The storm's track error was greatly affected after the boundary error propagated by advection to the storm region. The impact of the horizontal resolution on the forecast was such that the model with one degree resolution produced a fairly good track forecast up to 48 h, but failed to simulate some of the main structural features. In the experiments starting from the 0000 UTC September 25 initial field, the interior storm structure did not develop, and the storm exhibited too large a radius of maximum wind throughout the integration. However, the integrations starting from 1200 UTC September 22 developed a more intense storm, with a more realistic radius of maximum wind. These differences were due to the spinup time necessary for the storm to develop in the model when starting from a coarse resolution global analysis which did not adequately resolve the fine structure of the storm interior. This indicates the importance of proper specification of the storm in the initial field.

26. Kurihara, Yoshio, and Morris A Bender, 1989: On the structure of moving tropical cyclones In 18th Conference on Hurricanes and Tropical Meteorology, Boston, MA, American Meteorological Society, 188-189.
27. Kurihara, Yoshio, C Kerr, and Morris A Bender, 1989: An improved numerical scheme to treat the open lateral boundary of a regional model. Monthly Weather Review, 117(12), 2714-2722.
    [ Abstract PDF ]

    A numerical scheme proposed by Kurihara and Bender is modified so as to improve the behavior of open lateral boundaries of a regional model. In the new scheme, both the local values and the gradients of fields from a larger model are used to define the time-dependent reference values toward which the boundary gridpoint values of the regional model prediction are relaxed at each step of the model integration. Use of the gradients in the boundary forcing imposes constraints on the vorticity, divergence and baroclinicity fields for the regional model. The relaxation time of forcing is set to be short for the normal component of wind. For other variables, the relaxation time at a given boundary gridpoint depends on the wind direction at that gridpoint, with a minimum at a point of normal inflow and a maximum at a point of normal outflow. The forcing strength is reduced in the planetary boundary layer so that the boundary layer structure is determined mainly by the surface condition of the regional model. Also, a simple method to control the total mass in the regional model is described. Numerical results from 96-hour integrations with the improved scheme are compared with those from the previous scheme for the cases of the propagations of a wave and a vortex. The behavior of the model at the lateral boundary was noticeably improved with the use of the new scheme, while the solution in the interior domain was little affected by the scheme modification.

28. Bender, Morris A., and Yoshio Kurihara, 1987: A numerical study of the effect of the mountainous terrain of Japan on tropical cyclones In Short- and Medium-Range Numerical Weather Prediction, Collection of Papers Presented At WMO/IUGG NWP Symposium Tokyo, Japan, Aug. 4-8, 1986, Geneva, Switz, World Meteorological Organization, 651-663.
    [ Abstract PDF ]

    A triply-nested, movable mesh model was used to study the effects of mountainous terrain on the landfall of tropical cyclones onto the islands of Japan. The integration domain spanned 43° latitude and 47° longitude with finest resolution of 1/6°. Numerical experiments were separately performed for three cases. In each experiment a storm was embedded onto a stationary Haurwitz type wave at the initial time, and moved in a north-northeast direction at about 10 to 12 m s-1. In the first experiment, the tropical cyclone struck the southwest Izu Peninsula in eastern Japan. The second and third respectively made landfall on the Kii Peninsula in the central part of Japan and on the island of Kyushu in western Japan. In order to isolate some of the effects on the storm system resulting from its interaction with the mountainous terrain, these simulations were compared with supplemental experiments performed with a flat land condition. In all three cases it was found that the presence of the mountainous terrain greatly enhanced the storm decay after landfall. As the storm approached Tokyo Bay rapid weakening occurred as dry air from the mountain region to the west of Tokyo was advected into the eye and eyewall region of the storm. Upon leaving eastern Japan and again moving over open water, the storm never underwent reintensification. In the case of the storm leaving western Japan, reintensification over the Sea of Japan occurred very slowly as compared with the experiment run with a flat land distribution. Apparently, the above behavior was related to the structural change which occurred to the storm system during the passage over the mountainous islands. The precipitation pattern was also greatly affected by the presence of the mountainous terrain. As the storm made landfall over central Japan, the area of heaviest rainfall shifted to the right of the storm track, where strong upslope winds developed. This storm eventually travelled over the high mountains of east-central Japan and rapidly decayed by the end of the experiment. Although performed for an idealized experimental design, these experiments reveal some of the important effects the mountainous terrain may have on the behavior of tropical cyclones making landfall on Japan. Understanding these effects should prove useful in forecasting more accurately the behavior of the storms.

29. Bender, Morris A., Robert E Tuleya, and Yoshio Kurihara, 1987: A numerical study of the effect of island terrain on tropical cyclones. Monthly Weather Review, 115(1), 130-155.
    [ Abstract PDF ]

    A triply-nested. movable mesh model was used to study the behavior of tropical cyclones encountering island mountain ranges. The integration domain consisted of a 37° wide and 45° long channel, with an innermost mesh resolution of 1/6°. The storms used for this study were embedded in easterly flows of ~ 5 and ~ 10 m s^-1 initially. Realistic distributions of island topography at 1/6° resolution were inserted into the model domain for the region of the Caribbean, including the islands of Cuba, Hispaniola, and Puerto Rico; the island of Taiwan; and the region of Luzon in the northern Philippines. It was found that the islands affected the basic flow as well as the wind field directly associated with the storm system. The combination of these effects caused changes in the track and translational speed of the storm. In particular, in the case of the 5 m s^-1 easterly flow, the storm accelerated and veered to the north well before reaching Taiwan. For the other island distributions, the northward deflection of the track and the increase of translational speed occurred near and over the islands. After landfall, the surface pressure underwent rapid filling. As the tropical cyclone passed over Hispaniola, the surface low continued to move along with the upper level vortex as it transversed the mountain range, while over Luzon it became obscure before reforming on the lee side slope of the mountain. In case of Taiwan and the 10 m s^-1 easterly zonal flow, secondary surface lows developed behind the mountain range. The upper level vortex in this case became detached from the original surface low and eventually coupled with a secondary one. The intensity changes of the storm near and over the islands were strongly related to the latent energy supply and the vertical coherence of the storm system. Advection of dry air from near or above the mountain tops into the storm area caused significant weakening of all the storm moving with the weaker easterly flow. Storms leaving Hispaniola and moving over open sea quickly reintensified as their vertical structure remained coherent. On the other hand, storms leaving Luzon were disorganized and did not reintensify until several hours later when the vertical coherence of the systems was reestablished Although these experiments were performed for an idealized experimental design and basic flow, many observed storms have exhibited similar behavior in track deviation and decay. This implies that the effect of detailed topography should be considered if an accurate forecast of the storm direction and behavior is to be made.

30. Bender, Morris A., Robert E Tuleya, and Yoshio Kurihara, 1985: A numerical study of the effect of a mountain range on a landfalling tropical cyclone In 16th Conference on Hurricanes & Tropical Meteorology, Boston. MA, American Meteorological Society, 146-147.
31. Bender, Morris A., Robert E Tuleya, and Yoshio Kurihara, 1985: A numerical study of the effect of a mountain range on a landfalling tropical cyclone. Monthly Weather Review, 113(4), 567-582.
    [ Abstract PDF ]

    A triply-nested, movable mesh model was used to study the effects of a mountain range on a landfalling tropical cyclone embedded in an easterly flow of ~ 10 m s^-1. The integration domain consisted of a 37 degree wide and 45 degree long channel, with an innermost mesh resolution of 1/6 degrees. An idealized mountain range with maximum height of ~ 958 meters was placed parallel to the shoreline. The mountain range, which spanned 19 degrees in the north-south direction and 5 degrees in the east-west direction, was centered in the middle of the channel. Results obtained were compared with a previous landfall simulation, performed without the effect of the mountain range included. In particular, comparison was made of the total storm rainfall, maximum wind distribution and storm decay rate. It was found that the storm filled much more rapidly in the simulation run with the mountain included. The mountain range affected the decay rate through reduction in the supply of latent and kinetic energy into the storm circulation during, as well as after, passage of the storm over the mountain. It was found that a low-level, warm and dry region was produced where the storm winds descended the mountain slope. In order to better isolate the effect of the mountain on the basic easterly flow, a supplemental integration was performed for the flow without the storm. It revealed that the mountain range caused a significant change in the basic flow over the mountain as well as up to several hundred kilometers downstream and extending considerably above the mountain top. A low-level southerly jet was observed to the west of the mountain base.

32. Tuleya, Robert E., Morris A Bender, and Yoshio Kurihara, 1984: A simulation study of the landfall of tropical cyclones using a movable nested-mesh model. Monthly Weather Review, 112(1), 124-136.
    [ Abstract PDF ]

    By use of a triply nested, movable mesh model, several ideal simulations of tropical cyclone landfall were performed for a strong zonal flow of ~10 m s-1. The integration domain was a 37 X 45 degree channel with the innermost mesh having a 22 X 22 point resolution of 1/6 degree. General characteristics similar to observed landfalling tropical cyclones are obtained in the primary simulation experiment including an abrupt change in the low level (~68 m) winds at the coastline and a decay of the tropical cyclone as it moves inland. Additional interesting features subject to model and experimental limitations include: little noticeable track change of the model storm when compared to a control experiment with an ocean surface only; a possible temporary displacement of the center of the surface wind circulation from the surface pressure center at landfall; and a distinct decrease in kinetic energy generation and precipitation a few hours after landfall. The sensitivity to the specified land surface conditions was analyzed by performing additional experiments in which the land surface conditions including surface temperature, moisture, and distribution of surface roughness were changed. It was found that a reasonable change in some of these land conditions can make a considerable difference in behavior for a landfalling tropical cyclone. It was also shown that a small, less intense model storm fills less rapidly. This corresponds well with observations that many landfalling hurricanes decay to approximately the same asymptotic value one day after landfall.

33. Bender, Morris A., and Yoshio Kurihara, 1983: The energy budgets for the eye and eye wall of a numerically simulated tropical cyclone. Journal of the Meteorological Society of Japan, 61(2), 239-244.
    [ Abstract PDF ]

    Energy budgets are analyzed for a tropical cyclone simulated previously in a quadruply nested mesh model (Kurihara and Bender, 1982). It will be shown that the eddy kinetic energy within the eye is comparable in magnitude to that of the mean kinetic energy. It is supplied by import from the eye wall regions as well as by the conversion from total potential energy. At the same time it is converted to the kinetic energy of the mean flow and also lost by the dissipation. The influx of mean kinetic energy from the outer radii to the eye wall region and the export of potential energy both to the outer radii and to the eye region play important roles in the energetics of the eye wall region. Many obtained features agree well with those of a coarser resolution model (Tuleya and Kurihara, 1975) in which the eye of the vortex could not be resolved. This suggests that the eye structure has little impact on the energetics to the eye wall and outer regions of a tropical cyclone.

34. Kurihara, Yoshio, and Morris A Bender, 1983: A numerical scheme to treat the open lateral boundary of a limited area model. Monthly Weather Review, 111(3), 445-454.
    [ Abstract PDF ]

    A numerical scheme to treat the open lateral boundary of a limited-area primitive equation model was formulated. Although overspecification of the boundary condition is inevitable in the pointwise boundary setting, the scheme was designed to keep the overspecification to a minimum degree. To impose the boundary conditions, a damping technique was used. Special care was taken to deal with the boundary layer winds at the lateral boundary. The above scheme is most suitable when gravity waves do not prevail in the vicinity of the open boundary. The scheme was tested in the numerical integrations of prognostic equations for a Haurwitz-type wave. Experimental results are presented which indicate the utility of the proposed method.

35. Tuleya, Robert E., Morris A Bender, and Yoshio Kurihara, 1983: A numerical study of simulated hurricane landfall In Sixth Conference on Numerical Weather Prediction, Boston, MA, American Meteorological Society, 323-324.
36. Kurihara, Yoshio, and Morris A Bender, 1982: Structure and analysis of the eye of a numerically simulated tropical cyclone. Journal of the Meteorological Society of Japan, 60(1), 381-395.
    [ Abstract PDF ]

    A tropical cyclone has been simulated in a quadruply nested mesh model with finest grid resolution of about 5 km. At the center of the vortex, a compact eye was maintained. Azimuthal means as well as asymmetry of the eye and the eye wall structure are described. The asymmetric features within the eye wall moved cyclonically at a much smaller rotation rate than the cyclonic wind within the eye wall. Roles of the mean radial-vertical circulation, the eddy motions and the diffusion effect in the maintenance of the mean structure are analyzed. In the analysis, attention is given to the balance between the wind and pressure fields and also to the budgets of relative angular momentum, heat and water vapor. The eddy motions caused a cooling and moistening effect in the eye which counterbalanced a warming and drying effect due to the mean sinking motion.

37. Kurihara, Yoshio, and Morris A Bender, 1980: Use of a movable nested-mesh model for tracking a small vortex. Monthly Weather Review, 108(11), 1792-1809.
    [ Abstract PDF ]

    The mesh nesting strategy proposed by Kurihara et al. (1979) was used to construct a movable, nested-mesh, 11-level primitive equation model. The framework of the model is described in detail. With the use of a triply nested mesh system with 1 degree, 1/3 degree and 1/6 degree longitude-latitude resolution, a small intense dry vortex in a zonal flow of 10 m s^-1 was successfully advected for 48 h. The shape of the vortex was well preserved during the time integration which involved over 50 movements of the innermost mesh. The noise, which was excited when a mesh moved, was suppressed in ~4 minutes after the movement. For comparison, the results from similar experiments performed with reduced inner mesh resolutions are also presented.

38. Kurihara, Yoshio, and Morris A Bender, 1979: Supplementary note on "A scheme of dynamic initialization of the boundary layer in a primitive equation model". Monthly Weather Review, 107(9), 1219-1221.
    [ Abstract PDF ]

    A scheme is presented for improving the previously proposed method of dynamic initialization of the boundary layer in a primitive equation model (Kurihara and Tuleya, 1978). Performance of the revised scheme is shown for the case of a strong vortex superposed on a zonal flow.

39. Kurihara, Yoshio, G J Tripoli, and Morris A Bender, 1979: Design of a movable nested-mesh primitive equation model. Monthly Weather Review, 107(3), 239-249.
    [ Abstract PDF ]

    A numerical scheme to construct a two-way, movable, nested-mesh primitive equation model is proposed. Dynamical coupling in a two-way nesting system is performed at a dynamical interface which is separated from a mesh interface by two coarse-grid intervals. Dynamical interaction is achieved by a method which conserves mass, momentum and internal energy of the system. During the course of integration, the nested mesh moves so that the central position of the disturbance contained in the fine-mesh area never deviates from the center of the nest by more than one coarse-mesh interval. New grid data near the leading and trailing edges of the moving nest are obtained by an interpolation method which has a conservation property. The proposed methods of dynamical coupling and mesh movement were extensively tested by a one-dimensional shallow water equation model. Numerical results of these experiments are presented.

Direct link to page: http://www.gfdl.noaa.gov/bibliography/resultstest.php?author=1001