Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data segmentation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

Fluoropolymer Films Deposited by Argon Ion-Beam Sputtering of Polytetrafluoroethylene

The FT-IR, XPS and UV spectra of fluoropolymer films (SPTFE-I) deposited by argon ion-beam sputtering of polytetrafluoroethylene (PTFE) were obtained and compared with prior corresponding spectra of fluoropolymer films (SPTFE-P) deposited by argon rf plasma sputtering of PTFE. Although the F/C ratios for SPTFE-I and -P (1.63 and 1.51) were similar, their structures were quite different in that there was a much higher concentration of CF2 groups in SPTFE-I than in SPTFE-P, ca. 61 and 33% of the total carbon contents, respectively. The FT-IR spectra reflect that difference, that for SPTFE-I showing a distinct doublet at 1210 and 1150 per centimeter while that for SPTFE-P presents a broad, featureless band at ca. 1250 per centimeter. The absorbance of the 1210-per centimeter band in SPTFE-I was proportional to the thickness of the film, in the range of 50-400 nanometers. The SPTFE-I was more transparent in the UV than SPTFE-P at comparable thickness. The mechanism for SPTFE-I formation likely involves "chopping off" of oligomeric segments of PTFE as an accompaniment to "plasma" polymerization of TFE monomer or other fluorocarbon fragments generated in situ from PTFE on impact with energetic Ar ions. Data are presented for SPTFE-I deposits and the associated Ar(+) bombarded PTFE targets where a fresh target was used for each run or a single target was used for a sequence of runs.

Golub, Morton A.↗

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗

Aeroelastic Analysis of a Flexible Wing Wind Tunnel Model with Variable Camber Continuous Trailing Edge Flap Design

This paper presents data analysis of a flexible wing wind tunnel model with a variable camber continuous trailing edge flap (VCCTEF) design for drag minimization tested at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel test was designed to explore the relative merit of the VCCTEF concept for improved cruise efficiency through the use of low-cost aeroelastic model test techniques. The flexible wing model is a 10%-scale model of a typical transport wing and is constructed of woven fabric composites and foam core. The wing structural stiffness in bending is tailored to be half of the stiffness of a Boeing 757-era transport wing while the torsional stiffness is about the same. This stiffness reduction results in a wing tip deflection of about 10% of the wing semi-span. The VCCTEF is a multi-segment flap design having three chordwise camber segments and five spanwise flap sections for a total of 15 individual flap elements. The three chordwise camber segments can be positioned appropriately to create a desired trailing edge camber. Elastomeric material is used to cover the gaps in between the spanwise flap sections, thereby creating a continuous trailing edge. Wind tunnel data analysis conducted previously shows that the VCCTEF can achieve a drag reduction of up to 6.31% and an improvement in the lift-to-drag ratio (L=D) of up to 4.85%. A method for estimating the bending and torsional stiffnesses of the flexible wingUWAL wind tunnel model from static load test data is presented. The resulting estimation indicates that the stiffness of the flexible wing is significantly stiffer in torsion than in bending by as much as 9 to 1. The lift prediction for the flexible wing is computed by a coupled aerodynamic-structural model. The coupled model is developed by coupling a conceptual aerodynamic tool Vorlax with a finite-element model of the flexible wing via an automated geometry deformation tool. Based on the comparison of the lift curve slope, the lift prediction for the rigid wing is in good agreement with the estimated lift coefficients derived from the wind tunnel test data. Due to the movement of the VCCTEF during the wind tunnel test, uncertainty in the lift prediction due to the indicated variations of the VCCTEF deflection is studied. The results show a significant spread in the lift prediction which contradicts the consistency in the aerodynamic measurements, thus suggesting that the indicated variations as measured by the VICON system may not be reliable. The lift prediction of the flexible wing agrees very well with the measured lift curve for the baseline configuration. The computed bending deflection and wash-out twist of the flexible wing also match reasonably well with the aeroelastic deflection measurements. The results demonstrate the validity of the aerodynamic-structural tool for use to analyze aerodynamic performance of flexible wings.

Wind Tunnel↗

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

Augmented Method to Improve Thermal Data for the Figure Drift Thermal Distortion Predictions of the JWST OTIS Cryogenic Vacuum Test

The JWST Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and is scheduled for launch in 2018. The JWST OTE, including the 18 segment primary mirror, secondary mirror, and the Aft Optics Subsystem (AOS) are designed to be passively cooled and operate near 45K. These optical elements are supported by a complex composite backplane structure. As a part of the structural distortion model validation efforts, a series of tests are planned during the cryogenic vacuum test of the fully integrated flight hardware at NASA JSC Chamber A. The successful ends to the thermal-distortion phases are heavily dependent on the accurate temperature knowledge of the OTE structural members. However, the current temperature sensor allocations during the cryo-vac test may not have sufficient fidelity to provide accurate knowledge of the temperature distributions within the composite structure. A method based on an inverse distance relationship among the sensors and thermal model nodes was developed to improve the thermal data provided for the nanometer scale WaveFront Error (WFE) predictions. The Linear Distance Weighted Interpolation (LDWI) method was developed to augment the thermal model predictions based on the sparse sensor information. This paper will encompass the development of the LDWI method using the test data from the earlier pathfinder cryo-vac tests, and the results of the notional and as tested WFE predictions from the structural finite element model cases to characterize the accuracies of this LDWI method.

Thermal Modeling↗

A High Performance Computing Approach to Tree Cover Delineation in 1-m NAIP Imagery Using a Probabilistic Learning Framework

Tree cover delineation is a useful instrument in deriving Above Ground Biomass (AGB) density estimates from Very High Resolution (VHR) airborne imagery data. Numerous algorithms have been designed to address this problem, but most of them do not scale to these datasets, which are of the order of terabytes. In this paper, we present a semi-automated probabilistic framework for the segmentation and classification of 1-m National Agriculture Imagery Program (NAIP) for tree-cover delineation for the whole of Continental United States, using a High Performance Computing Architecture. Classification is performed using a multi-layer Feedforward Backpropagation Neural Network and segmentation is performed using a Statistical Region Merging algorithm. The results from the classification and segmentation algorithms are then consolidated into a structured prediction framework using a discriminative undirected probabilistic graphical model based on Conditional Random Field, which helps in capturing the higher order contextual dependencies between neighboring pixels. Once the final probability maps are generated, the framework is updated and re-trained by relabeling misclassified image patches. This leads to a significant improvement in the true positive rates and reduction in false positive rates. The tree cover maps were generated for the whole state of California, spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles. The framework produced true positive rates of around 88% for fragmented forests and 74% for urban tree cover areas, with false positive rates lower than 2% for both landscapes. Comparative studies with the National Land Cover Data (NLCD) algorithm and the LiDAR canopy height model (CHM) showed the effectiveness of our framework for generating accurate high-resolution tree-cover maps.

Segments↗

Unpaired image translation to mitigate domain shift in liquid argon time projection chamber detector responses

Deep learning algorithms often are developed and trained on a training dataset and deployed on test datasets. Any systematic difference between the training and a test dataset may severely degrade the final algorithm performance on the test dataset—what is known as the domain shift problem . This issue is prevalent in many scientific domains where algorithms are trained on simulated data but applied to real-world datasets. Typically, the domain shift problem is solved through various domain adaptation (DA) methods. However, these methods are often tailored for a specific downstream task, such as classification or semantic segmentation, and may not easily generalize to different tasks. This work explores the feasibility of using an alternative way to solve the domain shift problem that is not specific to any downstream algorithm. The proposed approach relies on modern Unpaired Image-to-Image (UI2I) translation techniques, designed to find translations between different image domains in a fully unsupervised fashion. In this study, the approach is applied to a domain shift problem commonly encountered in Liquid Argon Time Projection Chamber (LArTPC) detector research when seeking a way to translate samples between two differently distributed LArTPC detector datasets deterministically. This translation allows for mapping real-world data into the simulated data domain where the downstream algorithms can be run with much less domain-shift-related performance degradation. Conversely, using the translation from the simulated data to a real-world domain can increase the realism of the simulated dataset and reduce the magnitude of any systematic uncertainties. To evaluate the quality of the translations, we use both pixel-wise metrics and a downstream task to measure the effectiveness of UI2I methods for mitigating the domain shift problem. We adapted several popular UI2I translation algorithms to work on scientific data and demonstrated the viability of these techniques for solving the domain shift problem with LArTPC detector data. To facilitate further development of DA techniques for scientific datasets, the ‘Simple Liquid-Argon Track Samples’ dataset used in this study is also published.

97 MATHEMATICS AND COMPUTING↗

Time Domain Tool Validation Using ARES I-X Flight Data

The ARES I-X vehicle was launched from NASA's Kennedy Space Center (KSC) on October 28, 2009 at approximately 11:30 EDT. ARES I-X was the first test flight for NASA s ARES I launch vehicle, and it was the first non-Shuttle launch vehicle designed and flown by NASA since Saturn. The ARES I-X had a 4-segment solid rocket booster (SRB) first stage and a dummy upper stage (US) to emulate the properties of the ARES I US. During ARES I-X pre-flight modeling and analysis, six (6) independent time domain simulation tools were developed and cross validated. Each tool represents an independent implementation of a common set of models and parameters in a different simulation framework and architecture. Post flight data and reconstructed models provide the means to validate a subset of the simulations against actual flight data and to assess the accuracy of pre-flight dispersion analysis. Post flight data consists of telemetered Operational Flight Instrumentation (OFI) data primarily focused on flight computer outputs and sensor measurements as well as Best Estimated Trajectory (BET) data that estimates vehicle state information from all available measurement sources. While pre-flight models were found to provide a reasonable prediction of the vehicle flight, reconstructed models were generated to better represent and simulate the ARES I-X flight. Post flight reconstructed models include: SRB propulsion model, thrust vector bias models, mass properties, base aerodynamics, and Meteorological Estimated Trajectory (wind and atmospheric data). The result of the effort is a set of independently developed, high fidelity, time-domain simulation tools that have been cross validated and validated against flight data. This paper presents the process and results of high fidelity aerospace modeling, simulation, analysis and tool validation in the time domain.

Hough, Steven↗

Cloud and Precipitation Analyses using Merged Datasets from Two Airborne Microwave Radiometers Covering 10–183 GHz

Microwave radiometers provide valuable insight into the structure and characteristics of clouds and precipitation. In NASA’s airborne remote-sensing arsenal, the Advanced Microwave Precipitation Radiometer (AMPR) and the Conical Scanning Millimeter-wave Imaging Radiometer (CoSMIR) have been used extensively in field campaigns throughout the world. AMPR operates with four channels between 10.7 and 85.5 GHz, while CoSMIR operates with nine channels ranging from 50.3 to 183.31 GHz. Although these datasets provide key information when used separately, the combination of these radiometers covers virtually the full range of frequencies used by the Global Precipitation Measurement (GPM) Microwave Imager (GMI), enabling suborbital observations to compare with GPM spaceborne measurements. This presentation will detail the merger of AMPR and CoSMIR data during two NASA airborne field campaigns: the Integrated Precipitation and Hydrology Experiment (IPHEx) and the Olympic Mountains Experiment / Radar Definition Experiment (OLYMPEX/RADEX). Other field campaigns, such as the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS), may be included as well. Using the merged brightness temperature dataset, features of precipitating and non-precipitating clouds containing liquid and/or ice hydrometeors will be discussed from selected flight segments. Geophysical retrievals derived from these brightness temperatures using a one-dimensional variational (1DVAR) technique and/or multi-linear regression equations will also be employed in these analyses. Observed transitions between precipitating and non-precipitating systems will be explored in greater detail. Dropsonde data will be used to provide environmental contexts throughout each flight, and additional observations (e.g., from airborne and/or land-based radar) will be incorporated to supplement the radiometer-based results. The broader implications of these results and pathways for future work will also be discussed.

Corey G. Amiot↗

Wavefront Sensing of X-Ray Telescopes

Phase Retrieval analysis of off-axis or defocused focal-plane data from telescope optics has been proven effective in understanding misalignments and optical aberrations in normal incidence telescopes. The approach is used, e.g., in commissioning of the James Webb Space Telescope (JWST) segmented primary mirror. There is a similar need for evaluating low-order figure errors of grazing incidence mirrors and nested telescope assemblies. When implemented in these systems, phase retrieval does not depend on normal incidence access to each mirror (shell) surface and, therefore, provides an effective means for evaluating nested x-ray telescopes during integration and test. We have applied a well-known phase retrieval algorithm to grazing incidence telescopes. The algorithm uses the Levenberg-Marquardt optimization procedure to perform a non-linear least-squares fit of the telescope Point Spread Function (PSF). The algorithm can also retrieve low order figure errors at visible wavelengths where optical diffraction is the dominant defect in the PSF. In this paper we will present the analytical approach and its implementation for grazing incidence mirrors of the International X-Ray Observatory (IXO). We analyze the effects of low order axial surface errors individually, and in combination on the system PSF at 633 nanometers. We demonstrate via modeling that the wavefront sensing algorithm can recover axial errors (of the grazing incidence mirrors) to a small fraction of the known axial figure errors using simulated PSFs as input data to the algorithm.

Saha, Timo↗

Data from "A Bayesian Record Linkage Approach to Applications in Tree Demography Using Overlapping LiDAR Scans"

Processed LiDAR data and environmental covariates from 2015 and 2019 LiDAR scans in the Vicinity of Snodgrass Mountain (Western Colorado, USA), in a geographic subset used in primary analysis for the research paper.This package contains LiDAR-derived canopy height maps for 2015 and 2019, crown polygons derived from the height maps using a segmentation algorithm, and environmental covariates supporting the model of forest growth. Source datasets include August 2015 and August 2019 discrete-return LiDAR point clouds collected by Quantum Geospatial for terrain mapping purposes on behalf of the Colorado Hazard Mapping Program and the Colorado Water Conservation Board. Both datasets adhere to the USGS QL2 quality standard. The point cloud data were processed using the R package lidR to generate a canopy height model representing maximum vegetation height above the ground surface, using a pit-free algorithm.This dataset was compiled to assess how spatial patterns of tree growth in montane and subalpine forests are influenced by water and energy availability. Understanding these growth patterns can provide insight into forest dynamics in the Southern Rocky Mountains under changing climatic conditions.This dataset contains .tif, .csv, and .txt files. This dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗

Architecture of a catalytically active homotrimeric plant cellulose synthase complex

Cellulose is an essential plant cell wall component and represents the most abundant biopolymer on Earth. Supramolecular plant cellulose synthase complexes organize multiple linear glucose polymers into microfibrils as load-bearing wall components. We determined the structure of a poplar cellulose synthase CesA homotrimer that suggests a molecular basis for cellulose microfibril formation. This complex, stabilized by cytosolic plant-conserved regions and helical exchange within the transmembrane segments, forms three channels occupied by nascent cellulose polymers. Secretion steers the polymers toward a common exit point, which could facilitate protofibril formation. CesA’s N-terminal domains assemble into a cytosolic stalk that interacts with a microtubule-tethering protein and may thus be involved in CesA localization. Here, our data suggest how cellulose synthase complexes assemble and provide the molecular basis for plant cell wall engineering.

59 BASIC BIOLOGICAL SCIENCES↗

Automated Scoring of Morphological Changes in Images of Pentaerythritol Tetranitrate

Recent advances in characterization techniques that generate large datasets of material microstructure images require robust, automated image-processing. We applied an unsupervised anomaly detection method called feature anomaly detection system (FADS) to automatically detect and quantify microstructure changes in images of the explosive pentaerythritol tetranitrate (PETN) aged at various temperatures. We demonstrated the FADS approach on two-dimensional images extracted from computed tomography scans, but the same technique can be readily applied to other imaging modalities. FADS calculates anomaly scores on the basis of differences in filter activations of nominal and test data in pretrained convolutional neural networks. The FADS scores successfully differentiated between pristine PETN and PETN aged at a temperature where material coarsening occurred. Morphological metric analysis of segmented images verified observed trends in FADS scores as a function of aging temperature and aging time, specifically by calculating volume fractions, specific boundary lengths, two-point correlation functions, and local thicknesses. Here, the FADS technique has two important advantages compared to traditional morphological analysis: First, it uses grayscale images as input, rather than images that are segmented to separate the appropriate phases; and second, FADS scores capture any type of changes among image sets, rather than requiring prior knowledge or selection of a relevant set of metrics.

Accelerated aging↗

DECAL MDN Resolution Calibration

The simulated pixelated calorimeter uses 0.1 mm × 0.1 mm silicon pixels (0.0001 cm²) silicon pixels as active layers — far finer than current concepts like HGCAL (~0.5 cm²) — improving energy resolution through much finer segmentation. While the standard resolution formalism fits three numbers to a handful of discrete test-beam energies, this work is a first try at a more data-efficient alternative: learning the full response shape continuously in energy with a mixture density network. This continuous, differentiable surrogate is a natural building block for fast simulation of extremely complex calorimeters.

Wang, Peter [U. Chicago (main)]↗

Proposed methods for the reconstruction of the Pioneer/Venus probe entry trajectories

In connection with the Pioneer Venus Project, an orbiter and an array of atmospheric probes are sent to Venus in separate launches during 1978. The four probes, which are transported to the planet on an instrumented Bus are released approximately 20 days prior to encounters. There are three identical small probes and a large probe. Each is designed to make in-situ atmospheric measurements. One of the principal scientific objectives of the Pioneer Venus (P/V) multiprobe mission is the determination of the structure and composition of the Venusian atmosphere down to the planet surface. To enhance the science return it is of paramount importance to be able to accurately determine the spacecraft entry trajectories. In addition to the in-situ spacecraft measurements, earth-based Doppler data are available. The process felt to be most applicable for the P/V multiprobe entry trajectory reconstruction is considered, with emphasis placed on the highly dynamic upper entry segment of the trajectory.

Findlay, J. T.↗

Spatially Characterizing Effective Timber Supply

The structure of a computer-oriented cartographic model for assessing roundwood supply for generation of base load electricity is discussed. The model provides an analytical procedure for coupling spatial information of harvesting economics and owner willingness to sell stumpages. Supply is characterized in terms of standing timber; of accessibility considering various harvesting and hauling factors; and of availability as affected by ownership and residential patterns. Factors governing accessibility to timber include effective harvesting distance to haulic roads as modified by barriers and slopes. Haul distance is expressed in units that take into account the relative ease of travel along various road types to a central processing facility. Areas of accessible timber are grouped into spatial units, termed 'timbersheds', of common access to particular haul road segments that belong to unique 'transport zones'. Timber availability considerations include size of ownership parcels, housing density and excluded areas. The analysis techniques are demonstrated for a cartographic data base in western Massachusetts.

Berry, J. K.↗