Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

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 1,027 records · Page 57

Local Guided Wavefield Analysis for Characterization of Delaminations in Composites

Delaminations in composite laminates resulting from impact events may be accompanied by minimal indication of damage at the surface. As such, inspection techniques are required to ensure defects are within allowable limits. Conventional ultrasonic scanning techniques have been shown to effectively characterize the size and depth of delaminations but require physical contact with the structure. Alternatively, a noncontact scanning laser vibrometer may be used to measure guided wave propagation in the laminate structure. A local Fourier domain analysis method is presented for processing guided wavefield data to estimate spatially-dependent wavenumber values, which can be used to determine delamination depth. The technique is applied to simulated wavefields and results are analyzed to determine limitations of the technique with regards to determining defect size and depth. Finally, experimental wavefield data obtained in quasi-isotropic carbon fiber reinforced polymer (CFRP) laminates with impact damage is analyzed and wavenumber is measured to an accuracy of 8.5% in the region of shallow delaminations. Keywords: Ultrasonic wavefield imaging, Windowed Fourier transforms, Guided waves, Structural health monitoring, Nondestructive evaluation

Rogge, Matthew D.↗

Model-Based Estimation of Sampling-Caused Uncertainty in Aerosol Remote Sensing for Climate Research Applications

To evaluate the effect of sampling frequency on the global monthly mean aerosol optical thickness (AOT), we use 6 years of geographical coordinates of Moderate Resolution Imaging Spectroradiometer (MODIS) L2 aerosol data, daily global aerosol fields generated by the Goddard Institute for Space Studies General Circulation Model and the chemical transport models Global Ozone Chemistry Aerosol Radiation and Transport, Spectral Radiationtransport Model for Aerosol Species and Transport Model 5, at a spatial resolution between 1.125 deg × 1.125 deg and 2 deg × 3◦: the analysis is restricted to 60 deg S-60 deg N geographical latitude. We found that, in general, the MODIS coverage causes an underestimate of the global mean AOT over the ocean. The long-term mean absolute monthly difference between all and dark target (DT) pixels was 0.01-0.02 over the ocean and 0.03-0.09 over the land, depending on the model dataset. Negative DT biases peak during boreal summers, reaching 0.07-0.12 (30-45% of the global long-term mean AOT). Addition of the Deep Blue pixels tempers the seasonal dependence of the DT biases and reduces the mean AOT difference over land by 0.01-0.02. These results provide a quantitative measure of the effect the pixel exclusion due to cloud contamination, ocean sun-glint and land type has on the MODIS estimates of the global monthly mean AOT. We also simulate global monthly mean AOT estimates from measurements provided by pixel-wide along-track instruments such as the Aerosol Polarimetry Sensor and the Cloud-Aerosol LiDAR with Orthogonal Polarization. We estimate the probable range of the global AOT standard error for an along-track sensor to be 0.0005-0.0015 (ocean) and 0.0029-0.01 (land) or 0.5-1.2% and 1.1-4% of the corresponding global means. These estimates represent errors due to sampling only and do not include potential retrieval errors. They are smaller than or comparable to the published estimate of 0.01 as being a climatologically significant change in the global mean AOT, suggesting that sampling density is unlikely to limit the use of such instruments for climate applications at least on a global, monthly scale.

long-term variability↗

The Parallel System for Integrating Impact Models and Sectors (pSIMS)

We present a framework for massively parallel climate impact simulations: the parallel System for Integrating Impact Models and Sectors (pSIMS). This framework comprises a) tools for ingesting and converting large amounts of data to a versatile datatype based on a common geospatial grid; b) tools for translating this datatype into custom formats for site-based models; c) a scalable parallel framework for performing large ensemble simulations, using any one of a number of different impacts models, on clusters, supercomputers, distributed grids, or clouds; d) tools and data standards for reformatting outputs to common datatypes for analysis and visualization; and e) methodologies for aggregating these datatypes to arbitrary spatial scales such as administrative and environmental demarcations. By automating many time-consuming and error-prone aspects of large-scale climate impacts studies, pSIMS accelerates computational research, encourages model intercomparison, and enhances reproducibility of simulation results. We present the pSIMS design and use example assessments to demonstrate its multi-model, multi-scale, and multi-sector versatility.

crop modeling↗

Global Precipitation Measurement (GPM) Mission Products and Services at the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC)

The NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC) hosts and distributes GPM data within the NASA Earth Observation System Data Information System (EOSDIS). The GES DISC is also home to the data archive for the GPM predecessor, the Tropical Rainfall Measuring Mission (TRMM). Over the past 17 years, the GES DISC has served the scientific as well as other communities with TRMM data and user-friendly services. During the GPM era, the GES DISC will continue to provide user-friendly data services and customer support to users around the world. GPM products currently and to-be available: -Level-1 GPM Microwave Imager (GMI) and partner radiometer products, DPR products -Level-2 Goddard Profiling Algorithm (GPROF) GMI and partner products, DPR products -Level-3 daily and monthly products, DPR products -Integrated Multi-satellitE Retrievals for GPM (IMERG) products (early, late, and final) A dedicated Web portal (including user guides, etc.) has been developed for GPM data (http://disc.sci.gsfc.nasa.gov/gpm). Data services that are currently and to-be available include Google-like Mirador (http://mirador.gsfc.nasa.gov/) for data search and access; data access through various Web services (e.g., OPeNDAP, GDS, WMS, WCS); conversion into various formats (e.g., netCDF, HDF, KML (for Google Earth), ASCII); exploration, visualization, and statistical online analysis through Giovanni (http://giovanni.gsfc.nasa.gov); generation of value-added products; parameter and spatial subsetting; time aggregation; regridding; data version control and provenance; documentation; science support for proper data usage, FAQ, help desk; monitoring services (e.g. Current Conditions) for applications. The United User Interface (UUI) is the next step in the evolution of the GES DISC web site. It attempts to provide seamless access to data, information and services through a single interface without sending the user to different applications or URLs (e.g., search, access, subset, Giovanni, documents).

precipitation↗

PBL Height from AIRS, GPS RO, and MERRA-2 Products in NASA GES DISC and Their 10 Year Seasonal Mean Intercomparison

Within the planetary boundary layer (PBL), surface forcing response, drag, turbulence, and vertical mixing are important processes and play a more critical role here than in the overlying “free atmosphere”. The PBL Height (PBLH) is an important parameter in climate models, weather forecasts, and air quality prediction. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provides data processing, archiving, and distribution services for numerous Earth science products. PBLH is a parameter in three products served by GES DISC, which are from the Atmospheric Infrared Sounder (AIRS), the Global Positioning System (GPS) radio occultation (RO) experiment, and the NASA reanalysis product Modern-Era Retrospective analysis for Research and Applications – 2 (MERRA-2). These products have different spatial and temporal resolutions and coverages, and their PBLH definitions are also different. To better serve the PBL research community, we have summarized the specifications of these products. A ten-year seasonal mean intercomparison is also conducted to provide further guidance to users. The intercomparison results show that MERRA-2 has a much shallower PBL than AIRS and GPS RO. An experimental study indicates the different PBLH definition in MERRA-2 caused smaller values of PBLH. The improvement of the water vapor retrieval in AIRS version 7 over version 6 results in the version 7 PBLH agreeing better with GPS RO and MERRA-2 than version 6, especially near the equator and low latitudes.

Feng Ding↗

A Concept for the Deployment of a Large Lunar Crater Radio Telescope Using Teams of Tethered Robots

Kilometer-scale craters on the far side of the Moon have unique potential as future locations for large radio telescopes, which can observe the universe at wavelengths and frequencies (> 10 m, < 30 MHz) not possible with conventional Earth or orbital-based approaches. Distinct advantages of building a Lunar Crater Radio Telescope (LCRT) on the far side include i) isolation from radio noise due to the Earth’s ionosphere, orbiting satellites, and the Sun, ii) days of uninterrupted dark/cold sky viewing during lunar night, and iii) terrain geometry naturally suited for constructing the largest mesh antenna structure in the Solar System. A key challenge to constructing LCRT on the Moon is related to the complexity of deploying a 1-km diameter antenna and hanging receiver within a lunar crater whose diameter, depth, and slope are 3-5 km, 1 km, and ~30 degrees respectively. In this paper, we first evaluate the trade space for deploying a large, complex structure within a crater, and then provide a more detailed concept evaluation of our favored approach, which employs coordinated teams of tethered rovers to extract and suspend a folded antenna from a lander at the base of a crater. NASA’s Jet Propulsion Laboratory in collaboration with California Institute of Technology (Caltech) have developed a novel robotic system for accessing extremely steep terrains; the Axel rover is a two-wheeled rugged terrain vehicle that is supported by an electro-mechanical tether that provides power, data, and tensile support from a top-side anchor location. Recently, a pair of Axel robots have been used in a DuAxel configuration that allows for four-wheel driving and repeated passive anchoring at different locations. The DuAxel system has unique advantages for deploying an LCRT antenna, including the ability to deploy from a lander near a crater, drive a distance to the crater rim to deploy an Axel, and later, retract the deployed Axel in order to sequentially lift up sections of the antenna. Our proposed concept involves delivering a packaged antenna and receiver to the bottom-center of a crater floor on a lander, then later sending a team of multiple DuAxel rovers to retrieve guide wires from the lander, which are pulled to the top of the crater. We explore this concept in detail and provide some initial quantitative analysis to demonstrate the feasibility of our system with respect to the spatial and mass properties of the antenna as juxtaposed to DuAxel capabilities. Finally, we outline next steps towards validating our concept on the way to a future lunar deployment opportunity.

Hallinan, Greg↗

Towards a program of record of inland water quality: Exploiting present and heritage multispectral sensors for maximum information extraction

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This research exploits recent advancements in bio-optical modeling, cloud computing, and machine learning to enhance our capacity to leverage present and heritage satellite data. Recent research suggests that sensors with low spectral resolution, such as Sentinel 2 and Landsat missions, contain enough hidden spectral variation which can be exploited using data-driven approaches. The availability of three decades of archival imagery will open doors to discover global trends of eutrophication and increased cyanobacteria dominance and provide valuable insight to the development of predictive methodologies. Preliminary efforts in synthetic emulation of global natural inland waters will be discussed and contextualized against satellite radiometric measurement uncertainty, satellite data product uncertainties and causal signal ambiguity over the visible wavelength range, supported by high quality field and image data for selected inland aquatic sites. Insights on water quality estimation via data-driven machine learning models versus matrix inversions will be discussed, and how we can exploit spectral-spatial relationships in high spatial resolution data. A cross-sensor synergistic approach with detailed uncertainty analysis based on optical water types, will allow for unprecedented global snapshots of fine scale ecological dynamics of inland waters.

Inland↗

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

Soil Moisture Active Passive (SMAP) mission brightness temperature (T(b) ) observations are assimilated into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update simulations of surface and root-zone soil moisture. Different time series components of the T(b) observations are assimilated including anomalies, inter-annual variations, and high frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the uncertainties of modeled and observed T(b) is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. In a benchmark experiment, T(b) observations are assimilated using a spatially constant 4 Kelvin (K) observation uncertainty, as in the operational SMAP Level-4 algorithm. All DA experiments exhibit notable skill improvements in most regions. Improvements are largest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for surface soil moisture simulations (0.06). No notable difference in skill - neither improvement nor deterioration - is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. This may be explained by the presence of large observation operator errors, which are analytically shown to have the potential to render post-update uncertainty insensitive to inaccuracies in estimates of the Kalman gain. These results have important implications for the design of soil moisture DA systems, in particular for parameterizing model and observation uncertainties.

Hydrology↗

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic Density Structures (PDSs) are quasi-periodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (Lx) of tens to several thousand Mm, thus belonging to the class of “mesoscale structures”. Current interplanetary multi-spacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multi-spacecraft observations. For solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun-Earth line (Ly) ranged between ≈30 and ≈200 Earth Radii (RE). Here, we implemented a similar approach on the same parameters but focused on high density slow solar wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity obtaining two PDS groups based on their radial length scale: Lx1 ≈86RE and Lx2 ≈35RE. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the Ly order of magnitudes for the two PDS groups: Ly1 ≈340RE and Ly2 ≈187RE. Knowing the PDSs’ size scales is fundamental to constrain models aimed at reproducing these structures and is critical to better understand the PDS-magnetosphere coupling.

Interplanetary physics↗

Internal Consistency of the NVAP Water Vapor Dataset

The NVAP (NASA Water Vapor Project) dataset is a global dataset at 1 x 1 degree spatial resolution consisting of daily, pentad, and monthly atmospheric precipitable water (PW) products. The analysis blends measurements from the Television and Infrared Operational Satellite (TIROS) Operational Vertical Sounder (TOVS), the Special Sensor Microwave/Imager (SSM/I), and radiosonde observations into a daily collage of PW. The original dataset consisted of five years of data from 1988 to 1992. Recent updates have added three additional years (1993-1995) and incorporated procedural and algorithm changes from the original methodology. Since each of the PW sources (TOVS, SSM/I, and radiosonde) do not provide global coverage, each of these sources compliment one another by providing spatial coverage over regions and during times where the other is not available. For this type of spatial and temporal blending to be successful, each of the source components should have similar or compatible accuracies. If this is not the case, regional and time varying biases may be manifested in the NVAP dataset. This study examines the consistency of the NVAP source data by comparing daily collocated TOVS and SSM/I PW retrievals with collocated radiosonde PW observations. The daily PW intercomparisons are performed over the time period of the dataset and for various regions.

Suggs, Ronnie J.↗

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗

An Iterated Global Mascon Solution with Focus on Land Ice Mass Evolution

Land ice mass evolution is determined from a new GRACE global mascon solution. The solution is estimated directly from the reduction of the inter-satellite K-band range rate observations taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons are estimated with 10-day and 1-arc-degree equal area sampling, applying anisotropic constraints for enhanced temporal and spatial resolution of the recovered land ice signal. The details of the solution are presented including error and resolution analysis. An Ensemble Empirical Mode Decomposition (EEMD) adaptive filter is applied to the mascon solution time series to compute timing of balance seasons and annual mass balances. The details and causes of the spatial and temporal variability of the land ice regions studied are discussed.

Luthcke, S. B.↗

High Resolution Aerosol Data from MODIS Satellite for Urban Air Quality Studies

The Moderate Resolution Imaging Spectroradiometer (MODIS) provides daily global coverage, but the 10 km resolution of its aerosol optical depth (AOD) product is not suitable for studying spatial variability of aerosols in urban areas. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was developed for MODIS which provides AOD at 1 km resolution. Using MAIAC data, the relationship between MAIAC AOD and PM(sub 2.5) as measured by the 27 EPA ground monitoring stations was investigated. These results were also compared to conventional MODIS 10 km AOD retrievals (MOD04) for the same days and locations. The coefficients of determination for MOD04 and for MAIAC are R(exp 2) =0.45 and 0.50 respectively, suggested that AOD is a reasonably good proxy for PM(sub 2.5) ground concentrations. Finally, we studied the relationship between PM(sub 2.5) and AOD at the intra-urban scale (10 km) in Boston. The fine resolution results indicated spatial variability in particle concentration at a sub-10 kilometer scale. A local analysis for the Boston area showed that the AOD-PM(sub 2.5) relationship does not depend on relative humidity and air temperatures below approximately 7 C. The correlation improves for temperatures above 7 - 16 C. We found no dependence on the boundary layer height except when the former was in the range 250-500 m. Finally, we apply a mixed effects model approach to MAIAC aerosol optical depth (AOD) retrievals from MODIS to predict PM(sub 2.5) concentrations within the greater Boston area. With this approach we can control for the inherent day-to-day variability in the AOD-PM(sub 2.5) relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance. Our results show that the model-predicted PM(sub 2.5) mass concentrations are highly correlated with the actual observations (out-of-sample R(exp 2) of 0.86). Therefore, adjustment for the daily variability in the AOD-PM(sub 2.5) relationship provides a means for obtaining spatially-resolved PM(sub 2.5) concentrations.

remote sensing↗

Tracking Climate Effects on Plant-Pollinator Interaction Phenology with Satellites and Honey Bee Hives

Background/Question/Methods: The complexity of plant-pollinator interactions, the large number of species involved, and the lack of species response functions present challenges to understanding how these critical interactions may be impacted by climate and land cover change on large scales. Given the importance of this interaction for terrestrial ecosystems, it is desirable to develop new approaches. We monitor the daily weight change of honey bee (Apis mellifera) colonies to record the phenology of the Honey Bee Nectar Flow (HBNF) in a volunteer network (honeybeenet.gsfc.nasa.gov). The records document the successful interaction of a generalist pollinator with a variety of plant resources. We extract useful HBNF phenology metrics for three seasons. Sites currently exist in 35 states/provinces in North America, with a concentration in the Mid-Atlantic region. HBNF metrics are compared to standard phenology metrics derived from remotely sensed vegetation indices from NASA's MODIS sensor and published results from NOAA's A VHRR. At any given time the percentage of plants producing nectar is usually a sma11 fraction of the total satellite sensor signal. We are interested in determining how well the 'bulk' satellite vegetation parameters relate to the phenology of the HBNF, and how it varies spatially on landscape to continental scales. Results/Conclusions: We found the median and peak seasonal HBNF dates to be robust, with variation between replicate scale hives of only a few days. We developed quality assessment protocols to identify abnormal colony artifacts. Temporally, the peak and median of the HBNF in the Mid-Atlantic show a significant advance of 0.58 d/y beginning about 1970, very similar to that observed by the A VHRR since 1982 (0.57 d/y). Spatially, the HBNF metrics are highly correlated with elevation and winter minimum temperature distribution, and exhibit significant but regionally coherent inter-annual variation. The relationship between median of the spring HBNF with the "Green-up" metric from the 500 meter MODIS NDVI phenology product, for sites throughout the Eastern US 2000-2009, is well described by a single linear fit (r(exp 2) = 0.72). We conclude.that for the tree-dominated areas of the Eastern US at least the spring HBNF can be tracked very well by MODIS phenology. Analysis of other regions and seasons is presently underway but with more limited data. Spatial patterns in the eastern US and management implications will be presented and discussed.

Esaias, Wayne E.↗

Nominal 30-M Cropland Extent Map of Continental Africa by Integrating Pixel-Based and Object-Based Algorithms Using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, is a starting point to develop high-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based Geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of RF happened due to the noise of input training data, Support Vector Machine (SVM) was applied to compensate for such defects in specific areas. Third, the Recursive Hierarchical Segmentation (RHSeg) algorithm was employed to generate an object-oriented segmentation layer based on spectral and spatial properties from the same input data. This layer was merged with the pixel-based classification to improve segmentation accuracy. Accuracies of the merged 30-m crop extent product were computed using an error matrix approach in which 1754 independent validation samples were used. In addition, a comparison was performed with other available cropland maps as well as with LULC maps to show spatial similarity. Finally, the cropland area results derived from the map were compared with UN FAO statistics. The independent accuracy assessment showed a weighted overall accuracy of 94, with a producers accuracy of 85.9 (or omission error of 14.1), and users accuracy of 68.5 (commission error of 31.5) for the cropland class. The total net cropland area (TNCA) of Africa was estimated as 313 Mha for the nominal year 2015.

Cropland mapping; cropland areas; 30-m; Landsat-8;↗

Spectral analysis of gappy images with application to Sea WiFS data

For various reasons, such as the presence of clouds, satellite orbit configuration, data sampling pattern, remote sensing instruments yield irregularly spaced or gappy/patchy measurements. To obviate these problems, we use original, non-gridded satellite data and construct a spatio-temporal (or just spatial) autocorrelation function of the field's variation.

image analysis spectral analysis Sea WiFS↗

Chemical Observations of a Polar Vortex Intrusion

An intrusion of vortex edge air in D the interior of the Arctic polar vortex was observed on the January 31,2005 flight of the NASA DC-8 aircraft. This intrusion was identified as anomalously high values of ozone by the AROTAL and DIAL lidars. Our analysis shows that this intrusion formed when a blocking feature near Iceland collapsed, allowing edge air to sweep into the vortex interior. along the DC-8 flight track also shows the intrusion in both ozone and HNO3. Polar Stratospheric Clouds (PSCs) were observed by the DIAL lidar on the DC-8. The spatial variability of the PSCs can be explained using MLS HNO3 and H2O observations and meteorological analysis temperatures. We also estimate vortex denitrification using the relationship between N2O and HNO3. Reverse domain fill back trajectory calculations are used to focus on the features in the MLS data. The trajectory results improve the agreement between lidar measured ozone and MLS ozone and also improve the agreement between the HNO3 measurements PSC locations. The back trajectory calculations allow us to compute the local denitrification rate and reduction of HCl within the filament. We estimate a denitrification rate of about lO%/day after exposure to below PSC formation temperature. Analysis of Aura MLS observations made

Schoeberl, M. R.↗

Mid-Latitude Composition of Mars from Thermal and Epithermal Neutrons

Epithermal neutron data acquired by Mars Odyssey have been analyzed to determine global maps of water-equivalent hydrogen abundance. By assuming that hydrogen was distributed uniformly with depth within the surface, a map of minimum water abundance was obtained. The addition of thermal neutrons to this analysis could provide information needed to determine water stratigraphy. For example, thermal and epithermal neutrons have been used together to determine the depth and abundance of water-equivalent hydrogen of a buried layer in the south polar region. Because the emission of thermal neutrons from the Martian surface is sensitive to absorption by elements other than hydrogen, analysis of stratigraphy requires that the abundance of these elements be known. For example, recently published studies of the south polar region assumed that the Mars Pathfinder mean soil composition is representative of the regional soil composition. This assumption is partially motivated by the fact that Mars appears to have a well-mixed global dust cover and that the Pathfinder soil composition is representative of the mean composition of the Martian surface. In this study, we have analyzed thermal and epithermal neutron data measured by the neutron spectrometer subsystem of the gamma ray spectrometer to determine the spatial distribution of the composition of elements other than hydrogen. We have restricted our analysis to mid-latitude regions for which we have corrected the neutron counting data for variations in atmospheric thickness.

T H Prettyman↗