Search NASA⌕ Search

SEARCH · Search NASA

Results for “Topographic correction”

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.

124 records · Page 7

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean and at Landfall

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1' half-power width (two-way) across the aircraft ground track over a swath equal to 0.8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the off-nadir angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two dimensional FFT, and Doppler corrected. The open-ocean data were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving toward 330 deg at about 5 m/s. Individual waves up to 18 m height were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction. At some positions there were three different wave fields of comparable energy crossing each other. The NOAA aircraft spent over five hours within 180 km of the eye, and made five eye penetrations. On 26 August 1998, the NOAA aircraft flew at 2.2 km height when hurricane Bonnie was making landfall near Wilmington, NC, documenting the directional wave spectrum in the region between Charleston, SC and Cape Hatteras, NC. The aircraft flight lines included segments near and along the shoreline as well as far offshore. Animations of the directional wave spectrum spatial variation along the aircraft tracks on the two flights will be presented using a 100: 1 time compression.

Walsh, Edward J.↗

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean and at Landfall

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1 E half-power width (two-way) across the aircraft ground track over a swath equal to 0.8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the incidence angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two-dimensional FFT, and Doppler corrected. The open-ocean data were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving slowly to the north. Individual waves with heights up to 18 m were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction. At some positions there were three different wave fields of comparable energy crossing each other. The NOAA aircraft spent over five hours within 180 km of the hurricane Bonnie eye, and made five eye penetrations. A 3-minute animation of the directional wave spectrum spatial variation over this period will be shown as well as summary plots of the wave field spatial variation. On 26 August 1998, the NOAA aircraft flew at 2.2 km height when hurricane Bonnie was making landfall near Wilmington, NC, documenting the directional wave spectrum in the region between Charleston, SC and Cape Hatteras, NC. The aircraft ground track included both segments along the shoreline and Pamlico Sound as well as far offshore. An animation of the directional wave spectrum spatial variation at landfall will be presented and contrasted with the spatial variation when Bonnie was in the open ocean on 24 August 1998.

Walsh, E. J.↗

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1 deg half-power width (two-way) across the aircraft ground track over a swath equal to 0. 8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the incidence angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two-dimensional FFT, and Doppler corrected. The data presented were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving slowly to the north. Wave heights up to 18 m were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction and at times there were wave fields traveling at right angles to each other. The NOAA aircraft spent over five hours within 180 km of the hurricane Bonnie eye, and made five eye penetrations. A 2-minute animation of the directional wave spectrum spatial variation over this period will be shown.

Wright, C. W.↗

Geostatistical evaluation of satellite radar altimetry for high-resolution mapping of Lambert Glacier, Antarctica

The potential of satellite radar altimetry for high-resolution mapping of Antarctic ice streams is evaluated, using retracked and slope-corrected data from the Lambert Glacier and Amery Ice Shelf area, East Antarctica, acquired by Geosat during the Exact Repeat Mission (ERM), 1986-89. The map area includes lower Lambert Glacier north of 72.18 deg S, the southern Amery Ice Shelf, and the grounded inland ice sheet on both sides. The Geosat ERM altimetry is found to provide substantially more complete coverage than the 1978 Seasat altimetry, due to improved tracking. Variogram methods are used to estimate the noise levels in the data as a function of position throughout the map area. The spatial structure in the data is quantified by constructing experimental variograms using altimetry from the area of the grounding zone of Lambert Glacier, which is the area chiefly of interest in this topographically complex region. Kriging is employed to invert the along-track height measurements onto a fine-scale 3 km grid. The unsmoothed along-track Geosat ERM altimetry yields spatially continuous maps showing the main topographic features of lower Lambert Glacier, upper Amery Ice Shelf and the adjacent inland ice sheet. The probable position of the grounding line of Lambert Glacier is identified from a break in slope at the grounded ice/floating ice transition. The approximate standard error of the kriged map is inferred from the data noise levels.

Herzfeld, Ute C.↗

Thermomechanical evolution of oceanic fracture zones

To explain the topographic and gravitational signals of fracture zones (FZs), a realistic model is developed and compared with Seasat altimeter data across three major FZs. The model is based on the thermal, subsidence, and mechanical models of the oceanic lithosphere. The model is used to calculate the FZ topography and the evolution of the geoid signature across a FZ. When lithospheric flexure is included, the model explains the persistence of sharp FZ geoid steps on old seafloor as well as the asymmetry of the geoid steps across major FZs. The model does not explain the disappearance of the thermal part of the FZ geoid signature. The model profiles are used to access the various techniques for estimating the overall geoid offsets across FZs. It is shown that steps may be underestimated by as much as 50 percent. These predictions, if correct, will force reinterpretation of published geoid height versus age data derived from altimeter profiles across FZs.

Sandwell, D. T.↗

LANDSAT-D investigations in snow hydrology

Thematic mapper data for the southern Sierra Nevada area were registered to digital topographic data and compared to LANDSAT MSS and NOAA-7 AVHRR data of snow covered areas in order to determine the errors associated with using coarser resolution and to qualify the information lost when high resolution data are not available. Both the zenith and the azimuth variations in the radiative field are considered in an atmospheric radiative transfer model which deals with a plane parallel structured atmosphere composed of different layers, each assumed to be homogeneous in composition and to have a linear in tau temperature profile. Astronomical parameters for each layer are Earth-Sun distance and solar flux at the top of the atmosphere. Atmospheric parameters include pressure temperature, chemical composition of the air molecules, and the composition and size of the aerosol, water droplets, and ice crystals. Outputs of the model are the monochromatic radiance and irradiance at each level. The use of the model in atmospheric correction of remotely sensed data is discussed.

Dozier, J.↗

Depth and Distribution of CO2 Snow on Mars

The dynamic role of volatiles on the surface of Mars has been a subject of longstanding interest. In the pre-Viking era, much of the debate was necessarily addressed by theoretical considerations. A particularly influential treatment by Leighton and Murray put forth a simple model relying on solar energy balance, and led to the conclusion that the most prominent volatile exchanging with the atmosphere over seasonal cycles is carbon dioxide. Their model suggested that due to this exchange, atmospheric CO2 partial pressure is regulated by polar ice. While current thinking attributes a larger role to H2O ice than did the occasional thin polar coating this model predicted, the CO2 cycle appears to be essentially correct. There are a number of observational constraints on the seasonal exchange of surface volatiles with the atmosphere. The growth and retreat of polar CO2 frost is visible from Earth-based telescopes and from spacecraft in Mars orbit, both at visible wavelengths and in thermal IR properties of the surface. Recently, variations in Gamma ray and neutron fluxes have also been used to infer integrated changes in CO2 mass on the surface. Measurements made by Viking's Mars Atmospheric Water Detector experiment were sensitive to atmospheric H2O vapor abundance. Surface condensates and their transient nature were detected by the Viking landers. The study here is motivated by recent data collected by the Mars Global Surveyor, affording the opportunity to not only detect the lateral distribution of volatiles, but also to constrain the variable volumes of the reservoirs. We elaborate on a technique first employed by Smith et al. By examining averages of a large number of topographic measurements collected by the Mars Orbiter Laser Altimeter (MOLA), that study showed that the zonal pattern of deposition and sublimation of CO2 can be determined. In their first approach, reference surfaces were fit to all measurements in narrow latitude annuli, and the time dependent variations about those mean surfaces were examined. In their second approach, height measurements from pairs of tracks that cross on the surface were interpolated and differenced, forming a set of crossover residuals. These residuals were then examined as a function of time and latitude. The initial studies averaged over longitude to maximize signal and minimize noise in order to isolate the expected small signal. In this follow-up study we now attempt to extract the elevation change pattern also as a function of longitude, and we focus on the crossover approach.

Aharonson, Oded↗

Topographic slope from the SAR interferometric phase gradient

A new algorithm for the direct calculation of topographic slope maps from synthetic aperture radar (SAR) interferograms is presented. The algorithm derives slope maps without first requiring the creation of a digital elevational model (DEM) from the interferogram, thus obviating the need for high SNR in the interferogram and altitude calibration points for the scene. SAR data useful for interferometry has been collected by the Active Microwave Imager on board the ERS-1 satellite, when it was in a short period repeat orbit. Two passes of the radar sensor form a cross-track interferometric baseline. For a point target at some position (x,y,h), the interferometric phase difference phi is proportional to the difference in path lengths for the two sensor positions to the scatterer. Given the phase difference as measured in the complex interferogram and an accurate baseline geometry, the position of the scatterer, most significantly the height h, can be determined through triangulation. The interferometric phase measurement however is known only modulo 2-pi, and hence it is necessary to determine the correct multiple of 2-pi to add to the phase at each point to obtain an estimate of the actual phase with respect to an absolute datum. This phase unwrapping process is required for creating DEM's, and is difficult or impossible for regions of low SNR or SAR image layover. The new algorithm described here derives slope maps without requiring phase unwrapping.

Werner, Charles L.↗

Characterization of Martian near-subsurface materials by determination of cohesion and angle of internal friction

Back-analysis (reconstruction) of the stability of thirty avalanche chutes was performed in the very limited areas where high resolution imaging overlapped with available 1:500 K topographic map coverage. A new technique was developed to incorporate the third dimension (width) of an avalanche chute in stability back-analysis in order to yield unambiguous values of cohesion and angle of internal friction. The procedure is based upon extending the ordinary method of slices to three dimensions, in order to construct avalanche chute cross-sections whose widths and depths vary as a function of gradient, gravity, density of material, and phi and c. Applying the technique to the well documented slide at Lodalen, Norway as a test produces excellent correspondence with reality. Generally, the technique reveals that the width:depth ratio of any avalanche chute decreases with increasing contrast between the average slope angle and the angle of internal friction. Applying this technique to the martian avalanche chute yields results consistent with indications from earlier work, but with greater certainty. Values of cohesion and angle of internal friction identify the materials at the time of failure as moderately cohesive debris. If Sharp's identification of these features as avalanche chutes is correct, then the results here imply that weathering processes have had a significant effect to depths of tens of meters (where failure has occured) below the martian surface. It is also implied that on relatively steep slopes within Valles Marineris, sizable, unaltered, unmantled bedrock exposures for high resolution spectral and spatial scanning by Mars Observer may be scarce.

Sullivan, R. J.↗

Geocoding of AIRSAR/TOPSAR SAR Data

It has been demonstrated and recognized that radar interferometry is a promising method for the determination of digital elevation information and terrain slope from Synthetic Aperture Radar (SAR) data. An important application of Interferometric SAR (InSAR) data in areas with topographic variations is that the derived elevation and slope can be directly used for the absolute radiometric calibration of the amplitude SAR data as well as for scattering mechanisms analysis. On the other hand polarimetric SAR data has long been recognized as permitting a more complete inference of natural surfaces than a single channel radar system. In fact, imaging polarimetry provides the measurement of the amplitude and relative phase of all transmit and receive polarizations. On board the NASA DC-8 aircraft, NASA/JPL operates the multifrequency (P, L and C bands) multipolarimetric radar AIRSAR. The TOPSAR, a special mode of the AIRSAR system, is able to collect single-pass interferometric C- and/or L-band VV polarized data. A possible configuration of the AIRSAR/TOPSAR system is to acquire single-pass interferometric data at C-band VV polarization and polarimetric radar data at the two other lower frequencies. The advantage of this system configuration is to get digital topography information at the same time the radar data is collected. The digital elevation information can therefore be used to correctly calibrate the SAR data. This step is directly included in the new AIRSAR Integrated Processor. This processor uses a modification of the full motion compensation algorithm described by Madsen et al. (1993). However, the Digital Elevation Model (DEM) with the additional products such as local incidence angle map, and the SAR data are in a geometry which is not convenient, since especially DEMs must be referred to a specific cartographic reference system. Furthermore, geocoding of SAR data is important for multisensor and/or multitemporal purposes. In this paper, a procedure to geocode the new AIRSAR/TOPSAR data is presented. As an example an AIRSAR/TOPSAR image acquired in 1994 is geocoded and evaluated in terms of geometric accuracy.

Holecz, Francesco↗

Implementation of the NCAR Community Land Model (CLM) in the NASA/NCAR finite-volume Global Climate Model (fvGCM)

In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).

Radakovich, Jon D.↗

The Geographic Distribution of Boulder Halo Craters at Mid-to-High Latitudes on Mars

Extensive evidence exists for ground ice at mid-to-high latitudes on Mars, including results from neutron spectroscopy [1-3], thermal properties [4-5], geomorphology [e.g., 6-9], and the in situ observations of Mars Phoenix [10]. This ground ice has been hypothesized to be emplaced diffusively and fill pores [11], or to have accumulated by ice and dust deposition that draped or mantled the terrain [7, 12]. These two processes are not mutually exclusive; both potentially have occurred on Mars [5]. One of the landforms found in areas where ground ice is common on Mars are boulder halo craters [e.g., 13-15] (Figure 1), which are topographically muted impact craters that are filled by ice-rich regolith. They are outlined by boulders that trace a circular outline of the original crater rim. Boulder halos generally have distinctly higher boulder densities than the surrounding background plains and have few boulders in their interiors. The mechanism of boulder halo crater formation is somewhat uncertain. Our working model is that an impact event occurs with sufficient size to excavate to a depth greater than the boulder-poor, ice-rich soils. Excavated boulders are deposited around the crater's rim and in its proximal ejecta. Quite rapidly [14], the crater becomes infilled by icy soil. Rather than being buried, boulders in the halo remain at the surface, perhaps be-cause they 'float' relative to finer-grained materials [14, 16]. Regardless of the details of this process, the life-time of boulders at the surface is much greater than the timescale needed to remove most of the craters' topography. Physical weathering of rocks must be greatly out-paced by crater infilling (the opposite of what is typical, e.g., on the Moon [17]). The rapidity of this infilling is easiest to understand if icy mantling material is deposited and accumulates, rather than simply being added by pore filling of soils. If this model is correct, boulder halos only form when they excavate rock-producing materials from beneath the upper surface. Thus, the distribution and size of craters that result in boulders halos may provide in-sight into the thickness of the ice-rich surface layer in different locations. Note that this thickness is necessarily that of the ice-rich layer at the time of impact, not at present. This study is an initial survey of boulder halo crater locations in the 50deg to 80degN and 50deg to 80degS latitude bands on Mars.

Ice↗

Temporal and spatial variability in surface roughness and accumulation rate around 88° S from repeat airborne geophysical surveys

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88° S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability (<10 km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88° S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r(exp 2)=−0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r(exp 2)=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

Michael Studinger↗

Validation of ICESat-2 ATLAS Bathymetry and Analysis of ATLAS’s Bathymetric Mapping Performance

NASA’s Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) was launched in September, 2018. The satellite carries a single instrument, ATLAS (Advanced Topographic Laser Altimeter System), a green wavelength, photon-counting lidar, enabling global measurement and monitoring of elevation with a primary focus on the cryosphere. Although bathymetric mapping was not one of the design goals for ATLAS, pre-launch work by our research team showed the potential to map bathymetry with ICESat-2, using data from MABEL (Multiple Altimeter Beam Experimental Lidar), NASA’s high-altitude airborne ATLAS emulator, and adapting the laser-radar equation for ATLAS specific parameters. However, many of the sensor variables were only approximations, which limited a full assessment of the bathymetric mapping capabilities of ICESat-2 during pre-launch studies. Following the successful launch, preliminary analyses of the geolocated photon returns have been conducted for a number of coastal sites, revealing several salient examples of seafloor detection in water depths of up to ~40 m. The geolocated seafloor photon returns cannot be taken as bathymetric measurements, however, since the algorithm used to generate them is not designed to account for the refraction that occurs at the air–water interface or the corresponding change in the speed of light in the water column. This paper presents the first early on-orbit validation of ICESat-2 bathymetry and quantification of the bathymetric mapping performance of ATLAS using data acquired over St. Thomas, U.S. Virgin Islands. A refraction correction, developed and tested in this work, is applied, after which the ICESat-2 bathymetry is compared against high-accuracy airborne topo-bathymetric lidar reference data collected by the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA). The results show agreement to within 0.43—0.60 m root mean square error (RMSE) over 1 m grid resolution for these early on-orbit data. Refraction-corrected bottom return photons are then inspected for four coastal locations around the globe in relation to Visible Infrared Imaging Radiometer Suite (VIIRS) K(d)(490) data to empirically determine the maximum depth mapping capability of ATLAS as a function of water clarity. It is demonstrated that ATLAS has a maximum depth mapping capability of nearly 1 Secchi in depth for water depths up to 38 m and K(d)(490) in the range of 0.05–0.12 per m. Collectively, these results indicate the great potential for bathymetric mapping with ICESat-2, offering a promising new tool to assist in filling the global void in nearshore bathymetry.

ICESat-2↗

Discovery of a Plains Caldera Complex and Extinct Lava Lake in Arabia Terra, Mars: Implications for the Discovery of Additional Highland Volcanic Source Regions

Several irregularly shaped topographic depressions occur near the dichotomy boundary in northern Arabia Terra, Mars. The geomorphology of these features suggests that they formed by collapse, opposed to meteor impact. At least one depression (approx.55 by 85 km) displays geologic features indicating a complex, multi-stage collapse history. Features within and around the collapse structure indicate volcanic processes. The complex occurs within Hesperian ridged plains of likely volcanic origin and displays no crater rim or evidence for ejecta. Instead the depression consists of a series of circumferential graben and down-dropped blocks which also display upper surfaces similar to ridged plain lavas. Large blocks within the depression are tilted towards the crater center, and display graben that appear to have originally been linked with circumferential graben outside of the complex related to earlier collapse events. A nearly 700 m high mound exists along a graben within the complex that might be a vent. The deepest depression displays two sets of nearly continuous terraces, which we interpret as high-stands of a drained lava lake. These features appear similar to the black ledge described during the Kilauea Iki eruption in 1959. A lacustrine origin for the terraces seems unlikely because of the paucity of channels found in or around the depression that could be linked to aqueous surface processes. In addition, there is no obvious evidence for lacustrine sediments within the basin. Together with the presence of significant faulting that is indicative of collapse we conclude that this crater complex represents a large caldera formed in the Late Noachian to Early Hesperian. Other linear and irregular depressions in the region also might be linked to ancient volcanism. If that hypothesis is correct, it suggests that northern Arabia Terra could contain a large, previously unrecognized highland igneous province. Evacuation of magma via explosive and effusive activity produced localized collapse, might have contributed to nearby ridged plains, and pyroclastic materials erupted from these vents might have supplied sediments in fretted terrain and other deposits. The recognition of volcanoes within Arabia Terra expands the known extent of Noachian-Hesperian volcanism to cover much of the preserved martian highland crust.

Bleacher, Jacob↗

Orbit Determination for Mars Global Surveyor During Mapping

The Mars Global Surveyor (MGS) spacecraft reached a low-altitude circular orbit on February 4, 1999, after the termination of the second phase of aerobraking. The MGS spacecraft carries the Mars Orbiter Laser Altimeter (MOLA) whose primary goal is to derive a global, geodetically referenced 0.2 deg x 0.2 deg topographic grid of Mars with a vertical accuracy of better than 30 meters. During the interim science orbits in the' Hiatus mission phase (October - November 1997), and the Science Phasing Orbits (March - April, 1998, and June - July 1998) 208 passes of altimeter data were collected by the MOLA instrument. On March 1, 1999 the first ten orbits of MOLA altimeter data from the near-circular orbit were successfully returned from MGS by the Deep Space Network (DSN). Data will be collected from MOLA throughout the Mapping phase of the MCS mission, or for at least one Mars year (687 days). Whereas the interim orbits of Hiatus and SPO were highly eccentric, and altimeter data were only collected near periapsis when the spacecraft was below 785 km, the Mapping orbit of MGS is near circular, and altimeter data will be collected continuously at a rate of 10 Hz. The proper analysis of the altimeter data requires that the orbit of the MGS spacecraft be known to an accuracy comparable to that of the quality of the altimeter data. The altimeter has an ultimate precision of 30 cm on mostly flat surfaces, so ideally the orbits of the MGS spacecraft should be known to this level. This is a stringent requirement, and more realistic goals of orbit error for MGS are ten to thirty meters. In this paper we will discuss the force and measurement modelling required to achieve this objective. Issues in force modelling include the proper modelling of the gravity field of Mars, and the modelling of non-conservatives forces, including the development of a 'macro-model', in a similar fashion to TOPEX/POSEIDON and TDRSS. During Cruise and Aerobraking, the high gain antenna (HGA) was stowed on the +X face of the spacecraft. On March 29, 1999 the HGA will be deployed on a meter long boom which will remain Earth-pointing while the instrument panel (including the MOLA instrument) remains pointed at nadir. The tracking data must be corrected for the regular motion of the high gain antenna with respect to the center of mass, and the success of the MGS determination during Mapping will depend on correctly accounting for this offset in the measurement model.

Lemoine, F. G.↗