Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cloud parameterizations”

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.

720 records · Page 40

Mesoscale Simulations of a Florida Sea Breeze Using the PLACE Land Surface Model Coupled to a 1.5-Order Turbulence Parameterization

A sophisticated land-surface model, PLACE, the Parameterization for Land Atmospheric Convective Exchange, has been coupled to a 1.5-order turbulent kinetic energy (TKE) turbulence sub-model. Both have been incorporated into the Penn State/National Center for Atmospheric Research (PSU/NCAR) mesoscale model MM5. Such model improvements should have their greatest effect in conditions where surface contrasts dominate over dynamic processes, such as the simulation of warm-season, convective events. A validation study used the newly coupled model, MM5 TKE-PLACE, to simulate the evolution of Florida sea-breeze moist convection during the Convection and Precipitation Electrification Experiment (CaPE). Overall, eight simulations tested the sensitivity of the MM5 model to combinations of the new and default model physics, and initialization of soil moisture and temperature. The TKE-PLACE model produced more realistic surface sensible heat flux, lower biases for surface variables, more realistic rainfall, and cloud cover than the default model. Of the 8 simulations with different factors (i.e., model physics or initialization), TKE-PLACE compared very well when each simulation was ranked in terms of biases of the surface variables and rainfall, and percent and root mean square of cloud cover. A factor separation analysis showed that a successful simulation required the inclusion of a multi-layered, land surface soil vegetation model, realistic initial soil moisture, and higher order closure of the planetary boundary layer (PBL). These were needed to realistically model the effect of individual, joint, and synergistic contributions from the land surface and PBL on the CAPE sea-breeze, Lake Okeechobee lake breeze, and moist convection.

Lynn, Barry H.↗

Ocean-Atmosphere Interaction Over Agulhas Extension Meanders

Many years of high-resolution measurements by a number of space-based sensors and from Lagrangian drifters became available recently and are used to examine the persistent atmospheric imprints of the semi-permanent meanders of the Agulhas Extension Current (AEC), where strong surface current and temperature gradients are found. The sea surface temperature (SST) measured by the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and the chlorophyll concentration measured by the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) support the identification of the meanders and related ocean circulation by the drifters. The collocation of high and low magnitudes of equivalent neutral wind (ENW) measured by Quick Scatterometer (QuikSCAT), which is uniquely related to surface stress by definition, illustrates not only the stability dependence of turbulent mixing but also the unique stress measuring capability of the scatterometer. The observed rotation of ENW in opposition to the rotation of the surface current clearly demonstrates that the scatterometer measures stress rather than winds. The clear differences between the distributions of wind and stress and the possible inadequacy of turbulent parameterization affirm the need of surface stress vector measurements, which were not available before the scatterometers. The opposite sign of the stress vorticity to current vorticity implies that the atmosphere spins down the current rotation through momentum transport. Coincident high SST and ENW over the southern extension of the meander enhance evaporation and latent heat flux, which cools the ocean. The atmosphere is found to provide negative feedback to ocean current and temperature gradients. Distribution of ENW convergence implies ascending motion on the downwind side of local SST maxima and descending air on the upwind side and acceleration of surface wind stress over warm water (deceleration over cool water); the convection may escalate the contrast of ENW over warm and cool water set up by the dependence of turbulent mixing on stability; this relation exerts a positive feedback to the ENW-SST relation. The temperature sounding measured by the Atmospheric Infrared Sounder(AIRS) is consistent with the spatial coherence between the cloud-top temperature provided by the International Satellite Cloud Climatology Project (ISCCP) and SST. Thus ocean mesoscale SST anomalies associated with the persistent meanders may have a long-term effect well above the midlatitude atmospheric boundary layer, an observation not addressed in the past.

Atmospheric Infrared Sounder (AIRS)↗

Studying Diurnal Variations of Aerosols with NASA MERRA-2 Reanalysis Data

Aerosols play an important role in atmospheric dynamics, climate variations, and Earth's energy cycle by altering the radiation balance in the atmosphere through interaction with clouds, providing fertilizer for forests and canopy, and as a supply of iron to the ocean over long time periods. Studies suggest that much of the feedback between dust aerosols and dynamics is associated with diurnal and synoptic scale variability. However, the lack of sub-daily resolution of aerosols from satellite observations makes it difficult to study the diurnal characteristics, especially over tropical and subtropical regions. Investigation of this topic utilizes over 37 years of simulated global aerosol products from NASA atmospheric reanalysis, in the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) data set, available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). MERRA-2 covers the period 1980-present, and is continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using data from MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated using the MERRA-2 aerosol model, which interacts directly with radiation parameterization, and is radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Hourly, monthly, and monthly diurnal data are available at spatial resolution of 0.5o x 0.625o (latitude x longitude). By using MERRA-2 hourly and monthly diurnal products, different aerosol diurnal variabilities are observed over North America, Africa, Asia, and Australia, that may be due to different meteorological conditions and aerosol sources. The presentation will also provide an overview of MERRA-2 data services at GES DISC, such as how to find and download data, and how to quickly visualize and analyze data online with Giovanni.

aerosol↗

Constraints on Cumulus Parameterization from Simulations of Observed MJO Events

Two recent activities offer an opportunity to test general circulation model (GCM) convection and its interaction with large-scale dynamics for observed Madden-Julian oscillation (MJO) events. This study evaluates the sensitivity of the Goddard Institute for Space Studies (GISS) GCM to entrainment, rain evaporation, downdrafts, and cold pools. Single Column Model versions that restrict weakly entraining convection produce the most realistic dependence of convection depth on column water vapor (CWV) during the Atmospheric Radiation Measurement MJO Investigation Experiment at Gan Island. Differences among models are primarily at intermediate CWV where the transition from shallow to deeper convection occurs. GCM 20-day hindcasts during the Year of Tropical Convection that best capture the shallow–deep transition also produce strong MJOs, with significant predictability compared to Tropical Rainfall Measuring Mission data. The dry anomaly east of the disturbance on hindcast day 1 is a good predictor of MJO onset and evolution. Initial CWV there is near the shallow–deep transition point, implicating premature onset of deep convection as a predictor of a poor MJO simulation. Convection weakly moistens the dry region in good MJO simulations in the first week; weakening of large-scale subsidence over this time may also affect MJO onset. Longwave radiation anomalies are weakest in the worst model version, consistent with previous analyses of cloud/moisture greenhouse enhancement as the primary MJO energy source. The authors’ results suggest that both cloud-/moisture-radiative interactions and convection–moisture sensitivity are required to produce a successful MJO simulation.

Madden-Julian Oscillation↗

Clear-Sky Closure Studies of Tropospheric Aerosol and Water Vapor During ACE-2 Using Airborne Sunphotometer, Airborne In-Situ, Space-Borne, and Ground-Based Measurements

We report on clear-sky column closure experiments (CLEARCOLUMN) performed in the Canary Islands during the second Aerosol Characterization Experiment (ACE-2) in June/July 1997. We present CLEARCOLUMN results obtained by combining airborne sunphotometer and in-situ (a differential mobility analyzer, three optical particle counters, three nephelometers, and one absorption photometer) measurements taken aboard the Pelican aircraft, space-borne NOAA/AVHRR data and ground-based lidars. A wide range of aerosol types was encountered throughout the ACE-2 area, including background Atlantic marine, European pollution-derived, and (although less frequently than expected) African mineral dust. During the two days discussed here, vertical profiles flown in cloud free air masses revealed three distinctly different layers: a marine boundary layer (MBL) with varying pollution levels, an elevated dust layer, and a very clean layer between the MBL and the dust layer. Based on size-resolved composition information we have established an aerosol model that allows us to compute optical properties of the ambient aerosol using the optical particle counter results. In the dust, the agreement in layer AOD (lambda=380-1060 nm) is 3-8%. In the MBL there is a tendency for the in-situ results to be slightly lower than the sunphotometer measurements (10-17% at lambda=525 nm), but these differences are within the combined error bars of the measurements and computations. Aerosol size-distribudon closure based on in-situ size distributions and inverted sunphotometer extinction spectra has been achieved in the MBL (total surface area and volume agree within 0.2, and 7%, respectively) but not in the dust layer. The fact that the three nephelometers operated at three different relative humidities (RH) allowed to parameterize hygroscopic growth and to therefore estimate optical properties at ambient RH. The parameters derived for different aerosol types are themselves useful for the aerosol modeling community. The fact that the nephelometers and the absorption photometer sampled the aerosol through a cyclone make those measurements less useful for thee closure study carried out here. Large corrections (especially in the dust) had to be applied. Therefore, it is not surprising that closure with the sunphotometer was not always achieved. Agreement within 0.02 in AOD was achieved in the dust layer when the airorne sunphotometer extinction or AOD was compared to ground-based lidar measurements. We found that the presence of the elevated dust layers removes the good agreement between satellite and sunphotometer AOD usually found in the absence of the dust layer. We still compare the scattering phase functions used in the satellite retrieval with those obtained from the aerosol or the sunphotometer measurements.

Schmid, Beat↗

Goddard Cumulus Ensemble (GCE) Model: Application for Understanding Preciptation Processes

The global hydrological cycle is central to climate system interactions and the key to understanding their behavior. Rainfall and its associated precipitation processes are a key link in the hydrologic cycle. Fresh water provided by tropical rainfall and its variability can exert a large impact upon the structure of the upper ocean layer. In addition, approximately two-thirds of the global rain falls in the Tropics, while the associated latent heat release accounts for about three-fourths of the total heat energy for the Earth's atmosphere. Precipitation from convective cloud systems comprises a large portion of tropical heating and rainfall. Furthermore, the vertical distribution of convective latent-heat releases modulates large-scale tropical circulations (e.g., the 30-60-day intraseasonal oscillation), which, in turn, impacts midlatitude weather through teleconnection patterns such as those associated with El Nino. Shifts in these global circulations can result in prolonged periods of droughts and floods, thereby exerting a tremendous impact upon the biosphere and human habitation. And yet, monthly rainfall over the tropical oceans is still not known within a factor of two over large (5 degrees latitude by 5 degrees longitude) areas. Hence, the Tropical Rainfall Measuring Mission (TRMM), a joint U.S./Japan space project, can provide a more accurate measurement of rainfall as well as estimate the four-dimensional structure of diabatic heating over the global tropics. The distributions of rainfall and inferred heating can be used to advance our understanding of the global energy and water cycle. In addition, this information can be used for global circulation and climate models for testing and improving their parameterizations.

Tao, Wei-Kuo↗

The Fate of Saharan Dust Across the Atlantic and Implications for a Central American Dust Barrier

Saharan dust was observed over the Caribbean basin during the summer 2007 NASA Tropical Composition, Cloud, and Climate Coupling (TC4) field experiment. Airborne Cloud Physics Lidar (CPL) and satellite observations from MODIS suggest a barrier to dust transport across Central America into the eastern Pacific. We use the NASA GEOS-5 atmospheric transport model with online aerosol tracers to perform simulations of the TC4 time period in order to understand the nature of this barrier. Our simulations are driven by the Modem Era Retrospective-Analysis for Research and Applications (MERRA) meteorological analyses. We evaluate our baseline simulated dust distributions using MODIS and CALIOP satellite and ground-based AERONET sun photometer observations. GEOS-5 reproduces the observed location, magnitude, and timing of major dust events, but our baseline simulation does not develop as strong a barrier to dust transport across Central America as observations suggest. Analysis of the dust transport dynamics and lost processes suggest that while both mechanisms play a role in defining the dust transport barrier, loss processes by wet removal of dust are about twice as important as transport. Sensitivity analyses with our model showed that the dust barrier would not exist without convective scavenging over the Caribbean. The best agreement between our model and the observations was obtained when dust wet removal was parameterized to be more aggressive, treating the dust as we do hydrophilic aerosols.

Nowottnick, E.↗

Remote Sensing of Rain

The first problem addressed concerns passive-microwave rain retrievals. Most current approaches start by building off-line a cloud-model-derived database. Given data, the retrieval algorithms search the database for the microwave temperatures "closest" to the observed data, then after some fine-tuning (performed in different ways by different implementations) the rain is estimated to be that which corresponds to the selected (and fine-tuned) set of database temperatures. These approaches have three drawbacks: they cannot properly take into account the ambiguities which arise from the fact that several rain scenarios can produce the same observed temperatures; they are quite inefficient since they require manipulating a large database along with often complex "fine-tuning" procedures; and they cannot refine their estimates if additional data is available. This past year we have derived closed formulae relating observed microwave brightness temperatures, T(sub b), and the underlying rain rates, R: average T(sub b) =f (rain) and average rain = g (T(sub b)), along with the corresponding covariance matrices. These results are sufficient to describe the conditional probabilities p(R/T(sub b)) and p(T(sub b)/R) to second order. Progress has also been made towards deriving a robust description of the rain drop size distribution (DSD). The widespread approach consisting in parameterizing the DSD as a gamma-distribution in terms of the drop diameter D suffers from the facts that, in reality, the DSD is not a smooth function of D and that the largely arbitrary Gamma model imposes unintended behavior, which has implications on any quantities derived from the DSD model. We have therefore developed a non-parametric yet practical description of the DSD, which is particularly well-suited for use in remote-sensing applications. The diagram on the left shows a comparison between an actual DSD sample and the truncated non-parametric representation. One figure shows the relation between radar reflectivity and rain rate derived using this representation. Validation of the Tropical Rainfall Measuring Mission (TRMM) radar-radiometer combined R and DSD algorithm is underway. This algorithm was designed to make optimal use of the instantaneous reflectivity profiles measured by the TRMM radar and the microwave brightness temperatures measured by the TRMM passive radiometer. So far, it appears to be the most reliable TRMM rain algorithm.

Haddad, Ziad S.↗

Evaluating Greenland Surface-Mass-Balance and Firn-Densification Data Using ICESat-2 Altimetry

Surface-mass-balance (SMB) and firn-densification (FD) models are widely used in altimetry studies as a tool to separate atmospheric-driven from ice-dynamics-driven ice-sheet mass changes and to partition observed volume changes into ice-mass changes and firn-air-content changes. Until now, SMB models have been principally validated based on comparison with ice core and weather station data or comparison with widely separated flight radar-survey flight lines. Firn-densification models have been primarily validated based on their ability to match net densification over decades, as recorded in firn cores, and the short-term time-dependent component of densification has rarely been evaluated at all. The advent of systematic ice-sheet-wide repeated ice-surface-height measurements from ICESat-2 (the Ice Cloud, and land Elevation Satellite, 2) allows us to measure the net surface-height change of the Greenland ice sheet at quarterly resolution and compare the measured surface-height differences directly with those predicted by three FD–SMB models: MARv3.5.11 (Modèle Atmosphérique Régional version 3.5.11) and GSFCv1.1 and GSFCv1.2 (the Goddard Space Flight Center FD–SMB models version 1.1 and 1.2). By segregating the data by season and elevation, and based on the timing and magnitude of modelled processes in areas where we expect minimal ice-dynamics-driven height changes, we investigate the models' accuracy in predicting atmospherically driven height changes. We find that while all three models do well in predicting the large seasonal changes in the low-elevation parts of the ice sheet where melt rates are highest, two of the models (MARv3.5.11 and GSFCv1.1) systematically overpredict, by around a factor of 2, the magnitude of height changes in the high-elevation parts of the ice sheet, particularly those associated with melt events. This overprediction seems to be associated with the melt sensitivity of the models in the high-elevation part of the ice sheet. The third model, GSFCv1.2, which has an updated high-elevation melt parameterization, avoids this overprediction.

Benjamin E. Smith↗

Seasonal Distributions of Global Ocean Chlorophyll and Nutrients: Analysis with a Coupled Ocean General Circulation Biogeochemical, and Radiative Model

A coupled general ocean circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. The model is driven by climatological meteorological conditions, cloud cover, and sea surface temperature. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chorophytes, and picoplankton) and three nutrient groups (nitrate, ammonium, and silicate). Phytoplankton groups are initialized as homogeneous fields horizontally and vertically, and allowed to distribute themselves according to the prevailing conditions. Basin-scale model chlorophyll results are in very good agreement with CZCS pigments in virtually every global region. Seasonal variability observed in the CZCS is also well represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are also in good conformance, although occasional departures are apparent. Agreement of nitrate distributions with in situ data is even better, including seasonal dynamics, except for the equatorial Atlantic. The good agreement of the model with satellite and in situ data sources indicates that the model dynamics realistically simulate phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization, and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

Gregg, Watson W.↗

The classification of the Arctic Sea ice types and the determination of surface temperature using advanced very high resolution radiometer data

The accurate quantification of new ice and open water areas and surface temperatures within the sea ice packs is a key to the realistic parameterization of heat, moisture, and turbulence fluxes between ocean and atmosphere in the polar regions. Multispectral NOAA advanced very high resolution radiometer/2 (AVHRR/2) satellite images are analyzed to evaluate how effectively the data can be used to characterize sea ice in the Bering and Greenland seas, both in terms of surface type and physical temperature. The basis of the classification algorithm, which is developed using a late wintertime Bering Sea ice cover data, is that frequency distributions of 10.8- micrometers radiances provide four distinct peaks, represeting open water, new ice, young ice, and thick ice with a snow cover. The results are found to be spatially and temporally consistent. Possible sources of ambiguity, especially associated with wider temporal and spatial application of the technique, are discussed. An ice surface temperature algorithm is developed for the same study area by regressing thermal infrared data from 10.8- and 12.0- micrometers channels against station air temperatures, which are assumed to approximate the skin temperatures of adjacent snow and ice. The standard deviations of the results when compared with in situ data are about 0.5 K over leads and polynyas to about 0.5-1.5 K over thick ice. This study is based upon a set of in situ data limited in scope and coverage. Cloud masks are applied using a thresholding technique that utilizes 3.74- and 10.8- micrometers channel data. The temperature maps produced show coherence with surface features like new ice and leads, and consistency with corresponding surface type maps. Further studies are needed to better understand the effects of both the spatial and temporal variability in emissivity, aerosol and precipitable atmospheric ice particle distribution, and atmospheric temperature inversions.

Massom, Robert↗

An Automatic Light Rain Detection Algorithm on NASA MPLNET Lidar Observations in the Frame of WMO GALION Project

The water cycle strongly influences life on Earth. In particular, the precipitation modifies the atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially low-intensity precipitation) at global scale, bedsides improving our understanding of the hydrological cycle, it is crucial to reduce the associated uncertainty of the global climate models to correctly forecast future scenarios, i.e. to apply fast mitigation strategies. In this study we developed an algorithm to automatically detect precipitation from lidar measurements obtained by the National and Aeronautics Space Administration (NASA) Micropulse lidar network (MPLNET) permanent observational site in Goddard. The algorithm, once full operational, will deliver in Near Real Time (latency 1.5h) a new rain mask product that will be publicly available on MPLNET website as part of the new Version 3 Level 1.5 data. The methodology, based on an image processing technique, can detect only light precipitation events (defined by intensity and duration) as the morphological filters used through the detection process are applied on the lidar volume depolarization ratio range corrected composite images, i.e. heavy rain events are unusable as the lidar signal is completely extinguished after few meters in the precipitation or no signal detected because of the water accumulated on the receiver optics. Results from the algorithm, besides filling a gap in precipitation and virga detection by radars, are of particular interest for the scientific community because will help to better understand long-term aerosol-cloud interactions and aerosol atmospheric removal (scavenging effect) by rain as multi-year database being available for several MPLNET permanent observational sites across the globe. Moreover, we developed the automatic algorithm at Universitat Politecnica de Catalunya (UPC) Barcelona, the unique permanent observation station member of MPLNET and the European Aerosol Lidar Network (EARLINET) In the future the algorithm can be then easily applied to any other lidar and/or ceilometer network infrastructure in the frame of World Meteorological Organization (WMO) Global Aerosol Watch (GAW) aerosol lidar observation network (GALION)

Simone Lolli↗

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data: Part III. Using Combined PCA to Compare Spatiotemporal Variability of MODIS, MISR and OMI Aerosol Optical Depth

Satellite measurements of global aerosol properties are very useful in constraining aerosol parameterization in climate models. The reliability of different data sets in representing global and regional aerosol variability becomes an essential question. In this study, we present the results of a comparison using combined principal component analysis (CPCA), applied to monthly mean, mapped (Level 3) aerosol optical depth (AOD) product from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Ozone Monitoring Instrument (OMI). This technique effectively finds the common space-time variability in the multiple data sets by decomposing the combined AOD field. The results suggest that all of the sensors capture the globally important aerosol regimes, including dust, biomass burning, pollution, and mixed aerosol types. Nonetheless, differences are also noted. Specifically, compared with MISR and OMI, MODIS variability is significantly higher over South America, India, and the Sahel. MODIS deep blue AOD has a lower seasonal variability in North Africa, accompanied by a decreasing trend that is not found in either MISR or OMI AOD data. The narrow swath of MISR results in an underestimation of dust variability over the Taklamakan Desert. The MISR AOD data also exhibit overall lower variability in South America and the Sahel. OMI does not capture the Russian wild fire in 2010 nor the phase shift in biomass burning over East South America compared to Central South America, likely due to cloud contamination and the OMI row anomaly. OMI also indicates a much stronger (boreal) winter peak in South Africa compared with MODIS and MISR.

satellite observation↗

Overview of the New Version 3 NASA Micro-Pulse Lidar Network (MPLNET) Automatic Precipitation Detection Algorithm

Precipitation modifies atmospheric column thermodynamics through the process of evaporation and serves as a proxy for latent heat modulation. For this reason, a correct precipitation parameterization (especially for low-intensity precipitation) within global scale models is crucial. In addition to improving our modeling of the hydrological cycle, this will reduce the associated uncertainty of global climate models in correctly forecasting future scenarios, and will enable the application of mitigation strategies. In this manuscript we present a proof of concept algorithm to automatically detect precipitation from lidar measurements obtained from the National Aeronautics and Space Administration Micropulse lidar network (MPLNET). The algorithm, once tested and validated against other remote sensing instruments, will be operationally implemented into the network to deliver a near real time (latency <1.5 h) rain masking variable that will be publicly available on MPLNET website as part of the new Version 3 data products. The methodology, based on an image processing technique, detects only light precipitation events (defined by intensity and duration) such as light rain, drizzle, and virga. During heavy rain events, the lidar signal is completely extinguished after a few meters in the precipitation or it is unusable because of water accumulated on the receiver optics. Results from the algorithm, in addition to filling a gap in light rain, drizzle, and virga detection by radars, are of particular interest for the scientific community as they help to fully characterize the aerosol cycle, from emission to deposition, as precipitation is a crucial meteorological phenomenon accelerating atmospheric aerosol removal through the scavenging effect. Algorithm results will also help the understanding of long term aerosol–cloud interactions, exploiting the multi-year database from several MPLNET permanent observational sites across the globe. The algorithm is also applicable to other lidar and/or ceilometer network infrastructures in the framework of the Global Aerosol Watch (GAW) aerosol lidar observation network (GALION).

Simone Lolli↗

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

Data management and scientific integration within the Atmospheric Radiation Measurement Program

The Atmospheric Radiation Measurement (ARM) Program has been developed by the U.S. Department of Energy with the goal to improve the predictive capabilities of General Circulation Models (GCM's) in their treatment of clouds and radiative transfer effects. To achieve this goal, three experimental testbeds were designed for the deployment of instruments that will collect atmospheric data used to drive the GCM's. Each site, known as a Cloud and Radiation Testbed (CART), consists of a highly available, redundant data system for the collection of data from a variety of instrumentation. The first CART site was deployed in April 1992 in the Southern Great Plains (SGP), Lamont, Oklahoma, with the other two sites to follow in September 1995 in the Tropical Western Pacific and in 1997 on the North Slope of Alaska. Approximately 400 MB of data are transferred per day via the Internet from the SGP site to the ARM Experiment Center at Pacific Northwest Laboratory in Richland, Washington. The Experiment Center is central to the ARM data path and provides for the collection, processing, analysis, and delivery of ARM data. Data are received from the CART sites from a variety of instrumentation, observational systems, amd external data sources. The Experiment Center processes these data streams on a continuous basis to provide derived data products to the ARM Science Team in near real-time while providing a three-month running archive of data. A primary requirement of the ARM Program is to preserve and protect all data produced or acquired. This function is performed at Oak Ridge National Laboratory where leading edge technology is employed for the long-term storage of ARM data. The ARM Archive provides access to data for participation outside of the ARM Program. The ARM Program involves a collaborative effort by teams from various DOE National Laboratories, providing multi-disciplinary areas of expertise. This paper will discuss the collaborative methods in which the ARM teams translate the scientific goals of the Program into data products. By combining atmospheric scientists, systems engineers, and software engineers, the ARM Program has successfully designed and developed an environment where advances in understanding the parameterizations of GCM's can be made.

Gracio, Deborah K.↗

The Impact of Microphysics and Planetary Boundary Layer Physics on Model Simulation of U.S. Deep South Summer Convection

Inspection of output from various configurations of high-resolution, explicit convection forecast models such as the Weather Research and Forecasting (WRF) model indicates significant sensitivity to the choices of model physics pararneterizations employed. Some of the largest apparent sensitivities are related to the specifications of the cloud microphysics and planetary boundary layer physics packages. In addition, these sensitivities appear to be especially pronounced for the weakly-sheared, multicell modes of deep convection characteristic of the Deep South of the United States during the boreal summer. Possible ocean-land sensitivities also argue for further examination of the impacts of using unique ocean-land surface initialization datasets provided by the NASA Short-term Prediction Research and Transition (SPoRn Center to select NOAAlNWS weather forecast offices. To obtain better quantitative understanding of these sensitivities and also to determine the utility of the ocean-land initialization data, we have executed matrices of regional WRF forecasts for selected convective events near Mobile, AL (MOB), and Houston, TX (HGX). The matrices consist of identically initialized WRF 24-h forecasts using any of eight microphysics choices and any of three planetary boWldary layer choices. The resulting 24 simulations performed for each event within either the MOB or HGX regions are then compared to identify the sensitivities of various convective storm metrics to the physics choices. Particular emphasis is placed on sensitivities of precipitation timing, intensity, and coverage, as well as amount and coverage oflightuing activity diagnosed from storm kinematics and graupel in the mixed phase layer. The results confirm impressions gleaned from study of the behavior of variously configured WRF runs contained in the ensembles produced each spring at the Center for the Analysis and Prediction of Storms, but with the benefit of more straightforward control of the physics package choices. The design of the experiments thus allows for more direct interpretation of the sensitivities to each possible physics combination. The results should assist forecasters in their efforts to anticipate and correct for possible biases in simulated WRF convection patterns, and help the modeling community refine their model parameterizations.

McCaul, Eugene W., Jr.↗

Inter-annual and Long-term Temperature Variations in the Mesopause Region at High Latitudes Generated by the Stratospheric QBO

The Numerical Spectral Model (NSM) simulates the Quasi-biennial Oscillation (QBO) that dominates the zonal circulation of the lower stratosphere at low latitudes. In the model, the QBO is generated with parameterized small-scale gravity waves (GW), which are partially augmented in 3D with planetary waves owing to baroclinic instability. Due to GW filtering, the QBO extends into the upper mesosphere, evident in UARS zonal wind and TIMED temperature measurements. While the QBO zonal winds are confined to equatorial latitudes, even in simulations with latitude-independent wave source, the associated temperature variations extend to high latitudes. The meridional circulation redistributes some of the QBO energy to focus it partially onto the Polar Regions. The resulting QBO temperature variations away from the equator tend to increase at higher altitudes to produce inter-annual variations that can exceed 5 K in the polar mesopause region -- and our 3D model simulations show that the effect is variable from year to year and can produce large differences between the two hemispheres, presumably due to interactions involving the seasonal variations. Modeling studies with the NSM have shown that long-term variations can also be generated by the QBO interacting with the seasonal cycles through OW node-filtering. A 30-month QBO, optimally synchronized by the 6-month Semi-Annual Oscillation (SAO), thus produces a 5-year or semi-decadal (SD) oscillation -- and observational evidence for that has been provided by a recent analysis of stratospheric NCEP data. In a simulation with the 2D version of the NSM, this SD oscillation extends into the upper mesosphere, and we present results to show that the related temperature variations could contribute significantly to the long-term variations of the polar mesopause region. Quasi-decadal variations could furthermore arise from the modeled solar cycle modulations of the QBO and 12-month annual oscillation. Our numerical results are discussed in the context of the observed low summer temperatures reproduced by the model, to demonstrate that the above interannual and long-term variations could contribute significantly to the climatology of Polar Mesospheric Clouds (PMC) investigated by the Aeronomy of Ice in the Mesosphere (AIM) mission.

Mayr, Hans G.↗