SEARCH · Engineering Papers

Results for “global shape 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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

NASA NTRS2014

Range-Specific High-Resolution Mesoscale Model Setup: Data Assimilation

Mesoscale weather conditions can have an adverse effect on space launch, landing, and ground processing at the Eastern Range (ER) in Florida and Wallops Flight Facility (WFF) in Virginia. During summer, land-sea interactions across Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) lead to sea breeze front formation, which can spawn deep convection that can hinder operations and endanger personnel and resources. Many other weak locally driven low-level boundaries and their interactions with the sea breeze front and each other can also initiate deep convection in the KSC/CCAFS area. Some of these other boundaries include the Indian River breeze front, Banana River breeze front, outflows from previous convection, horizontal convective rolls, convergence lines from other inland bodies of water such as Lake Okeechobee, the trailing convergence line from convergence of sea breeze fronts due to the shape of Cape Canaveral, frictional convergence lines from the islands in the Bahamas, convergence lines from soil moisture differences, convergence lines from cloud shading, and others. All these subtle weak boundary interactions often make forecasting of operationally important weather very difficult at KSC/CCAFS during the convective season (May-Oct). These convective processes often build quickly, last a short time (60 minutes or less), and occur over small distances, all of which also poses a significant challenge to the local forecasters who are responsible for issuing weather advisories, watches, and warnings. Surface winds during the transition seasons of spring and fall pose the most difficulties for the forecasters at WFF. They also encounter problems forecasting convective activity and temperature during those seasons. Therefore, accurate mesoscale model forecasts are needed to aid in their decision making. Both the ER and WFF would benefit greatly from high-resolution mesoscale model output to better forecast a variety of unique weather phenomena. Global and national scale models cannot properly resolve important local-scale weather features at each location due to their horizontal resolutions being much too coarse. Therefore, a properly tuned model at a high resolution is needed to provide improved capability. This task is a multi-year effort in which the Applied Meteorology Unit (AMU) will tune the Weather Research and Forecasting (WRF) model individually for each range. The goal of the first year, the results of which are in this report, was to tune the WRF model based on the best model resolution and run time while using reasonable computing capabilities. To accomplish this, the ER and WFF supported the tasking of the AMU to perform a number of sensitivity tests in order to determine the best model configuration for operational use at each of the ranges to best predict winds, precipitation, and temperature (Watson 2013). This task is a continuation of that work and will provide a recommended local data assimilation (DA) and numerical forecast model design optimized for the ER and WFF to support space launch activities. The model will be optimized for local weather challenges at both ranges.

Global and national scale models↗

NASA NTRS2014

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = −0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution.

Albedo algorithm↗

NASA NTRS1995

The structure of the distanct geomagnetic tail during long periods of northward IMF

We have used a newly developed, parallelized, global MHD magnetosphere - ionosphere simulation model with a 400 R(sub E) long tail to study the evolution, structure, and dynamics of the distant magnetotail during extended periods of northward interplanetary magnetic field (IMF). We find that the tail evolves to a nearly time stationary structure about one solar wind transit time after the IMF turns northward. Four regions of different magnetic topology can be distinguished which extend at least to the end of the simulation box at 400 R(sub E). Besides lobe field lines and open solar wind field lines tailward of an X-line, there is a broad boundary layer of closed field lines which we call the tail flank boundary layer (TFBL). Just inside the TFBL there is a region of closed field loops. Besides the X-line we find two O-lines which are enclosed by the closed field loops and are roughly aligned with the tail axis. Together they form a U shaped separator between the northward and the southward plasma sheet fields.

Raeder, J.↗

NASA NTRS2003

MODIS Retrieval of Dust Aerosol

The MODerate resolution Imaging Spectroradiometer (MODIS) currently aboard both the Terra and Aqua satellites produces a suite of products designed to characterize global aerosol distribution, optical thickness and particle size. Never before has a space-borne instrument been able to provide such detailed information, operationally, on a nearly global basis every day. The three years of Terra-MODIS data have been validated by comparing with co-located AERONET observations of aerosol optical thickness and derivations of aerosol size parameters. Some 8000 comparison points located at 133 AERONET sites around the globe show that the MODIS aerosol optical thickness retrievals are accurate to within the pre-launch expectations. However, the validation in regions dominated by desert dust is less accurate than in regions dominated by fine mode aerosol or background marine sea salt. The discrepancy is most apparent in retrievals of aerosol size parameters over ocean. In dust situations, the MODIS algorithm tends to under predict particle size because the reflectances at top of atmosphere measured by MODIS exhibit the stronger spectral signature expected by smaller particles. This pattern is consistent with the angular and spectral signature of non-spherical particles. All possible aerosol models in the MODIS Look-Up Tables were constructed from Mie theory, assuming a spherical shape. Using a combination of MODIS and AERONET observations, in regimes dominated by desert dust, we construct phase functions, empirically, with no assumption of particle shape. These new phase functions are introduced into the MODIS algorithm, in lieu of the original options for large dust-like particles. The results will be analyzed and examined.

Remer, Lorraine A.↗

NASA NTRS2012

Rapid Process to Generate Beam Envelopes for Optical System Analysis

The task of evaluating obstructions in the optical throughput of an optical system requires the use of two disciplines, and hence, two models: optical models for the details of optical propagation, and mechanical models for determining the actual structure that exists in the optical system. Previous analysis methods for creating beam envelopes (or cones of light) for use in this obstruction analysis were found to be cumbersome to calculate and take significant time and resources to complete. A new process was developed that takes less time to complete beam envelope analysis, is more accurate and less dependent upon manual node tracking to create the beam envelopes, and eases the burden on the mechanical CAD (computer-aided design) designers to form the beam solids. This algorithm allows rapid generation of beam envelopes for optical system obstruction analysis. Ray trace information is taken from optical design software and used to generate CAD objects that represent the boundary of the beam envelopes for detailed analysis in mechanical CAD software. Matlab is used to call ray trace data from the optical model for all fields and entrance pupil points of interest. These are chosen to be the edge of each space, so that these rays produce the bounding volume for the beam. The x and y global coordinate data is collected on the surface planes of interest, typically an image of the field and entrance pupil internal of the optical system. This x and y coordinate data is then evaluated using a convex hull algorithm, which removes any internal points, which are unnecessary to produce the bounding volume of interest. At this point, tolerances can be applied to expand the size of either the field or aperture, depending on the allocations. Once this minimum set of coordinates on the pupil and field is obtained, a new set of rays is generated between the field plane and aperture plane (or vice-versa). These rays are then evaluated at planes between the aperture and field, at a desired number of steps perceived necessary to build up the bounding volume or cone shape. At each plane, the ray coordinates are again evaluated using the convex hull algorithm to reduce the data to a minimal set. When all of the coordinates of interest are obtained for every plane of the propagation, the data is formatted into an xyz file suitable for FRED optical analysis software to import and create a STEP file of the data. This results in a spiral-like structure that is easily imported by mechanical CAD users who can then use an automated algorithm to wrap a skin around it and create a solid that represents the beam.

Howard, Joseph↗

NASA NTRS2014

A Preliminary Analysis of Precipitation Properties and Processes during NASA GPM IFloodS

The Iowa Flood Studies (IFloodS) is a NASA Global Precipitation Measurement (GPM) ground measurement campaign, which took place in eastern Iowa from May 1 to June 15, 2013. The goals of the field campaign were to collect detailed measurements of surface precipitation using ground instruments and advanced weather radars while simultaneously collecting data from satellites passing overhead. Data collected by the radars and other ground instruments, such as disdrometers and rain gauges, will be used to characterize precipitation properties throughout the vertical column, including the precipitation type (e.g., rain, graupel, hail, aggregates, ice crystals), precipitation amounts (e.g., rain rate), and the size and shape of raindrops. The impact of physical processes, such as aggregation, melting, breakup and coalescence on the measured liquid and ice precipitation properties will be investigated. These ground observations will ultimately be used to improve rainfall estimates from satellites and in particular the algorithms that interpret raw data for the upcoming GPM mission's Core Observatory satellite, which launches in 2014. The various precipitation data collected will eventually be used as input to flood forecasting models in an effort to improve capabilities and test the utility and limitations of satellite precipitation data for flood forecasting. In this preliminary study, the focus will be on analysis of NASA NPOL (S‐band, polarimetric) radar (e.g., radar reflectivity, differential reflectivity, differential phase, correlation coefficient) and NASA 2D Video Disdrometers (2DVDs) measurements. Quality control and processing of the radar and disdrometer data sets will be outlined. In analyzing preliminary cases, particular emphasis will be placed on 1) documenting the evolution of the rain drop size distribution (DSD) as a function of column melting processes and 2) assessing the impact of range on ground‐based polarimetric radar estimates of DSD properties.

Carey, Lawrence↗

NASA NTRS2024

Enabling the Assimilation of CrIS Shortwave Infrared Observations in Global NWP at NOAA. Part II: OSEs and Results

The assimilation of data from hyperspectral infrared sounders in global data assimilation systems has historically been focused on observations in the longwave infrared (LWIR) region of the spectrum (650-1100 cm−1), despite the often concurrent availability of measurements from the shortwave infrared (SWIR) region of the spectrum (2150-2550 cm−1), because issues (like solar effects) have generally prevented the assimilation of SWIR observations. Recent advances in radiative transfer models have worked to address some of the previous challenges in simulating SWIR observations, and the assimilation of SWIR data (e.g. from potential future small-satellites) is now a feasible prospect. Still, a better understanding of how these observations perform in a data assimilation system and impact resulting analyses and NWP forecasts is necessary. In this study, the value of SWIR observations in global NWP is assessed by assimilating SWIR observations from the Cross-track Infrared Sounder (CrIS) in NOAA’s Global Data Assimilation System (GDAS). The methodologies used to enable the assimilation of these observations, including the implementation of a scene-dependent observation error and the enhancement of quality control procedures, are discussed, as are the results of Observing System Experiments (OSEs) conducted to evaluate the impact of assimilating SWIR observations on forecast skill. The overall results show that SWIR assimilation produces similar forecast impacts to LWIR assimilation. The ability to demonstrate that the assimilation or SWIR observations in NWP is a realistic prospect may help to shape future constellations of small-satellites to serve as a beneficial complement to the current constellation if hyperspectral IR sounders.

Erin Jones↗

NASA NTRS2001

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗

NASA NTRS2012

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Cole, Benjamin H.↗

NASA NTRS2004

On Study of Air/Space-borne Dual-Wavelength Radar for Estimates of Rain Profiles

In this study, a framework is discussed to apply air/space-borne dual-wavelength radar for the estimation of characteristic parameters of hydrometeors. The focus of our study is on the Global Precipitation Measurements (GPM) precipitation radar, a dual-wavelength radar that operates at Ku (13.8 GHz) and Ka (35 GHz) bands. As the droplet size distributions (DSD) of rain are expressed as the Gamma function, a procedure is described to derive the median volume diameter (D(sub 0)) and particle number concentration (N(sub T)) of rain. The correspondences of an important quantity of dual-wavelength radar, defined as deferential frequency ratio (DFR), to the D(sub 0) in the melting region are given as a function of the distance from the 0 C isotherm. A self-consistent iterative algorithm that shows a promising to account for rain attenuation of radar and infer the DSD without use of surface reference technique (SRT) is examined by applying it to the apparent radar reflectivity profiles simulated from the DSD model and then comparing the estimates with the model (true) results. For light to moderate rain the self-consistent rain profiling approach converges to unique and correct solutions only if the same shape factors of Gamma functions are used both to generate and retrieve the rain profiles, but does not converges to the true solutions if the DSD form is not chosen correctly. To further examine the dual-wavelength techniques, the self-consistent algorithm, along with forward and backward rain profiling algorithms, is then applied to the measurements taken from the 2nd generation Precipitation Radar (PR-2) built by Jet Propulsion Laboratory. It is found that rain profiles estimated from the forward and backward approaches are not sensitive to shape factor of DSD Gamma distribution, but the self-consistent method is.

Liao, Liang↗

NASA NTRS2007

Effects of RAM Exposure on a Low Earth Orbit BroadBand Radiometer (BBR): CERES Experience and Implications for EarthCARE

In order to best detect real changes in the Earth's climate system it is estimated that space based instrumentation measuring the global Earth Radiation Budget (ERB) must remain calibrated with a stability of 0.3% per decade. This level of stability is beyond the specified accuracy of existing ERB programs such as the Clouds and the Earth's Radiant Energy System (CERES, using three broadband radiometric scanning channels: the shortwave (SW 0.3 - 5?m), total (0.3 - >100 micron), and window (8 - 12 micron)). When in low earth orbit, it has been shown that optical response to blue-UV radiance can be reduced significantly (> 3%) due to UV hardened contaminants deposited on the surface of the optics. Evidence suggests that exposure of telescope optics to the forward looking ram direction is the primary cause of this contamination build up. With typical onboard calibration lamps emitting very low energy in the blue-UV region, this darkening is not directly measurable using standard internal calibration techniques. This paper describes a study using a model of ram exposure induced contaminant deposition and darkening, in conjunction standard established in-flight vicarious and internal calibration techniques to derive the spectral shape of the darkening to which a broadband instrument is subjected. The results of the model when applied to the CERES instruments are shown. These findings are of great importance to the EarthCARE project, whose BBR uses one broadband telescope permanently looking forward at 45 degrees, with continual exposure to the ram direction. Specific attention may therefore be needed in the design of BBR optics and on-board calibration in order to prevent or compensate for the spectral darkening seen in the CERES project.

Matthews, Grant↗

NASA NTRSDate not supplied

Composition Reanalyses at NASA’s Global Modeling and Assimilation Office

Scientific Assessments of Ozone Depletion 2018 and 2022 report the recovery of Antarctic ozone, consistent with diminishing concentrations of halogens in the stratosphere. Both the size and the depth of the springtime ozone holes exhibit downward trends, particularly in September. Superimposed on the recovery trend are significant year-to-year fluctuations arising from dynamical variability potentially influenced by climate forcing and volcanic eruptions. Since 2018 the Antarctic stratosphere has experienced a wild parade of austral spring seasons with a sudden stratospheric warming and exceptionally high ozone in 2019 followed by large and long-lasting ozone holes in the following years, prompting questions about the speed and detectability of Antarctic ozone recovery and its potential connections to human-induced climate change. From the Australian New Years’ wildfires of 2020 to volcanic aerosols to the unprecedented injection of water vapor by the Hunga volcanic eruption, recent scientific literature on the subject investigates many factors that may or may not have contributed to this remarkable interannual variability. This presentation summarizes dynamical influences on springtime Antarctic stratospheric ozone in recent years and discusses these results in the context of long-term trends. In particular, this work seeks to isolate dynamical and chemical contributions to the Antarctic springs since 2019 using a consistent methodological framework. Our analysis uses constituent and meteorological fields from the MERRA-2 Stratospheric Composition Reanalysis of Aura MLS (M2-SCREAM) in combination with data from the Atmospheric Chemistry Experiment – Fourier Transform Spectrometer and from ozone sondes. Using the high spatial and temporal resolution of M2-SCREAM, we analyze the daily evolution of southern polar ozone, HCl, N2O, and water vapor in dynamical flow-following coordinates relative to the location, size, and shape of the polar vortex. For vertically integrated stratospheric ozone we emphasize the importance of the wave-influenced geometry of the polar vortex and its vertical alignment. The ultimate purpose of this project is to contribute to the quantification of the speed of the Antarctic ozone recovery and associated uncertainties.

Krzysztof Wargan↗

NASA NTRS1993

Resolving topographic detail on Venus by modeling complex Magellan altimetry echoes

Magellan's altimeter is providing some of the finest resolution topography of Venus achieved to date. Nevertheless, efforts continue to improve the topographic resolution whenever possible. One effort to this end is stereoscopic imaging, which provides topography at scales similar to that of the synthetic aperture radar (SAR). However, this technique requires two SAR images of the same site to be obtained and limits the utility of this method. In this paper, we present another method to resolve topographic features at scales smaller than that of an altimeter footprint, which is more globally applicable than the stereoscopic approach. Each pulse which is transmitted by Magellan's altimeter scatters from the planet and echoes to the receiver, delayed based on the distance between the spacecraft and each surface element. As resolved in time, each element of an altimetry echo represents the sum of all points on the surface which are equidistant from the spacecraft. Thus, individual returns, as a function of time, create an echo profile which may be used to derive properties of the surface, such as the scattering law or, in this case, the topography within the footprint. The Magellan project has derived some of this information by fitting model templates to radar echo profiles. The templates are calculated based on Hagfor's Law, which assumes a smooth, gently undulating surface. In most regions these templates provide a reasonable fit to the observed echo profile; however, in some cases the surface departs from these simple assumptions and more complex profiles are observed. Specifically, we note that sub-footprint topographic relief apparently has a strong effect on the shape of the echo profile. To demonstrate the effects of sub-resolution relief on echo profiles, we have calculated the echo shapes from a wide range of simple topographic models. At this point, our topographic models have emphasized surfaces where only two dominant elevations are contained within a footprint, such as graben, ridges, crater rims, and central features in impact craters.

Lovell, Amy J.↗

NASA NTRS2019

Mid-Latitude Freshwater Availability Reduced by Projected Vegetation Responses to Climate Change

Plants are expected to generate more global-scale runoff under increasing atmospheric carbon dioxide concentrations through their influence on surface resistance to evapotranspiration. Recent studies using Earth System Models from phase 5 of the Coupled Model Intercomparison Project ostensibly reaffirm this result, further suggesting that plants will ameliorate the dire reductions in water availability projected by other studies that use aridity metrics. Here we complicate this narrative by analysing the change in precipitation partitioning to plants, runoff and storage in multiple Earth system models under both high carbon dioxide concentrations and warming. We show that projected plant responses directly reduce future runoff across vast swaths of North America, Europe and Asia because bulk canopy water demands increase with additional vegetation growth and longer and warmer growing seasons. These runoff declines occur despite increased surface resistance to evapotranspiration and vegetation total water use efficiency, even in regions with increasing or unchanging precipitation. We demonstrate that constraining the large uncertainty in the multimodel ensemble with regional-scale observations of evapotranspiration partitioning strengthens these results. We conclude that terrestrial vegetation plays a large and unresolved role in shaping future regional freshwater availability, one that will not ubiquitously ameliorate future warming-driven surface drying.

Mankin, Justin S.↗

NASA NTRS1994

Pinatubo Aerosol Evolution: Using Composite Data Sets to Build the Global- to Micro-Scale Picture and Assess Consistency of Different Measurements

This paper brings together experimental evidence required to build realistic models of the global evolution of physical, chemical, and optical properties of the aerosol resulting from the 1991 Pinatubo volcanic eruption. Such models are needed to compute the effects of the aerosol on atmospheric chemistry, dynamics, radiation, and temperature. Whereas there is now a large and growing body of post-Pinatubo measurements by a variety of techniques, some results are in conflict, and a self-consistent, unified picture is needed, along with an assessment of remaining uncertainties. This paper examines data from photometers, radiometers, impactors, optical counters/sizers, and lidars operated on the ground, aircraft, balloons, and spacecraft. Example data sources include: - Tracking sunphotometers and lidars at Mauna Loa Observatory (MLO) and on the DC-8 - Particle spectrometers and wire impactors on the ER-2 and DC-8 - Dustsondes (particle counters/sizers on balloons) - SAGE II, SAM II, AVHRR, CLAES, and ISAMS sensors on a variety of satellites. We assess the mutual consistency of these disparate data sets and recommend 'consensus" properties and uncertainties in the process of developing a composite data set. Recommended properties include the spatial and temporal evolution of particle chemical composition, shape, wavelength and temperature-dependent refractive index, size distribution, and optical depth spectra. Supporting references are cited and representative data shown.

Russell, P. B.↗

NASA NTRS2003

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The May 2-7, 1998, Storm

Complete description of a self-consistent model for magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves, and back on waves, are considered self-consistently by solving both equations on a global magnetospheric scale under non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗

NASA NTRS2000

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗