Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sparse approximation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Some methods of encoding simple visual images for use with a sparse distributed memory, with applications to character recognition

To study the problems of encoding visual images for use with a Sparse Distributed Memory (SDM), I consider a specific class of images- those that consist of several pieces, each of which is a line segment or an arc of a circle. This class includes line drawings of characters such as letters of the alphabet. I give a method of representing a segment of an arc by five numbers in a continuous way; that is, similar arcs have similar representations. I also give methods for encoding these numbers as bit strings in an approximately continuous way. The set of possible segments and arcs may be viewed as a five-dimensional manifold M, whose structure is like a Mobious strip. An image, considered to be an unordered set of segments and arcs, is therefore represented by a set of points in M - one for each piece. I then discuss the problem of constructing a preprocessor to find the segments and arcs in these images, although a preprocessor has not been developed. I also describe a possible extension of the representation.

Jaeckel, Louis A.↗
NASA NTRS1988

Using data tagging to improve the performance of Kanerva's sparse distributed memory

The standard formulation of Kanerva's sparse distributed memory (SDM) involves the selection of a large number of data storage locations, followed by averaging the data contained in those locations to reconstruct the stored data. A variant of this model is discussed, in which the predominant pattern is the focus of reconstruction. First, one architecture is proposed which returns the predominant pattern rather than the average pattern. However, this model will require too much storage for most uses. Next, a hybrid model is proposed, called tagged SDM, which approximates the results of the predominant pattern machine, but is nearly as efficient as Kanerva's original formulation. Finally, some experimental results are shown which confirm that significant improvements in the recall capability of SDM can be achieved using the tagged architecture.

Rogers, David↗
NASA NTRS1988

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗
NASA NTRS2014

The Uncertainty of Biomass Estimates from Modeled ICESat-2 Returns Across a Boreal Forest Gradient

The Forest Light (FLIGHT) radiative transfer model was used to examine the uncertainty of vegetation structure measurements from NASA's planned ICESat-2 photon counting light detection and ranging (LiDAR) instrument across a synthetic Larix forest gradient in the taiga-tundra ecotone. The simulations demonstrate how measurements from the planned spaceborne mission, which differ from those of previous LiDAR systems, may perform across a boreal forest to non-forest structure gradient in globally important ecological region of northern Siberia. We used a modified version of FLIGHT to simulate the acquisition parameters of ICESat-2. Modeled returns were analyzed from collections of sequential footprints along LiDAR tracks (link-scales) of lengths ranging from 20 m-90 m. These link-scales traversed synthetic forest stands that were initialized with parameters drawn from field surveys in Siberian Larix forests. LiDAR returns from vegetation were compiled for 100 simulated LiDAR collections for each 10 Mg · ha(exp -1) interval in the 0-100 Mg · ha(exp -1) above-ground biomass density (AGB) forest gradient. Canopy height metrics were computed and AGB was inferred from empirical models. The root mean square error (RMSE) and RMSE uncertainty associated with the distribution of inferred AGB within each AGB interval across the gradient was examined. Simulation results of the bright daylight and low vegetation reflectivity conditions for collecting photon counting LiDAR with no topographic relief show that 1-2 photons are returned for 79%-88% of LiDAR shots. Signal photons account for approximately 67% of all LiDAR returns, while approximately 50% of shots result in 1 signal photon returned. The proportion of these signal photon returns do not differ significantly (p greater than 0.05) for AGB intervals greater than 20 Mg · ha(exp -1). The 50m link-scale approximates the finest horizontal resolution (length) at which photon counting LiDAR collection provides strong model fits and minimizes forest structure uncertainty in the synthetic Larix stands. At this link-scale AGB greater than 20 Mg · ha(exp -1) has AGB error from 20-50% at the 95% confidence level. These results suggest that the theoretical sensitivity of ICESat-2 photon counting LiDAR measurements alone lack the ability to consistently discern differences in inferred AGB at 10 Mg · ha(exp -1) intervals in sparse forests characteristic of the taiga-tundra ecotone.

LiDAR↗
NASA NTRS1991

Implicit solvers for unstructured meshes

Implicit methods were developed and tested for unstructured mesh computations. The approximate system which arises from the Newton linearization of the nonlinear evolution operator is solved by using the preconditioned GMRES (Generalized Minimum Residual) technique. Three different preconditioners were studied, namely, the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over relaxation (SSOR). The preconditioners were optimized to have good vectorization properties. SSOR and ILU were also studied as iterative schemes. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also studied. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗
NASA NTRS1983

Atmospheric effects on radiometry from zenith of a plane with dark vertical protrusions

Effects of an optically thin plane-parallel scattering atmosphere on radiometric imaging from the zenith of a specific surface-type are analyzed. The surface model was previously developed to describe arid steppe, where the sparse vegetation forms dark vertical protrusions from the bright soil-plane. The analysis is in terms of the surface reflectivity to the zenith r sub p for the direct beam, which is formulated as r sub p = r sub i exp (-s tan theta sub 0), where v sub i is the Lambert law reflectivity of the soil, the protrusions parameters s is the projection on a vertical plane of protrusions per unit area and theta sub 0 is the zenith angle. The surface reflectivity r sub p is approximately equal to that for the global irradiance (which is directly measured in the field) only for a narrow range of the solar zenith angles. The effects of the atmosphere when imaging large uniform areas of this type are comparable to those in imaging a Lambert surface with a reflectivity r sub p. Thus, the effects can be approximated by those in the case of a dark Lambert surface (analyzed previously), inasmuch as r sub p is smaller than the soil reflectivity r sub i for any off-zenith illumination. The surface becomes effectively darker with increasing solar zenith angle. Adjacency effects of a reflection from one area and scattering in the instantaneous field of view (object pixel) are analyzed as cross radiance and cross irradiance.

Otterman, J.↗
NASA NTRSDate not supplied

Development of the Ames Global Hyperspectral Synthetic Dataset

This study develops the surface BRDF (bidirectional reflectance distribution function) product of the Ames Global Hyperspectral Synthetic Dataset (AGHSD), based on the corresponding MODIS products, to support the NASA Surface Biology and Geology mission development. A main challenge in deriving a hyperspectral dataset from the multi-band satellite products is how to identify a succinct yet robust algorithm that allow us to infer BRDF at unobserved wavelengths based on the few observed bands. Using the theories of radiative transfer in vegetation canopies, we arrive at a simple equation that accurately approximates hyperspectral surface BRDF as the weighted sum of components from the soil and the vegetation. Each of the components is modeled by the product of the spectrally-dependent optical properties of a surface element (the spectra of the soil surface reflectance, the leaf single albedo, or the canopy scattering coefficient) and a spectrally-independent bidirectional scattering function. The optical properties of the soil and the vegetation can be obtained from existing spectral libraries or model simulations. The bidirectional scattering functions are represented by the Ross-Thick-Li-Sparse BRDF model, where the linear coefficients are estimated with regression analysis from the multi-band MODIS data. We validate the algorithm with simulations by Monte Carlo Ray Tracing model experiments, and the results are highly consistent with the theoretic derivation. We apply the algorithm to generate the AGHSD BRDF product at 1km and 8-day resolutions for the year of 2019. The results are biogeochemically and physically coherent and consistent, and thus serve the goal to support the science and application development of the SBG community.

Hyperspectral↗
NASA NTRS2000

The DYNAMO Orbiter Project: High Resolution Mapping of Gravity/Magnetic Fields and In Situ Investigation of Mars Atmospheric Escape

Dynamo is a small Mars orbiter planned to be launched in 2005 or 2007, in the frame of the NASA/CNES Mars exploration program. It is aimed at improving gravity and magnetic field resolution, in order to better understand the magnetic, geologic and thermal history of Mars, and at characterizing current atmospheric escape, which is still poorly constrained. These objectives are achieved by using a low periapsis orbit, similar to the one used by the Mars Global Surveyor spacecraft during its aerobraking phases. The proposed periapsis altitude for Dynamo of 120-130 km, coupled with the global distribution of periapses to be obtained during one Martian year of operation, through about 5000 low passes, will produce a magnetic/gravity field data set with approximately five times the spatial resolution of MGS. Low periapsis provides a unique opportunity to investigate the chemical and dynamical properties of the deep ionosphere, thermosphere, and the interaction between the atmosphere and the solar wind, therefore atmospheric escape, which may have played a crucial role in removing atmosphere, and water, from the planet. There is much room for debate on the importance of current atmosphere escape processes in the evolution of the Martian atmosphere, as early "exotic" processes including hydrodynamic escape and impact erosion are traditionally invoked to explain the apparent sparse inventory of present-day volatiles. Yet, the combination of low surface gravity and the absence of a substantial internally generated magnetic field have undeniable effects on what we observe today. In addition to the current losses in the forms of Jeans and photochemical escape of neutrals, there are solar wind interaction-related erosion mechanisms because the upper atmosphere is directly exposed to the solar wind. The solar wind related loss rates, while now comparable to those of a modest comet, nonetheless occur continuously, with the intriguing possibility of important cumulative and/or enhanced effects over the several billion years of the solar system's life. If the detailed history of the Martian internal field could be traced back, and the current escape processes could be understood well enough to model the expected stronger losses under early Sun conditions, one could go a long way toward constraining this part of the mysterious history of Mars' atmosphere.

Smrekar, S.↗
NASA NTRS1992

The geophysical signal of the martian global dichotomy

A first-order tectonic question for Mars is the origin and nature of the global dichotomy (GD) separating approximately the northern and southern hemispheres of the planet. It is appropriate to focus on hypotheses for the origin of the GD as well as on geophysical models of related internal structure that are constrained by present-day observations. There are basic planetary scale observations that relate to the GD: (1) the dichotomy boundary separates two fundamentally different elevations on the planet, as the terrain to the north is lower by an average of about 3 km; (2) the boundary separates terrain of regionally distinct crater ages, heavily cratered (older) in the south and sparsely cratered (younger) in the north; (3) the amount of ancient crust apparently removed from north of the dichotomy boundary cannot be accounted for by simple surface erosion and deposition in the south, and the constraint becomes particularly severe if isostatic adjustment is presumed to have accompanied this process. This last point leads to the supposition that some type of interior process must have been responsible for the creation of the GD. An obvious way to create the observed elevation difference between the two hemispheres is with a thinner crust in the north, although a denser cust would also work. Hypotheses for producing a thinner northern crust include preferential sub-crustal erosion, a giant impact, and simply invoking the crustal thickness difference as a primordial feature of the planet. If the GD represents a fundamental change in the crustal thickness of Mars, then there should be geophysical evidence of this. The center-of-figure to center-of-mass offset of the planet may be related to the GD, but the Tharsis topography must certainly also contribute. If the Tharsis and GD effects can be separated, then a crustal thickness model can be tested, though the results will not be unique. The gravity field provides another geophysical constraint on GD models. A simple change in thickness of an isostatically compensated crust should show a characteristic gravity signal across the dichotomy boundary. There is a gravity anomaly that is clearly associated with the boundary in regions where the gravity signal is not cluttered by contributions from other features such as Tharsis and Elysium. The gravity signal has a complex spatial relationship with the boundary, and efforts are presently underway to model this anomaly.

Phillips, Roger J.↗
NASA NTRS1994

Improvements to TOVS retrievals over sea ice and applications to estimating Arctic energy fluxes

Modeling studies suggest that polar regions play a major role in modulating the Earth's climate and that they may be more sensitive than lower latitudes to climate change. Until recently, however, data from meteorological stations poleward of 70 degs have been sparse, and consequently, our understanding of air-sea-ice interaction processes is relatively poor. Satellite-borne sensors now offer a promising opportunity to observe polar regions and ultimately to improve parameterizations of energy transfer processes in climate models. This study focuses on the application of the TIROS-N operational vertical sounder (TOVS) to sea-ice-covered regions in the nonmelt season. TOVS radiances are processed with the improved initialization inversion ('3I') algorithm, providng estimates of layer-average temperature and moisture, cloud conditions, and surface characteristics at a horizontal resolution of approximately 100 km x 100 km. Although TOVS has flown continuously on polar-orbiting satellites since 1978, its potential has not been realized in high latitudes because the quality of retrievals is often significantly lower over sea ice and snow than over the surfaces. The recent availability of three Arctic data sets has provided an opportunity to validate TOVS retrievals: the first from the Coordinated Eastern Arctic Experiment (CEAREX) in winter 1988/1989, the second from the LeadEx field program in spring 1992, and the third from Russian drifting ice stations. Comparisons with these data reveal deficiencies in TOVS retrievals over sea ice during the cold season; e.g., ice surface temperature is often 5 to 15 K too warm, microwave emissivity is approximately 15% too low at large view angles, clear/cloudy scenes are sometimes misidentified, and low-level inversions are often not captured. In this study, methods to reduce these errors are investigated. Improvements to the ice surface temperature retrieval have reduced rms errors from approximately 7 K to 3 K; correction of microwaves (50 GHz) brightness temperatures for view angle dependence now allows the surface type (open water versus sea ice) to be determined; modifications to the clear/cloud tests have improved cloud detection over sea ice, especially where it is inhomogeneous; and improved surface temperature estimates have resulted in more successful diagnoses of low-level stratification. Applications of improved retrievals to studies of the Arctic energy budget are encouraging. Preliminary calculations of air-ice stress vectors, 10-m wind speed, and horizontal advection of heat and moisture are presented.

Francis, Jennifer A.↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗
NASA NTRS2017

Assimilation of AERONET and MODIS AOT Observations Using Variational and Ensemble Data Assimilation Methods and Its Impact on Aerosol Forecasting Skill

Data assimilation of Aerosol Robotic Network (AERONET) and Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol optical thickness (AOT) for aerosol forecasting was tested within the Navy Aerosol Analysis Prediction System (NAAPS) framework, using variational and ensemble data assimilation methods. Navy aerosol forecasting currently makes use of a deterministic NAAPS simulation coupled to Navy Variational Data Assimilation System for aerosol optical depth, a two-dimensional variational data assimilation system, for MODIS AOT assimilation. An ensemble version of NAAPS (ENAAPS) coupled to an ensemble adjustment Kalman filter (EAKF) from the Data Assimilation Research Testbed was recently developed, allowing for a range of data assimilation and forecasting experiments to be run with deterministic NAAPS and ENAAPS. The main findings are that the EAKF, with its flow-dependent error covariances, makes better use of sparse observations such as AERONET AOT. Assimilating individual AERONET observations in the two-dimensional variational system can increase the analysis errors when observations are located in high AOT gradient regions. By including AERONET with MODIS AOT assimilation, the magnitudes of peak aerosol events (AOT> 1) were better captured with improved temporal variability, especially in India and Asia where aerosol prediction is a challenge. Assimilating AERONET AOT with MODIS had little impact on the 24 h forecast skill compared to MODIS assimilation only, but differences were found downwind of AERONET sites. The 24 h forecast skill was approximately the same for forecasts initialized with analyses from AERONET AOT assimilation alone compared to MODIS assimilation, particularly in regions where the AERONET network is dense; including the United States and Europe, indicating that AERONET could serve as a backup observation network for over-land synoptic-scale aerosol events.

Rubin, Juli I.↗
NASA NTRS1989

Totally parallel multilevel algorithms for sparse elliptic systems

The fastest known algorithms for the solution of a large elliptic boundary value problem on a massively parallel hypercube all require O(log(n)) floating point operations and O(log(n)) distance-1 communications, if massively parallel is defined to mean a number of processors proportional to the size n of the problem. The Totally Parallel Multilevel Algorithm (TPMA) that has, as special cases, four of these fast algorithms is described. These four algorithms are Parallel Superconvergent Multigrid (PSMG), Robust Multigrid, the Fast Fourier Transformation (FFT) based Spectral Algorithm, and Parallel Cyclic Reduction. The algorithm TPMA, when described recursively, has four steps: (1) project to a collection of interlaced, coarser problems at the next lower level; (2) apply TPMA, recursively, to each of these lower level problems, solving directly at the lowest level; (3) interpolate these approximate solutions to the finer grid, and to verage them to form an approximate solution on this grid; and (4) refine this approximate solution with a defect-correction step, using a local approximate inverse. Choice of the projection operator (P), the interpolation operator (Q), and the smoother (S) determines the class of problems on which TPMA is most effective. There are special cases in which the first three steps produce an exact solution, and the smoother is not needed (e.g., constant coefficient operators).

Frederickson, Paul O.↗
NASA NTRS1995

Linear and nonlinear trending and prediction for AVHRR time series data

The variability of AVHRR calibration coefficient in time was analyzed using algorithms of linear and non-linear time series analysis. Specifically we have used the spline trend modeling, autoregressive process analysis, incremental neural network learning algorithm and redundancy functional testing. The analysis performed on available AVHRR data sets revealed that (1) the calibration data have nonlinear dependencies, (2) the calibration data depend strongly on the target temperature, (3) both calibration coefficients and the temperature time series can be modeled, in the first approximation, as autonomous dynamical systems, (4) the high frequency residuals of the analyzed data sets can be best modeled as an autoregressive process of the 10th degree. We have dealt with a nonlinear identification problem and the problem of noise filtering (data smoothing). The system identification and filtering are significant problems for AVHRR data sets. The algorithms outlined in this study can be used for the future EOS missions. Prediction and smoothing algorithms for time series of calibration data provide a functional characterization of the data. Those algorithms can be particularly useful when calibration data are incomplete or sparse.

Smid, J.↗
NASA NTRS1981

Genesis of recent silicic magmatism in the Medicine Lake Highland, California - Evidence from cognate inclusions found at Little Glass Mountain

Sparse, granular inclusions of early-formed minerals found within the Little Glass Mountain rhyolite flows in northern California are shown to provide a means of characterizing the physical conditions, at depth, beneath the Medicine Lake Highland during the latest phase of volcanic activity. Mineral compositions, in combination with thermodynamic calculations and experiments, suggest crystalization at a pressure of 5,200 bars within a 966-836 C temperature range; implying that mineral segregation and equilibration occurred at a depth of 15-18 km beneath the surface. In addition, mass balance calculations indicate that the Medicine Lake flow is a close approximation to the parental magma for the latest silicic lavas.

Mertzman, S. A.↗
NASA NTRS2002

Image Dis-Integration for Improved Plasmasphere Visualization

Traditionally, study of the plasmasphere has involved terrestrial observation of local characteristics. Global modeling of the plasmasphere in such an observation regime made use of an ensemble of (sparse) local measurements. Recently, sensors aboard the IMAGE (Imager for Magnetopause-to-Aurora Global Exploration) satellite (in particular, the EUV (Extreme Ultra Violet) Imager) have created the potential for truly global study of the plasmasphere. IMAGE was launched in spring of 2000 in an orbit with apogee altitude 7.2 R E (Earth radii) and perigee altitude 1000 km. IMAGE's EUV sensor allows an external view of the distribution of cold plasma in the plasmasphere to be acquired. EUV is designed to image light emission at 30.4 nanometers, which is the emission wavelength of the He + ion in the presence of solar radiation. He + makes up approximately 15-20% of the plasma in the plasmasphere, thus imaging of He + enables determination of plasma distribution. The EUV instrument provides a 90° by 84° field of view which is imaged as an equally spaced 150x140 pixel array on a spherical imaging surface. The EUV produces an image approximately every 10 minutes when the sensor is operating. Since EUV images contain line-of-sight integrations of plasma distributions, they do not directly express equatorial plane density (which would enable comparison of observed plasma distributions with predictions from models). Furthermore, the plasma density at any point in three-space is not known. The goal of our work was development of a technique that can enable plasma density to be determined throughout three-space. Our approach to creation of a three-space representation of the plasma distribution involves disintegrating the EUV lines of sight to form a volumetric map of plasma densities.

Timothy S. Newman↗
NASA NTRS1993

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗
NASA NTRS1985

Blended SST fields combined SMMR/ship data (summary of technique)

All remote sensors have bias problems which vary in space and time. The infrared measurements are very sparse in cloudy areas: the SMMR has coverage limitations due to the interference of land areas and only half time operation. If one simply averaged the data from these various sources, the combination of the biases and varying sampling would introduce artifacts into the analyzed product which would reflect the sampling and which could be quite misleading. The remote sensors provide reasonable estimates of gradients of SST over length scales greater than the sensor resolution but smaller than hemispheric. Ships can be considered unbiased but are poorly distributed on small scales except in the densest shipping lanes. An analysis scheme based on a suggestion by Holl was developed which exploits just this approximately complementarity to provide an improved SST analysis. The resulting accuracy appears to be in the required approx. 0.5 C range for the particular case of using SMMR and ship data to produce two degree latitude by two degrees longitude monthly analyses.

Wilheit, T. T.↗