Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit modeling”

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

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

190 records · Page 11

Confronting Models with Data: The GEWEX Cloud Systems Study

The GEWEX Cloud System Study (GCSS; GEWEX is the Global Energy and Water Cycle Experiment) was organized to promote development of improved parameterizations of cloud systems for use in climate and numerical weather prediction models, with an emphasis on the climate applications. The strategy of GCSS is to use two distinct kinds of models to analyze and understand observations of the behavior of several different types of clouds systems. Cloud-system-resolving models (CSRMs) have high enough spatial and temporal resolutions to represent individual cloud elements, but cover a wide enough range of space and time scales to permit statistical analysis of simulated cloud systems. Results from CSRMs are compared with detailed observations, representing specific cases based on field experiments, and also with statistical composites obtained from satellite and meteorological analyses. Single-column models (SCMs) are the surgically extracted column physics of atmospheric general circulation models. SCMs are used to test cloud parameterizations in an un-coupled mode, by comparison with field data and statistical composites. In the original GCSS strategy, data is collected in various field programs and provided to the CSRM Community, which uses the data to "certify" the CSRMs as reliable tools for the simulation of particular cloud regimes, and then uses the CSRMs to develop parameterizations, which are provided to the GCM Community. We report here the results of a re-thinking of the scientific strategy of GCSS, which takes into account the practical issues that arise in confronting models with data. The main elements of the proposed new strategy are a more active role for the large-scale modeling community, and an explicit recognition of the importance of data integration.

Randall, David↗

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗

A Meteorological Model's Dependence on Radiation Update Frequency

Numerical weather models are used to simulate circulations in the atmosphere including clouds and precipitation by applying a set of mathematical equations over a three-dimensional grid. The grid is composed of discrete points at which the meteorological variables are defined. As computing power continues to rise these models are being used at finer grid spacing, but they must still cover a wide range of scales. Some of the physics that must be accounted for in the model cannot be explicitly resolved, and their effects, therefore, must be estimated or "parameterized". Some of these parameterizations are computationally expensive. To alleviate the problem, they are not always updated at the time resolution of the model with the assumption being that the impact will be small. In this study, a coupled land-atmosphere model is used to assess the impact of less frequent updates of the computationally expensive radiation physics for a case on June 6, 2002, that occurred during a field experiment over the central plains known as International H20 Project (IHOP). The model was tested using both the original conditions, which were dry, and with modified conditions wherein moisture was added to the lower part of the atmosphere to produce clouds and precipitation (i.e., a wet case). For each of the conditions (i.e., dry and wet), four set of experiments were conducted wherein the model was run for a period of 24 hours and the radiation fields (including both incoming solar and outgoing longwave) were updated every 1, 3, 10, and 100 time steps. Statistical tests indicated that average quantities of surface variables for both the dry and wet cases were the same for the various update frequencies. However, spatially the results could be quite different especially in the wet case after it began to rain. The near-surface wind field was found to be different most of the time even for the dry case. In the wet case, rain intensities and average vertical profiles of heating associated with cloudy areas were found to differ for the various radiation update frequencies. The latter implies that the mean state of the model could be different as a result of not updating the radiation fields every time step and has important implications for longer term climate studies

Eastman, Joseph L.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Tropospheric Column Ozone Response to ENSO in GEOS-5 Assimilation of OMI and MLS Ozone Data

We use GEOS-5 analyses of Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) ozone observations to investigate the magnitude and spatial distribution of the El Nino Southern Oscillation (ENSO) influence on tropospheric column ozone (TCO) into the middle latitudes. This study provides the first explicit spatially resolved characterization of the ENSO influence and demonstrates coherent patterns and teleconnections impacting the TCO in the extratropics. The response is evaluated and characterized by both the variance explained and sensitivity of TCO to the Nino 3.4 index. The tropospheric response in the tropics agrees well with previous studies and verifies the analyses. A two-lobed response symmetric about the Equator in the western Pacific/Indonesian region seen in some prior studies and not in others is confirmed here. This two-lobed response is consistent with the large-scale vertical transport. We also find that the large-scale transport in the tropics dominates the response compared to the small-scale convective transport. The ozone response is weaker in the middle latitudes, but a significant explained variance of the TCO is found over several small regions, including the central United States. However, the sensitivity of TCO to the Nino 3.4 index is statistically significant over a large area of the middle latitudes. The sensitivity maxima and minima coincide with anomalous anti-cyclonic and cyclonic circulations where the associated vertical transport is consistent with the sign of the sensitivity. Also, ENSO related changes to the mean tropopause height can contribute significantly to the midlatitude response. Comparisons to a 22-year chemical transport model simulation demonstrate that these results from the 9- year assimilation are representative of the longer term. This investigation brings insight to several seemingly disparate prior studies of the El Nino influence on tropospheric ozone in the middle latitudes.

ML↗

Updraft Dynamics and Microphysics: On the Added Value of The Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

aerosols↗

Updraft Dynamics and Microphysics: On the Added Value of the Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether or under what conditions a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using two different reference frames: static cloudy-updraft grid cells versus tracked cumulus thermals. In both frameworks we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy-updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy-updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy-updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

Daniel Hernandez-Deckers↗

Can Aerosol Direct Radiative Effects Account for Analysis Increments of Temperature in the Tropical Atlantic?

In the late 1990's, prior to the launch of the Terra satellite, atmospheric general circulation models (GCMs) did not include aerosol processes because aerosols were not properly monitored on a global scale and their spatial distributions were not known well enough for their incorporation in operational GCMs. At the time of the first GEOS Reanalysis (Schubert et al. 1993), long time series of analysis increments (the corrections to the atmospheric state by all available meteorological observations) became readily available, enabling detailed analysis of the GEOS-1 errors on a global scale. Such analysis revealed that temperature biases were particularly pronounced in the Tropical Atlantic region, with patterns depicting a remarkable similarity to dust plumes emanating from the African continent as evidenced by TOMS aerosol index maps. Yoram Kaufman was instrumental encouraging us to pursue this issue further, resulting in the study reported in Alpert et al. (1998) where we attempted to assess aerosol forcing by studying the errors of a the GEOS-1 GCM without aerosol physics within a data assimilation system. Based on this analysis, Alpert et al. (1998) put forward that dust aerosols are an important source of inaccuracies in numerical weather-prediction models in the Tropical Atlantic region, although a direct verification of this hypothesis was not possible back then. Nearly 20 years later, numerical prediction models have increased in resolution and complexity of physical parameterizations, including the representation of aerosols and their interactions with the circulation. Moreover, with the advent of NASA's EOS program and subsequent satellites, atmospheric aerosols are now monitored globally on a routine basis, and their assimilation in global models are becoming well established. In this talk we will reexamine the Alpert et al. (1998) hypothesis using the most recent version of the GEOS-5 Data Assimilation System with assimilation of aerosols. We will explicitly calculate the impact of aerosols on the temperature analysis increments in the tropical Atlantic and assess the extent to which inclusion of atmospheric aerosols have reduced these increments.

General Circulation Models↗

Regional Impacts of Woodland Expansion on Nitrogen Oxide Emissions from Texas Savannahs: Combining Field, Modeling and Remote Sensing Approaches

Woody encroachment has contributed to documented changes world-wide and locally in the southwestern U.S. Specifically, in North Texas rangelands encroaching mesquite (Prosopis glandulosa var. glandulosa) a known N-fixing species has caused changes in aboveground biomass. While measurements of aboveground plant production are relatively common, measures of soil N availability are scarce and vary widely. N trace gas emissions (nitric and nitrous oxide) flom soils reflect patterns in current N cycling rates and availability as they are stimulated by inputs of organic and inorganic N. Quantification of N oxide emissions from savanna soils may depend upon the spatial distribution of woody plant canopies, and specifically upon the changes in N availability and cycling and subsequent N trace gas production as influenced by the shift from herbaceous to woody vegetation type. The main goal of this research was to determine whether remotely sensible parameters of vegetation structure and soil type could be used to quantify biogeochemical changes in N at local, landscape and regional scales. To accomplish this goal, field-based measurements of N trace gases were carried out between 2000-2001, encompassing the acquisition of imaging spectrometer data from the NASA Airborne Visible and Infrared Imaging Spectrometer (AVIRIS) on September 29, 2001. Both biotic (vegetation type and soil organic N) and abiotic (soil type, soil pH, temperature, soil moisture, and soil inorganic N) controls were analyzed for their contributions to observed spatial and temporal variation in soil N gas fluxes. These plot level studies were used to develop relationships between spatially extensive, field-based measurements of N oxide fluxes and remotely sensible aboveground vegetation and soil properties, and to evaluate the short-term controls over N oxide emissions through intensive field wetting experiments. The relationship between N oxide emissions, remotely-sensed parameters (vegetation cover, and soil type), and physical controls (soil moisture, and temperature) permitted the regional scale quantification of soil N oxides emissions. Landscape scale analysis linking N oxide emissions with cover change revealed an alleviation from N limitation following mesquite invasion. This study demonstrated the advantage of using N trace gases as a measure of ecosystem N availability combined with remote sensing to characterize the spatial heterogeneity in ecosystem parameters at a scale commensurate with field-based measurements of these properties. Woody vegetation encroachment provided an opportunity to capitalize on detection of the remotely-sensible parameter of woody cover as it relates to belowground biogeochemical processes that determine N trace gas production. The first spatially-explicit estimates of NO flux were calculated based on Prosopis fractional cover derived from high resolution remote sensing estimates of fractional woody cover (< 4 m) for a 120 sq km region of North Texas. An assessment of both N stocks and fluxes from the study revealed an alleviation of N limitation at this site experiencing recent woody encroachment. Many arid and semi-arid regions of the world are experiencing woody invasions, often of N-fixing species. The issue of woody encroachment is in the center of an ecological and political debate. Improving the links between biogeochemical processes and remote sensing of ecosystem properties will improve our understanding of biogeochemical processes at the regional scale, thus providing a means to address issues of land-use and land-cover change.

Asner, Gregory P.↗

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves, Precipitating Ring Current Ions, and Thermal Electron Heating - 2

This paper is dedicated to further presentations and discussions of the results from our new global self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2006; here referred to as Paper 1]. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation [for details see Paper 1]. To demonstrate the effects of the EMIC wave propagation and refraction on the RC proton precipitations and heating of the thermal plasmaspheric electrons we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. Firstly, the wave induced precipitations have a quite fine structure, and are highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 (raised dot) 10(exp 6) [(cm (raised dot) s (raised dot) sr)(sup -l)] are observed during the main and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not simply connected to the most intense EMIC waves. The character of the EMIC wave power spectral density distribution over the equatorial wave normal angle is an extremely crucial for the effectiveness of the RC ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from other ring current model [Kozyra et al., 1997] reveals that although we observe a qualitative agreement between localizations of the wave induced fluxes in the models, there is no quantitative agreement between the magnitudes of these fluxes. These differences are mainly due to a qualitative difference between the characters of the EMIC wave power spectral density distributions over the equatorial wave normal angle. Finally, the two energy sources to the plasmaspheric electrons are considered; (i) the heat fluxes caused by the EMIC wave energy absorption due to Landau resonance, and (ii) the heat fluxes due to Coulomb energy degradation of the RC o(+) ions. The heat fluxes caused by the EMIC wave energy absorption due to Landau resonance are observed in the postnoon-premidnight MLT sector, and maximize at the magnitude of 10l1 (eV/(cm(sup 2)(raised dot) s) at L=3.25, MLT=22 at 3400 UT after 1 May, 0000 UT. The greatest Coulomb energy deposition rates are about 2 (raised dot) 10(sup 10)(eV/(cm(sup 2)(raised dot) s) and observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the RC with plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, the RC, and the topside ionosphere [Gurgiolo et al., 20051.

Khazanov, G. V.↗