A Multiyear Ensemble Simulation of the U.S. Climate with a Stretched-Grid GCMFox-Rabinovitz, Michael S.; Berbery, Ernesto Hugo; Takacs, Lawrence L.; Govindaraju, Ravi C.
doi: 10.1175/MWR2956.1pmid: N/A
Multiyear (1987–97) limited ensemble integrations using a stretched-grid GCM, previously developed and experimented with by the authors, are employed for U.S. regional climate simulations. The ensemble members (six in total) are produced at two different regional resolutions: three members with 60-km and the other three members with 10-km regional resolution. The use of these two finer and coarser regional resolution ensemble members allows one to examine the impact of resolution on the overall quality of the simulated regional fields. For the multiyear ensemble simulations, an efficient regional downscaling to realistic mesoscales has been obtained. The ensemble means of the midtroposphere prognostic variables (height and meridional wind) show an overall good resemblance to the global reanalysis, especially for summer. Low-level features like the warm season Great Plains low-level jet are well represented in the simulations. During winter the 100-km simulations develop a southward wind east of the Rockies that is present neither in the reanalyses nor in the 60-km simulations. The analysis of the annual mean precipitation and its variance reveals that the ensemble simulations reproduce many of the observed features of a high-resolution rain gauge dataset analyzed on a 0.5° × 0.5° grid. Signal-to-noise ratios are larger than 1.5 s over a major part of the United States, especially over the Midwest and also over the mountainous regions like the Rockies and the Appalachians, suggesting that the orographic forcing is contributing to a larger signal. The ratios are smaller toward the eastern and western U.S. coastlines. This result could be attributed, at least in part, to limits in the representation of the land–sea contrasts. For comparison purposes, an additional simulation has been performed using a global uniform 2° × 2.5° grid with the same number of global grid points as those of the above stretched grids. The stretched-grid GCM ensemble means show, overall, a better regional depiction of features than those of the uniform-grid GCM. The results of the study show that even using limited ensemble integrations with a state-of-the-art stretched-grid GCM is beneficial for reducing the uncertainty of the multiyear regional climate simulation, especially when using finer 60-km regional resolution.
Terdiurnal Surface-Pressure Oscillations over the Continental United StatesRay, Richard D.; Poulose, Susan
doi: 10.1175/MWR2988.1pmid: N/A
The small terdiurnal pressure oscillation S 3 ( p ) is determined over the conterminous United States by analyzing long time series of hourly barometer data from 180 stations. Spectral analysis of these time series reveals that the terdiurnal band is dominated by three or four spectral peaks, separated in frequency by 1 cpy. The central peak at 3 cpd is invariably smaller than the two immediate side peaks, indicative of extraordinarily strong seasonal variations in the tide. The largest terdiurnal tide occurs over the south-central United States in winter where amplitudes exceed 300 μ bar. Summertime amplitudes are roughly one-half as large. Summer and winter tides are almost completely out of phase, with rapid 180° shifts occurring in the equinox seasons when amplitudes are very small.
Comparison between Tornadic and Nontornadic Mesocyclones Using the Vorticity (Pseudovorticity) Line TechniqueCai, Huaqing
doi: 10.1175/MWR2990.1pmid: N/A
Comparisons between tornadic and nontornadic mesocyclones using the concept of fractal geometry are presented. Both the maximum vertical vorticity ( ζ max ) and pseudovorticity ( ζ pv ) associated with a mesocyclone at low levels are found to be scaling with the horizontal grid spacing ( ε ) according to a power-law relationship. The linear least square best fitting of ln ( ζ max ) or ln ( ζ pv ) versus ln ( ε ) for different scales can be obtained for each mesocyclone at a certain time, and it is named the vorticity (pseudovorticity) line of a mesocyclone. Different mesocyclones have different vorticity (pseudovorticity) line slopes that are closely related to the fractal dimension of vorticity (pseudovorticity) of a mesocyclone as a possible fractal structure. Various factors that may affect the accurate estimate of the vorticity (pseudovorticity) line of a mesocyclone are also discussed in detail. Differences between tornadic and nontornadic mesocyclones are found in terms of the slope of vorticity (pseudovorticity) lines based on three tornadic and two nontornadic mesocyclones. A possible reason why previous studies were not able to identify the difference(s) between tornadic and nontornadic mesocyclones is discussed. Self-similarity (scale invariance), which is a basic characteristic of a fractal structure, seems to be valid between tornado and mesocyclone scales based on the analysis of the vorticity (pseudovorticity) line of the tornadic Kellerville, Texas, mesocyclone. It is hypothesized that a steeper slope of the vorticity (pseudovorticity) line may be indicative of a tornadic mesocyclone.
Flow over Heated Terrain. Part I: Linear Theory and Idealized Numerical SimulationsCrook, N. Andrew; Tucker, Donna F.
doi: 10.1175/MWR2964.1pmid: N/A
The flow past heated topography is examined with both linear and nonlinear models. It is first shown that the forcing of an obstacle with horizontally homogenous surface heating can be approximated by the forcing of an obstacle with surface heating isolated over the obstacle. The small-amplitude flow past an obstacle with isolated heating is then examined with a linear model. Under the linear approximation, the flow response to heated topography is simply the addition of the separate responses to thermal and orographic forcing. These separate responses are first considered individually and then the combined response is examined. Nondimensional parameters are developed that measure the relative importance of thermal and orographic forcing. Nonaxisymmetric forcing is then considered by examining the flow along and across a heated elliptically shaped obstacle. It is shown that the low-level lifting is maximized when the flow is along the major axis of the obstacle. The linear solutions are then tested in a nonlinear anelastic model. The response to a heat source and orography are first examined separately. Good agreement is found between nonlinear and linear models for the individual responses to thermal and orographic forcing. The case of uniformly heated flow past an obstacle is then examined. In these simulations, the thermal response is isolated by subtracting the orographic-only response from the full thermal–orographic response. The numerical simulations are able to capture the main features of the thermal response. Finally, numerical simulations of the flow along and across an elliptically shaped heated obstacle are examined, where it is verified that the lifting is maximized when the flow is along the major axis of the obstacle. These results are extended in Part II of this study to examine the moist convective response to flow over both idealized terrain and the complex terrain of the Rocky Mountains of the United States.
Flow over Heated Terrain. Part II: Generation of Convective PrecipitationTucker, Donna F.; Crook, N. Andrew
doi: 10.1175/MWR2965.1pmid: N/A
Previous studies have shown that thunderstorms in the Rocky Mountain region have preferred areas in which to form. There has been some indication that these areas depend on the midtropospheric wind direction. A nonhydrostatic model with a terrain-following horizontal grid is employed to investigate the initiation of precipitating convection over heated topography. Horizontally homogeneous meteorological conditions with no directional shear in the vertical wind profile are used. The numerical simulations indicate that precipitating convection was more likely to be generated downwind of ridges than upwind of them. Initiation of these storms was more likely downwind of ridges with their long axis parallel to the wind direction than downwind of ridges with their long axis perpendicular to the wind direction. In Part I of this study it was shown that heating-induced convergence is larger downwind of a ridge with its longer axis parallel to the wind direction. For the orographic configuration of the Rocky Mountains, total precipitation is maximized for southerly and northwesterly winds. Slower wind speeds are more likely and faster wind speeds are less likely to produce convective storms. Soundings with larger instability are more likely to produce convection. The soundings with a greater temperature lapse rate produce more initiation locations, and soundings with greater moisture produce greater amounts of precipitation. Even though a number of assumptions were made for this study, the authors believe the results explain a significant amount of the observed variability in the initiation locations of precipitating convection in the Rocky Mountains during the summer. Because of the theoretical basis for this work, detailed in Part I of this study, the authors believe it should explain convective initiation in other mountainous areas that are subject to strong solar heating.
The Electrical Structure of Two Supercell Storms during STEPSMacGorman, Donald R.; Rust, W. David; Krehbiel, Paul; Rison, William; Bruning, Eric; Wiens, Kyle
doi: 10.1175/MWR2994.1pmid: N/A
Balloon soundings were made through two supercell storms during the Severe Thunderstorm Electrification and Precipitation Study (STEPS) in summer 2000. Instruments measured the vector electric field, temperature, pressure, relative humidity, and balloon location. For the first time, soundings penetrated both the strong updraft and the rainy downdraft region of the same supercell storm. In both storms, the strong updraft had fewer vertically separated charge regions than found near the rainy downdraft, and the updraft’s lowest charge was elevated higher, its bottom being near the 40-dB Z boundary of the weak-echo vault. The simpler, elevated charge structure is consistent with the noninductive graupel–ice mechanism dominating charge generation in updrafts. In the weak-echo vault, the amount of frozen precipitation and the time for particle interactions are too small for significant charging. Inductive charging mechanisms and lightning may contribute to the additional charge regions found at lower altitudes outside the updraft. Lightning mapping showed that the in-cloud channels of a positive ground flash could be in any one of the three vertically separated positive charge regions found outside the updraft, but were in the middle region, at 6–8 km MSL, for most positive ground flashes. The observations are consistent with the electrical structure of these storms having been inverted in polarity from that of most storms elsewhere. It is hypothesized that the observed inverted-polarity cloud flashes and positive ground flashes were caused by inverted-polarity storm structure, possibly due to a larger than usual rime accretion rate for graupel in a strong updraft.
The “Owl Horn” Radar Signature in Developing Southern Plains SupercellsKramar, Matthew R.; Bluestein, Howard B.; Pazmany, Andrew L.; Tuttle, John D.
doi: 10.1175/MWR2992.1pmid: N/A
During spring 2001 in the Southern Plains, a recurring, hitherto undocumented reflectivity signature that the authors have called the “Owl Horn” signature (because the radar reflectivity pattern resembles the profile of the Great Horned Owl) was observed on a mobile, X-band radar display. The reflectivity signature was always located at the rear side of a developing supercell, spanned the entire rear side of the storm, and was always seen on low-level plan position indicator (PPI) scans. It lasted on the order of only 5–10 min and was not an artifact of the radar. A study of the Owl Horn signature was undertaken using the Tracking Radar Echoes by Correlation technique (TREC) to estimate the wind field. TREC has previously been applied to clear-air and hurricane environments, and to the internal motions of severe storms, but not to their evolution. The characteristics of the signature are presented, and then, through the application of TREC to the radar reflectivity data (Doppler wind data were not available in 2001) collected during May and June 2001, the horizontal wind field was estimated around and in the Owl Horn signature. Instances of the Owl Horn in numerical model storm simulations were investigated. The numerical simulations were used to identify conditions under which the signature occurs, the process by which it is created is discussed, and its dependence upon the environmental wind shear is examined. Results indicate that the hodograph shape and magnitude influence the production of the Owl Horn signature. Supercell-magnitude shear is required, and some curvature—particularly low-level curvature—is essential to the production of the feature. The Owl Horn signature is formed when horizontal vorticity is tilted into the vertical by expanding outflow through a positive feedback mechanism with the outflow.
The Interaction of Supertyphoon Maemi (2003) with a Warm Ocean EddyLin, I-I.; Wu, Chun-Chieh; Emanuel, Kerry A.; Lee, I-Huan; Wu, Chau-Ron; Pun, Iam-Fei
doi: 10.1175/MWR3005.1pmid: N/A
Understanding the interaction of ocean eddies with tropical cyclones is critical for improving the understanding and prediction of the tropical cyclone intensity change. Here an investigation is presented of the interaction between Supertyphoon Maemi, the most intense tropical cyclone in 2003, and a warm ocean eddy in the western North Pacific. In September 2003, Maemi passed directly over a prominent (700 km × 500 km) warm ocean eddy when passing over the 22°N eddy-rich zone in the northwest Pacific Ocean. Analyses of satellite altimetry and the best-track data from the Joint Typhoon Warning Center show that during the 36 h of the Maemi–eddy encounter, Maemi’s intensity (in 1-min sustained wind) shot up from 41 m s −1 to its peak of 77 m s −1 . Maemi subsequently devastated the southern Korean peninsula. Based on results from the Coupled Hurricane Intensity Prediction System and satellite microwave sea surface temperature observations, it is suggested that the warm eddies act as an effective insulator between typhoons and the deeper ocean cold water. The typhoon’s self-induced sea surface temperature cooling is suppressed owing to the presence of the thicker upper-ocean mixed layer in the warm eddy, which prevents the deeper cold water from being entrained into the upper-ocean mixed layer. As simulated using the Coupled Hurricane Intensity Prediction System, the incorporation of the eddy information yields an evident improvement on Maemi’s intensity evolution, with its peak intensity increased by one category and maintained at category-5 strength for a longer period (36 h) of time. Without the presence of the warm ocean eddy, the intensification is less rapid. This study can serve as a starting point in the largely speculative and unexplored field of typhoon–warm ocean eddy interaction in the western North Pacific. Given the abundance of ocean eddies and intense typhoons in the western North Pacific, these results highlight the importance of a systematic and in-depth investigation of the interaction between typhoons and western North Pacific eddies.
A Refractive Index Mapping Operator for Assimilation of Occultation DataSyndergaard, Stig; Kursinski, E. Robert; Herman, Benjamin M.; Lane, Emily M.; Flittner, David E.
doi: 10.1175/MWR3001.1pmid: N/A
This paper describes the details of a fast, linear, forward-inverse refractive index mapping operator that can be used for assimilation of occultation data of various kinds into NWP models. Basically, the mapping consists of the integration of the refractive index along finite straight lines, mimicking the observational geometry as well as the subsequent retrieval of a refractive index profile, assuming spherical symmetry. Line integrals are discretized such that the refractivity is evaluated along the horizontal at fixed levels that can be chosen to coincide with the pressure levels of an NWP model. Integration of the hydrostatic equation at a large number of locations is thereby avoided. The mapping operator is tested using an idealized model of a weather front with large horizontal gradients. Mapped refractivity profiles are compared with retrieved refractivity profiles obtained via accurate 3D ray tracing simulations of GPS radio occultation events with ray path tangent points near the weather front. The simulations indicate that the mapping is a good representation of occultation measurements, including the influence large horizontal gradients have on retrieved refractivity profiles. To further the results, a simple ad hoc modification is introduced to approximately account for the ray path bending near the tangent points. The forward-inverse mapping allows for the near cancellation of otherwise crude approximations—for example, straight-line propagation—and the general concept could perhaps be adapted for the development of fast and accurate observation operators for the assimilation of other types of remote sensing data.
Environmental Distinctions between Cellular and Slabular Convective LinesJames, Richard P.; Fritsch, J. Michael; Markowski, Paul M.
doi: 10.1175/MWR3002.1pmid: N/A
The organizational mode of quasi-linear convective systems often falls within a spectrum of modes described by a line of discrete cells on one end (“cellular”) and an unbroken two-dimensional swath of ascent on the other (“slabular”). Convective events exhibiting distinctly cellular or slabular characteristics over the continental United States were compiled, and composite soundings of the respective inflow environments were constructed. The most notable difference between the environments of slabs and cells occurred in the wind profiles; lines organized as slabs existed in much stronger low-level line-relative inflow and stronger low-level shear. A compressible model with high resolution (Δ x = 500 m) was used to investigate the effects of varying environmental conditions on the nature of the convective overturning. The numerical results show that highly cellular convective lines are favored when the environmental conditions and initiation procedure allow the convectively generated cold pools to remain separate from one another. The transition to a continuous along-line cold pool and gust front leads to the generation of a more “solid” line of convection, as dynamic pressure forcing above the downshear edge of the cold outflow creates a swath of quasi-two-dimensional ascent. Using both full-physics simulations and a simplified cold-pool model, it is demonstrated that the magnitude of the two-dimensional ascent in slabular convective systems is closely related to the integrated cold-pool strength. It is concluded that slabular organization tends to occur under conditions that favor the development of a strong, contiguous cold pool. The tendency to produce slabular convection is therefore enhanced by environmental conditions such as large CAPE, weak convective inhibition, strong along-line winds, and moderately strong cross-line wind shear.