Search NASA⌕ Search

SEARCH · Search NASA

Results for “Principal Component 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 163 records · Page 9

The Global Oscillation Network Group site survey. 1: Data collection and analysis methods

The Global Oscillation Network Group (GONG) Project is planning to place a set of instruments around the world to observe solar oscillations as continuously as possible for at least three years. The Project has now chosen the sites that will comprise the network. This paper describes the methods of data collection and analysis that were used to make this decision. Solar irradiance data were collected with a one-minute cadence at fifteen sites around the world and analyzed to produce statistics of cloud cover, atmospheric extinction, and transparency power spectra at the individual sites. Nearly 200 reasonable six-site networks were assembled from the individual stations, and a set of statistical measures of the performance of the networks was analyzed using a principal component analysis. An accompanying paper presents the results of the survey.

Hill, Frank↗

NASA Tech Briefs, January 2005

Topics covered include: Fiber-Optic Sensor Would Monitor Growth of Polymer Film; Sensors for Pointing Moving Instruments Toward Each Other; Pd/CeO2/SiC Chemical Sensors; Microparticle Flow Sensor; Scattering-Type Surface-Plasmon-Resonance Biosensors; Diode-Laser-Based Spectrometer for Sensing Gases; Improved Cathode Structure for a Direct Methanol Fuel Cell; X-Band, 17-Watt Solid-State Power Amplifier; Improved Anode for a Direct Methanol Fuel Cell; Tools for Designing and Analyzing Structures; Interactive Display of Scenes with Annotations; Solving Common Mathematical Problems; Tools for Basic Statistical Analysis; Program Calculates Forces in Bolted Structural Joints; Integrated Structural Analysis and Test Program; Molybdate Coatings for Protecting Aluminum Against Corrosion; Synthesizing Diamond from Liquid Feedstock; Modifying Silicates for Better Dispersion in Nanocomposites; Powder-Collection System for Ultrasonic/Sonic Drill/Corer; Semiautomated, Reproducible Batch Processing of Soy; Hydrogen Peroxide Enhances Removal of NOx from Flue Gases; Subsurface Ice Probe; Real-Time Simulation of Aeroheating of the Hyper-X Airplane; Using Laser-Induced Incandescence To Measure Soot in Exhaust; Method of Real-Time Principal-Component Analysis; Insect-Inspired Flight Control for Unmanned Aerial Vehicles; Domain Compilation for Embedded Real-Time Planning; Semantic Metrics for Analysis of Software; Simulation of Laser Cooling and Trapping in Engineering Applications; Large Fluvial Fans and Exploration for Hydrocarbons; Doping-Induced Interband Gain in InAs/AlSb Quantum Wells; Development of Software for a Lidar-Altimeter Processor; Upgrading the Space Shuttle Caution and Warning System; and Fractal Reference Signals in Pulse-Width Modulation.

Source record↗

Real Time Intelligent Target Detection and Analysis with Machine Vision

We present an algorithm for detecting a specified set of targets for an Automatic Target Recognition (ATR) application. ATR involves processing images for detecting, classifying, and tracking targets embedded in a background scene. We address the problem of discriminating between targets and nontarget objects in a scene by evaluating 40x40 image blocks belonging to an image. Each image block is first projected onto a set of templates specifically designed to separate images of targets embedded in a typical background scene from those background images without targets. These filters are found using directed principal component analysis which maximally separates the two groups. The projected images are then clustered into one of n classes based on a minimum distance to a set of n cluster prototypes. These cluster prototypes have previously been identified using a modified clustering algorithm based on prior sensed data. Each projected image pattern is then fed into the associated cluster's trained neural network for classification. A detailed description of our algorithm will be given in this paper. We outline our methodology for designing the templates, describe our modified clustering algorithm, and provide details on the neural network classifiers. Evaluation of the overall algorithm demonstrates that our detection rates approach 96% with a false positive rate of less than 0.03%.

Howard, Ayanna↗

The Global Oscillation Network Group site survey, 2: Results

The Global Oscillation Network Group (GONG) Project will place a network of instruments around the world to observe solar oscillations as continuously as possible for three years. The Project has now chosen the six network sites based on analysis of survey data from fifteen sites around the world. The chosen sites are: Big Bear Solar Observatory, California; Mauna Loa Solar Observatory, Hawaii; Learmonth Solar Observatory, Australia; Udaipur Solar Observatory, India; Observatorio del Teide, Tenerife; and Cerro Tololo Interamerican Observatory, Chile. Total solar intensity at each site yields information on local cloud cover, extinction coefficient, and transparency fluctuations. In addition, the performance of 192 reasonable networks assembled from the individual site records is compared using a statistical principal components analysis. An accompanying paper descibes the analysis methods in detail; here we present the results of both the network and individual site analyses. The selected network has a duty cycle of 93.3%, in good agreement with numerical simulations. The power spectrum of the network observing window shows a first diurnal sidelobe height of 3 x 10(exp -4) with respect to the central component, an improvement of a factor of 1300 over a single site. The background level of the network spectrum is lower by a factor of 50 compared to a single-site spectrum.

Hill, Frank↗

Anomaly Detection for the Roman Space Telescope Wide Field Instrument’s Science Data Processing Pipeline

The Roman Space Telescope (RST) Wide Field Instrument (WFI) will be utilizing a preliminary Science Data Processing (SDP) pipeline during its Integration and Test, and to some extent during Operations, to track basic statistics and identify known features such as cosmic rays, snowballs as well as possible anomalies in raw detector data. In our detectors, these anomalies appear as jumps in the ramp of a readout and are classified as cosmic rays if they appear as a streak or snowballs if they’re more circular. The WFI employs an array of 18 H4RG-10 detectors that collect image samples. Each set of raw frames within a non-destructive exposure is packaged by the SDP pipeline into image cubes for each detector. Each cube is a time series of 4096 × 4096 accumulating pixel frames. The preliminary analysis pipeline is used to locate anomalies in these time-series accumulation frames and identify the type of anomaly, either natural phenomena or detector characteristic. To compare different methods, we’ve implemented both heuristic-based and data-driven methods to identify anomalies. For the heuristic-based approach, we identify snowballs and cosmic rays by the size and shape of outlier pixel clusters between consecutive frames. For data driven methods, we evaluated a Convolutional Neural Network (CNN) model, and more traditional methods like Principal Component Analysis (PCA). CNN is a supervised learning/classification method. Thus, we used a labeled dataset of anomalies to perform segmentation of the image and identify anomalies. We used previously identified cosmic rays and snowballs to measure the accuracy and efficiency of the mentioned approaches. In evaluating these methods, we aim to pick the best fit for the SDP pipeline’s anomaly detection in terms of both performance and runtime.

Paul Horton↗

Statistical analysis of astronomical data containing upper bounds - General methods and examples drawn from X-ray astronomy

Statistical procedures taken from the field of survival analysis have been adapted to astronomical usage and have been applied to a sample of stars in the B-V color range between 0.1 and 0.5 with measured soft X-ray luminosities and projected equatorial velocities. The two-sample problem and linear regression problem with arbitrarily censored data were studied. A new method for determining the linear regression between two random variables in the presence of arbitrary censoring has been developed which can also be used for a likelihood-ratio test for the independence of two random variables and for principal-component analysis in the presence of arbitrary censoring. The required numerical computations can be carried out straightforwardly and rapidly.

Schmitt, J. H. M. M.↗

Visualization of Global Sensitivity Analysis Results Based on a Combination of Linearly Dependent and Independent Directions

A useful technique for the validation and verification of complex flight systems is Monte Carlo Filtering -- a global sensitivity analysis that tries to find the inputs and ranges that are most likely to lead to a subset of the outputs. A thorough exploration of the parameter space for complex integrated systems may require thousands of experiments and hundreds of controlled and measured variables. Tools for analyzing this space often have limitations caused by the numerical problems associated with high dimensionality and caused by the assumption of independence of all of the dimensions. To combat both of these limitations, we propose a technique that uses a combination of the original variables with the derived variables obtained during a principal component analysis.

Davies, Misty D.↗

IMERG and GPCP Seasonality and Response to Climate and Weather Variability

The Integrated Multi-satellitE Retrievals for GPM (IMERG) and the The Global Precipitation Climatology Project (GPCP) are two of the most popular precipitation products. IMERG is a relatively new dataset that targets the needs primarily of the hydrological community by resolving the hydrological cycle of precipitation at fine temporal (30-minutes) and spatial (10-km) scales. IMERG only recently exceeded the user base of the highly successful, but discontinued in 2019, TRMM Multi-satellite Precipitation Analysis (TMPA). GPCP, on the other hand, has traditionally been strong in the climate research community, and recently has been revised under the framework of NASA's Making Earth System Data records for Use in Research Environments (MEaSUREs) program. Both IMERG and GPCP are similar in the underlying approaches to achieve global coverage, in particular using satellite microwave and infrared observations, and adjusting the precipitation retrieval with rain gauge information. While there is a tendency to use both datasets interchangeably, differences remain and some of them limit the usage of IMERG as a climate data record at this point. By applying Principal Component Analysis, we identify the differences and similarities mode-by-mode, and further the guidance on suggested usage of IMERG as a climate data record. As an example in the attached figure, IMERG vs GPSP differences in the explained variance by the two leading seasonal modes can be identified, mainly in the extreme southern latitudes, and over western boundary currents (Gulfstream and Kuroshio). At the preparation time for this presentation, the new version "07"of IMERG was in the works that may resolve the issues presented here. Nevertheless, our analysis can help to gauge the uncertainties of the studies already done using the currently existing IMERG version "06", and evaluate the improvements in the upcoming version "07".

Andrey Savtchenko↗

Incorporating Endmember Variability into Spectral Mixture Analysis Through Endmember Bundles

Variation in canopy structure and biochemistry induces a concomitant variation in the top-of-canopy spectral reflectance of a vegetation type. Hence, the use of a single endmember spectrum to track the fractional abundance of a given vegetation cover in a hyperspectral image may result in fractions with considerable error. One solution to the problem of endmember variability is to increase the number of endmembers used in a spectral mixture analysis of the image. For example, there could be several tree endmembers in the analysis because of differences in leaf area index (LAI) and multiple scatterings between leaves and stems. However, it is often difficult in terms of computer or human interaction time to select more than six or seven endmembers and any non-removable noise, as well as the number of uncorrelated bands in the image, limits the number of endmembers that can be discriminated. Moreover, as endmembers proliferate, their interpretation becomes increasingly difficult and often applications simply need the aerial fractions of a few land cover components which comprise most of the scene. In order to incorporate endmember variability into spectral mixture analysis, we propose representing a landscape component type not with one endmember spectrum but with a set or bundle of spectra, each of which is feasible as the spectrum of an instance of the component (e.g., in the case of a tree component, each spectrum could reasonably be the spectral reflectance of a tree canopy). These endmember bundles can be used with nonlinear optimization algorithms to find upper and lower bounds on endmember fractions. This approach to endmember variability naturally evolved from previous work in deriving endmembers from the data itself by fitting a triangle, tetrahedron or, more generally, a simplex to the data cloud reduced in dimension by a principal component analysis. Conceptually, endmember variability could make it difficult to find a simplex that both surrounds the data cloud and has vertices that are realistic endmember spectra with reflectances between 0 and 1. In this paper, we create endmember bundles and bounding fraction images for an AVIRIS subscene simulated with a plant canopy radiative transfer model. The simulated subscene is spatially patterned after a subscene from the AVIRIS image acquired August, 1993 over La Copita, Texas. In addition, for comparison, we performed a traditional unmixing with image endmembers.

Bateson, C. Ann↗

An Initial Analysis of LANDSAT-4 Thematic Mapper Data for the Discrimination of Agricultural, Forested Wetlands, and Urban Land Cover

The capabilities of TM data for discriminating land covers within three particular cultural and ecological realms was assessed. The agricultural investigation in Poinsett County, Arkansas illustrates that TM data can successfully be used to discriminate a variety of crop cover types within the study area. The single-date TM classification produced results that were significantly better than those developed from multitemporal MSS data. For the Reelfoot Lake area of Tennessee TM data, processed using unsupervised signature development techniques, produced a detailed classification of forested wetlands with excellent accuracy. Even in a small city of approximately 15,000 people (Union City, Tennessee). TM data can successfully be used to spectrally distinguish specific urban classes. Furthermore, the principal components analysis evaluation of the data shows that through photointerpretation, it is possible to distinguish individual buildings and roof responses with the TM.

Quattrochi, D. A.↗

Study of LANDSAT-D thematic mapper performance as applied to hydrocarbon exploration

Analysis of the tapes of the Detroit, Michigan scene, which were received in fully processed format with geometric and radiometric correction, shows evidence of an along line data slip every sixteenth line in TM channel 2. Very large scale products were therefore generated in false color using channels 1, 3, and 4. Subjective evaluation of these enhanced scenes indicates that they are acceptable for interpretation at scales up to 1:50,000 and should be useful for change mapping probably up to 1:24,000 scale. The significant striping visible in water bodies for both the natural color and false color products indicates that the detector calibration is probably performing below the preflight specification. Variance-covariance matrices were computed and principal component analysis were performed for a set of 512 x 512 windows within the Arkansas scene. Initial analysis shows the shortwave infrared channels (TM 5 and 6) are a highly significant data source. The thermal channel (TM 7) shows negative correlation with TM 1 through 4.

Source record↗

An initial analysis of LANDSAT 4 Thematic Mapper data for the classification of agricultural, forested wetland, and urban land covers

An initial analysis of LANDSAT 4 thematic mapper (TM) data for the delineation and classification of agricultural, forested wetland, and urban land covers was conducted. A study area in Poinsett County, Arkansas was used to evaluate a classification of agricultural lands derived from multitemporal LANDSAT multispectral scanner (MSS) data in comparison with a classification of TM data for the same area. Data over Reelfoot Lake in northwestern Tennessee were utilized to evaluate the TM for delineating forested wetland species. A classification of the study area was assessed for accuracy in discriminating five forested wetland categories. Finally, the TM data were used to identify urban features within a small city. A computer generated classification of Union City, Tennessee was analyzed for accuracy in delineating urban land covers. An evaluation of digitally enhanced TM data using principal components analysis to facilitate photointerpretation of urban features was also performed.

Quattrochi, D. A.↗

Extension and Statistical Analysis of the GACP Aerosol Optical Thickness Record.

The primary product of the Global Aerosol Climatology Project (GACP) is a continuous record of the aerosol optical thickness (AOT) over the oceans. It is based on channel-1 and -2 radiance data from the Advanced Very High Resolution Radiometer (AVHRR) instruments flown on successive National Oceanic and Atmospheric Administration (NOAA) platforms. We extend the previous GACP dataset by four years through the end of 2009 using NOAA-17 and -18 AVHRR radiances recalibrated against MODerate resolution Imaging Spectroradiometer (MODIS) radiance data, thereby making the GACP record almost three decades long. The temporal overlap of over three years of the new NOAA-17 and the previous NOAA-16 record reveals an excellent agreement of the corresponding global monthly mean AOT values, thereby confirming the robustness of the vicarious radiance calibration used in the original GACP product. The temporal overlap of the NOAA-17 and -18 instruments is used to introduce a small additive adjustment to the channel-2 calibration of the latter resulting in a consistent record with increased data density. The Principal Component Analysis (PCA) of the newly extended GACP record shows that most of the volcanic AOT variability can be isolated into one mode responsible for ~12% of the total variance. This conclusion is confirmed by a combined PCA analysis of the GACP, MODIS, andMulti-angle Imaging SpectroRadiometer (MISR) AOTs during the volcano-free period fromFebruary 2000 to December 2009.We show that the modes responsible for the tropospheric AOT variability in the three datasets agree well in terms of correlation and spatial patterns. A previously identified negative AOT trend which started in the late 1980s and continued into the early 2000s is confirmed. Its magnitude and duration indicate that it was caused by changes in tropospheric aerosols. The latest multi-satellite segment of the GACP record shows that this trend tapered off, with no noticeable AOT change after 2002. This result is consistent with the MODIS andMISR AOT records as well as with the recent gradual reversal frombrightening to dimming revealed by surface flux measurements in many aerosol producing regions. Thus the robustness of the GACP record is confirmed, increasing our confidence in the validity of the negative trend. Although the nominal negative GACP AOT trend could partially be an artifact of increasing aerosol absorption, we argue that the time dependence of the GACP record, including the latest flat period, is more consistent with the actual decrease in the tropospheric AOT.

aerosols↗

Neural network uncertainty assessment using Bayesian statistics: a remote sensing application

Neural network (NN) techniques have proved successful for many regression problems, in particular for remote sensing; however, uncertainty estimates are rarely provided. In this article, a Bayesian technique to evaluate uncertainties of the NN parameters (i.e., synaptic weights) is first presented. In contrast to more traditional approaches based on point estimation of the NN weights, we assess uncertainties on such estimates to monitor the robustness of the NN model. These theoretical developments are illustrated by applying them to the problem of retrieving surface skin temperature, microwave surface emissivities, and integrated water vapor content from a combined analysis of satellite microwave and infrared observations over land. The weight uncertainty estimates are then used to compute analytically the uncertainties in the network outputs (i.e., error bars and correlation structure of these errors). Such quantities are very important for evaluating any application of an NN model. The uncertainties on the NN Jacobians are then considered in the third part of this article. Used for regression fitting, NN models can be used effectively to represent highly nonlinear, multivariate functions. In this situation, most emphasis is put on estimating the output errors, but almost no attention has been given to errors associated with the internal structure of the regression model. The complex structure of dependency inside the NN is the essence of the model, and assessing its quality, coherency, and physical character makes all the difference between a blackbox model with small output errors and a reliable, robust, and physically coherent model. Such dependency structures are described to the first order by the NN Jacobians: they indicate the sensitivity of one output with respect to the inputs of the model for given input data. We use a Monte Carlo integration procedure to estimate the robustness of the NN Jacobians. A regularization strategy based on principal component analysis is proposed to suppress the multicollinearities in order to make these Jacobians robust and physically meaningful.

Neural Networks (Computer)↗

Passive microwave observations of the Wedell Sea during austral winter and early spring

The results of multispectral passive microwave observations (6.7 to 90-GHz) are presented from the cruises of the FS Polarstern in the Weddell Sea from July to December 1986. This paper includes primarily the analysis of radiometric observations taken at ice station sites. Averaged emissivity spectra for first-year (FY) ice were relatively constant throughout the experiment and were not statistically different from FY ice signatures in the Arctic. Detailed ice characterization was carried out at each site to compare the microwave signatures of the ice with the physical properties. Absorption optical depths of FY ice were found to be sufficiently high that only the structure in the upper portions of the ice contributed significantly to interstation emissivity variations. The emissivities at 90-GHz, e(90), had the greatest variance. Both e(90) at vertical polarization and GR(sub e)(90, 18.7)(defined as (e(sub V)(90)-e(sub V)(18.7))/e(sub V)(90 + e(sub V)(18.7)) depended on the scattering optical depth which is a function of the snow grain diameter and layer thickness. The variance showed a latitude dependence and is probably due to an increase in the strength of snow metamorphism nearer the northern edge of the ice pack. The contribution of variations of near-surface brine volume to the emissivity was not significant over the range of values encountered at the station sites. Emissivity spectra are presented for a range of thin ice types. Unsupervised principal component analysis produced three significant eigenvectors and showed a separation among four different surface types: open water, thin ice, FY ice, and FY ice with a thick snow cover. A comparison with SMMR satellite data showed that average ice concentrations derived from the ship's ice watch log were consistent with the satellite concentrations. The surface based emissivities for FY ice were also compared with emissivities calculated from scanning multichannel microwave radiometer (SMMR) satellite radiances. Best agreement was found at 6.7 and 10-GHz, while at 18 and 37-GHz, SMMR emissivities were slightly lower than surface based results. For the three lower frequencies agreement was found within a confidence limit of 95% and for 37-GHz within about 90%.

Grenfell, T. C.↗

Predicting the Seawater Chemistry of an Ocean World Using Machine Learning on Isotopic Measurements of Volatile CO2

Introduction: Given the long time intervals required for data transmission to and from ocean worlds targets, low bandwidth for data transmission, time required for data processing and analysis, and potentially extreme radiation environments (e.g., Europa), it is clear that ocean worlds missions will need more autonomous flight instruments and software in order to achieve established science goals. Protracted time intervals for data analysis (e.g., Europa Lander) strongly motivates the development of rapid, consistent and streamlined methods for interpreting data from flight mass spectrometers to e.g., determine how mass spectra from a plume or surface liquid/ice relates to the surface/subsurface. Since mass spectrometry also has the potential to correctly identify biosignatures[1], it is imperative that such methods for interpreting data are consistent and accurate. We used 848 isotope ratio mass spectra from laboratory analyses of CO2 that interacted with ocean worlds-relevant seawaters as a ‘training’ dataset for ‘unsupervised’ machine learning. In unsupervised learning, characteristics of the data are not labeled or linked, and any similarities found only result from the neural network. CO2 isotopologues analyzed for this dataset mimic the remote measurements of CO2 by a flight mass spectrometer, and are detailed in Theiling [2]. From this dataset, we used measured features of the spectra, such as retention time, intensity, and (isotopologue) mass ratios as inputs for our autoencoder neural network. Our neural network was trained to find similarities in these and other spectral features for seawaters of a particular composition and amount of initial CO2. Successful training then created an output of these similarities for various seawaters, which included MgSO4, Na2SO4, NaCl, MgCl2, KCl, and NaHCO3, and combinations of these salts. We then applied dimensionality reduction techniques such as Principal Component Analysis (PCA), T-Distributed Stochastic Neighbor Embedding (TSNE), and Uniform Manifold Approximation and Projection (UMAP) to demonstrate latent data features as a two-dimensional projection in a unitless, high-dimensional space. In this projection, a data point represents the combined effect of spectral features such as intensity, retention time, and isotope ratio. Our initial UMAP demonstrates data clustering (organization of the data by the neural network) based on the amount of CO2 that had initially interacted with each seawater. Further training using more ‘supervised’ learning techniques demonstrate strong clustering of preliminary data based on initial CO2 concentration, seawater chemical composition, and ionic strength (salinity). Our preliminary work therefore suggests that machine learning has the potential to identify compositional variants of an ocean world seawater based on mass spectra from volatile CO2 measurements. Acknowledgments: This work was funded through a Strategic Task Group at NASA Goddard Space Flight Center. The training dataset was collected through funding from the Oklahoma Space Grant Consortium. References: [1] Pappalardo, R. et al. (2013) Astrobiology, 13, 740–773. [2] Theiling (2020) Icarus, 114216.

Europa↗

Analysis of Salinity Intrusion in the San Francisco Bay-Delta using a GA- Optimized Neural Net, and Application of the Model to Prediction in the Elkhorn Slough Habitat

The San Francisco Bay Delta is a large hydrodynamic complex that incorporates the Sacramento and San Joaquin Estuaries, the Burman Marsh, and the San Francisco Bay proper. Competition exists for the use of this extensive water system both from the fisheries industry, the agricultural industry, and from the marine and estuarine animal species within the Delta. As tidal fluctuations occur, more saline water pushes upstream allowing fish to migrate beyond the Burman Marsh for breeding and habitat occupation. However, the agriculture industry does not want extensive salinity intrusion to impact water quality for human and plant consumption. The balance is regulated by pumping stations located alone the estuaries and reservoirs whereby flushing of fresh water keeps the saline intrusion at bay. The pumping schedule is driven by data collected at various locations within the Bay Delta and by numerical models that predict the salinity intrusion as part of a larger model of the system. The Interagency Ecological Program (IEP) for the San Francisco Bay/Sacramento-San Joaquin Estuary collects, monitors, and archives the data, and the Department of Water Resources provides a numerical model simulation (DSM2) from which predictions are made that drive the pumping schedule. A problem with this procedure is that the numerical simulation takes roughly 16 hours to complete a C:~ prediction. We have created a neural net, optimized with a genetic algorithm, that takes as input the archived data from multiple stations and predicts stage, salinity, and flow at the Carquinez Straits (at the downstream end of the Burman Marsh). This model seems to be robust in its predictions and operates much faster than the current numerical DSM2 model. Because the system is strongly tidal driven, we used both Principal Component Analysis and Fast Fourier Transforms to discover dominant features within the IEP data. We then filtered out the dominant tidal forcing to discover non-primary tidal effects, and used this to enhance the neural network by mapping input-output relationships in a more efficient manner. Furthermore, the neural network implicitly incorporates both the hydrodynamic and water quality models into a single predictive system. Although our model has not yet been enhanced to demonstrate improve pumping schedules, it has the possibility to support better decision-making procedures that may then be implemented by State agencies if desired. Our intention is now to use this model in the smaller Elkhorn Slough complex near Monterey Bay where no such hydrodynamic model currently exists. At the Elkhorn Slough, we are fusing the neural net model of tidally-driven flow with in situ flow data and airborne and satellite remote sensation data. These further constrain the behavior of the model in predicting the longer-term health and future of this vital estuary.

Thompson, David E.↗

Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.

Liu, Antony K.↗