Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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 523 records · Page 29

Advances in Land Data Assimilation at the NASA Goddard Space Flight Center

Research in land surface data assimilation has grown rapidly over the last decade. In this presentation we provide a brief overview of key research contributions by the NASA Goddard Space Flight Center (GSFC). The GSFC contributions to land assimilation primarily include the continued development and application of the Land Information System (US) and the ensemble Kalman filter (EnKF). In particular, we have developed a method to generate perturbation fields that are correlated in space, time, and across variables and that permit the flexible modeling of errors in land surface models and observations, along with an adaptive filtering approach that estimates observation and model error input parameters. A percentile-based scaling method that addresses soil moisture biases in model and observational estimates opened the path to the successful application of land data assimilation to satellite retrievals of surface soil moisture. Assimilation of AMSR-E surface soil moisture retrievals into the NASA Catchment model provided superior surface and root zone assimilation products (when validated against in situ measurements and compared to the model estimates or satellite observations alone). The multi-model capabilities of US were used to investigate the role of subsurface physics in the assimilation of surface soil moisture observations. Results indicate that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Building on this experience, GSFC leads the development of the Level 4 Surface and Root-Zone Soil Moisture (L4_SM) product for the planned NASA Soil-Moisture-Active-Passive (SMAP) mission. A key milestone was the design and execution of an Observing System Simulation Experiment that quantified the contribution of soil moisture retrievals to land data assimilation products as a function of retrieval and land model skill and yielded an estimate of the error budget for the SMAP L4_SM product. Terrestrial water storage observations from GRACE satellite system were also successfully assimilated into the NASA Catchment model and provided improved estimates of groundwater variability when compared to the model estimates alone. Moreover, satellite-based land surface temperature (LST) observations from the ISCCP archive were assimilated using a bias estimation module that was specifically designed for LST assimilation. As with soil moisture, LST assimilation provides modest yet statistically significant improvements when compared to the model or satellite observations alone. To achieve the improvement, however, the LST assimilation algorithm must be adapted to the specific formulation of LST in the land model. An improved method for the assimilation of snow cover observations was also developed. Finally, the coupling of LIS to the mesoscale Weather Research and Forecasting (WRF) model enabled investigations into how the sensitivity of land-atmosphere interactions to the specific choice of planetary boundary layer scheme and land surface model varies across surface moisture regimes, and how it can be quantified and evaluated against observations. The on-going development and integration of land assimilation modules into the Land Information System will enable the use of GSFC software with a variety of land models and make it accessible to the research community.

Reichle, Rolf↗

Evaluation of Enhanced High Resolution MODIS/AMSR-E SSTs and the Impact on Regional Weather Forecast

Over the last few years, the NASA Short-term Prediction Research and Transition (SPoRT) Center has been generating a 1-km sea surface temperature (SST) composite derived from retrievals of the Moderate Resolution Imaging Spectroradiometer (MODIS) for use in operational diagnostics and regional model initialization. With the assumption that the day-to-day variation in the SST is nominal, individual MODIS passes aboard the Earth Observing System (EOS) Aqua and Terra satellites are used to create and update four composite SST products each day at 0400, 0700, 1600, and 1900 UTC, valid over the western Atlantic and Caribbean waters. A six month study from February to August 2007 over the marine areas surrounding southern Florida was conducted to compare the use of the MODIS SST composite versus the Real-Time Global SST analysis to initialize the Weather Research and Forecasting (WRF) model. Substantial changes in the forecast heat fluxes were seen at times in the marine boundary layer, but relatively little overall improvement was measured in the sensible weather elements. The limited improvement in the WRF model forecasts could be attributed to the diurnal changes in SST seen in the MODIS SST composites but not accounted for by the model. Furthermore, cloud contamination caused extended periods when individual passes of MODIS were unable to update the SSTs, leading to substantial SST latency and a cool bias during the early summer months. In order to alleviate the latency problems, the SPoRT Center recently enhanced its MODIS SST composite by incorporating information from the Advanced Microwave Scanning Radiometer-EOS (AMSR-E) instruments as well as the Operational Sea Surface Temperature and Sea Ice Analysis. These enhancements substantially decreased the latency due to cloud cover and improved the bias and correlation of the composites at available marine point observations. While these enhancements improved upon the modeled cold bias using the original MODIS SSTs, the discernable impacts on the WRF model were still somewhat limited. This paper explores several factors that may have contributed to this result. First, the original methodology to initialize the model used the most recent SST composite available in a hypothetical real ]time configuration, often matching the forecast initial time with an SST field that was 5-8 hours offset. To minimize the differences that result from the diurnal variations in SST, the previous day fs SST composite is incorporated at a time closest to the model initialization hour (e.g. 1600 UTC composite at 1500 UTC model initialization). Second, the diurnal change seen in the MODIS SST composites was not represented by the WRF model in previous simulations, since the SSTs were held constant throughout the model integration. To address this issue, we explore the use of a water skin-temperature diurnal cycle prediction capability within v3.1 of the WRF model to better represent fluctuations in marine surface forcing. Finally, the verification of the WRF model is limited to very few over-water sites, many of which are located near the coastlines. In order to measure the open ocean improvements from the AMSR-E, we could use an independent 2-dimensional, satellite-derived data set to validate the forecast model by applying an object-based verification method. Such a validation technique could aid in better understanding the benefits of the mesoscale SST spatial structure to regional models applications.

Schiferl, Luke D.↗

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

simulation↗

Optically thin cirrus clouds - Radiative impact on the warm pool

The role of the radiative effects of thin cirrus clouds in the energy balance of the 'warm pool' region is examined with reference to data obtained by the Infrared Interferometer Spectrometer (IRIS) flown on the Nimbus 4 satellite in 1970. First, the IRIS observations and a method for deriving the thin cirrus information are briefly discussed. A simple radiative energy balance model is then developed which is applicable to the mesoscale convective systems in the 'warm pool'. The radiative energy balance model, which does not explicitly account for the transports, is used to demonstrate the net radiative heating produced by the thin cirrus clouds.

Prabhakara, C.↗

Proceedings of Shuttle Environmental Effects Program Review

Measurements of Titan exhaust cloud effluents are documented and compared, mesoscale and microphysical acid rain models are described, and a submesoscale model is proposed. Various instruments and facilities for measuring ice nuclei and other constituents of solid rocket motor exhaust effluents are discussed. Regional air quality monitoring and rain collection systems are described, and the ecological impact of solid rocket motor exhaust effluents is examined. The potential effect of space shuttle launches is estimated where data are adequate.

Potter, A. E.↗

Numerical simulations of a transverse indirect circulation and low-level jet in the exit region of an upper-level jet

A numerical study was performed of a severe weather event (tornado) which occurred on May 10, 1973 in the Ohio region. The situation was modeled with a primitive equation mesoscale dynamic formulation. Account was taken of precipitation, the planetary boundary layer parameters as bulk quantities, the vertical pressure gradient, and lateral boundary conditions based on radiosonde data. Two 12-hr simulations, adiabatic and nondivergent, respectively, were analyzed for relationships between upper and lower level jets. In the adiabatic formulation, a transverse circulation with a low level jet formed at the exit region of the upper level jet. The nondivergent situation led to similar, but weaker, phenomena. Both forms suggest that indirect circulation in the exit zone of an upper level jet is strongly influenced by the initial structure of the jet.

Brill, K. F.↗

Numerical simulations of the subsynoptic features associated with the AVE-SESAME I case. I - The preconvective environment

The extensive diagnostic calculations made possible by the AVE-SESAME I database are used in combination with numerical simulations from the Mesoscale Atmospheric Simulation System (MASS) model to examine the dynamics of the meso-alpha-scale features during the preconvective period from 1130 to 2030 UTC on April 10, 1979. The version of the MASS model used in this investigation is presented, and an overview of the general synoptic conditions present at the time of model initialization is presented along with the data used to initialize the model. The dynamical processes present in the numerical simulations are presented and compared with analyses of the observational data from this and other investigations of this case. The relative importance of the adiabatic and diabatic processes in creating and then initiating the release of the convective instability is discussed.

Zack, John W.↗

A mesoscale numerical simulation of an index for use in the forecasting of the potential severity of wet microbursts

This paper demonstrates the feasibility of using the wet microburst potential intensity index developed by Proctor (1989) for predicting the onset of a wet microburst. After describing the index and the principles involved in its derivation, the version 4.0 of the Mesoscale Atmospheric Simulation System (MASS) model described by Karyampudi et al. (1988) was used to illustrate how the index would perform as a forecast product. An analysis is presented of the use of the microburst potential intensity index in the numerical simulation of the DFR aircraft accident on August 2, 1985 in a downburst environment.

Zack, John W.↗

The role of the sea-surface temperature distribution on numerically simulated cyclogenesis during ERICA

The goal was to quantify the extent to which a sea surface temperature (SST) front can influence cyclogenesis. The approach was to use the Drexel Limited-Area Mesoscale Prediction System (LAMPS) dynamical model to simulate cyclogenesis over various SST fields. Research during the past year focused on the development and testing of a four dimensional data assimilation (FDDA) technique within LAMPS. The technique is a continuous dynamical assimilation where forcing terms are added to the governing model equations to gradually nudge the model solution toward a gridded analysis. Here, the nudging is used as a dynamic initialization tool during a 12 hour preforecast to generate model balanced initial conditions for a subsequent 24 hour numerical prediction. Tests were performed to determine which variables to nudge and how to specify the four dimensional weighting function used to scale the nudging terms. To date, optimal results were obtained by nudging the u and v components of the wind along with the potential temperature. The weighting function ranged from 0 to 1 and varies in time as a quadratic polynomial. It was initialized at 0, reached its maximum at 9 hours into the preforecast, and fell back at 0 to 12 hours. The nudging terms are included in the model equations for all grid points except those within the model predicted oceanic boundary layer. This design attempts to confine changes imposed by the specified SST field to the oceanic boundary layer during the preforecast period.

Lapenta, William M.↗

Application of a linear spectral model to the study of Amazonian squall lines during GTE/ABLE 2B

A linear nonhydrostatic spectral model is run with the basic state, or large scale, vertical profiles of temperature and wind observed prior to convective development along the northern coast of South America during the GTE/ABLE 2B. The model produces unstable modes with mesoscale wavelength and propagation speed comparable to observed Amazonian squall lines. Several tests with different vertical profiles of low-level winds lead to the conclusion that a shallow and/or weak low-level jet either does not produce a scale selection or, if it does, the selected mode is stationary, indicating the absence of a propagating disturbance. A 700-mbar jet of 13 m/s, with a 600-mbar wind speed greater or equal to 10 m/s, is enough to produce unstable modes with propagating features resembling those of observed Amazonian squall lines. However, a deep layer of moderate winds (about 10 m/s) may produce similar results even in the absence of a low-level wind maximum. The implications in terms of short-term weather forecasting are discussed.

Silva Dias, Maria A. F.↗

Evaluation of cloud shading effects on the generation and modification of mesoscale circulations

This study presents scale analysis and numerical model evaluations of the impact of cloud shading on (1) the development of sea breeze and thermally induced upslope flows and (2) the generation of mesoscale circulations between cloudy areas adjacent to clear areas. Based on the assumption of modifications of solar and longwave radiation which are typical for some overcast conditions, it was found that the reduction in the first type of circulation is significant and most noticeable in the vertical velocities. In the second case, thermally induced circulations in favorable conditions appear to approach the typical intensity of the sea breeze.

Segal, M.↗

Evaluation of the 29-km Eta Model for Weather Support to the United States Space Program

The Applied Meteorology Unit (AMU) conducted a year-long evaluation of NCEP's 29-km mesoscale Eta (meso-eta) weather prediction model in order to identify added value to forecast operations in support of the United States space program. The evaluation was stratified over warm and cool seasons and considered both objective and subjective verification methodologies. Objective verification results generally indicate that meso-eta model point forecasts at selected stations exhibit minimal error growth in terms of RMS errors and are reasonably unbiased. Conversely, results from the subjective verification demonstrate that model forecasts of developing weather events such as thunderstorms, sea breezes, and cold fronts, are not always as accurate as implied by the seasonal error statistics. Sea-breeze case studies reveal that the model generates a dynamically-consistent thermally direct circulation over the Florida peninsula, although at a larger scale than observed. Thunderstorm verification reveals that the meso-eta model is capable of predicting areas of organized convection, particularly during the late afternoon hours but is not capable of forecasting individual thunderstorms. Verification of cold fronts during the cool season reveals that the model is capable of forecasting a majority of cold frontal passages through east central Florida to within +1-h of observed frontal passage.

Manobianco, John↗

Evaluation of Cloud-Resolving and Limited Area Model Intercomparison Simulations Using TWP-ICE Observations: Precipitation Microphysics - Part 2

Ten 3-D cloud-resolving model (CRM) simulations and four 3-D limited area model (LAM) simulations of an intense mesoscale convective system observed on 23-24 January 2006 during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) are compared with each other and with observations and retrievals from a scanning polarimetric radar, colocated UHF and VHF vertical profilers, and a Joss-Waldvogel disdrometer in an attempt to explain a low bias in simulated stratiform rainfall. Despite different forcing methodologies, similar precipitation microphysics errors appear in CRMs and LAMs with differences that depend on the details of the bulk microphysics scheme used. One-moment schemes produce too many small raindrops, which biases Doppler velocities low, but produces rainwater contents (RWCs) that are similar to observed. Two-moment rain schemes with a gamma shape parameter (mu) of 0 produce excessive size sorting, which leads to larger Doppler velocities than those produced in one-moment schemes but lower RWCs. Two-moment schemes also produce a convective median volume diameter distribution that is too broad relative to observations and, thus, may have issues balancing raindrop formation, collision-coalescence, and raindrop breakup. Assuming a mu of 2.5 rather than 0 for the raindrop size distribution improves one-moment scheme biases, and allowing mu to have values greater than 0 may improve excessive size sorting in two-moment schemes. Underpredicted stratiform rain rates are associated with underpredicted ice water contents at the melting level rather than excessive rain evaporation, in turn likely associated with convective detrainment that is too high in the troposphere and mesoscale circulations that are too weak. A limited domain size also prevents a large, well-developed stratiform region like the one observed from developing in CRMs, although LAMs also fail to produce such a region.

three dimensional models↗

Ongoing Breakthroughs in Convective Parameterization

While the increase of computer power mobilizes a part of the community towards models with explicit convection or based on machine learning, we review the part of the literature dedicated to convective parameterization development for large-scale forecast and climate models. Recent findings: Many developments are underway to overcome endemic limitations of traditional convective parameterizations, either in unified or multi-object frameworks: scale-aware and stochastic approaches, new prognostic equations or representations of new components such as cold pools. Understanding their impact on the emergent properties of a model remains challenging, due to subsequent tuning of parameters and the limited understanding given by traditional metrics. Summary: Further effort still needs to be dedicated to the representation of the life cycle of convective systems, in particular their mesoscale organization and associated cloud cover. The development of more process-oriented metrics based on new observations is also needed to help quantify model improvement and better understand the mechanisms of climate change.

parameterizations for large-scale models↗

Configuring the HYSPLIT Model for National Weather Service Forecast Office and Spaceflight Meteorology Group Applications

For expedience in delivering dispersion guidance in the diversity of operational situations, National Weather Service Melbourne (MLB) and Spaceflight Meteorology Group (SMG) are becoming increasingly reliant on the PC-based version of the HYSPLIT model run through a graphical user interface (GUI). While the GUI offers unique advantages when compared to traditional methods, it is difficult for forecasters to run and manage in an operational environment. To alleviate the difficulty in providing scheduled real-time trajectory and concentration guidance, the Applied Meteorology Unit (AMU) configured a Linux version of the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) (HYSPLIT) model that ingests the National Centers for Environmental Prediction (NCEP) guidance, such as the North American Mesoscale (NAM) and the Rapid Update Cycle (RUC) models. The AMU configured the HYSPLIT system to automatically download the NCEP model products, convert the meteorological grids into HYSPLIT binary format, run the model from several pre-selected latitude/longitude sites, and post-process the data to create output graphics. In addition, the AMU configured several software programs to convert local Weather Research and Forecast (WRF) model output into HYSPLIT format.

Dreher, Joseph G.↗

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↗

Continuous assimilation of simulated Geosat altimetric sea level into an eddy-resolving numerical ocean model. I - Sea level differences. II - Referenced sea level differences

The optimal interpolation method of Lorenc (1981) was used to conduct continuous assimilation of altimetric sea level differences from the simulated Geosat exact repeat mission (ERM) into a three-layer quasi-geostrophic eddy-resolving numerical ocean box model that simulates the statistics of mesoscale eddy activity in the western North Pacific. Assimilation was conducted continuously as the Geosat tracks appeared in simulated real time/space, with each track repeating every 17 days, but occurring at different times and locations within the 17-day period, as would have occurred in a realistic nowcast situation. This interpolation method was also used to conduct the assimilation of referenced altimetric sea level differences into the same model, performing the referencing of altimetric sea sevel differences by using the simulated sea level. The results of this dynamical interpolation procedure are compared with those of a statistical (i.e., optimum) interpolation procedure.

White, Warren B.↗

Mesoscale simulations of the November 25-26 and December 5-6 cirrus cases using the RAMS model

The Regional Atmospheric Modeling System (RAMS), developed at Colorado State University, was used during the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) 2 (13 Nov. through 6 Dec. 1991) to provide real time forecasts of cirrus clouds. Forecasts were run once a day, initializing with the 0000 UTC dataset provided by NOAA (Forecast Systems Laboratory (FSL) Mesoscale Analysis and Prediction System (MAPS)). In order to obtain better agreement with observations, a second set of simulations were done for the FIRE 2 cases that occurred on 25-26 Nov. and 5-6 Dec. In this set of simulations, a more complex radiation scheme was used, the Chen/Cotton radiation scheme, along with the nucleation of ice occurring at ice supersaturations as opposed to nucleation occurring at water supersaturations that was done in the actual forecast version. The runs using these more complex schemes took longer wall clock time (7-9 hours for the actual forecasts as compared to 12-14 hrs for the runs using the more complex schemes) however, the final results of the simulations were definitely improved upon. Comparisons between these two sets of simulations are given. Now underway are simulations of these cases using a closed analytical solution for the auto-conversion of ice from a pristine ice class (sizes less than about 50 microns in effective diameter) to a snow class (effective diameters on the order of several hundred microns). This solution is employed along with a new scheme for the nucleation of ice crystals due to Meyers et al and Demott et al. The scheme is derived assuming complete gamma distributions for both the pristine and snow classes. The time rate of change of the number concentration and mass mixing-ratio of each distribution is found by calculating either the flux of crystals that grow beyond a certain critical diameter by vapor deposition in an ice supersaturated regime or by calculating the flux of crystals that evaporate to sizes below that same critical effective diameter.

Harrington, J. L.↗