Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mapping variables”

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 775 records · Page 43

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spacetime quantum mechanics for bosonic and fermionic systems

We provide a Hilbert space approach to quantum mechanics where space and time are treated on an equal footing. Our approach replaces the standard dependence on an external classical time parameter with a spacetime-symmetric algebraic structure, thereby unifying the axioms that traditionally distinguish the treatment of spacelike and timelike separations. Standard quantum evolution can be recovered from timelike correlators, defined by means of a quantum action operator, a quantum version of the action of classical mechanics. The corresponding map also provides an alternative perspective on the path integral formulation that, in the case of fermions, does not require the use of Grassmann variables. In addition, the formalism can be interpreted in terms of generalized quantum states, codifying both the conventional information of a quantum system at a given time and its evolution. We show that these states are solutions to a quantum principle of stationary action grounded in timelike correlations and pseudo-entropies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

RatXcan: A framework for cross-species integration of genome-wide association and gene expression data

Genome-wide association studies (GWAS) have implicated specific alleles and genes as risk factors for numerous complex traits. However, translating GWAS results into biologically and therapeutically meaningful discoveries remains extremely challenging. Most GWAS results identify noncoding regions of the genome, suggesting that differences in gene regulation are the major driver of trait variability. To better integrate GWAS results with gene regulatory polymorphisms, we previously developed PrediXcan (also known as “transcriptome-wide association studies” orTWAS), which maps SNPs to predicted gene expression using GWAS data. In this study, we developed RatXcan, a framework that extends this methodology to outbred heterogeneous stock (HS) rats. RatXcan accounts for the close familial relationships among HS rats by modeling the relatedness with a random effect that encodes the genetic relatedness. RatXcan also corrects for polygenic-driven inflation because of the equivalence between a relatedness random effect and the infinitesimal polygenic model. To develop RatXcan, we trained transcript predictors for 8,934 genes using reference genotype and expression data from five rat brain regions. We found that the cis genetic architecture of gene expression in both rats and humans was sparse and similar across brain tissues. We tested the association between predicted expression in rats and two example traits (body length and BMI) using phenotype and genotype data from 5,401 densely genotyped HS rats and identified a significant enrichment between the genes associated with rat and human body length and BMI. Thus, RatXcan represents a valuable tool for identifying the relationship between gene expression and phenotypes across species and paves the way to explore shared biological mechanisms of complex traits.

Genetics & Heredity↗

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Subsatellite measurements of the lunar magnetic field

The Apollo 15 subsatellite magnetometer data have been used to map the lunar magnetic field over a narrow band of the lunar surface. Within this band the magnetic field is generally stronger and more variable over the farside highlands than the nearside maria. The correspondence between the strong variable lunar field regions and the source regions for limb compressions suggests that limb compressions arise as the result of the deflection of the solar wind just upstream of the terminator by the lunar magnetic field. Using this apparent relationship between field strength and limb compression source regions, it is deduced that the field strength in the northern farside highlands is not as strong as in the southern hemisphere at similar longitudes. Simultaneous measurements of the interplanetary magnetic field obtained by Explorer 35 and the Apollo 15 subsatellite above the dayside hemisphere are essentially identical. Thus, both instruments are measuring the undisturbed interplanetary field.

Russell, C. T.↗

Estimation of thermal flux and emissivity of the land surface from multispectral aircraft data

In order to evaluate the importance of surface thermal flux and emissivity variations on surface and boundary layer processes, a technique that uses thermal data from an airborne multispectral scanner to determine the surface skin temperature and thermal emissivity over a regional area has been developed. These values are used to estimate the total flux density emanating from the surface and at the top of the atmosphere. Data from the multispectral atmospheric mapping sensor (MAMS) collected during the First ISLSCP Field Experiment (FIFE) are used to develop the technique, and to show the time and space variability of the flux values. The ground truth data available during FIFE provide a unique resource to evaluate this technique.

Jedlovec, Gary J.↗

Calibration, navigation, and registration of MAMS data for FIFE

The International Satellite Land Surface Climatology Project (ISLSCP) was conducted to study the interaction of the atmosphere with the land surface and the research problems associated with the interpretation of satellite data over the Earth's land surface. The experimental objectives of the First ISLSCP Field Experiment (FIFE) were the simultaneous acquisition of satellite, atmospheric, and surface data and to use these data to understand the processes controlling energy/mass exchange at the surface. The experiment site is a 15 x 15 km area southeast of Manhattan, Kansas, intersected by Interstate 70 and Kansas highway 177. The Konza Prairie portion is 5 x 5 km and is a controlled experiment site consisting primarily of native tall grass prairie vegetation. The remainder of the site is grazing and farm land with trees along creek beds that are scattered over the area. Airborne multispectral imagery from the Multispectral Atmospheric Mapping Sensor (MAMS) was collected over this region on two days during Intensive Field Campaign-1 (1FC-1) to study the time and space variability of remotely-sensed geophysical parameters. These datasets consist of multiple overflights covering about a 60-min period during late morning on June 4, 1987 and shortly after dark on the following day. Image data from each overpass were calibrated and Earth located with respect to each other using aircraft inertial navigation system parameters and ground control points. These were the first MAMS flights made with 10-bit thermal data.

Jedlovec, G. J.↗

Diffuse 0.511 MeV line emission and the distribution of positron annihilation in the galaxy

We model the distribution in the Galaxy of 0.511 MeV line emission and the implied positron annihilation using plausible distributions based on observations at other photon energies. We use 0.511 MeV line observations from the Galactic center and from directions away from the Galactic center to normalize the two-dimensional 0.511 MeV sky maps. We find that in order to understand all of the available data it is necessary to invoke the presence of a time-variable component of 0.511 MeV line emission in addition to an underlying diffuse component. The large 0.511 MeV fluxes observed with broad field of view detectors, such as the SMM gamma-ray spectrometer, can be reconciled with the small fluxes seen with narrow field-of-view instruments (GRIS and OSSE) from directions away from the Galactic center if the diffuse emission follows the distribution of Galactic novae for which the recently detected 0.511 MeV line emission from the Galactic center with OSSE is predominantly of diffuse origin. For future GRIS and HEXAGONE observations from the direction of the Galactic center we predict a minimum 0.511 MeV line flux of (6.0 +/- 0.9) x 10 exp -4 photons/sq cm/s.

Skibo, J. G.↗

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗

Bulk Insolation Models as Predictors for Locations for High Lunar Hydrogen Concentrations

In this study we consider the bulk effects of surface illumination on topography (insolation) and the possible thermodynamic effects on the Moon's hydrogen budget. Insolation is important as one of the dominant loss processes governing distributions of hydrogen volatiles on the Earth, Mars and most recently Mercury. We evaluated three types of high latitude > 65 deg., illumination models that were derived from the Lunar Observing Laser Altimetry (LOLA) digital elevation models (DEM)'s. These models reflect varying accounts of solar flux interactions with the Moon's near-surface. We correlate these models with orbital collimated epithermal neutron measurements made by the Lunar Exploration Neutron Detector (LEND). LEND's measurements derive the Moon's spatial distributions of hydrogen concentration. To perform this analysis we transformed the topographic model into an insolation model described by two variables as each pixels 1) slope and 2) slope angular orientation with respect to the pole. We then decomposed the illumination models and epithermal maps as a function of the insolation model and correlate the datasets.

Hydrogen Concentrations↗

Preliminary assessment of the accuracy and precision of TOPEX/POSEIDON altimeter data with respect to the large-scale ocean circulation

TOPEX/POSEIDON sea surface height measurements are examined for quantitative consistency with known elements of the oceanic general circulation and its variability. Project-provided corrections were accepted but are at tested as part of the overall results. The ocean was treated as static over each 10-day repeat cycle and maps constructed of the absolute sea surface topography from simple averages in 2 deg x 2 deg bins. A hybrid geoid model formed from a combination of the recent Joint Gravity Model-2 and the project-provided Ohio State University geoid was used to estimate the absolute topography in each 10-day period. Results are examined in terms of the annual average, seasonal average, seasonal variations, and variations near the repeat period. Conclusion are as follows: the orbit error is now difficult to observe, having been reduced to a level at or below the level of other error sources; the geoid dominates the error budget of the estimates of the absolute topography; the estimated seasonal cycle is consistent with prior estimates; shorter-period variability is dominated on the largest scales by an oscillation near 50 days in spherical harmonics Y(sup m)(sub 1)(theta, lambda) with an amplitude near 10 cm, close to the simplest alias of the M(sub 2) tide. This spectral peak and others visible in the periodograms support the hypothesis that the largest remaining time-dependent errors lie in the tidal models. Though discrepancies attribute to the geoid are within the formal uncertainties of the good estimates, removal of them is urgent for circulation studies. Current gross accuracy of the TOPEX/POSEIDON mission is in the range of 5-10 cm, distributed overbroad band of frequencies and wavenumbers. In finite bands, accuracies approach the 1-cm level, and expected improvements arising from extended mission duration should reduce these numbers by nearly an order of magnitude.

Wunsch, Carl↗

ALMA Spectral Imaging of Titan Contemporaneous with Cassini's Grand Finale

The Cassini mission performed 127 targeted ybys of Titan during its 13 yr mission to Saturn, culminating in the Grand Finale between 2017 April and September. Here we demonstrate the use of the Atacama Large Millimeter/ submillimeter Array (ALMA) to continue Cassini's legacy for chemical and climatological studies of Titan's atmosphere. Whole-hemisphere, interferometric spectral maps of HCN, HNC, HC3N, CH3CN, C2H3CN, C2H5CN,and C3H8 were obtained using ALMA in 2017 May at moderate (0 2, or 1300 km) spatial resolution, revealing the effects of seasonally variable chemistry and dynamics on the distribution of each species. The ALMA submillimeter observations of HCN and HC3N are consistent with Cassini infrared data on these species, obtained in the same month. Chemical/dynamical lifetimes of a few years are inferred for C2H3CN and C2H5CN, in reasonably close agreement with the latest chemical models incorporating the sticking of C2H5CN to stratospheric aerosol particles. ALMA radial limb ux proles provide column density information as a function of altitude, revealing maximum abundances in the thermosphere (above 600 km) for HCN, HNC, HC3N, and C2H5CN. This constitutes the rst detailed measurement of the spatial distribution of HNC, which is found to be conned predominantly to altitudes above 730 ± 60 km. The HNC emission map shows an east-west hemisphericasymmetry of 13% ± 3%. These results are consistent with very rapid production (and loss) of HNC in Titan's uppermost atmosphere, making this molecule an effective probe of short-timescale (diurnal) ionospheric processes.

Cordiner, M. A.↗

Canopy Cover Predictions using Ground Observations and Remotely Sensed Data

Maps of vegetation status are needed at many scales, from the field level to monitor ecosystem condition to the global level to understand the carbon cycle. Status is quantified by such variables as leaf area index, biomass, and fraction of canopy cover. Current methods of predicting vegetation variables use remote sensing data to provide a spatially exhaustive data source. In a study in western Montana, several hundred ground observations made by the US Forest Service on tenth-acre conifer plots were used to develop aspatial regression and geostatistical prediction models. Normalized Difference Vegetation Index (NDVI) values from Landsat Thematic Mapper images were used as ancillary data. These models were then used to predict canopy cover at unsampled locations in a 97 square kilometer region on the boundary of the Flathead National Forest and the Bob Marshall Wilderness. Independent data from two dates six years apart were used for validation. Given the assumption that actual canopy cover remained relatively unchanged within this time period, partial validation can be achieved by measuring the correspondence of the two maps. This criterion results in ranking the aspatial regression maps as less accurate than the geostatistically generated maps. The geostatistical approach emphasizes ground measurements more heavily than does aspatid regression. Geostatistical simulations of canopy cover also provide a means of describing uncertainty about the patterns of canopy cover.

Dungan, Jennifer L.↗

Multivariate Bayesian Optimization of CoO Nanoparticles for CO 2 Hydrogenation Catalysis

The hydrogenation of CO 2 holds promise for transforming the production of renewable fuels and chemicals. However, the challenge lies in developing robust and selective catalysts for this process. Transition metal oxide catalysts, particularly cobalt oxide, have shown potential for CO 2 hydrogenation, with performance heavily reliant on crystal phase and morphology. Achieving precise control over these catalyst attributes through colloidal nanoparticle synthesis could pave the way for catalyst and process advancement. Yet, navigating the complexities of colloidal nanoparticle syntheses, governed by numerous input variables, poses a significant challenge in systematically controlling resultant catalyst features. We present a multivariate Bayesian optimization, coupled with a data-driven classifier, to map the synthetic design space for colloidal CoO nanoparticles and simultaneously optimize them for multiple catalytically relevant features within a target crystalline phase. The optimized experimental conditions yielded small, phase-pure rock salt CoO nanoparticles of uniform size and shape. These optimized nanoparticles were then supported on SiO 2 and assessed for thermocatalytic CO 2 hydrogenation against larger, polydisperse CoO nanoparticles on SiO 2 and a conventionally prepared catalyst. The optimized CoO/SiO 2 catalyst consistently exhibited higher activity and CH 4 selectivity (ca. 98%) across various pretreatment reduction temperatures as compared to the other catalysts. This remarkable performance was attributed to particle stability and consistent H* surface coverage, even after undergoing the highest temperature reduction, achieving a more stable catalytic species that resists sintering and carbon occlusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Luxury iron uptake and storage in pennate diatoms from the equatorial Pacific Ocean

Abstract Iron is a key micronutrient for ocean phytoplankton, and the availability of iron controls primary production and community composition in large regions of the ocean. Pennate diatoms, a phytoplankton group that responds to iron additions in low-iron areas, can have highly variable iron contents, and some groups such as Pseudo-nitzschia, are known to use ferritin to store iron for later use. We quantified and mapped the intracellular accumulation of iron by a natural population of Pseudo-nitzschia from the Fe-limited equatorial Pacific Ocean. A total of 48 h after iron addition, nearly half of the accumulated iron was localized in storage bodies adjacent to chloroplasts believed to represent ferritin. Over the subsequent 48 h, stored iron was distributed to the rest of the cell through subsequent growth and division, partially supporting the iron contents of the daughter cells. This study provides the first quantitative view into the cellular trafficking of iron in a globally relevant phytoplankton group and demonstrates the unique capabilities of synchrotron-based element imaging approaches.

59 BASIC BIOLOGICAL SCIENCES↗

Nested active learning for efficient model contextualization and parameterization: pathway to generating simulated populations using multi-scale computational models

There is increasing interest in the use of mechanism-based multi-scale computational models (such as agent-based models (ABMs)) to generate simulated clinical populations in order to discover and evaluate potential diagnostic and therapeutic modalities. The description of the environment in which a biomedical simulation operates (model context) and parameterization of internal model rules (model content) requires the optimization of a large number of free parameters. In this work, we utilize a nested active learning (AL) workflow to efficiently parameterize and contextualize an ABM of systemic inflammation used to examine sepsis. Contextual parameter space was examined using four parameters external to the model’s rule set. The model’s internal parameterization, which represents gene expression and associated cellular behaviors, was explored through the augmentation or inhibition of signaling pathways for 12 signaling mediators associated with inflammation and wound healing. We have implemented a nested AL approach in which the clinically relevant (CR) model environment space for a given internal model parameterization is mapped using a small Artificial Neural Network (ANN). The outer AL level workflow is a larger ANN that uses AL to efficiently regress the volume and centroid location of the CR space given by a single internal parameterization. We have reduced the number of simulations required to efficiently map the CR parameter space of this model by approximately 99%. In addition, we have shown that more complex models with a larger number of variables may expect further improvements in efficiency.

97 MATHEMATICS AND COMPUTING↗

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗

Experimental philosophy leading to a small scale digital data base of the conterminous United States for designing experiments with remotely sensed data

Research using satellite remotely sensed data, even within any single scientific discipline, often lacked a unifying principle or strategy with which to plan or integrate studies conducted over an area so large that exhaustive examination is infeasible, e.g., the U.S.A. However, such a series of studies would seem to be at the heart of what makes satellite remote sensing unique, that is the ability to select for study from among remotely sensed data sets distributed widely over the U.S., over time, where the resources do not exist to examine all of them. Using this philosophical underpinning and the concept of a unifying principle, an operational procedure for developing a sampling strategy and formal testable hypotheses was constructed. The procedure is applicable across disciplines, when the investigator restates the research question in symbolic form, i.e., quantifies it. The procedure is set within the statistical framework of general linear models. The dependent variable is any arbitrary function of remotely sensed data and the independent variables are values or levels of factors which represent regional climatic conditions and/or properties of the Earth's surface. These factors are operationally defined as maps from the U.S. National Atlas (U.S.G.S., 1970). Eighty-five maps from the National Atlas, representing climatic and surface attributes, were automated by point counting at an effective resolution of one observation every 17.6 km (11 miles) yielding 22,505 observations per map. The maps were registered to one another in a two step procedure producing a coarse, then fine scale registration. After registration, the maps were iteratively checked for errors using manual and automated procedures. The error free maps were annotated with identification and legend information and then stored as card images, one map to a file. A sampling design will be accomplished through a regionalization analysis of the National Atlas data base (presently being conducted). From this analysis a map of homogeneous regions of the U.S.A. will be created and samples (LANDSAT scenes) assigned by region.

Labovitz, M. L.↗