Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatially distributed model evaluation”

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 19 records

Application of spatial Poisson process models to air mass thunderstorm rainfall

Eight years of summer storm rainfall observations from 93 stations in and around the 154 sq km Walnut Gulch catchment of the Agricultural Research Service, U.S. Department of Agriculture, in Arizona are processed to yield the total station depths of 428 storms. Statistical analysis of these random fields yields the first two moments, the spatial correlation and variance functions, and the spatial distribution of total rainfall for each storm. The absolute and relative worth of three Poisson models are evaluated by comparing their prediction of the spatial distribution of storm rainfall with observations from the second half of the sample. The effect of interstorm parameter variation is examined.

Eagleson, P. S.↗

Cometary atmospheres: Modeling the spatial distribution of observed neutral radicals

Progress on modeling the spatial distributions of cometary radicals is described. The Monte Carlo particle-trajectory model was generalized to include the full time dependencies of initial comet expansion velocities, nucleus vaporization rates, photochemical lifetimes and photon emission rates which enter the problem through the comet's changing heliocentric distance and velocity. The effect of multiple collisions in the transition zone from collisional coupling to true free flow were also included. Currently available observations of the spatial distributions of the neutral radicals, as well as the latest available photochemical data were re-evaluated. Preliminary exploratory model results testing the effects of various processes on observable spatial distributions are also discussed.

Combi, M. R.↗

Evaluating Global Emission Inventories of Biogenic Bromocarbons

Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidizing capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF (National Science Foundation) HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38 %) to 0.78 (115 %) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24 %) to 1.25 (167 %) ppt. We also use aircraft observations made during the 2011 Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) campaign, in the tropical western Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (Br(VSLS/y)). Our simulations show Br(VSLS/y) ranges from approximately 4.0 to 8.0 ppt depending on the inventory. We report an optimized estimate at the lower end of this range (approximately 4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.

bromocarbons↗

Assessment of global and regional tectonic models for faulting in the ancient terrains of Mars

The structural mapping of tectonic landforms on Mars has provided a great deal of information regarding the planet's tectonic history. The present paper has the objective to test models of regional and global fracturing by evaluating the spatial distribution and orientation of the ancient heavily cratered terrain (HCT) structural landforms. Attention is given to background information concerning earlier studies, landform selection and analysis, rose diagrams, tests for Tharsis generation, and a calculation of thermal stress in the early Martian lithosphere. The contribution of classical Tharsis-centered tectonism to the production of HCT structural landforms was found to be very small. Few structural landforms investigated are either perpendicular to calculated principal stress trajectories or radial to Tharsis.

Schultz, R. A.↗

Spatial quantification of maps or images - Cell size or pixel size implications

The paper discusses spatial quantification of maps or images. Spatial accuracies evaluated by mapping and inventory measures were observed with a range of cell sizes for a 10.36 sq km soil survey segment. The distribution of interboundary distances in a map was evaluated as unique, characterizing information; experimentally observed relationships between accuracy and cell size and interboundary distance distribution parameters and cell size led to a model of the process of spatial quantization. Grid positioning effects were evaluated and were significant only at the level of individual mapping units; two generations of 'universal' process models were derived, implemented, and evaluated.

Wehde, M. E.↗

Feeding Ten Billion People Is Possible Within Four Terrestrial Planetary Boundaries

Global agriculture puts heavy pressure on planetary boundaries, posing the challenge to achieve future food security without compromising Earth system resilience. On the basis of process-detailed, spatially explicit representation of four interlinked planetary boundaries (biosphere integrity, land-system change, freshwater use, nitrogen flows) and agricultural systems in an internally consistent model framework, we here show that almost half of current global food production depends on planetary boundary transgressions. Hotspot regions, mainly in Asia, even face simultaneous transgression of multiple underlying local boundaries. If these boundaries were strictly respected, the present food system could provide a balanced diet (2,355 kcal per capita per day) for 3.4 billion people only. However, as we also demonstrate, transformation towards more sustainable production and consumption patterns could support 10.2 billion people within the planetary boundaries analysed. Key prerequisites are spatially redistributed cropland, improved water–nutrient management, food waste reduction and dietary changes. Adoption of the Sustainable Development Goals by all nations in 2015 is the first ever commitment to a world development path that safeguards the stability of the Earth system as a prerequisite for meeting universal human standards1. The longstanding challenge of achieving food security through sustainable agriculture is particularly acute in this context as world agriculture is a leading cause for the current transgressions of multiple planetary boundaries (PBs) globally and regionally2–5. The PB framework is a comprehensive scientific attempt to synoptically define our planet’s biogeophysical limits to anthropogenic interference. It suggests bounds to nine interacting processes that together delineate a Holocene-like Earth system state. The Holocene is chosen as the reference state as it is the only period known to provide a safe operating space for a world population of several billion people, and according to a precautionary principle, the PBs are set in sufficient distance from processes that may critically undermine Earth system resilience and global sustainability. A challenging question, thus, is whether human development goals such as food security can be met while maintaining multiple PBs along with their subglobal manifestations. Further PB transgressions could jeopardize the chances of providing sufficient food for a world population projected to be wealthier and reach >9 billion by 2050. This conundrum portrays a tradeoff between Earth’s biophysical carrying capacity and humankind’s rising food demand, calling in response for radical rethinking of food production and consumption patterns6–9. Yield gap closures, avoidance of excessive input use, shifts towards less resource-demanding diets, food waste reductions and efficient international trade are crucial options for sustainably increasing the food supply10–15. For example, enhancing water-use efficiency on irrigated and rain-fed farms can triple or quadruple crop yields in low-performing systems, suggesting possible global gains of >20% (ref. 16). Even higher gains appear feasible through globally optimized configurations of the land-use pattern17, and cutting food losses by half could generate food for another billion people18. Thus, collective large-scale implementation of such options could sustain food for a further growing world population19. Yet achieving this within a safe operating space as defined by PBs requires not only a halt to but actually a reversal of existing PB transgressions. Previous studies suggest that such a reconciliation might be possible, but these were based on aggregate representations of PBs (not accounting for the spatial patterns of limits, transgressions and interactions) or considered only one boundary in isolation17,20–23. Here, we systematically quantify to what extent current food production depends on local to global transgressions of the PBs for biosphere integrity, land-system change, freshwater use and nitrogen (N) flows, along with the potential of a range of solutions to avoid these transgressions and still increase food supply (Table 1). To this end, we configured an internally consistent process-based model of the terrestrial biosphere including agriculture (LPJmL) with multiple spatially distributed PBs and their interactions. LPJmL is among the longest-established and best-evaluated biosphere models, showing robust performance regarding simulation of, for example, carbon, water and crop yield dynamics (Supplementary Figs. 1 and 2 and Supplementary Table 1; see ref. 24 for a comprehensive benchmarking and Supplementary Methods for more detail on model evaluations). In principle following established definitions4, we refine the computation of some PBs with respect to their regional patterns and interactions (Methods), providing globally gridded precautionary limits to human interference with the Earth system at a level of great detail. In particular, we account for the evidence that many PBs need to be represented spatially explicitly4 to cover their

Gerten, Dieter↗

The Distribution of Snow Black Carbon observed in the Arctic and Compared to the GISS-PUCCINI Model

In this study, we evaluate the ability of the latest NASA GISS composition-climate model, GISS-E2- PUCCINI, to simulate the spatial distribution of snow BC (sBC) in the Arctic relative to present-day observations. Radiative forcing due to BC deposition onto Arctic snow and sea ice is also estimated. Two sets of model simulations are analyzed, where meteorology is linearly relaxed towards National Centers for Environmental Prediction (NCEP) and towards NASA Modern Era Reanalysis for Research and Applications (MERRA) reanalyses. Results indicate that the modeled concentrations of sBC are comparable with presentday observations in and around the Arctic Ocean, except for apparent underestimation at a few sites in the Russian Arctic. That said, the model has some biases in its simulated spatial distribution of BC deposition to the Arctic. The simulations from the two model runs are roughly equal, indicating that discrepancies between model and observations come from other sources. Underestimation of biomass burning emissions in Northern Eurasia may be the main cause of the low biases in the Russian Arctic. Comparisons of modeled aerosol BC (aBC) with long-term surface observations at Barrow, Alert, Zeppelin and Nord stations show significant underestimation in winter and spring concentrations in the Arctic (most significant in Alaska), although the simulated seasonality of aBC has been greatly improved relative to earlier model versions. This is consistent with simulated biases in vertical profiles of aBC, with underestimation in the lower and middle troposphere but overestimation in the upper troposphere and lower stratosphere, suggesting that the wet removal processes in the current model may be too weak or that vertical transport is too rapid, although the simulated BC lifetime seems reasonable. The combination of observations and modeling provides a comprehensive distribution of sBC over the Arctic. On the basis of this distribution, we estimate the decrease in snow and sea ice albedo and the resulting radiative forcing. We suggest that the albedo reduction due to BC deposition presents significant space-time variations, with highest mean reductions of 1.25% in the Russian Arctic, which are much larger than those in other Arctic regions (0.39% to 0.64 %). The averaged value over the Arctic north of 66degN is 0.4-0.6% during spring, leading to regional surface radiative forcings of 0.7, 1.1 and 1.0Wm(exp−2) in spring 2007, 2008 and 2009, respectively.

Dou, T.↗

Impacts of aerosol size-dependent below-cloud scavenging on tropospheric aerosol in the NASA GEOS model

Cloud scavenging is the major sink for a suite of tropospheric aerosols, and thus largely affects the global distribution and lifetime of aerosols as well as the formation of clouds. Simplified model parameterizations of cloud scavenging of aerosols substantially contribute to large uncertainties in the simulated aerosol burdens, spatial distributions, lifetimes, as well as aerosol radiative forcing. Here we evaluate the wet scavenging scheme in the current NASA GEOS model, develop more physically-based parameterizations of below-cloud scavenging, and assess the impacts of the new scavenging parameterizations on the GEOS simulations of aerosols. The current below-cloud scavenging (BCS; washout) scheme uses a first-order removal with a constant scavenging coefficient. We replaced the current scheme with a size-dependent parameterization of Croft et al.(Atmos. Chem. Phys., 2009; hereafter referred to as C09). The C09 BCS parameterizations calculate aerosol below-cloud scavenging coefficients based on a size-dependent collision efficiency between aerosol particles and raindrops (or snow crystals). A look-up table of compiled collision efficiencies provides values for sixty aerosol sizes and nine rain flux rates. Snow-aerosol collision efficiencies are a function of the sixty aerosol sizes only. We assume a log-normal size distribution for aerosols (bulk) in GOCART. The collision efficiencies are then calculated online using a bilinear interpolation from the look-up table for the given aerosol sizes and precipitation flux rates in the model.

aerosols↗

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗

Spectrally Resolved Flux Derived from Collocated AIRS and CERES Observations and its Application in Model Validation: Clear-Sky Over the Tropic Oceans - Part I

Spectrally resolved outgoing IR flux, the integrand of the outgoing longwave radiation (OLR), has its unique value in evaluating model simulations. Here we describe an algorithm of deriving such clear-sky outgoing spectral flux through the whole IR region from the collocated Atmospheric Infrared Sounder (AIRS) and the Clouds & the Earth's Radiant Energy System (CERES) measurements over the tropical oceans. Based on the scene types and corresponding angular distribution models (ADMs) used in the CERES Single Satellite Footprint (SSF) dataset, spectrally-dependent ADMs are developed and used to estimate the spectral flux at each AIRS channel. A multivariate linear prediction scheme is then used to estimate spectral fluxes at frequencies not covered by the AIRS instrument. The whole algorithm is validated using synthetic spectra as well as the CERES OLR measurements. Using the GFDL AM2 model simulation as a case study, the application of the derived clear-sky outgoing spectral flux in model evaluation is illustrated. By comparing the observed and simulated spectral flux in 2004, compensating errors in the simulated OLR from different absorption bands can be revealed, so does the errors from frequencies within a given absorption band. Discrepancies between the simulated and observed spatial distributions and seasonal evolutions of the spectral fluxes at different spectral ranges are further discussed. The methodology described in this study can be applied to other surface types as well as cloudy-sky observations and corresponding model evaluations.

Huang, Xianglei↗

A simulator for evaluating methods for the detection of lesion-deficit associations

Although much has been learned about the functional organization of the human brain through lesion-deficit analysis, the variety of statistical and image-processing methods developed for this purpose precludes a closed-form analysis of the statistical power of these systems. Therefore, we developed a lesion-deficit simulator (LDS), which generates artificial subjects, each of which consists of a set of functional deficits, and a brain image with lesions; the deficits and lesions conform to predefined distributions. We used probability distributions to model the number, sizes, and spatial distribution of lesions, to model the structure-function associations, and to model registration error. We used the LDS to evaluate, as examples, the effects of the complexities and strengths of lesion-deficit associations, and of registration error, on the power of lesion-deficit analysis. We measured the numbers of recovered associations from these simulated data, as a function of the number of subjects analyzed, the strengths and number of associations in the statistical model, the number of structures associated with a particular function, and the prior probabilities of structures being abnormal. The number of subjects required to recover the simulated lesion-deficit associations was found to have an inverse relationship to the strength of associations, and to the smallest probability in the structure-function model. The number of structures associated with a particular function (i.e., the complexity of associations) had a much greater effect on the performance of the analysis method than did the total number of associations. We also found that registration error of 5 mm or less reduces the number of associations discovered by approximately 13% compared to perfect registration. The LDS provides a flexible framework for evaluating many aspects of lesion-deficit analysis.

NASA Program Biomedical Research and Countermeasur↗

Ridge belts - Evidence for regional- and local-scale deformation on the surface of Venus

Models for the formation of the ridge belt fan assemblage on Venus are evaluated through consideration of the orientation, spatial distribution, topographic expression, and wavelengths of observed tectonic surface features. A compressional mechanism is favored for long wavelength deformation corresponding to the spacing of ridge belts (300 - 400 km). However, short wavelength ridges and grooves (10 - 20 km) that are contained within belts and trend parallel to them are likely to be compressional and extensional, and to reflect both regional and local stress fields. It is hypothesized that large-scale (much greater than ridge spacing) early-stage mantle downwelling is the source of regional compression responsible for the establishment of the long wavelength of deformation.

Zuber, M. T.↗

Evaluating Secondary Inorganic Aerosols in Three Dimensions

The spatial distribution of aerosols and their chemical composition dictates whether aerosols have a cooling or a warming effect on the climate system. Hence, properly modeling the three-dimensional distribution of aerosols is a crucial step for coherent climate simulations. Since surface measurement networks only give 2-D data, and most satellites supply integrated column information, it is thus important to integrate aircraft measurements in climate model evaluations. In this study, the vertical distribution of secondary inorganic aerosol (i.e., sulfate, ammonium, and nitrate) is evaluated against a collection of 14 AMS flight campaigns and surface measurements from 2000 to 2010 in the USA and Europe. GISS ModelE2 is used with multiple aerosol microphysics (MATRIX, OMA) and thermodynamic (ISORROPIA II, EQSAM) configurations. Our results show that the MATRIX microphysical scheme improves the model performance for sulfate, but that there is a systematic underestimation of ammonium and nitrate over the USA and Europe in all model configurations. In terms of gaseous precursors, nitric acid concentrations are largely underestimated at the surface while overestimated in the higher levels of the model. Heterogeneous reactions on dust surfaces are an important sink for nitric acid, even high in the troposphere. At high altitudes, nitrate formation is calculated to be ammonia limited. The underestimation of ammonium and nitrate in polluted regions is most likely caused by a too simplified treatment of the NH3/NH4(+) partitioning which affects the HNO3/NO3(-) partitioning.

aerosols↗

Modeling the smoky troposphere of the southeast Atlantic: a comparison to ORACLES airborne observations from September of 2016

In the southeast Atlantic, well-defined smoke plumes from Africa advect over marine boundary layer cloud decks; both are most extensive around September, when most of the smoke resides in the free troposphere. A framework is put forth for evaluating the performance of a range of global and regional atmospheric composition models against observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) airborne mission in September 2016. A strength of the comparison is a focus on the spatial distribution of a wider range of aerosol composition and optical properties than has been done previously. The sparse airborne observations are aggregated into approximately 2° grid boxes and into three vertical layers: 3–6 km, the layer from cloud top to 3 km, and the cloud-topped marine boundary layer. Simulated aerosol extensive properties suggest that the flight-day observations are reasonably representative of the regional monthly average, with systematic deviations of 30 % or less. Evaluation against observations indicates that all models have strengths and weaknesses, and there is no single model that is superior to all the others in all metrics evaluated. Whereas all six models typically place the top of the smoke layer within 0–500 m of the airborne lidar observations, the models tend to place the smoke layer bottom 300–1400 m lower than the observations. A spatial pattern emerges, in which most models underestimate the mean of most smoke quantities (black carbon, extinction, carbon monoxide) on the diagonal corridor between 16° S, 6° E, and 10° S, 0° E, in the 3–6 km layer, and overestimate them further south, closer to the coast, where less aerosol is present. Model representations of the above-cloud aerosol optical depth differ more widely. Most models overestimate the organic aerosol mass concentrations relative to those of black carbon, and with less skill, indicating model uncertainties in secondary organic aerosol processes. Regional-mean free-tropospheric model ambient single scattering albedos vary widely, between 0.83 and 0.93 compared with in situ dry measurements centered at 0.86, despite minimal impact of humidification on particulate scattering. The modeled ratios of the particulate extinction to the sum of the black carbon and organic aerosol mass concentrations (a mass extinction efficiency proxy) are typically too low and vary too little spatially, with significant inter-model differences. Most models overestimate the carbonaceous mass within the offshore boundary layer. Overall, the diversity in the model biases suggests that different model processes are responsible. The wide range of model optical properties requires further scrutiny because of their importance for radiative effect estimates.

Yohei Shinozuka↗

Scientific Objectives, Measurement Needs, and Challenges Motivating the PARAGON Aerosol Initiative

Aerosols are involved in a complex set of processes that operate across many spatial and temporal scales. Understanding these processes, and ensuring their accurate representation in models of transport, radiation transfer, and climate, requires knowledge of aerosol physical, chemical, and optical properties and the distributions of these properties in space and time. To derive aerosol climate forcing, aerosol optical and microphysical properties and their spatial and temporal distributions, and aerosol interactions with clouds, need to be understood. Such data are also required in conjunction with size-resolved chemical composition in order to evaluate chemical transport models and to distinguish natural and anthropogenic forcing. Other basic parameters needed for modeling the radiative influences of aerosols are surface reflectivity and three-dimensional cloud fields. This large suite of parameters mandates an integrated observing and modeling system of commensurate scope. The Progressive Aerosol Retrieval and Assimilation Global Observing Network (PARAGON) concept, designed to meet this requirement, is motivated by the need to understand climate system sensitivity to changes in atmospheric constituents, to reduce climate model uncertainties, and to analyze diverse collections of data pertaining to aerosols. This paper highlights several challenges resulting from the complexity of the problem. Approaches for dealing with them are offered in the set of companion papers.

optical depth↗

The Analysis of the Patterns of Radiation-Induced DNA Damage Foci by a Stochastic Monte Carlo Model of DNA Double Strand Breaks Induction by Heavy Ions and Image Segmentation Software

To create a generalized mechanistic model of DNA damage in human cells that will generate analytical and image data corresponding to experimentally observed DNA damage foci and will help to improve the experimental foci yields by simulating spatial foci patterns and resolving problems with quantitative image analysis. Material and Methods: The analysis of patterns of RIFs (radiation-induced foci) produced by low- and high-LET (linear energy transfer) radiation was conducted by using a Monte Carlo model that combines the heavy ion track structure with characteristics of the human genome on the level of chromosomes. The foci patterns were also simulated in the maximum projection plane for flat nuclei. Some data analysis was done with the help of image segmentation software that identifies individual classes of RIFs and colocolized RIFs, which is of importance to some experimental assays that assign DNA damage a dual phosphorescent signal. Results: The model predicts the spatial and genomic distributions of DNA DSBs (double strand breaks) and associated RIFs in a human cell nucleus for a particular dose of either low- or high-LET radiation. We used the model to do analyses for different irradiation scenarios. In the beam-parallel-to-the-disk-of-a-flattened-nucleus scenario we found that the foci appeared to be merged due to their high density, while, in the perpendicular-beam scenario, the foci appeared as one bright spot per hit. The statistics and spatial distribution of regions of densely arranged foci, termed DNA foci chains, were predicted numerically using this model. Another analysis was done to evaluate the number of ion hits per nucleus, which were visible from streaks of closely located foci. In another analysis, our image segmentaiton software determined foci yields directly from images with single-class or colocolized foci. Conclusions: We showed that DSB clustering needs to be taken into account to determine the true DNA damage foci yield, which helps to determine the DSB yield. Using the model analysis, a researcher can refine the DSB yield per nucleus per particle. We showed that purely geometric artifacts, present in the experimental images, can be analytically resolved with the model, and that the quantization of track hits and DSB yields can be provided to the experimentalists who use enumeration of radiation-induced foci in immunofluorescence experiments using proteins that detect DNA damage. An automated image segmentaiton software can prove useful in a faster and more precise object counting for colocolized foci images.

Ponomarev, Artem↗

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗