Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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,171 records · Page 65

The Intensity Distribution of Faint Gamma-Ray Bursts Detected with BATSE

We have recently completed a search of 6 years of archival BATSE data for gamma-ray bursts (GRBS) that were too faint to activate the real-time burst detection system running on board the spacecraft. These "nontriggered" bursts can be combined with the "triggered" bursts detected on board to produce a GRB intensity distribution that reaches peak fluxes a factor of approximately 2 lower than could be studied previously. The value of the statistic (in Euclidean space) for the bursts we detect is 0.177 plus or minus 0.006. This surprisingly low value is obtained because we detected very few bursts on the 4.096 s and 8.192 s timescales (where most bursts have their highest signal-to-noise ratio) that were not already detected on the 1.024 s timescale. If allowance is made for a power-law distribution of intrinsic peak luminosities, the extended peak flux distribution is consistent with models in which the redshift distribution of the gamma-ray burst rate approximately traces the star formation history of the universe. We argue that this class of models is preferred over those in which the burst rate is independent of redshift. We use the peak flux distribution to derive a limit of 10% (99% confidence) on the fraction of the total burst rate that could be contributed by a spatially homogeneous (in Euclidean space) subpopulation of burst sources, such as type lb/c supernovae. These results lend support to the conclusions of previous studies predicting that relatively few faint "classical" GRBs will be found below the BATSE onboard detection threshold.

Kommers, Jefferson M.↗

RGS Spectroscopy of the Cygnus Loop XA Knot

The observations were performed at the end of April 2002, and the data were received in July 2002. Unfortunately, the observations were badly compromised by high levels of background radiation; one of the three observations was lost entirely. Two replacement observations were scheduled for November 2002, and were only made available in January of 2003. A preliminary analysis of the RGS data has been performed. Analysis of these data is complicated by the extended nature of the source, and the limited statistics. Nevertheless, examination of the 'm lambda' vs. 'xdsp-corr' data shows that the bright knot is visible in C VI Ly-alpha and Ly-beta, the O VII triplet and He-beta, and O VIII Ly-alpha. The dispersed C VI Ly-alpha image is the first image of the remnant at that energy. The F, I, R components of the O VII triplet can be seen, and the 0 VII He-beta is distinct from the O VIII Ly-alpha (although the dispersed diffuse emission provides a pseudo-continuum). Trial-and-error adjustment of the rgsproc parameters allowed the extraction of a narrow (approx. 1 arcmin wide) region centered on the knot; corresponding background spectra were extracted from the sky background data sets. The xspec rgsxsrc model was used to convolve the response with a Chandra spatial image. Because the spatial distribution of the emission is energy-dependent, the use of a single convolution kernel is only approximate.

Mushotzky, Richard F.↗

Gamma-Ray Observations of the Supernova Remnant RX J0852.0-4622 with the Fermi Large Area Telescope

We report on gamma-ray observations of the supernova remnant (SNR) RX J0852.04622 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. In the Fermi-LAT data, we find a spatially extended source at the location of the SNR. The extension is consistent with the SNR size seen in other wavelengths such as X-rays and TeV gamma rays, leading to the identification of the gamma-ray source with the SNR. The spectrum is well described as a power law with a photon index of = 1.85 0.06 (stat)+0.18 0.19 (sys), which smoothly connects to the H.E.S.S. spectrum in the TeV energy band. We discuss the gamma-ray emission mechanism based on multiwavelength data. The broadband data can be fit well by a model in which the gamma rays are of hadronic origin. We also consider a scenario with inverse Compton scattering of electrons as the emission mechanism of the gamma rays. Although the leptonic model predicts a harder spectrum in the Fermi-LAT energy range, the model can fit the data considering the statistical and systematic errors.

Tanaka, T.↗

Hyperspectral Mapping of the Invasive Species Pepperweed and the Development of a Habitat Suitability Model

Mapping and predicting the spatial distribution of invasive plant species is central to habitat management, however difficult to implement at landscape and regional scales. Remote sensing techniques can reduce the impact field campaigns have on these ecologically sensitive areas and can provide a regional and multi-temporal view of invasive species spread. Invasive perennial pepperweed (Lepidium latifolium) is now widespread in fragmented estuaries of the South San Francisco Bay, and is shown to degrade native vegetation in estuaries and adjacent habitats, thereby reducing forage and shelter for wildlife. The purpose of this study is to map the present distribution of pepperweed in estuarine areas of the South San Francisco Bay Salt Pond Restoration Project (Alviso, CA), and create a habitat suitability model to predict future spread. Pepperweed reflectance data were collected in-situ with a GER 1500 spectroradiometer along with 88 corresponding pepperweed presence and absence points used for building the statistical models. The spectral angle mapper (SAM) classification algorithm was used to distinguish the reflectance spectrum of pepperweed and map its distribution using an image from EO-1 Hyperion. To map pepperweed, we performed a supervised classification on an ASTER image with a resulting classification accuracy of 71.8%. We generated a weighted overlay analysis model within a geographic information system (GIS) framework to predict areas in the study site most susceptible to pepperweed colonization. Variables for the model included propensity for disturbance, status of pond restoration, proximity to water channels, and terrain curvature. A Generalized Additive Model (GAM) was also used to generate a probability map and investigate the statistical probability that each variable contributed to predict pepperweed spread. Results from the GAM revealed distance to channels, distance to ponds and curvature were statistically significant (p < 0.01) in determining the locations of suitable pepperweed habitats.

Pepperweed↗

End-to-end Modeling of Lyot Coronagraph for LUVOIR: Sensitivity to Segment-Level Aberrations

Direct imaging and spectroscopic characterization of atmospheres of Earth-like exoplanets require starlight suppression to detect a faint planetary signal at small angular separations. Future space telescopes will require segmentation to survey a statistically significant population of nearby stars systems by increasing aperture size beyond the limits set by traditional monolithic apertures; however, segmentation also represents an additional challenge due to the introduction of temporal pupil discontinuities in terms of both phase and amplitude. For this study, a baseline Lyot coronagraph is used with the Large UV / Optical /Infrared (LUVOIR) mission’s 15-m on-axis, segmented optical telescope assembly (OTA). We present a high-fidelity, end-to-end modeling approach developed as part of the Segmented Coronagraph Design and Analysis (SCDA) group that includes diffractive optical propagations at the 1e-10 contrast level across the entire optical train including the OTA and intermediate re-imaging optics inside the coronagraph instrument. We perform a tolerancing analysis to assess the contrast sensitivity to wavefront perturbations for a set of key telescope terms including segment-level low-order and mid-spatial frequency aberrations. As an input to the end-to-end telescope model, we use two industry-provided segment-level wavefront aberration datasets pertaining: (1) thermal gradients and accelerations (courtesy of Ball-Aerospace and L3-Harris), and (2) telescope dynamics including line-of-sight and segment jitter (courtesy of Lockheed-Martin). Contrast performance outputs of the coronagraph are fed into a mission design reference simulator to assess the scientific yield sensitivity. End-to-end modeling will incorporate performance improvements due to active wavefront sensing and control and/or segment-level metrology to improve system performance.

Optics↗

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

Using Clustering to Establish Climate Regimes from PCM Output

A multivariate statistical clustering technique--based on the k-means algorithm of Hartigan has been used to extract patterns of climatological significance from 200 years of general circulation model (GCM) output. Originally developed and implemented on a Beowulf-style parallel computer constructed by Hoffman and Hargrove from surplus commodity desktop PCs, the high performance parallel clustering algorithm was previously applied to the derivation of ecoregions from map stacks of 9 and 25 geophysical conditions or variables for the conterminous U.S. at a resolution of 1 sq km. Now applied both across space and through time, the clustering technique yields temporally-varying climate regimes predicted by transient runs of the Parallel Climate Model (PCM). Using a business-as-usual (BAU) scenario and clustering four fields of significance to the global water cycle (surface temperature, precipitation, soil moisture, and snow depth) from 1871 through 2098, the authors' analysis shows an increase in spatial area occupied by the cluster or climate regime which typifies desert regions (i.e., an increase in desertification) and a decrease in the spatial area occupied by the climate regime typifying winter-time high latitude perma-frost regions. The patterns of cluster changes have been analyzed to understand the predicted variability in the water cycle on global and continental scales. In addition, representative climate regimes were determined by taking three 10-year averages of the fields 100 years apart for northern hemisphere winter (December, January, and February) and summer (June, July, and August). The result is global maps of typical seasonal climate regimes for 100 years in the past, for the present, and for 100 years into the future. Using three-dimensional data or phase space representations of these climate regimes (i.e., the cluster centroids), the authors demonstrate the portion of this phase space occupied by the land surface at all points in space and time. Any single spot on the globe will exist in one of these climate regimes at any single point in time. By incrementing time, that same spot will trace out a trajectory or orbit between and among these climate regimes (or atmospheric states) in phase (or state) space. When a geographic region enters a state it never previously visited, a climatic change is said to have occurred. Tracing out the entire trajectory of a single spot on the globe yields a 'manifold' in state space representing the shape of its predicted climate occupancy. This sort of analysis enables a researcher to more easily grasp the multivariate behavior of the climate system.

Oglesby, Robert↗

Regression Bias in using Solar Wind Measurements

Simultaneous solar wind measurements from the solar wind monitors, WIND and ACE, differ due to the spatial and temporal structure of the solar wind. Correlation studies that use these measurements as input may infer an incorrect correlation due to uncertainties arising from this spatial and temporal structure, especially at extreme and rare solar wind values. In particular, regression analysis will lead to a regression function whose slope is biased towards the mean value of the measurement parameter. This article demonstrates this regression bias by comparing simultaneous ACE and WIND solar wind measurements. A non-linear regression analysis between them leads to a perception of underestimation of extreme values of one measurement on average over the other. Using numerical experiments, we show that popular regression analysis techniques such as linear least-squares, orthogonal least-squares, and non-linear regression are not immune to this bias. Hence while using solar wind parameters as an independent variable in a correlation or regression analysis, random uncertainty in the independent variable can create unintended biases in the response of the dependent variable. More generally, the regression to the mean effect can impact both event-based, statistical studies of magnetospheric response to solar wind forcing.

magnetosphere↗

Monitoring Forest Regrowth Using a Multi-Platform Time Series

Over the past 50 years, the forests of western Washington and Oregon have been extensively harvested for timber. This has resulted in a heterogeneous mosaic of remaining mature forests, clear-cuts, new plantations, and second-growth stands that now occur in areas that formerly were dominated by extensive old-growth forests and younger forests resulting from fire disturbance. Traditionally, determination of seral stage and stand condition have been made using aerial photography and spot field observations, a methodology that is not only time- and resource-intensive, but falls short of providing current information on a regional scale. These limitations may be solved, in part, through the use of multispectral images which can cover large areas at spatial resolutions in the order of tens of meters. The use of multiple images comprising a time series potentially can be used to monitor land use (e.g. cutting and replanting), and to observe natural processes such as regeneration, maturation and phenologic change. These processes are more likely to be spectrally observed in a time series composed of images taken during different seasons over a long period of time. Therefore, for many areas, it may be necessary to use a variety of images taken with different imaging systems. A common framework for interpretation is needed that reduces topographic, atmospheric, instrumental, effects as well as differences in lighting geometry between images. The present state of remote-sensing technology in general use does not realize the full potential of the multispectral data in areas of high topographic relief. For example, the primary method for analyzing images of forested landscapes in the Northwest has been with statistical classifiers (e.g. parallelepiped, nearest-neighbor, maximum likelihood, etc.), often applied to uncalibrated multispectral data. Although this approach has produced useful information from individual images in some areas, landcover classes defined by these techniques typically are not consistent for the same scene imaged under different illumination conditions, especially in the mountainous regions. In addition, it is difficult to correct for atmospheric and instrumental differences between multiple scenes in a time series. In this paper, we present an approach for monitoring forest cutting/regrowth in a semi-mountainous portion of the southern Gifford Pinchot National Forest using a multisensor-time series composed of MSS, TM, and AVIRIS images.

Sabol, Donald E., Jr.↗

The Impact of Subsampling on MODIS Level-3 Statistics of Cloud Optical Thickness and Effective Radius

The MODIS Level-3 optical thickness and effective radius cloud product is a gridded l deg. x 1 deg. dataset that is derived from aggregation and subsampling at 5 km of 1 km, resolution Level-2 orbital swath data (Level-2 granules). This study examines the impact of the 5 km subsampling on the mean, standard deviation and inhomogeneity parameter statistics of optical thickness and effective radius. The methodology is simple and consists of estimating mean errors for a large collection of Terra and Aqua Level-2 granules by taking the difference of the statistics at the original and subsampled resolutions. It is shown that the Level-3 sampling does not affect the various quantities investigated to the same degree, with second order moments suffering greater subsampling errors, as expected. Mean errors drop dramatically when averages over a sufficient number of regions (e.g., monthly and/or latitudinal averages) are taken, pointing to a dominance of errors that are of random nature. When histograms built from subsampled data with the same binning rules as in the Level-3 dataset are used to reconstruct the quantities of interest, the mean errors do not deteriorate significantly. The results in this paper provide guidance to users of MODIS Level-3 optical thickness and effective radius cloud products on the range of errors due to subsampling they should expect and perhaps account for, in scientific work with this dataset. In general, subsampling errors should not be a serious concern when moderate temporal and/or spatial averaging is performed.

Oreopoulos, Lazaros↗

Validation of the AMSU-B Bias Corrections Based on Satellite Measurements from SSM/T-2

The NOAA-15 Advanced Microwave Sounding Unit-B (AMSU-B) was designed in the same spirit as the Special Sensor Microwave Water Vapor Profiler (SSM/T-2) on board the DMSP F11-14 satellites, to perform remote sensing of spatial and temporal variations in mid and upper troposphere humidity. While the SSM/T-2 instruments have a 48 km spatial resolution at nadir and 28 beam positions per scan, AMSU-B provides an improvement with a 16 km spatial resolution at nadir and 90 beam positions per scan. The AMSU-B instrument, though, has been experiencing radio frequency interference (RFI) contamination from the NOAA-15 transmitters whose effect is dependent upon channel, geographic location, and current spacecraft antenna configuration. This has lead to large cross-track biases reaching as high as 100 Kelvin for channel 17 (150 GHz) and 50 Kelvin for channel 19 (183 +/-3 GHz). NOAA-NESDIS has recently provided a series of bias corrections for AMSU-B data starting from March, 1999. These corrections are available for each of the five channels, for every third field of view, and for three cycles within an eight second period. There is also a quality indicator in each data record to indicate whether or not the bias corrections should be applied. As a precursor to performing retrievals of mid and upper troposphere humidity, a validation study is performed by statistically analyzing the differences between the F14 SSM/T-2 and the bias corrected AMSU-B brightness temperatures for three months in the spring of 1999.

Kolodner, Marc A.↗

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↗

A Method for Obtaining High Time and Spatial Resolution Convective Cloud Top Data for the TTL

A method for obtaining high time and spatial resolution convective cloud top data for the TTL Leonhard Pfister, Eric Jensen, Rei Ueyama, Eliot Atlas, and Maria Navarro Convective systems in the tropics have a maximum in the cloud top altitude distribution of about 13.5 km. However, there is a significant tail to this distribution -- a few percent reach the cold point tropopause (CPT) at 16.5 km, and there has been clear evidence of convective mass deposited as high as 19 km in the tropics. The region between 13.5 km and the cold point tropopause is transitional, between the free tropical troposphere where convective mixing dominates, and the stratosphere where slow upward ascent dominates. In this region (the Tropical Tropopause Layer), convective injection, slow ascent, and mixing from midlatitudes all have similar time scales. So, even though only a few percent of convective systems reach the CPT, convection is important. Space Based Lidar and cloud radar measurements have yielded information about long term average statistical distributions of cloud altitude as a function of location. However, we also need time-dependent cloud top altitude and cloud top potential temperature information, primarily to understand the water vapor and TTL cloud distributions. This is because the effect of convection depends on the local temperature, and on the subsequent temperature history. Time dependent cloud top information is also needed to understand short-lived tracers because cross-isentropic flow is time and space dependent. This paper presents a method of obtaining time and space dependent convective cloud top theta (and altitude) information using 3-hourly geostationary brightness temperature data, coupled with global 3 -hourly rainfall estimates and temperature analyses. We explore different mixing algorithms to obtain the most reasonable agreement with near-simultaneous observations by cloudsat and calipso. Observations of short-lived tracers from ATTREX, coupled with short-term trajectories are used to test the method's accuracy. An important caveat is the ambiguity of evaluating convective cloud top altitudes under from combined cloudsat and calipso measurements.

cold point tropopause (CPT)↗

The Role of Spatial Disorientation in Fatal General Aviation Accidents

In-flight Spatial Disorientation (SD) in pilots is a serious threat to aviation safety. Indeed, SD may play a much larger role in aviation accidents than the approximate 6-8% reported by the National Transportation Safety Board (NTSB) each year, because some accidents coded by the NTSB as aircraft control-not maintained (ACNM) may actually result from SD. The purpose of this study is to determine whether SD is underestimated as a cause of fatal general aviation (GA) accidents in the NTSB database. Fatal GA airplane accidents occurring between January 1995 and December 1999 were reviewed from the NTSB aviation accident database. Cases coded as ACNM or SD as the probable cause were selected for review by a panel of aerospace medicine specialists. Using a rating scale, each rater was instructed to determine if SD was the probable cause of the accident. Agreement between the raters and agreement between the raters and the NTSB were evaluated by Kappa statistics. The raters agreed that 11 out of 20 (55%) accidents coded by the NTSB as ACNM were probably caused by SD (p less than 0.05). Agreement between the raters and the NTSB did not reach significance (p greater than 0.05). The 95% C.I. for the sampling population estimated that between 33-77% of cases that the NTSB identified as ACNM could be identified by aerospace medicine experts as SD. Aerospace medicine specialists agreed that some cases coded by the NTSB as ACNM were probably caused by SD. Consequently, a larger number of accidents may be caused by the pilot succumbing to SD than indicated in the NTSB database. This new information should encourage regulating agencies to insure that pilots receive SD recognition training, enabling them to take appropriate corrective actions during flight. This could lead to new training standards, ultimately saving lives among GA airplane pilots.

Scheuring, RIchard↗

The Collection 6 'dark-target' MODIS Aerosol Products

Aerosol retrieval algorithms are applied to Moderate resolution Imaging Spectroradiometer (MODIS) sensors on both Terra and Aqua, creating two streams of decade-plus aerosol information. Products of aerosol optical depth (AOD) and aerosol size are used for many applications, but the primary concern is that these global products are comprehensive and consistent enough for use in climate studies. One of our major customers is the international modeling comparison study known as AEROCOM, which relies on the MODIS data as a benchmark. In order to keep up with the needs of AEROCOM and other MODIS data users, while utilizing new science and tools, we have improved the algorithms and products. The code, and the associated products, will be known as Collection 6 (C6). While not a major overhaul from the previous Collection 5 (C5) version, there are enough changes that there are significant impacts to the products and their interpretation. In its entirety, the C6 algorithm is comprised of three sub-algorithms for retrieving aerosol properties over different surfaces: These include the dark-target DT algorithms to retrieve over (1) ocean and (2) vegetated-dark-soiled land, plus the (3) Deep Blue (DB) algorithm, originally developed to retrieve over desert-arid land. Focusing on the two DT algorithms, we have updated assumptions for central wavelengths, Rayleigh optical depths and gas (H2O, O3, CO2, etc.) absorption corrections, while relaxing the solar zenith angle limit (up to 84) to increase pole-ward coverage. For DT-land, we have updated the cloud mask to allow heavy smoke retrievals, fine-tuned the assignments for aerosol type as function of season location, corrected bugs in the Quality Assurance (QA) logic, and added diagnostic parameters such as topographic altitude. For DT-ocean, improvements include a revised cloud mask for thin-cirrus detection, inclusion of wind speed dependence in the retrieval, updates to logic of QA Confidence flag (QAC) assignment, and additions of important diagnostic information. At the same time as we have introduced algorithm changes, we have also accounted for upstream changes including: new instrument calibration, revised land-sea masking, and changed cloud masking. Upstream changes also impact the coverage and global statistics of the retrieved AOD. Although our responsibility is to the DT code and products, we have also added a product that merges DT and DB product over semi-arid land surfaces to provide a more gap-free dataset, primarily for visualization purposes. Preliminary validation shows that compared to surface-based sunphotometer data, the C6, Level 2 (along swath) DT-products compare at least as well as those from C5. C6 will include new diagnostic information about clouds in the aerosol field, including an aerosol cloud mask at 500 m resolution, and calculations of the distance to the nearest cloud from clear pixels. Finally, we have revised the strategy for aggregating and averaging the Level 2 (swath) data to become Level 3 (gridded) data. All together, the changes to the DT algorithms will result in reduced global AOD (by 0.02) over ocean and increased AOD (by 0.02) over land, along with changes in spatial coverage. Changes in calibration will have more impact to Terras time series, especially over land. This will result in a significant reduction in artificial differences in the Terra and Aqua datasets, and will stabilize the MODIS data as a target for AEROCOM studie

Aerosol retrieval algorithms↗

Advanced Water Vapor Lidar Detection System

In the present water vapor lidar system, the detected signal is sent over long cables to a waveform digitizer in a CAMAC crate. This has the disadvantage of transmitting analog signals for a relatively long distance, which is subjected to pickup noise, leading to a decrease in the signal to noise ratio. Generally, errors in the measurement of water vapor with the DIAL method arise from both random and systematic sources. Systematic errors in DIAL measurements are caused by both atmospheric and instrumentation effects. The selection of the on-line alexandrite laser with a narrow linewidth, suitable intensity and high spectral purity, and its operation at the center of the water vapor lines, ensures minimum influence in the DIAL measurement that are caused by the laser spectral distribution and avoid system overloads. Random errors are caused by noise in the detected signal. Variability of the photon statistics in the lidar return signal, noise resulting from detector dark current, and noise in the background signal are the main sources of random error. This type of error can be minimized by maximizing the signal to noise ratio. The increase in the signal to noise ratio can be achieved by several ways. One way is to increase the laser pulse energy, by increasing its amplitude or the pulse repetition rate. Another way, is to use a detector system with higher quantum efficiency and lower noise, on the other hand, the selection of a narrow band optical filter that rejects most of the day background light and retains high optical efficiency is an important issue. Following acquisition of the lidar data, we minimize random errors in the DIAL measurement by averaging the data, but this will result in the reduction of the vertical and horizontal resolutions. Thus, a trade off is necessary to achieve a balance between the spatial resolution and the measurement precision. Therefore, the main goal of this research effort is to increase the signal to noise ratio by a factor of 10 over the current system, using a newly evaluated, very low noise avalanche photo diode detector and constructing a 10 MHz waveform digitizer which will replace the current CAMAC system.

Elsayed-Ali, Hani↗

Bulk Compositional Trends in Meteorites: A Guide for Analysis and Interpretation of NEAR XGRS Data from Asteroid 433 EROS

The Near Earth Asteroid Rendezvous (NEAR) spacecraft is to orbit the S-class asteroid 433 Eros for about one year beginning on February 14, 2000. The X-ray/gamma-ray O, Mg, Al, Si, Fe, and K; possibly H, Ca, S, Ti, and Th) of Eros with a spatial resolution ranging from a few km for X-rays to approx. 25% of the asteroid's surface for gamma-rays. The major scientific goals for the NEAR XGRS are to relate the composition of Eros to known classes of meteorites, to assess compositional heterogeneity and to identify geological processes that have occurred on the asteroid. Comparing remote-sensing data from asteroids to laboratory data from meteorites requires that the latter be well determined and understood. How well particular classes of meteorites can be identified as analogues of Eros depends not only on the error of the XGRS measurement, but also on the spread in abundances observed among different members of a given meteorite class. To prepare for the return of XGRS data from Eros, we have compiled a large database of bulk elemental compositions of meteorites, using data from a wide variety of published and unpublished sources. Custom software was developed to easily extract statistical information and make plots of data from different meteorite classes. Here, we use the meteorite compositional database to investigate which abundances and abundance ratios, of those measureable by the NEAR XGRS, are most diagnostic for distinguishing meteorite classes and identifying geological processes that have occurred on the samples' parent asteroids.

Nittler, Larry R.↗

Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology

Changing trends in ecosystem productivity can be quantified using satellite observations of Normalized Difference Vegetation Index (NDVI). However, the estimation of trends from NDVI time series differs substantially depending on analyzed satellite dataset, the corresponding spatiotemporal resolution, and the applied statistical method. Here we compare the performance of a wide range of trend estimation methods and demonstrate that performance decreases with increasing inter-annual variability in the NDVI time series. Trend slope estimates based on annual aggregated time series or based on a seasonal-trend model show better performances than methods that remove the seasonal cycle of the time series. A breakpoint detection analysis reveals that an overestimation of breakpoints in NDVI trends can result in wrong or even opposite trend estimates. Based on our results, we give practical recommendations for the application of trend methods on long-term NDVI time series. Particularly, we apply and compare different methods on NDVI time series in Alaska, where both greening and browning trends have been previously observed. Here, the multi-method uncertainty of NDVI trends is quantified through the application of the different trend estimation methods. Our results indicate that greening NDVI trends in Alaska are more spatially and temporally prevalent than browning trends. We also show that detected breakpoints in NDVI trends tend to coincide with large fires. Overall, our analyses demonstrate that seasonal trend methods need to be improved against inter-annual variability to quantify changing trends in ecosystem productivity with higher accuracy.

Inter-Annual↗