Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonhydrostatic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Examples of data assimilation in mesoscale models

The keynote address was the problem of physical initialization of mesoscale models. The classic purpose of physical or diabatic initialization is to reduce or eliminate the spin-up error caused by the lack, at the initial time, of the fully developed vertical circulations required to support regions of large rainfall rates. However, even if a model has no spin-up problem, imposition of observed moisture and heating rate information during assimilation can improve quantitative precipitation forecasts, especially early in the forecast. The two key issues in physical initialization are the choice of assimilating technique and sources of hydrologic/hydrometeor data. Another example of data assimilation in mesoscale models was presented in a series of meso-beta scale model experiments with and 11 km version of the MASS model designed to investigate the sensitivity of convective initiation forced by thermally direct circulations resulting from differential surface heating to four dimensional assimilation of surface and radar data. The results of these simulations underscore the need to accurately initialize and simulate grid and sub-grid scale clouds in meso- beta scale models. The status of the application of the CSU-RAMS mesoscale model by the NOAA Forecast Systems Lab for producing real-time forecasts with 10-60 km mesh resolutions over (4000 km)(exp 2) domains for use by the aviation community was reported. Either MAPS or LAPS model data are used to initialize the RAMS model on a 12-h cycle. The use of MAPS (Mesoscale Analysis and Prediction System) model was discussed. Also discussed was the mesobeta-scale data assimilation using a triply-nested nonhydrostatic version of the MM5 model.

Carr, Fred↗

Session on validation of coupled models

The session on validation of coupled models is reviewed. The current use of a mesoscale model with a grid size of 20-km during STORM-FEST in 1992 has proven to be extremely valuable. The availability of forecast products at a much higher temporal and spatial resolution was very helpful for mesoscale forecasting, mission planning, and the guidance of research aircraft. Recent numerical simulation of ocean cyclones and mesoscsle convective systems using nonhydrostatic cloud/mesoscale models with a grid size as small as 2-km have demonstrated the potential of these models for predicting mesoscale convective systems, squall lines, hurricane rainbands, mesoscale gravity waves, and mesoscale frontal structures embedded within an extratropical cyclone. Although mesoscale/cloud scale models have demonstrated strong potential for use in operational forecasting, very limited quantitative evaluation (and verification) of these models were performed. As a result, the accuracy, the systematic biases, and the useful forecasts limits were not properly defined for these models. Also, no serious attempts were made to use these models for operational prediction of mesoscale convective systems.

Kuo, Bill↗

Case study of a low-reflectivity pulsating microburst - Numerical simulation of the Denver, 8 July 1989, storm

Results from a 3D simulation of a very strong microburst, detected by the Low-Level Windshear Alert System, are presented. The simulation was performed using the Terminal Area Simulation System (TASS) model (Proctor, 1987) which is a 3D nonhydrostatic cloud model that includes parameterizations for both liquid and ice-phase microphysics. Emphasis is placed om the pulsating characteristic and the very-low radar reflectivity of the event. The simulation revealed that the dissipating cloud induces a type-1 pulsating microburst event, with at least three distinct microburst pulses over a period of 20 min.

Proctor, Fred H.↗

Case study of a low-reflectivity pulsating microburst: Numerical simulation of the Denver, 8 July 1989, storm

On 8 July 1989, a very strong microburst was detected by the Low-Level Windshear Alert system (LLWAS), within the approach corridor just north of Denver Stapleton Airport. The microburst was encountered by a Boeing 737-200 in a 'go-around' configuration which was reported to have lost considerable air speed and altitude during penetration. Data from LLWAS revealed a pulsating microburst with an estimated peak velocity change of 48 m/s. Wilson et al. reported that the microburst was accompanied by no apparent visible clues such as rain or virga, although blowing dust was present. Weather service hourly reports indicated virga in all quadrants near the time of the event. A National Center for Atmospheric Research (NCAR) research Doppler radar was operating; but according to Wilson et al., meaningful velocity could not be measured within the microburst due to low radar-reflectivity factor and poor siting for windshear detection at Stapleton. This paper presents results from the three-dimensional numerical simulation of this event, using the Terminal Area Simulation System (TASS) model. The TASS model is a three-dimensional nonhydrostatic cloud model that includes parameterizations for both liquid and ice phase microphysics, and has been used in investigations of both wet and dry microburst case studies. The focus of this paper is the pulsating characteristic and the very-low radar reflectivity of this event. Most of the surface outflow contained no precipitation. Such an event may be difficult to detect by radar.

Proctor, Fred H.↗

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗

Numerical simulations of upstream blocking, columnar disturbances, and bores in stably stratified shear flows over an obstacle

A two-dimensional, nonhydrostatic, elastic numerical model has been used to study the generation of gravity waves for a stably stratified shear flow over an obstacle. When a low-level wind shear is included in the simulation, we find that the predictions for noticeable upstream effects based on Froude number for a uniform flow are no longer accurate. Upstream effects are encountered in the form of upstream propagating columnar disturbances and internal bores away from the obstacle. The limited parameter space studies conducted in this study suggest that the ratio of the shear depth to the obstacle heigh (d/H), the obstacle aspect ratio (H/L), and the Froude number (U/NH) are instrumental in determining the strength and the existence of these upstream disturbances. Thus, the present theoretical and empirical understanding of the importance of the Froude number for determining the nature of upstream effects should be modified substantially to include additional nondimensional parameters when shear is present.

Chen, Chaing↗

Gravity waves in Titan's atmosphere

Scintillations (high frequency variations) observed in the radio signal during the occultation of Voyager 1 by Titan (Hinson and Tyler, 1983) provide information concerning neutral atmospheric density fluctuations on scales on hundreds of meters to a few kilometers. Those seen at altitudes higher than 25 km above the surface were interpreted by Hinson and Tyler as being caused by linear, freely propagating (energy-conserving) gravity waves, but this interpretation was found to be inconsistent with the scintillation data below the 25-km altitude level. Here an attempt is made to interpret the entire scintillation profile between the surface and the 90-km altitude level in terms of gravity waves generated at the surface. Numerical calculations of the density fluctuations caused by two-dimensional, nonhydrostatic, finite-amplitude gravity waves propagating vertically through Titan's atmosphere are performed to produce synthetic scintillation profiles for comparison with the observations. The numerical model accurately treats the effects of wave transience, nonlinearity, and breakdown due to convective instability in the overturned part of the wave. The high-altitude scintillation data were accurately recovered with a freely propagating wave solution, confirming the analytic model of Hinson and Tyler. It is found that the low-altitude scintillation data can be fit by a model where a component of the gravity waves becomes convectively unstable and breaks near the 15 km level. The large-scale structure of the observed scintillation profile in the entire altitude range between 5 and 85 km can be simulated by a model where the freely propagating and breaking waves are forced at the surface simultaneously. Further analysis of the Voyager 1 Titan low-altitude scintillation data, using inversion theory appropriate for strong scattering, could potentially remove some of the ambiguities remaining in this analysis and allow a better determination of the strength and source of the waves.

Friedson, A. James↗

Evaluation of Bogus Vortex Techniques with Four-Dimensional Variational Data Assimilation

The effectiveness of techniques for creating "bogus" vortices in numerical simulations of hurricanes is examined by using the Penn State/NCAR nonhydrostatic mesoscale model (MM5) and its adjoint system. A series of four-dimensional variational data assimilation (4-D VAR) experiments is conducted to generate an initial vortex for Hurricane Georges (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and surface wind information into the mesoscale numerical model. Several different strategies are tested for improving the vortex representation. The initial vortices produced by the 4-D VAR technique are able to reproduce many of the structural features of mature hurricanes. The vortices also result in significant improvements to the hurricane forecasts in terms of both intensity and track. In particular, with assimilation of only bogus sea-level pressure information, the response in the wind field is contained largely within the divergent component, with strong convergence leading to strong upward motion near the center. Although the intensity of the initial vortex seems to be well represented, a dramatic spin down of the storm occurs within the first 6 h of the forecast. With assimilation of bogus surface wind data only, an expected dominance of the rotational component of the wind field is generated, but the minimum pressure is adjusted inadequately compared to the actual hurricane minimum pressure. Only when both the bogus surface pressure and wind information are assimilated together does the model produce a vortex that represents the actual intensity of the hurricane and results in significant improvements to forecasts of both hurricane intensity and track.

Pu, Zhao-Xia↗

Microphysics, Radiation and Surface Processes in the Goddard Cumulus Ensemble (GCE) Model

The response of cloud systems to their environment is an important link in a chain of processes responsible for monsoons, frontal depression, El Nino Southern Oscillation (ENSO) episodes and other climate variations (e.g., 30-60 day intra-seasonal oscillations). Numerical models of cloud properties provide essential insights into the interactions of clouds with each other, with their surroundings, and with land and ocean surfaces. Significant advances are currently being made in the modeling of rainfall and rain-related cloud processes, ranging in scales from the very small up to the simulation of an extensive population of raining cumulus clouds in a tropical- or midlatitude-storm environment. The Goddard Cumulus Ensemble (GCE) model is a multi-dimensional nonhydrostatic dynamic/microphysical cloud resolving model. It has been used to simulate many different mesoscale convective systems that occurred in various geographic locations. In this paper, recent GCE model improvements (microphysics, radiation and surface processes) will be described as well as their impact on the development of precipitation events from various geographic locations. The performance of these new physical processes will be examined by comparing the model results with observations. In addition, the explicit interactive processes between cloud, radiation and surface processes will be discussed.

Tao, Wei-Kuo↗

An Interconnected Network of Core-Forming Melts Produced by Shear Deformation

The formation mechanism of terrestrial planetary is still poorly understood, and has been the subject of numerous experimental studies. Several mechanisms have been proposed by which metal-mainly iron with some nickel-could have been extracted from a silicate mantle to form the core. Most recent models involve gravitational sinking of molten metal or metal sulphide through a partially or fully molten mantle that is often referred to as a'magma ocean. Alternative models invoke percolation of molten metal along an interconnected network (that is, porous flow) through a solid silicate matrix. But experimental studies performed at high pressures have shown that, under hydrostatic conditions, these melts do not form an interconnected network, leading to the widespread assumption that formation of metallic cores requires a magma ocean. In contrast, here we present experiments which demonstrate that shear deformation to large strains can interconnect a significant fraction of initially isolated pockets of metal and metal sulphide melts in a solid matrix of polycrystalline olivine. Therefore, in a dynamic (nonhydrostatic) environment, percolation remains a viable mechanism for the segregation and migration of core-forming melts in a solid silicate mantle.

Bruhn, D.↗

Microphysics, Meteorology, Microwave and Modeling of Mediterranean Storms: The M(sup 5) Problem

Comprehensive understanding of the microphysical nature of Mediterranean storms requires a combination of in situ meteorological data analysis and radar-passive microwave data analysis, effectively integrated with numerical modeling studies at various scales, particularly from synoptic scale down to mesoscale. The microphysical properties of and their controls on severe storms are intrinsically related to meteorological processes under which storms have evolved, processes which eventually select and control the dominant microphysical properties themselves. Insofar as hazardous Mediterranean storms, highlighted by the September 25-28/1992 Genova flood event, the October 5-7/1998 Friuli flood event, and the October 13-15/2000 Piemonte flood event (all taking place in northern Italy), developing a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution and the heterogeneous nature of precipitation fields within the storm domains. This involves convective development, stratiform transition and decay, orographic lifting, and sloped frontal lifting proc esses. This also involves vertical motions and thermodynamical instabilities governing physical processes that determine details of the liquid/ice water contents, size distributions, and fall rates of the various modes of hydrometeors found within the storm environments. This paper presents detailed 4-dimensional analyses of the microphysical elements of the three severe Mediterranean storms identified above, investigated with the aid of SSM/I and TRMM satellite measurements (and other remote sensing measurements). The analyses are guided by nonhydrostatic mesoscale model simulations at high resolution of the intense rain producing portions of the storm environments. The results emphasize how meteorological controls taking place at the large scale, coupled with localized terrain controls, ultimately determine the most salient features of the bulk microphysical properties of the storms. These results have bearing on precipitation remote sensing from space, and the role of modeling in designing precipitation retrieval algorithms.

Smith, Eric A.↗

The Structure of A Pacific Narrow Cold Frontal Rainband

A NOAA P-3 instrumented aircraft observed an intense, fast-moving narrow cold frontal Farmhand as it approached the Pacific Northwest coast on 19 February 2001 during the Pacific Coastal Jets Experiment. Pseudo-dual-Doppler analyses performed on the airborne Doppler radar data while the frontal system was well offshore indicated that a narrow ribbon of very high radar reflectively convective cores characterized the Farmhand at low levels with echo tops to approximately 4-5 km. The NCFR exhibited gaps in its narrow ribbon of high reflectively, probably as a result of hydrodynamic instability all no its advancing cold pool leading edge. In contrast to some earlier studies of cold frontal rainbands, density current theory described well the motion of the overall front. The character of the updraft structure associated with the heavy rainfall at its leading edge varied across the gap region. The vertical shear of the cross-frontal low-level ambient flow exerted a strong influence on the updraft character, consistent with theoretical arguments developed for squall lines describing the balance of vorticity at the leading edge. In short regions south of the gaps the vertical wind shear was strongest with the updrafts and rain shafts more intense, narrower, and more erect or even downshear tilted. North of the gaps the wind shear weakened with less intense Dihedrals which tilted upshear with a broader band of rainfall. Simulations using a nonhydrostatic mesoscale nested grid model are used to investigate the gap regions, particularly the balance of cold pool induced to pre-frontal ambient shears at the leading edge. Observations confirm the model results that the updraft character depends on the balance of vorticity at the leading edge. Downshear-tilted updrafts imply that convection south of the gap regions would weaken with time relative to the frontal segments north of the gaps since inflow air would be affected by passage through the heavy rain region before ascent, suggesting a mechanism for gap filling.

Jorgensen, David P.↗

GMS-based"Future Time" Rainfall Data Assimilation for Mesoscale Weather Prediction over Korean Peninsula and Future Prospects with GPM Satellite Measurements

This study examines the use of satellite-derived nowcasted (short-term forecasted) rainfall over 3-hour time periods to gain an equivalent time increment in initializing a nonhydrostatic mesoscale model used for predicting convective rainfall events over the Korean peninsula. Infrared (IR) window measurements from the Japanese Geostationary Meteorological Satellite (GMS) are used to specify latent heating for a spinup period of the model - but in future time -- thus initializing in advance of actual time in the framework of a prediction scenario. The main scientific objective of the study is to investigate the strengths and weaknesses of this approach insofar as data assimilation, in which the nowcasted assimilation data are derived independently of the prognostic model itself. Although there have been various recent improvements in formulating the dynamics, thermodynamics, and microphysics of mesoscale models, as well as computer advances which allow the use of high resolution cloud-resolving grids and explicit latent heating over regional domains, spinup remains at the forefront of unresolved mesoscale modeling problems. In general, non-realistic spinup limits the skill in predicting the spatial-temporal distribution of convection and precipitation, primarily in the early hours of a. forecast, stemming from standard prognostic variables not representing the initial diabatic heating field produced by the ambient convection and cloud fields. The long-term goal of this research is to improve short-range (12-hour) quantitative precipitation forecasting (QPF) over the Korean peninsula through the use of innovative data assimilation methods based on geosynchronous satellite measurements. As a step in th~s direction, a non-standard data assimilation experiment in conjunction with GMS-retrieved nowcasted rainfall information introduced to the mesoscale model is conducted. The 3-hourly precipitation forecast information is assimilated through nudging the associated diabatic heating during the early stages of a forecast period. This procedure is expected to enhance details in the moisture field during model integration, and thus improve spinup performance, assuming the errors in the future time latent heating data ate less than intrinsic model background errors.

Smith, Eric A.↗

Comparisons and Contrasts Between Mediterranean Flood-Producing Orographic Storms and Tropical Cyclones Modified by Complex Orography

Often, late-year, flood-producing storms in the Mediterranean Alps region arise from baroclinic disturbances developing on westerly longwaves, with the incipient storm circulations appearing well west of the flood zone. Regardless of the genesis locations and before the storms begin their flow interactions with the orography of the Mediterranean Alps, a process which can produce intense long-lived rainfall over the mountainous terrain of Italy and France, most of their moisture is drawn from warm Mediterranean waters while important elements of their thermodynamic structure are derived from airmass properties originating in north Africa and the Mediterranean basin itself. In fact, because of strong thermodynamic contrasts between the baroclinic zone producing the incipient disturbance and high potential temperature air from the south advecting into and wrapping around the storm center, a Mediterranean storm whose origins are from an evolving shortwave baroclinic perturbation, can become a quasi-warm core hybrid in term of its energetics composition, a process somewhat akin to the development of polar lows. This process begs the question of whether there are similarities between the autumn-winter season flood producing storms of the Mediterranean Alps and warm core tropical cyclones which become altered dynamically and thermodynamically by landfall over steep orographic terrain. This study investigates, from a comparison and contrast perspective, how the similarities and differences between the two types of storms manifest themselves over their respective life cycles. Before the orographically-induced precipitation-release stage in which both types of storm have evolved to the same hydrometeorlogical conclusion, the two storm categories are transforming from entirely distinctive origins, i.e., the classic eastward traveling, baroclinic mid-latitude cyclone forming in western Europe, and the classic eastward-traveling, warm core tropical cyclone developing in the western Pacific. The analyses use high resolution, nonhydrostatic model simulations, satellite data, and various types of ground data to conduct the study. The flood events of Piemonte- 2000 in northern Italy and Typhoon Rusa-2002 in South Korea are selected to focus the analysis on well documented meteorological events.

Smith, Eric A.↗

An Investigation of the Characterization of Cloud Contamination in Hyperspectral Radiances

In regions lacking direct observations, the assimilation of radiances from infrared and microwave sounders is the primary method for characterizing the atmosphere in the analysis process. In recent years, technological advances have led to the launching of more advanced sounders, particularly in the thermal infrared spectrum. With the advent of these hyperspectral sounders, the amount of data available for the analysis process has and will continue to be dramatically increased. However, the utilization of infrared radiances in variational assimilation can be problematic in the presence of clouds; specifically the assessment of the presence of clouds in an instantaneous field of view (IFOV) and the contamination in the individual channels within the IFOV. Various techniques have been developed to determine if a channel is contaminated by clouds. The work presented in this paper and subsequent presentation will investigate traditional techniques and compare them to a new technique, the C02 sorting technique, which utilizes the high spectral resolution of the Atmospheric Infrared Sounder (AIRS) within the framework of the Gridpoint Statistical Interpolation (GSI) 3DVAR system. Ultimately, this work is done in preparation for the assessment of short-term forecast impacts with the regional assimilation of AIRS radiances within the analysis fields of the Weather Research and Forecast Nonhydrostatic Mesoscale Model (WRF-NMM) at the NASA Short-term Prediction Research and Transition (SPORT) Center.

McCarty, William↗

Internal Structure of Rhea

We model the interior of Rhea on the basis of observational constraints and the results from geodynamical models available in the literature. Ten main types of models are defined, depending on the presence or absence of a high-pressure ice layer (ice II), and the extent of separation of the rock component from the volatiles. The degree-two gravity coefficients are computed for each of these models in order to assess which properties of the interior are likely to be inferred from Cassini radio science measurements scheduled on 26 November 2005. C22 greater than 2.5 x 10(exp -4) indicates that the satellite is undifferentiated, except for a slight increase in density with depth resulting from material self-compression. C22 between 1.67 x 10(exp -4)(lower bound) and 1.90 x 10(exp -4) indicates the presence of a rocky core, whose radius can be determined from the satellite's mass and ices densities, for a given temperature profile. For other values, most of the ten models cannot be distinguished from each other. However, assumptions on the density of the rock phase, presence or absence of ice II, and the degree of differentiation could allow a unique model to be determined in many cases. While the calculation presented in this work assumes that Rhea is in hydrostatic equilibrium, it is likely that Rhea' gravity field is partly affected by nonhydrostatic anomalies.

gravity science↗

Impact of MODIS High-Resolution Sea-Surface Temperatures on WRF Forecasts at NWS Miami, FL

Over the past few years,studies at the Short-term Prediction Research and Transition (SPoRT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) composite sea-surface temperature (SST) products in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. The recent paper by LaCasse et al. (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPoRT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The scientific hypothesis being tested is: More accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running the WRF system in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software; The EMS is a standalone modeling system capable of downloading the necessary daily datasets, and initializing, running and displaying WRF forecasts in the NWS Advanced Weather Interactive Processing System (AWIPS) with little intervention required by forecasters. Twenty-seven hour forecasts are run daily with start times of 0300,0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and the far western portions of the Bahamas, the Florida Keys, the Straights of Florida, and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS, invoking the diabatic. "hot-start" capability. In this WRF model "hot-start", the LAPS-analyzed cloud and precipitation features are converted into model microphysics fields with enhanced vertical velocity profiles, effectively reducing the model spin-up time required to predict precipitation systems. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at l/12 degree resolution (approx. 9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPoRT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA in every respect except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. The MODIS SST composites for initializing the SPoRT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST composites into the SPoRTWRF runs is staggered such that the 0400UTC composite initializes the 0900 UTC WRF, the 0700 UTC composite initializes the 1500 UTC WRF, the 1600 UTC composite initializes the 2100 UTC WRF, and the 1900 UTC composite initializes the 0300 UTC WRF. A comparison of the SPoRT and Miami forecasts is underway in 2007, and includes quantitative verification of near-surface temperature, dewpoint, and wind forecasts at surface observation locations. In addition, particular days of interest are being analyzed to determine the impact of the MODIS SST data on the development and evolution of predicted sea/land-breeze circulations, clouds, and precipitation. This paper will present verification results comparing the NWS MIA forecasts the SPoRT experimental WRF forecasts, and highlight any substantial differences noted in the predicted mesoscale phenomena.

Case, Jonathan L.↗

A Multi-Season Study of the Effects of MODIS Sea-Surface Temperatures on Operational WRF Forecasts at NWS Miami, FL

Studies at the Short-term Prediction Research and Transition (SPORT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) sea-surface temperature (SST) composites in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. Recent work by LaCasse et al (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPORT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The project's goal is to determine whether more accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running WRF in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software. Twenty-seven hour forecasts are run dally initialized at 0300, 0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at 1/12deg resolution (approx.9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPORT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water, The MODIS SST composites for initializing the SPORT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST data into the SPORT WRF runs is staggered such that SSTs are updated with a new composite every six hours in each of the WRF runs. From mid-February to July 2007, over 500 parallel WRF simulations have been collected for analysis and verification. This paper will present verification results comparing the NWS MIA operational WRF runs to the SPORT experimental runs, and highlight any substantial differences noted in the predicted mesoscale phenomena for specific cases.

Case, Jonathan L.↗