ALBERT

Description: Background information and exercises are provided to: (1) establish or expand understanding of the concepts, methods, and terminology of computer processing of image producing data; (2) develop insight into the advantages of computer based image processing compared with the photointerpretation approach for processing, classifying, interpreting, and applying remote sensing data; (3) foster a broad perspective on the principal of the main techniques for image enhancement, pattern recognition, and thematic classification; (4) appreciate the pros and cons of batch and interactive modes of image analysis; (5) examine and evaluate some specific computer generated products for subscenes in Pennsylvania and New Jersey; and (6) interrelate these particular examples of output with more theoretical explanations of computer processing strategies and procedures.

Description: Images characteristics of geographic regions other than the northeastern part of the United States are presented for interpretation. Pre- and post-eruption imagery of Mt. St. Helens volcano serves to demonstrate the advantages of thermal infrared sensing, and the potential for developing a timely, decision oriented thematic map to be used in solving drought-related problems in Upper Volta is examined to show the applicability of satellite remote sensing in all geographic areas.

Description: The geostationary and polar satellites comprising the current operational system are discussed. The data acquisition capabilities of both satellite types and their complementary functions are reviewed. The advanced very high resolution radiometer on the TIROS N satellites is particularly addressed along with the imaging and atmospheric sounding instrumentation aboard the GOES satellites. The dissemination of the satellite data to the prospective users is also discussed.

Description: The 1980 eruptions of Mount St. Helens in southeast Washington resulted in a pronounced effect on the surface and ground water resources of the state. In response to the volcanic activity, the U.S. Geological Survey intensified statewide surface and ground water sampling programs to determine the nature and magnitude of the volcanic-induced variations. Streams to the east of Mount St. Helens received the major ash fallout. Chemical effects were best noted in smaller streams sampled 60 to 70 miles northeast of Mount St. Helens. The chemical variations observed were pronounced but short lived. Sulfate and chloride increases in anionic composition were prevalent immediately following the eruption; however, the original bicarbonate predominance was again attained within several days. Suspended iron and aluminum concentrations were similarly elevated during the period of greatest ash deposition (highest turbidity); however, the dissolved concentrations remained relatively constant. Depressions of pH were minor and short lived. Streams draining to the south, tributaries to the Columbia river, showed little observable changes in water chemistry. Streams draining to the west (Toutle river and its tributaries) were compositionally affected by the various volcanic activities. Chloride and sulfate anion percentage exceeded the bicarbonate percentage up to one month following the eruption period. Streams and lakes sampled in the immediate vicinity of Mount St. Helens, in addition to trace metals, contained organic compounds derived from decomposing wood buried in the debris deposits. This organic material may constitute a significant source of organic compounds to surface and ground water for some time to come.

Description: The benefits-to-cost ratio of satellite remote sensing, both as a substitute for conventional methods of monitoring and assessing resources, and as a supplement to these methods is examined using a model which analyzes the cost of aerial photography versus satellite scanner for producing and interpreting an image of the Earth's surface sized to LANDSAT dimensions. Examples of cost savings are tabulated for ground surveys, aerial photos, and LANDSAT. Possible additional benefits from LANDSAT D are assessed. The way in which satellites fit into more comprehensive models for resources management is discussed. It is shown that remote sensing is but one essential component in a complex system that aggregates technical. Socioeconomics, political, cultural, and other factors in the human decision process.

Description: Computer processing facilitates extraction of information from every pixel by executing a variety of functional operations, called processed algorithms, in general or specialized routines. The best results are obtained when data from more than one multispectral band are used together. Multivariate tatistical analysis, computer tape characteristics, processing modes, and a choice of systems (batch or interactive) are discussed. The major operations in computer processing elaborated include: preprocessing, enhancement, effects of rationing, and classification. Techniques for multisource data correlation are considered with emphasis on geobased systems.

Description: Learning objectives include: (1) developing an understanding of the implications of the term "near surface observations"; (2) associating the appearance of large ground features as seen in satellite imagery with their appearance as seen from the ground; (3) grasping criteria and procedures for selecting training sites on the ground for use in supervised classification; (4) running through an example of training site selection; (5) becoming familiar with several methods of accuracy assessment; (6) becoming aware of the approach and value of making supporting measurements of the spectral and physical properties of materials on the ground and from aircraft; (7) taking note of the different types of instruments used in making specific ground measurements; and (8) appreciating the rationale underlying laboratory and field studies on or near the Earth's surface for the purpose of developing new sensor systems.

Description: Learning objectives include: (1) developing a facility for applying conventional techniques of photointerpretation to small scale (satellite) imager; (2) promoting the ability to locate, identify, and interpret small natural and man made surface features in a LANDSAT image; (3) using supporting imagery, such as aerial and space photography, to conduct specific applications analyses; (4) learning to apply change detection techniques to recognize and explain transient and temporal events in individual or seasonal imagery; (5) producing photointerpretation maps that define major surface units, themes, or classes; (6) classifying or analyzing a scene for specific discipline applications in geology, agriculture, forestry, hyrology, coastal wetlands, and environmental pollution; and (7) evaluating both advantages and shortcomings in relying on the photointerpretive approach (rather than computer based analytical approach) for extracting information from LANDSAT data.

Description: Activities are provided for: (1) developing insight into the way in which the LANDSAT MSS produces multispectral data; (2) promoting understanding of what a "pixel" means in a LANDSAT image and the implications of the term "mixed pixel"; (3) explaining the concept of spectral signatures; (4) deriving a simple signature for a class or feature by analysis: of the four band images; (5) understanding the production of false color composites; (6) appreciating the use of color additive techniques; (7) preparing Diazo images; and (8) making quick visual identifications of major land cover types by their characteristic gray tones or colors in LANDSAT images.

Description: Learning objectives of the activities provided include: (1) reading the annotation of a LANDSAT image; (2) becoming acquainted with the characteristics of 1:1,000,000 scale transparencies and prints of MSS images; (3) noting the general information visible in LANDSAT photo products; (4) observing changes of appearance of any ground feature or class in the black and white images made from the four MSS bands and the characteristic color of each class in color composites; (5) determining the degree to which a LANDSAT image meets map accuracy standards and can be fitted to map projections; (6) assessing the effects of LANDSAT enlargements and scale changes and of the limitations of satellite resolution relative to aerial photos; (7) observing the influence of time of acquisition (season) on a scene; (8) getting a feel for image quality as dependent on processing and photoreproduction; (9) appreciating the characteristics of the RBV and thermal band imagery obtained from LANDSAT-3; and (10) becoming familiar with certain attributes of adjacent LANDSAT images which permit them to be joined in mosaics and to be viewed in stereo.

Description: Programs for the development and operation of meteorological satellites from the TIROS 1 satellite and the establishment of NASA through the 1960's are described. The technical problems confronted in the development of the early satellite systems are discussed in addition to issues in international involvement and program support. The TIROS and Nimbus series satellites are primarily addressed.

Description: The primary LANDSAT mission and the system requirements are summarized and pertinent parameters of the spacecraft, its orbit, and payload are tabulated. The history acquisition to entry into the archives for storage and product generation and dissemination is recounted. The LANDSAT D data handling plan is discussed showing requirements for both the MSS CCT and the thematic mapper CCT.

Description: The term remote sensing is defined as well as ideas such as class, pattern, feature, pattern recognition, feature extraction, and theme. The electromagnetic spectrum is examined especially those wavelength regions available to remote sensing. Relevant energy and wave propagation laws are discussed and the characteristics of emitted and reflected radiation and their detection are investigated. The identification of classes by their spectral signatures, the multispectral approach, and the principal types of sensors and platforms used in remote sensing are also considered.

Description: By using the most complete available records of direct beam radiation and volcanic eruptions, an historical analysis of the role of the latter in modulating the former was made. A very simple fallout and dispersion model was applied to the historical chronology of explosive eruptions. The resulting time series explains about 77 percent of the radiation variance, as well as suggests that tropical and subpolar eruptions are more important than mid-latitude eruptions in their impact on the stratospheric aerosol optical depth. The simpler climatic models indicate that past hemispheric temperature can be stimulated very well with volcanic and CO2 inputs and suggest that climate forecasting will also require volcano forecasting. There is some evidence that this is possible some years in advance.

Description: Visible and infrared pictures from two Geostationary Operational Environmental Satellite Systems satellites, in circular orbits at about 19,000 nautical miles, are available continuously at approximately 30 minute intervals. Still pictures and film loops from this system vividly depict the events associated with the May 18, 1980 eruption of Mount St. Helens. The initial explosion, shock wave, and visible horizontal dust distribution during the following week are readily apparent. Meteorological wind and height fields permit the inference of the vertical distribution of volcanic dust as well as explain the atmospheric behavior which caused the visible and nonvisible dust distribution.

Description: Examples of the use of geostationary satellites in meteorology are given. Studies of the rate of change of cumulus clouds and cloud systems and wind parameter determination from cloud motions are reviewed. Computer processed imagery products are also discussed.

Description: The Defense Meteorological Satellite Program is described with particular emphasis on the military applications of METSAT data. Satellite operational support, data processing and image quality requirements are discussed.

Description: A historical overview of the pioneer projects for the development of meteorological satellites is given. In addition, the parallel development of the responsible space agencies and panels is addressed. The Vanguard 2 satellite, the first Earth radiation experiment, and the vidicon equipped TIROS-1 satellite are discussed.

Description: A brief review of the effects of climate and weather on the 1980 Mount St. Helens eruptions and the subsequent dispersion of ash and gases and the reciprocal influences of the eruptions on climate and climatology is presented. The effects of mesoscale destruction of snow fields and vegetation, a revised mountain profile, and ash deposits are addressed along with impacts on hemispheric climate and disruption of normal climatological observations, in the areas directly affected by the explosions and ashfall. Environmental and economic consequences are also considered.

Description: Geographic coverage frequency and geographic shot density for a satellite borne Doppler lidar wind velocity measuring system are measured. The equations of motion of the light path on the ground were derived and a computer program devised to compute shot density and coverage frequency by latitude-longitude sections. The equations for the coverage boundaries were derived and a computer program developed to plot these boundaries, thus making it possible, after an application of a map coloring algorithm, to actually see the areas of multiple coverage. A theoretical cross-swath shot density function that gives close approximations in certain cases was also derived. This information should aid in the design of an efficient data-processing system for the Doppler lidar.

Description: Calculation of eletric field vectors within and near an isolated thundercloud (which has a given volume charge distribution) is envisaged. The maximum field strength within the thundercloud is calculated. The effects of screening layers, both above and below the thundercloud may be investigated, as well time-dependent potential problem. The study can lead to a better understanding of the charge distributions in a thundercloud, which in turn may shed some light on the actual mechanism of charging of a thundercloud.

Description: The history of remote sensing is reviewed and the scope and versatility of the several remote sensing systems already in orbit are discussed, especially those with sensors operating in other EM spectral modes. The multisensor approach is examined by interrelating LANDSAT observations with data from other satellite systems. The basic principles and practices underlying the use of thermal infrared and radar sensors are explored and the types of observations and interpretations emanating from the Nimbus, Heat Capacity Mapping Mission, and SEASAT programs are examined. Approved or proposed Earth resources oriented missions for the 1980's previewed include LANDSAT D, Stereosat, Gravsat, the French satellite SPOT-1, and multimission modular spacecraft launched from space shuttle. The pushbroom imager, the linear array pushbroom radiometer, the multispectral linear array, and the operational LANDSAT observing system, to be designed the LANDSAT-E series are also envisioned for this decade.

Description: The use of computer aided classification of LANDSAT data in developing water quality plans for New Jersey watersheds is used to exemplify how a state natural resource management program benefits from satellite imagery. The transition of a research and development system into an operational remote sensing system to help decision makers is demonstrated. Nontechnial issues that can assist (or hinder) an agency in adopting a new technology are examined. The progress of LANDSAT use by state government from the earliest stage of curiosity through to incorporation in actual state planning methods is charted. Potential applications of LANDSAT data to real information needs and solutions to management problems are examined. The problems and mistakes that occurred in using LANDSAT data in the past are discussed as well as the ways by which these problems were overcome.

Description: Information and activities are provided to: (1) enhance the ability to distinguish between a Geographic Information System (GIS) and a data management system; (2) develop understanding of spatial data handling by conventional methods versus the automated approach; (3) promote awareness of GIS design and capabilities; (4) foster understanding of the concepts and problems of data base development and management; (5) facilitate recognition of how a computerized GIS can model conditions in the present "real world" to project conditions in the future; and (6) appreciate the utility of integrating LANDSAT and other remotely sensed data into the GIS.

Description: Future developments in satellite meteorology are proposed and examined in the light of policy and funding changes.

Description: Atmospheric and surface fields produced from the TIROS N high resolution infrared sounder/microwave sounding unit data are discussed. The data were analyzed by direct physical inversion of the multispectral radiative transfer equation. Sea ice mapping and sea surface temperature determination are addressed.

Description: Some of the research opportunities that are available from geostationary satellite data, especially from MONEX during the Global Weather Experiment are described. In addition, the International Cloud Climatology Program is described and the use of satellite imagery in studying mesoscale convection complexes is discussed. Finally, information derived from satellite data on ocean winds is examined.

Description: The evidence that volcanic eruptions affect climate is reviewed. Single explosive volcanic eruptions cool the surface by about 0.3 C and warm the stratosphere by several degrees. Although these changes are of small magnitude, there have been several years in which these hemispheric average temperature changes were accompanied by severely abnormal weather. An example is 1816, the "year without summer" which followed the 1815 eruption of Tambora. In addition to statistical correlations between volcanoes and climate, a good theoretical understanding exists. The magnitude of the climatic changes anticipated following volcanic explosions agrees well with the observations. Volcanoes affect climate because volcanic particles in the atmosphere upset the balance between solar energy absorbed by the Earth and infrared energy emitted by the Earth. These interactions can be observed. The most important ejecta from volcanoes is not volcanic ash but sulfur dioxide which converts into sulfuric acid droplets in the stratosphere. For an eruption with its explosive magnitude, Mount St. Helens injected surprisingly little sulfur into the stratosphere. The amount of sulfuric acid formed is much smaller than that observed following significant eruptions and is too small to create major climatic shifts. However, the Mount St. Helens eruption has provided an opportunity to measure many properties of volcanic debris not previously measured and has therefore been of significant value in improving our knowledge of the relations between volcanic activity and climate.

Description: The spectral and spatial requirements for remote sensing in the next decade are presented. The requirements presented were obtained through extensive literature research and discussions with leading members of the various remote sensing research communities. In the 0.35-2.5 micron region of the spectrum, numerous bands will be needed at bandwidths as narrow as 10-20 nanometers. There is also growing interest in the thermal infrared (8-14 microns). Spatial resolution (instantaneous field of view) of 5 to 10 meters will be of great benefit to many fields of remote sensing.

Description: The presence of positive serial correlation (autocorrelation) in remotely sensed data results in an underestimate of the variance-covariance matrix when calculated using contiguous pixels. This underestimate produces an inflation in F statistics. For a set of Thematic Mapper Simulator data (TMS), used to test the ability to discriminate a known geobotanical anomaly from its background, the inflation in F statistics related to serial correlation is between 7 and 70 times. This means that significance tests of means of the spectral bands initially appear to suggest that the anomalous site is very different in spectral reflectance and emittance from its background sites. However, this difference often disappears and is always dramatically reduced when compared to frequency distributions of test statistics produced by the comparison of simulated training sets possessing equal means, but which are composed of autocorrelated observations. Previously announced in STAR as N82-25602

Description: A complete automatic unsupervised classification approach has been developed to estimate the proportion of corn, soybeans, and other crops in a Landsat 5-by-6-nanometer segment. The technique involves extracting crop growth variables that are predictable from agromet models, it provides a signature extension not possible before. The technique has been applied in 3 years (1978-1980) for 56 segments and shows neither a year-to-year effect nor any geographic effect. The results to date are considered far better than those obtainable from any other currently available technique. Nontheless, the technique can be improved to reduce the bias.

Description: Since NASA's program of research concerning remote sensing was initiated in the 1960s, one of its major objectives has been to advance the state-of-the-art in machine processing of satellite acquired multispectral data. Possibilities have been studied regarding a use of these data to identify type, to monitor condition, and to estimate the ontogenetic stage of cultural vegetation. The present investigation provides a review of the state-of-the-art of the technology used to make remote sensing crop production estimates in foreign regions. Attention is given to Landsat data acquisition, aspects of registration and preprocessing, questions of data transformation, data modeling, proportion estimation, labeling, development stage models, crop condition models, and an outlook regarding future developments.

Description: Automated information extraction procedures for analysis of multitemporal Landsat data in non-U.S. crop inventory and monitoring are reviewed. Experiments to develop and evaluate crop area estimation technologies for spring small grains, summer crops, corn, and soybeans are discussed. Previously announced in STAR as N82-32793

Description: Agriculture, which occupies a vital position in the economy of the State of California, depends crucially on the available water. The California Department of Water Resources (DWR) is, therefore, greatly concerned with the total water requirements for agricultural applications. In view of the limitations of an area-limited, single-date survey system, the DWR has been cooperating with NASA and the University of California in a study of the applicability of Landsat imagery and digital data as an aid in making decisions concerning the management of water resources. Attention is given to a statewide inventory of irrigated land, computer-assisted estimation and mapping of irrigated land, and a crop type analysis using Landsat digital data.

Description: The California Integrated Remote Sensing System's San Bernardino County Project was devised to study the utilization of a data base at a number of jurisdictional levels. The present paper discusses the implementation of change-detection and masking techniques in the updating of Landsat-derived land-cover maps. A baseline landcover classification was first created from a 1976 image, then the adjusted 1976 image was compared with a 1979 scene by the techniques of (1) multidate image classification, (2) difference image-distribution tails thresholding, (3) difference image classification, and (4) multi-dimensional chi-square analysis of a difference image. The union of the results of methods 1, 3 and 4 was used to create a mask of possible change areas between 1976 and 1979, which served to limit analysis of the update image and reduce comparison errors in unchanged areas. The techniques of spatial smoothing of change-detection products, and of combining results of difference change-detection algorithms are also shown to improve Landsat change-detection accuracies.

Description: A technique is examined for improving the comparability of Landsat multisppectral scanner (MSS) data acquired on different dates. The technique involves conversion of digital brightness counts to relative radiance values measured in energy units (milliwatts per square centimeter-steradian). The statistical data of signature from 23 land cover (or biomass) classifications derived from all three Landsats were compared before and after the radiance normalization. Significant convergence occurred among the data sets for mean spectral values and the variances associated with each of seven major land cover types for MSS bands 4, 5, and 7. Overall, the variance attributed to the sensor component was reduced from 5.39 to 2.69 percent, with the largest decrease occurring in band 4 (14.4 percent to 3.7 percent).

Description: The problem of estimating the parameters which determine a mixture density is reviewed as well as maximum likelihood estimation for it. A particular iterative procedure for numerically approximating maximum likelihood estimates for mixture density problems is considered. This EM algorithm, is a specialization to the mixture density context of a general algorithm of the same name used to approximate maximum likelihood estimates for incomplete data problems. The formulation and theoretical and practical properties of the EM algorithm for mixture densities are discussed focussing in particular on mixtures of densities from exponential families.

Description: An important question in cluster analysis and pattern recognition is the determination of the number of clusters into which a given population should be divided. Frequently, particularly when certain specific clustering methods are being used, the number of clusters is taken to be equal to the number of modes, or local maxima, in the probability density function underlying the given data set. The use of kernal density estimates in mode estimation is discussed. The test statistic to be used is defined and a bootstrap technique for assessing significance is given. An illustrative application is followed by an examination of the asymptotic behavior of the test statistic.

Description: Experience gained in applying nonparametric maximum likelihood techniques of density estimation to judge the comparative quality of various estimators is reported. Two invariate data sets of one hundered samples (one Cauchy, one natural normal) are considered as well as studies in the multivariate case.

Description: The problem of estimating a density f on R sup d from a sample Xz(1),...,X(n) of independent identically distributed random vectors is critically examined, and some recent results in the field are reviewed. The following statements are qualified: (1) For any sequence of density estimates f(n), any arbitrary slow rate of convergence to 0 is possible for E(integral/f(n)-fl); (2) In theoretical comparisons of density estimates, integral/f(n)-f/ should be used and not integral/f(n)-f/sup p, p 1; and (3) For most reasonable nonparametric density estimates, either there is convergence of integral/f(n)-f/ (and then the convergence is in the strongest possible sense for all f), or there is no convergence (even in the weakest possible sense for a single f). There is no intermediate situation.

Description: A method of determining three-dimensional motion and structure from two image frames is presented. The method requires eight point correspondences between the two frames, from which motion and structure parameters are determined by solving a set of eight linear equations and a singular value decomposition of a 3x3 matrix. It is shown that the solution thus obtained is unique.

Description: Methods for the segmentation and recognition of individual targets sensed with forward looking infrared detectors are discussed. Particular attention is given to an adaptive multi-scenario target screener.

Description: The structural analysis of remotely sensed imagery is defined and basic techniques for implementing the process are described. Structural analysis uses knowledge of the properties of an entity, its parts and their relationships, and the relationships in which it participates at a higher level to locate and recognize objects in a visual scene. The representation of structural knowledge, the development of algorithms for using the knowledge to help analyze an image, and techniques for storage and retrieval of relational models are addressed.

Description: Computer-based techniques for automatic photointerpretation based upon information derived from texture and morphology analysis of images are discussed. By automatic photointerpretation, is meant the determination of semantic descriptions of the content of the images by computer. To perform semantic analysis of morphology, a heirarchical structure of knowledge representation was developed. The simplest elements in a morphology are strokes, which are used to form alphabets. The alphabets are the elements for generating words, which are used to describe the function or property of an object or a region. The words are the elements for constructing sentences, which are used for semantic description of the content of the image. Photointerpretation based upon morphology is then augmented by textural information. Textural analysis is performed using a pixel-vector approach.

Description: Methods used to measure the geometrical properties of regions in a segmented image are discussed including the use of centroids, moments, and principle axes. In addition, statistical picture properties, particularly those which describe visual texture, are discussed. Gray level statistics, local property statistics, and autocorrelation and power spectrum are addressed.

Description: Research into the role of cooperative, or relaxation, processes in image segmentation is surveyed. Cooperative processes can be employed at several levels of the segmentation process as a preprocessing enhancement step, during supervised or unsupervised pixel classification and, finally, for the interpretation of image segments based on segment properties and relations.

Description: The methodologies and capabilities of image segmentation techniques are reviewed. Single linkage schemes, hybrid linkage schemes, centroid linkage schemes, histogram mode seeking, spatial clustering, and split and merge schemes are addressed.

Description: Given N distinct points (X sub i, Y sub i) and N real numbers Z sub i, BSPLASH constructs a function G (x, y) that satisfies G (x sub i, y sub i) = Z sub i for i = 1,..., N. This C(2) interpolant consists of a bicubic spline approximation and Shepard's bivariate interpolant.

Description: The physical and geometric interpretation of a surface fitting technique, multiquadric equations (MQ), is discussed. The relation of this method and its reciprocal form to biharmonic and harmonic techniques respectively, is pointed out. In both forms the solutions may be viewed as being numerical approximations of an integral equation in which an unknown density function is the physical source for disturbing potential (reciprocal-MQ) or elastic displacement (MQ).

Description: An algorithm is described for constructing a smooth computable function, f, defined over the surface of a sphere and interpolating a set of n data values, u sub i, associated with n locations, P sub i, on the surface of the sphere. The interpolation function, f, will be continuous and have continuous first partial derivatives. The locations, p sub i, are not required to lie on any type of regular grid.

Description: A variety of numerical methods for fitting a function to data given at a set of points scattered throughout a domain in the plane are surveyed. Four classes of methods are discussed: (1) global interpolation; (2) local interpolation; (3) global approximation; and (4) local approximation. Also, two-stage methods and contouring are discussed. The surfaces constructed include polynomials, spline functions, and rational functions, among others.

Description: Some mathematical/statistical problems within the AgRISTARS program amendable to investigations involving the use of surface fitting techniques are overviewed. The Bayes and maximum likelihood rules, bias determination, regression estimators, parameter estimation, and classifier design are addressed.

Description: A simplified, time-dependent energy balance climate model is run at the latitude belt 40-50 deg N. The model solves for the temperatures of the land, air, and 12 vertical oceanic layers, and it includes a wind stirred mixed layer. A change in model ocean optical turbidity from relatively clear (Jerlov 1) to particle rich (Jerlov III) conditions decreases the effective thickness of the oceanic layer in which heat is stored seasonally and increases the seasonal variation of sea surface temperature by 2-3 C. A decrease in the liquid water content of clouds by a factor of 4 warms the model climate and increases the seasonal variation of sea surface temperatures by 2-3 C. A 2-3 C change in the seasonal variation of sea surface temperature is also obtained by varying oceanic mixing through a factor of 2 change in the surface wind speed.

Description: Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

Description: A study has been conducted to determine the cause of a major synoptic difference in the 72 h GLAS model forecasts from 0000 GMT 19 February 1976 that resulted from the inclusion of satellite data. The prognostic differences that resulted in diverging cyclone paths between the forecast that included satellite temperature soundings (SAT) and the forecast that excluded satellite sounding data (NOSAT) have been traced to initial state differences in the upper level wind and temperature patterns. These modifications enhanced the variation of thermal vorticity and thermal advection across the cyclone center and in the SAT case, gave a greater initial rate of movement of the upper-level vorticity maximum associated with the surface cyclone.

Description: Storms on 20 May 1977 generated a vast cirrus deck. Disturbed areas at storm top had equivalent black-body temperatures (T sub BB) much lower than the tropopause temperature, indicative of overshooting tops. The area of T sub BB not greater than -71 C represents the area of convective activity penetrating 2 km above the tropopause. This area was relatively large after cloud tops and radar reflectivities reached their maximum heights. It became much smaller during tornadoes when reflectivities were decreasing. T sub BB was at a minimum at the time of mesocyclone formation. The Del City storm had two periods of growth, as indicated both by reflectivities and the T sub BB areas. The mesocyclone was first detected during the second less intense period of growth; the tornado occurred during decreasing reflectivities. The maintenance of large areas of relatively low T sub BB after tornado dissipation is ascribed to continued convection on the flanks of the storm and to residual updrafts in a thick anvil cloud.

Description: The effect of horizontal model resolution on satellite data impact has been studied for two versions of the GLAS second-order general circulation model: the C-model with a 4-deg latitude by 5-deg longitude resolution and the F-model with a 2.5-deg latitude and 3-deg longitude resolution. It is found that the 48-72 h forecast skill of the GLAS model was significantly improved by the increased resolution. Initial state differences between the SAT and NOSAT cycles using the F-model were on the average smaller than the corresponding differences with the C-model. However, the F-model cycle differences exhibited a smaller scale structure and, in some cases, larger gradients.

Description: Digital time series data at 16 heights within two storms were collected at vertical incidence with a 10-cm Doppler radar. On several occasions during data collection, lightning echoes were observed as increased reflectivity on an oscilloscope display. Simultaneously, lightning signals from nearby electric field change antennas were recorded on an analog recorder together with the radar echoes. Reflectivity, mean velocity, and Doppler spectra were examined by means of time series analysis for times during and after lightning discharges. Spectra from locations where lightning occurred show peaks, due to the motion of the lightning channel at the air speed. These peaks are considerably narrower than the ones due to precipitation. Besides indicating the vertical air velocity that can then be used to estimate hydrometeor-size distribution, the lightning spectra provide a convenient means to estimate the radar cross section of the channel. Subsequent to one discharge, we deduce that a rapid change in the orientation of hydrometeors occurred within the resolution volume.

Description: Image processing and enhancement techniques for improving the geologic utility of digital satellite radar images are reviewed. Preprocessing techniques such as mean and variance correction on a range or azimuth line by line basis to provide uniformly illuminated swaths, median value filtering for four-look imagery to eliminate speckle, and geometric rectification using a priori elevation data. Examples are presented of application of preprocessing methods to Seasat and Landsat data, and Seasat SAR imagery was coregistered with Landsat imagery to form composite scenes. A polynomial was developed to distort the radar picture to fit the Landsat image of a 90 x 90 km sq grid, using Landsat color ratios with Seasat intensities. Subsequent linear discrimination analysis was employed to discriminate rock types from known areas. Seasat additions to the Landsat data improved rock identification by 7%.

Description: An examination of the possibilities of using Landsat data to simulate NOAA-6 Advanced Very High Resolution Radiometer (AVHRR) data on two channels, as well as using actual NOAA-6 imagery, for large-scale hydrological studies is presented. A running average was obtained of 18 consecutive pixels of 1 km resolution taken by the Landsat scanners were scaled up to 8-bit data and investigated for different gray levels. AVHRR data comprising five channels of 10-bit, band-interleaved information covering 10 deg latitude were analyzed and a suitable pixel grid was chosen for comparison with the Landsat data in a supervised classification format, an unsupervised mode, and with ground truth. Landcover delineation was explored by removing snow, water, and cloud features from the cluster analysis, and resulted in less than 10% difference. Low resolution large-scale data was determined useful for characterizing some landcover features if weekly and/or monthly updates are maintained.

Description: To investigate the possibility that significant amounts of tropical tropospheric air may be convectively introduced into the stratosphere, aerosol samplings over Panama were made at various altitudes using a wire impactor collector. The percentage of particle sizes less than the mean mode decreases with height above the tropopause, suggesting depletion of small particles, possibly due to coagulation. Larger aerosols (greater than 0.3 micron in diam.) are more abundant farther above the tropopause, indicating growth, mainly by condensation. The total particle concentration decreases with increasing height above the tropopause, and also with increasing temperature. Aerosols containing smaller-size particles are thus found closer to the tropopause, and larger-size, more-evolved aerosols occur at higher altitudes. These data indicate that convective activity at the Intertropical Convergence Zone may be a source mechanism for stratospheric aerosols.

Description: Although mean circulations are generally credited with dehydration of the earth's stratosphere, convective instability in the tropics converts mean circulations to small residuals of local convective circulations. The effects of large cumulonimbus which penetrate the stratosphere and form huge anvils in the lower stratosphere are discussed with respect to hydration and dehydration of the stratosphere. Radiative heating at anvil base combined with cooling at anvil top drives a dehydration engine considered essential to explain the dry stratosphere. Seasonal and longitudinal variations in dehydration potentials are examined with maximum potential attributed to Micronesian area during winter and early spring.

Description: The sensitivity of a Goddard Laboratory for Atmospheric Sciences global objective analysis cycle to the addition of FGGE level II-b data is assessed. The GOAS system comprises a predictive continuity provided by a model first-guess forecast integrated from a previous forecast and updated by data gathered in the interim. FGGE data originated in the Jan.-Mar. 1979 period and were acquired by rawinsondes, pilot balloons, surface stations, satellites, ships, and drifting buoys deployed during SOP-1. Focussing on 2-5 and 8-day forecasts, comparisons were made of the 6 hr forecast error at the 300 mb height in three experiments using all, no-satellite (NOSAT), and without rawinsondes or pilot balloons modes. Larger errors occurred in the case of NOSAT, while significant corrections to the GOAS predictions were noted using all the FGGE data. It was concluded that all forecasts were improved by inclusion of full FGGE data sets, including forecasting beyond one week.

Description: Observations of increases of stratospheric condensation nuclei suggest a photo-initiated sulphuric acid vapour formation process in spring in polar regions. It is proposed that the sulphuric acid rapidly forms condensation nuclei through attachment to negatively charged multi-ion complexes and that the process may be modulated through variations in solar activity.

Description: Orographically-induced lee-wave clouds were observed over New Mexico by a multichannel scanning radiometer on Skylab during December 1973. Channels centered at 0.83, 1.61 and 2.125 microns were used to determine the cloud optical thickness, thermodynamic phase and effective particle size. An additional channel centered at 11.4 microns was used to determine cloud-top temperature, which was corroborated through comparison with the stereographically determined cloud top altitudes and conventional temperature soundings. Analysis of the measured near-infrared reflection functions at 1.61 and 2.125 microns are most easily interpreted as indicating the presence of liquid-phase water droplets. This interpretation is not conclusive even after considerable effort to understand possible sources for misinterpretation. However, if accepted the resulting phase determination is considered anomalous due to the inferred cloud-top temperatures being in the -32 to -47 C range. Theory for the homogeneous nucleation of pure supercooled liquid water droplets predicts very short lifetimes for the liquid phase at these cold temperatures. A possible explanation for the observations is that the wave-clouds are composed of solution droplets. Impurities in the cloud droplets could decrease the homogeneous freezing rate for these droplets, permitting them to exist for a longer time in the liquid phase, at the cold temperatures found.

Description: When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Description: (Previously announced in STAR as N81-19529)

Description: Problems associated with the use of Landsat multispectral scanner (MSS) imagery for the detection of urban growth and land use patterns are discussed. The presence of vegetation, either original or added between scanning periods, has been found to dramatically effect the range of signatures in a given area. Different land use developmental stages have been successfully identified by means of 1:50,000 scale panchromatic aerial photography, a resolution only considered possible by spaceborne instrumentation with the advent of the Landsat D satellite. Textural information generated through the grey-tone spatial-dependency matrix for the Landsat band 5 data is compared for different years and a change detection algorithm is described. It is found that the addition of vegetation during development after the removal of natural vegetation resulted in error of omission in the single band data, which must therefore only be used in concert with other data sources.

Description: Procedures necessary for the development of a generalized three-dimensional perspective software capability in support of graphic, topographic, and color mapping of Landsat data are reviewed. The NASA Earth Resources Laboratory developed the procedures in order to facilitate the processing and analysis of disparate, geographically oriented base maps from aircraft and satellite sensors. Perspective displays are obtained through a translation of the space-viewed object to a vantage point coordinate system, followed by a rotation through two angles for alignment along the vantage line of sight, and finally a perspective transformation to yield two-dimensional displays with no hidden lines. Matrix equations for the transformations are reviewed, including scaling, and block diagrams are provided of the data and perspective software systems. The classification data plane may be mapped onto a topographic elevation data plane.

Description: A modified Kubelka-Munk model has been utilized to derive useful equations for the analysis of apparent canopy reflectance. Based on the solution to the model simple working equations were formulated by employing reflectance characteristic parameters. The relationships derived show the asymptotic nature of reflectance data that is typically observed in remote sensing studies of plant biomass. They also establish the range of expected apparent canopy reflectance values for specific plant canopy types. The usefulness of the simplified equations was demonstrated by the exceptionally close fit of the theoretical curves to two separately acquired data sets for alfalfa and shortgrass prairie canopies.

Description: Meteorological satellites and their capabilities are described. Future satellite configurations and instrumentation are discussed in the light of future user needs. In addition to the continuation of existing baseline products and services, the goals for improvement of the geosynchronous system through the 1990's will be: increasing spacial resolution in the visible and infrared channels; increasing vertical mean layer temperature resolution; adding the ability to image the solar disk; and upgrading the ground systems. Other improvements are discussed.

Description: The meteorological satellite had its substantive origin in the analytical process that helped initiate America's military satellite program. Its impetus lay in the desire to acquire current meteorological information in large areas for which normal meteorological observational data were not available on a day-to-day basis. Serious consideration was given to the feasibility of reconnaissance from meteorological satellites. The conceptualization of a meteorological satellite is discussed along with the early research which gave substance to that concept.

Description: Originally conceived to be the future operational weather satellite, the Nimbus series of seven satellites became the work horse and experimental backbone of the NASA/NOAA space research program. Early problems in developing the attitude control system and proving its operational characteristics in the simulated space environment are discussed. The program proved itself many times over; contributing substantially to the scientific knowledge of atmospherics and weather. As application/research vehicles, the seven Nimbus spacecraft were used for the development, test and application of a variety of new and advanced meteorological and geophysical remote sensing instruments and associated data transmission and processing techniques.

Description: The following working objectives were adopted: (1) define the current state of knowledge concerning the role of multispectral imaging science in hydrology; (2) identify critical areas where gaps in our knowledge limit opportunities for significant improvements in our understanding of the hydrologic processes; (3) evaluate the potential of multispectral imaging sciences as tools to close these gaps in knowledge; and (4) develop guidelines for a series of remote-sensing-based experiments that would help close these gaps in knowledge and, thereby, provide man with the improved scientific base necessary for better utilization of the world's water resource. The resulting documentation is intended to provide guidance for multispectral imaging programs in the hydrologic sciences with special emphasis on the visible and infrared (IR) wavelengths.

Description: Specific activity in information extraction science (taken to include data handling) is needed to: help identify the bounds of practical missions; identify potential data handling and analysis scenarios; identify the required enabling technology; and identify the requirements for a design data base to be used by the disciplines in determining potential parameters for future missions. It was defined that specific analysis topics were a function of the discipline involved, and therefore no attempt was made to define any specific analysis developments required. Rather, it was recognized that a number of generic data handling requirements exist whose solutions cannot be typically supported by the disciplines. The areas of concern were therefore defined as: data handling aspects of system design considerations; enabling technology for data handling, with specific attention to rectification and registration; and enabling technology for analysis. Within each of these areas, the following topics were addressed: state of the art (current status and contributing factors); critical issues; and recommendations for research and/or development.

Description: The use of multispectral imaging systems in land use, geomorphological, and cartographic applications is examined. Spatial, spectral, and geometric resolution requirements for photo interpretation and multispectral pattern recognition are discussed. The potential contributions of these systems in specific experiments is also included.

Description: The role that TM and MSS data play in the analysis of spatial patterns for land use/land cover, geomorphology studies, and the development of cartographic products is discussed. Tables listing geographic science data gaps, mission data requirements, and possible future remote sensing missions are included.

Description: One of the most important objectives of the International Satellite Cloud Climatology Project (ISCCP) is the quantification of cloud cover and associated radiation parameters for climate research. The time-varying nature of cloud cover requires a measurement system that obtains cloud radiance at sufficient time intervals for determining accurate values of cloudiness for the appropriate averaging period. Current plans for the ISCCP call for such day to be acquired every 3 hours during both data and night to ensure that the diurnal cycle of cloud cover is adequately sampled. The satellite system proposed is an array of geostationary and polar-orbiting satellites. The error in mean cloud cover estimates for observation systems which do not include geostationary satellites are quantified. Candidates for such a system include two NOAA Sun-synchronous satellites with daylight equatorial-crossing times of 0730 and 1430, and a USSR satellite in an 81 deg. inclination, 900-km altitude orbit.

Description: The diurnal variations of soil temperature and moisture content were simulated for a bare agricultural field in the San Joaquin Valley in California. The simulation pertained to the first 72 hours of drying, from saturation, of a sandy, clay loam soil. The results were compared with measurements of soil temperature and moisture content made at the field. Calculated and measured values of soil temperature trends agreed in general, but model results of moisture trends did not replicate observed diurnal effects evident at depths 4 centimeters or more below the surface.

Description: A mathematical model for the use of simultaneous Doppler derived correlated ranges in the geometric mode is presented. The model was tested with data taken during the EDOC-2 campaign with different integration intervals. The results of this adjustment are compared with the EDOC-2 adopted solution and those from an uncorrelated model used earlier to provide more economical calculations. It is shown that the correlated mode is superior to the uncorrelated one when the optimum integration interval of 23 seconds is used.

Description: (Previously announced in STAR as N82-22587)

Description: The degree of multiple scattering encountered by nadir-directed lidar in clear atmosphere conditions was assessed by the Monte Carlo method based on a model of the vertical distribution of aerosol scattering in the atmosphere. The lowest 3 km of the atmosphere were regarded as containing 90% of the aerosol optical thickness, with the unity normalized value for the forward peak of the aerosol phase function being 4.38 at a 694.3 nm laser wavelength. Results were obtained for two lidar receiver heights for three receiver field of view (FOV) angular halfwidths. An increase in receiver height was determined to cause a significant increase in the amount of angular scattering of the lidar signal. A factor of 20 change in receiver height produced an order of magnitude change in the single to multiple scattering ratio.

Description: The state of the art in remote sensing of the earth and the planets was discussed in terms of sensor performance, signal processing, and data interpretation. Particular attention was given to lidar for characterizing atmospheric particulates, the modulation of short waves by long ocean gravity waves, and runoff modeling for snow-covered areas. The use of NOAA-6 spacecraft AVHRR data to explore hydrologic land surface features, the effects of soil moisture and vegetation canopies on microwave and thermal microwave emissions, and regional scale evapotranspiration rate determination through satellite IR data are examined. A Shuttle experiment to demonstrate high accuracy global time and frequency transfer is described, along with features of the proposed Gravsat, radar image processing for rock-type discrimination, and passive microwave sensing of temperature and salinity in coastal zones.

Description: (Previously announced in STAR as N81-24494)

Description: Applications of Doppler radar to detection of storm hazards are reviewed. Normal radar sweeps reveal data on reflectivity fields of rain drops, ionized lightning paths, and irregularities in humidity and temperature. Doppler radar permits identification of the targets' speed toward or away from the transmitter through interpretation of the shifts in the microwave frequency. Wind velocity fields can be characterized in three dimensions by the use of two radar units, with a Nyquist limit on the highest wind speeds that may be recorded. Comparisons with models numerically derived from Doppler radar data show substantial agreement in storm formation predictions based on information gathered before the storm. Examples are provided of tornado observations with expanded Nyquist limits, gust fronts, turbulence, lightning and storm structures. Obtaining vertical velocities from reflectivity spectra is discussed.

Description: Geometric characteristics of the spherical earth are shown to be responsible for the increase of variance with latitude of zonally averaged meteorological statistics. An analytic model is constructed to display the effect of a spherical geometry on zonal averages, employing a sphere labeled with radial unit vectors in a real, stochastic field expanded in complex spherical harmonics. The variance of a zonally averaged field is found to be expressible in terms of the spectrum of the vector field of the spherical harmonics. A maximum variance is then located at the poles, and the ratio of the variance to the zonally averaged grid-point variance, weighted by the cosine of the latitude, yields the zonal correlation typical of the latitude. An example is provided for the 500 mb level in the Northern Hemisphere compared to 15 years of data. Variance is determined to increase north of 60 deg latitude.

Description: The decay curves of various triboluminescent-excited materials were obtained, including well-crystallized and poorly crystallized kaolin, bentonite, quartz, sodium chloride, and chalk calcite. A qualitative increase in triboluminescence was observed for kaolin dipped in water or tryptophan solution compared to dry kaolin, and for frozen kaolin and montmorillonite pastes. Theoretical explanations for the tryptophan effect are discussed.

Description: (Previously announced in STAR as N82-17715)

Description: Definitions of the data acquisition and interpretation goals which were sought in the Florida Area Cumulus Experiment (FACE) phase-II experiment as confirmation of results from FACE-I are presented. The FACE-1 trials in cloud seeding to increase rainfall were conducted in 1970, '71, '73, '75, and '76 and phase II proceeded from June 1978-August 1980. The data from the second set were to be treated with a block randomization scheme using the mean vector wind speed and direction and superposition of rainfall data, considering fixed and floating target clouds. A first level of probability was configured to eliminate cells from further investigation, and then two other levels are included in the analyses, each of continually increasing levels of likelihood. The instances of rainfall with seeding cover periods of either 60 or more or 60 or less seeding flares being ignited.

Description: Observations of cloud top height, backscattering, and signal depolarization have been obtained by a lidar system operating onboard a high-altitude research aircraft. The transmitter for the cloud lidar system is a doubled Nd:YAG laser operating at 5 Hz. The system functions as a fully automated sensor under microprocessor control and operates from a nominal 19-km altitude. Measurements have been acquired over a wide variety of cloud cover in conjunction with passive visible and infrared measurements. Initial observation results are reported

Description: The details of an iterative radiative transfer code for computing the intensity and degree of polarization of diffuse radiation in models of the ocean-atmosphere system are described. The present code neglects the upwelling radiation from below the ocean surface and as such can be applied to the part of the spectrum where the absorption by water is strong. To establish the reliability of the numerical scheme and the computer code, the results are compared with those of Fraser and Walker (1968), Dave (1972), and Mullamaa (1964); they are found to be in excellent agreement. The computations also show that both the intensity and the degree of polarization of the upwelling diffuse radiation at the top of the atmosphere vary significantly when the rough ocean at the base of the atmosphere is replaced by a Lambertian surface that reflects the same energy as the rough ocean.

Description: The goals and methods for global precipitation measurements from space are examined. Records currently exist only for visible and IR scans of cloud properties, and are applied from GEO for detecting diurnal variations in precipitation. Microwave radiometry is noted to be a suitable method for supplementing the visible and IR data for measuring stratiform oceanic precipitation, and when used at up to 3 microns can detect areas, if not amounts, of precipitation from GEO. Applications of radar altimeters are proposed in terms of modifications to the Seasat-type 2.2 cm radar, the use of surface target attentuation radar, of frequency agile rain radar, or of adaptive pointing radar. Soil moisture sensing is available with passive microwave radiometry in the 20-50 cm bands, or active sensing in the 5-8 cm bands. The utilization of GARP ground truth data is explored, along with statistical methods for treating the data samples.

Description: Scatterometer data from Seasat of the Amazon rain forest were examined to determine if the region is suitable to use as a reference for calibration of radars. The consistency of Skylab data viewing the Amazon region prompted the analyses, and the Seasat-A scatterometer system (SASS) gathered data of the same region at varying angles. The instrument employed a 100 W 14.6 GHz signal with the reflected power sampled 61 times during each 1.89 sec measurement period. Doppler filters were used in 15 parallel channels of reception, and represented areas 20 km by 50-70 km. Tests were made of regional and temporal stability of the Amazon area, with five measurement angles averaged at different incidence angles to find the mean deviation, which was found to be less than 0.5 dB. Diurnal effects were confined to early morning, and further tests are recommended to obtain results for an entire year, to develop screening methods for thick clouds and rain, and to determine the deviation more precisely.

Description: The effects of adopting new definitive precession and equinox corrections on the terrestrial reference frame was investigated. It is noted that: (1) the effect on polar motion is a diurnal periodic term with an amplitude increasing linearly in time whole on UT1 it is a linear term; (2) general principles are given to determine the effects of small rotations of the frame of a conventional inertial reference system (CIS) on the frame of the conventional terrestrial reference system (CTS); (3) seven CTS options are presented, one of which is necessary to accommodate such rotation. Accommodating possible future changes in the astronomical nutation is discussed. The effects of differences which may exist between the various CTS's and CIS's on Earth rotation parameters (ERP) and how these differences can be determined are examined. It is shown that the CTS differences can be determined from observations made at the same site, while the CIS differences by comparing the ERP's determined by the different techniques during the same time period.

Description: Data requirements, tasks to be accomplished, and the technical approaches to be used in identifying the characteristics of rice for crop inventories are listed as well as methods for estimating crop development and assessing its conditions.

Description: Directional reflectance factors of an alfalfa canopy were determined and related to canopy structure, agronomic variables, and irradiance conditions at four periods during a cutting cycle. Nadir and off-nadir reflectance factors decreased with increasing biomass in Thematic Mapper band 3(0.63-0.69 micrometer) and increased with increasing biomass in band 4(0.76-0.90 micrometer). The sensor view angle had less impact on perceived reflectance as the alfalfa progressed from an erectophile canopy of stems after harvest to a near planophile canopy of leaves at maturity. Studies of directional reflectance are needed for testing and upgrading vegetation canopy models and to aid in the complex interpretation problems presented by aircraft scanners and pointable satellites where illumination and viewing geometries may vary widely. Distinct changes in the patterns of radiance observed by a sensor as structural and biomass changes occur are keys to monitoring the growth and condition of crops.

Description: The balance of potential enstrophy and its relationship to vacillation cycles and the sudden warming is studied for a beta-channel model of the stratosphere. It is shown that the mean flow cannot be steady in the presence of large-amplitude quasi-geostrophic waves (approximately 1-0.25 geopotential kilometers /gpkm/) when any dissipation is present, and the maximum wave amplitude allowed is approximately 2 gpkm. If wave forcing (transience plus dissipation) is artificially maintained, the mean flow decelerates slowly at first then explosively as the potential vorticity gradient of the basic state is wiped out over the channel. This process is called wave saturation. The initial phase of the explosive deceleration resembles both the observed and modeled mean flow evolution during a sudden stratospheric warming. A simple vacillation model based upon these ideas shows remarkable similarity to the results of Holton and Mass (1976) and Davies (1981).