SEARCH · Engineering Papers
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.
In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.
We demonstrate the practicality and effectiveness of using a Green's functions estimation approach for adjusting uncertain parameters in an Earth system model (ESM). This estimation approach has previously been applied to an intermediate-complexity climate model and to individual ESM components, e.g., ocean, sea ice, or carbon cycle components. Here, the Green's functions approach is applied to a state-of-the-art ESM that comprises a global atmosphere/land configuration of the Goddard Earth Observing System (GEOS) coupled to an ocean and sea ice configuration of the Massachusetts Institute of Technology general circulation model (MITgcm). Horizontal grid spacing is approximately 110 km for GEOS and 37–110 km for MITgcm. In addition to the reference GEOS-MITgcm simulation, we carried out a series of model sensitivity experiments, in which 20 uncertain parameters are perturbed. These “control” parameters can be used to adjust sea ice, microphysics, turbulence, radiation, and surface schemes in the coupled simulation. We defined eight observational targets: sea ice fraction, net surface shortwave radiation, downward longwave radiation, near-surface temperature, sea surface temperature, sea surface salinity, and ocean temperature and salinity at 300 m. We applied the Green's functions approach to optimize the values of the 20 control parameters so as to minimize a weighted least-squares distance between the model and the eight observational targets. The new experiment with the optimized parameters resulted in a total cost reduction of 9 % relative to a simulation that had already been adjusted using other methods. The optimized experiment attained a balanced cost reduction over most of the observational targets. We also report on results from a set of sensitivity experiments that are not used in the final optimized simulation but helped explore options and guided the optimization process. These experiments include an assessment of sensitivity to the number of control parameters and to the selection of observational targets and weights in the cost function. Based on these sensitivity experiments, we selected a specific definition for the cost function. The sensitivity experiments also revealed a decreasing overall cost as the number of control variables was increased. In summary, we recommend using the Green's functions estimation approach as an additional fine-tuning step in the model development process. The method is not a replacement for modelers' experience in choosing and adjusting sensitive model parameters. Instead, it is an additional practical and effective tool for carrying out final adjustments of uncertain ESM parameters.
During the past decade, numerical weather and global non-hydrostatic models have started using more complex microphysical schemes originally developed for high resolution cloud resolving models (CRMs) with 1-2 km or less horizontal resolutions. These microphysical schemes affect the dynamic through the release of latent heat (buoyancy loading and pressure gradient) the radiation through the cloud coverage (vertical distribution of cloud species), and surface processes through rainfall (both amount and intensity). Recently, several major improvements of ice microphysical processes (or schemes) have been developed for cloud-resolving model (Goddard Cumulus Ensemble, GCE, model) and regional scale (Weather Research and Forecast, WRF) model. These improvements include an improved 3-ICE (cloud ice, snow and graupel) scheme (Lang et al. 2010); a 4-ICE (cloud ice, snow, graupel and hail) scheme and a spectral bin microphysics scheme and two different two-moment microphysics schemes. The performance of these schemes has been evaluated by using observational data from TRMM and other major field campaigns. In this talk, we will present the high-resolution (1 km) GeE and WRF model simulations and compared the simulated model results with observation from recent field campaigns [i.e., midlatitude continental spring season (MC3E; 2010), high latitude cold-season (C3VP, 2007; GCPEx, 2012), and tropical oceanic (TWP-ICE, 2006)].
The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.
The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central U.S. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined. This will be mainly achieved by performing sensitivity tests using various horizontal (from 1- to 5-kilometers) and vertical resolutions (from 20- to 200-meters in the lower troposphere to 200- to 500-m in the middle and upper troposphere). Sensitivity tests using various microphysical schemes (warm rain only, and three ice with either graupel or hail) will also be performed. The PBL and diurnal variation of precipitation processes will also be evaluated. The budgets will be calculated for different regions (i.e., convective and stratiform regions).
The Goddard Cumulus Ensemble (GCE) model examines the impact of various microphysical schemes, and vertical and horizontal resolution in the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing rations of all water phases (vapor, liquid and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes and with surface processes. For idealized hurricanes, an axisymmetric version of the GCE model was developed and used to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was use to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central U.S. FOr the cloud ensemble, the GCE model was integrated for several days in order to have a good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolutio will be examined. This will be mainly achieved by performing sensitivity tests using various horizontal (from 1-to 5-km) and vertical resolutions (from 20- to 200-m in the lower troposphere to 200- to 500m in the middle and upper troposphere). Sensitivity test using various microphysical schemes (warm rain only, and three ice with either graupel or hail) will also be performed. The thermodynamic and water budget associated with various types of precipitation systems will also be evaluated. The budgets will be calculated for different regions (i.e., convective and stratiform regions).
The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution ion the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble f clouds. The model variables include horizontal and vertical velocities, potential temperatures, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulated squall lines that developed in the western Pacific, eastern Atlantic and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined. This will be mainly achieved by performing sensitivity tests using various horizontal (from 1- to 5-km) and vertical resolutions (from 20- to 200-m in the lower troposphere to 200- to 500-m in the middle and upper troposphere). Sensitivity tests using various microphysical schemes (warm rain only, and three ice with either graupel or hail) will also be performed. The thermodynamic and water budget associated with various types of precipitation systems will also be evaluated. The budgets will be calculated for different regions (i.e., convective and stratiform regions).
Several different ice nucleation parameterizations in two different General Circulation Models (GCMs) are used to understand the effects of ice nucleation on the mean climate state, and the Aerosol Indirect Effects (AIE) of cirrus clouds on climate. Simulations have a range of ice microphysical states that are consistent with the spread of observations, but many simulations have higher present-day ice crystal number concentrations than in-situ observations. These different states result from different parameterizations of ice cloud nucleation processes, and feature different balances of homogeneous and heterogeneous nucleation. Black carbon aerosols have a small (0.06 Wm(exp-2) and not statistically significant AIE when included as ice nuclei, for nucleation efficiencies within the range of laboratory measurements. Indirect effects of anthropogenic aerosols on cirrus clouds occur as a consequence of increasing anthropogenic sulfur emissions with different mechanisms important in different models. In one model this is due to increases in homogeneous nucleation fraction, and in the other due to increases in heterogeneous nucleation with coated dust. The magnitude of the effect is the same however. The resulting ice AIE does not seem strongly dependent on the balance between homogeneous and heterogeneous ice nucleation. Regional effects can reach several Wm2. Indirect effects are slightly larger for those states with less homogeneous nucleation and lower ice number concentration in the base state. The total ice AIE is estimated at 0.27 +/- 0.10 Wm(exp-2) (1 sigma uncertainty). This represents a 20% offset of the simulated total shortwave AIE for ice and liquid clouds of 1.6 Wm(sup-2).
NASA's Marshall Space Flight Center and the University of Alabama in Huntsville (UAHuntsville) are collaborating with the 45th Weather Squadron (45WS) at Cape Canaveral Air Force Station (CCAFS) to enable improved nowcasting of lightning cessation. This project centers on use of dual-polarimetric radar capabilities, and in particular, the new C-band dual polarimetric weather radar acquired by the 45WS. Special emphasis is placed on the development of a physically-based operational algorithm to predict lightning cessation. While previous studies have developed statistically based lightning cessation algorithms driven primarily by trending in the actual total lightning flash rate, we believe that dual polarimetric radar variables offer the possibility to improve existing algorithms through the inclusion of physically meaningful trends reflecting interactions between in-cloud electric fields and ice-microphysics. Specifically, decades of polarimetric radar research using propagation differential phase has demonstrated the presence of distinct phase and ice crystal alignment signatures in the presence of strong electric fields associated with lightning. One question yet to be addressed is: To what extent can propagation phase-based ice-crystal alignment signatures be used to nowcast the cessation of lightning activity in a given storm? Accordingly, data from the UAHuntsville Advanced Radar for Meteorological and Operational Research (ARMOR) along with the NASA-MSFC North Alabama Lightning Mapping Array are used in this study to investigate the radar signatures present before and after lightning cessation. Thus far our case study results suggest that the negative differential phase shift signature weakens and disappears after the analyzed storms ceased lightning production (i.e., after the last lightning flash occurred). This is a key observation because it suggests that while strong electric fields may still have been present, the lightning cessation signature was encompassed in the period of the polarimetric negative phase shift signature. To the extent this behavior is repeatable in other cases, even if only in a substantial fraction of those cases, the analysis suggests that differential propagation phase may prove to be a useful parameter for future lightning cessation algorithms. Indeed, a preliminary analysis of 15+ cases has shown additional indications of the weakening and disappearance of this ice alignment signature with lightning cessation. A summary of these case-study results is presented.
Previous studies of tropical lightning activity using Tropical Rainfall Measurement Mission (TRMM) Lightning Imaging Sensor (LIS) data performed analyses of lightning behavior over mesoscale "feature" scales or over uniform grids. In order to study lightning and the governing ice microphysics intrinsic to thunderstorms at a more process-specific scale (i.e., the scale over which electrification processes and lightning occur in a "unit" thunderstorm), a new convective cell-scale database was developed by analyzing and refining the University of Utah's Precipitation Features database and retaining precipitation data parameters computed from the TRMM precipitation radar (PR), microwave imager (TMI) and LIS instruments. The resulting data base was to conduct a limited four-year study of tropical continental convection occurring over the Amazon Basin, Congo, Maritime Continent and the western Pacific Ocean. The analysis reveals expected strong correlations between lightning flash counts per cell and ice proxies, such as ice water path, minimum and average 85GHz brightness temperatures, and 18dBz echo top heights above the freezing level in all regimes, as well as regime-specific relationships between lighting flash counts and PR-derived surface rainfall rates. Additionally, radar CFADs were used to partition the 3D structure of cells in each regime at different flash counts. The resulting cell-scale analyses are compared to previous mesoscale feature and gridded studies wherever possible.
Recently launched cloud observing satellites provide information about the vertical structure of deep convection and its microphysical characteristics. In this study, CloudSat reflectivity data is stratified by cloud type, and the contoured frequency by altitude diagrams reveal a doublearc structure in deep convective cores (DCCs) above 8 km. This suggests two distinct hydrometeor modes (snow versus hail/graupel) controlling variability in reflectivity profiles. The day–night contrast in the double arcs is about four times larger than the wet–dry season contrast. Using QuickBeam, the vertical reflectivity structure of DCCs is analyzed in two versions of the Superparameterized Community Atmospheric Model (SP-CAM) with single-moment (no graupel) and double-moment (with graupel) microphysics. Doublemoment microphysics shows better agreement with observed reflectivity profiles; however, neither model variant captures the double-arc structure. Ultimately, the results show that simulating realistic DCC vertical structure and its variability requires accurate representation of ice microphysics, in particular the hail/graupel modes, though this alone is insufficient.
During this grant, covering the period from September 1998 to December 2001, we continued the investigation of the role of turbulent mixing in the wake of subsonic aircraft initiated in 1994 for NASA's Atmospheric Effects of Aviation Project. The goal of the research has been to provide sufficient understanding and quantitative analytical capability to assess the dynamical, chemical, and microphysical interactions in the near-field wake that have the greatest potential to influence the global atmospheric impact of the projected fleet of subsonic aircraft. Through large-eddy simulations we have shown that turbulence in the early wake dynamics can have a strong effect on both the ice microphysics of contrail evolution and on wake chemistry. The wake vortex dynamics are the primary determinant of the vertical extent of the contrail; this together with the local wind shear largely determines the horizontal extent. The fraction of the initial ice crystals surviving the wake vortex dynamics, their spatial distribution, and the ice mass distribution are all sensitive to the aircraft type, assumed initial ice crystal number, and ambient humidity and turbulence conditions. Our model indicates that there is a significant range of conditions for which a smaller aircraft such as a B737 produces as significant a persistent contrail as a larger aircraft such as a B747, even though the latter consumes almost five times as much fuel. Large-eddy simulations of the near wake of a B757 provided a fine-grained chemical-dynamical representation of simplified NOx - HOx chemistry in wakes of ages from a few seconds to several minutes. By sampling the simulated data in a manner similar to that of in situ aircraft measurements it was possible to provide a likely explanation for a puzzle uncovered in the 1996 SUCCESS flight measurements of OH and HO2 The results illustrate the importance of considering fluid dynamics effects in interpreting chemistry results when mixing rates and species fluctuations are large, and demonstrate the feasibility of using 3D unsteady LES with coupled chemistry to study such phenomena.
The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have a good sampling of cloud statistics.
The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, eastern Atlantic and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined.
Profiles of the microphysical properties of clouds and raincells are essential in many areas of atmospheric research and operational meteorology. In order to enhance the understanding of the nonlinear and underconstrained relationships between cloud and hydrometeor microphysical profiles and passive microwave brightness temperatures, estimations of cloud profiles for an anvil, a convective, and an updraft region of an oceanic squall were performed. The estimations relied on comparisons between radiative transfer calculations of incrementally estimated microphysical profiles and concurrent dual-altitude wideband brightness temperatures from the 22 February 1993 flight during TOGA-COARE. The wideband observations (10--220 GHz) are necessary for estimating cloud profiles reaching up to 20 km. The low frequencies enhance the rain and cloud water profiles, while the high frequencies are required to detail the higher altitude ice microphysics. A microphysical profile was estimated for each of the three regions of the storm. Each of the three estimated profiles produced calculated brightness temperatures within approximately 10 K of the observations. A majority, of the total iterative adjustment were to the estimated profile's frozen hydrometeor characteristics and were necessary to match the high frequency calculations with the observations. This indicates a need to validate cloud resolving models using high frequencies. Some difficulties matching the 37 GHz observation channels on the DC-8 and ER-2 aircrafts with the calculations simulated at the two aircraft heights (approximately 11 km and 20 km, respectively) were noted and potential causes presented.
The variability of cirrus ice microphysical properties is investigated using observations obtained during the Small Particles in Cirrus (SPARTICUS) campaign. An existing approach that represents a size distribution (SD) as a single gamma function using an ellipsoid of equally realizable solutions in (N(sub 0), lambda, mu) phase space is modified to automatically identify multiple modes in SDs and characterize each mode by such an ellipsoid. The modified approach is applied to ice crystals with maximum dimension D greater than15 micrometers collected by the 2-D stereo and 2-D precipitation probes on the Stratton Park Engineering Company Learjet. The dependencies of N(sub 0), mu, and lambda from each mode, total number concentration, bulk extinction, ice water content (IWC), and mass median maximum dimension D(sub mm) as a function of temperature T and cirrus type are then analyzed. The changes in the observed codependencies between N(sub 0), mu, and lambda, bulk extinction, IWC, and D(sub mm) with environmental conditions indicate that particles were larger at higher T during SPARTICUS. At most two modes were observed in any SD during SPARTICUS, with the average boundary between them at 115 micrometers, similar to past studies not using probes with shatter mitigating tips and artifact removal algorithms. The bimodality of the SDs increased with T. This and the differences in N(sub 0), mu, and lambda between the modes suggest that particles with smaller D nucleated more recently than particles with larger D, which grew via vapor deposition and aggregation. Because smaller crystals, whose concentrations are uncertain, make marginal contributions to higher order moments, the use of higher moments for evaluating model fields is suggested.
One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and in-coming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a CRM, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. The GCE model has been used to understand the following: 1) water and energy cycles and their roles in the tropical climate system; 2) the vertical redistribution of ozone and trace constituents by individual clouds and well organized convective systems over various spatial scales; 3) the relationship between the vertical distribution of latent heating (phase change of water) and the large-scale (pre-storm) environment; 4) the validity of assumptions used in the representation of cloud processes in climate and global circulation models; and 5) the representation of cloud microphysical processes and their interaction with radiative forcing over tropical and midlatitude regions. Four-dimensional cloud and latent heating fields simulated from the GCE model have been provided to the TRMM Science Data and Information System (TSDIS) to develop and improve algorithms for retrieving rainfall and latent heating rates for TRMM and the NASA Earth Observing System (EOS). More than 90 referred papers using the GCE model have been published in the last two decades. Also, more than 10 national and international universities are currently using the GCE model for research and teaching. In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.
We advocate for several modeling approaches suitable to prove understanding of and capability to model aerosol behavior in nucleating droplets and ice and to model regional cloud systems at the mesoscale. The full complement of data that is ideal for particular study types is increasingly well understood; the ASR ISDAC Science Plan provides an excellent example of listing the data required for each study type targeted by that campaign. Whereas discussion here focused on the microscale and mesoscale, other studies should obviously focus on the regional, seasonal and global scale. When considering ice microphysics alone, laser focus on unconstrained model components could be usefully intensified (e.g., collision-coalescence kernels). Similarly, a focus on boundary layer structure would be deeply valuable to a wide range of model types. The expense of observations motivates efficient planning for observational data sets to support specific study approaches; the concept of closure is often useful (e.g., the surface energy budget).