Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth profiling”

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 397 records · Page 22

Hydrographic data from R/V endeavor cruise #90

The final cruise of the NSF sponsored Warm Core Rings Program studied a Warm Core Ring (WCR) in the Fall of 1982 as it formed from a large northward meander of the Gulf Stream. This ring, known as 82-H or the eighth ring identified in 1982, formed over the New England Seamounts near 39.5 deg N, 65 deg W. Surveys using Expendable Bathythermographs, Conductivity-Temperature-Depth-Oxygen stations and Doppler Current Profiling provide a look at the genesis of a WCR. These measurements reveal that WCR 82-H separated from the Gulf Stream sometime between October 2-5. This ring was a typical WCR with a diameter of about 200 km and speeds in the high velocity core of the 175 cm/sec. Satellite imagery of 82-H following the cruise showed that it drifted WSW in the Slope Water region at almost 9 km/day, had at least one interaction with the Gulf Stream and was last observed on February 8, 1983 at 39 deg N, 72 deg W.

Stalcup, M. D.↗

Graben and extension in northern Tharsis, Mars

A morphometric analysis of north-south trending graben at 35 deg N in northern Tharsis is used to provide the first quantitative measurements of crustal extension in and strain for areas of the Martian crust that have undergone both extreme and mild deformation. These estimates also demonstrate the viability of the Viking data for making such quantitative assessment. Graben in the region accommodate about 8 km of post Early Amazonian east-west extension. This extension corresponds to a net regional strain of 0.45-0.61 percent; local strains vary from less than 1 percent to 5 percent. Extension is nonuniform and localized near 110 deg (south of Alba Patera and north of Ceraunius Fossae) because of reactivation of preexisting fractures in underlying basement and possibly due to superposition of stress fields from Tharsis and Alba Patera. A major mechanical discontinuity in the shallow crust, probably representing the base of or a strength contrast within the megaregolith, is observed. It occurs at a depth of 2 km on the profile flanks and is depressed to about 7.5 km near 110 deg due to loading of the crust by the volcano.

Plescia, J. B.↗

A synoptic of H-alpha line profile in the T Tauri star SU Aurigae

We present a catalog of 106 high spectral resolution observations of the H-alpha line profile in the T Tauri star SU Aurigae, obtained during the period from 1986 October through 1990 November. The spectra were acquired during joint synoptic programs to observe selected T Tauri stars using the Hamilton Echelle Spectrometer of the Lick Observatory and the solar-stellar spectrograph at the McMath telescope of the National Solar Observatory on Kitt Peak. A restricted set of Mg II h and k line profiles was also obtained in a coordinated program involving the International Ultraviolet Explorer (IUE) satellite observatory and the McMath solar-stellar facility. Striking variability is evident on a nightly basis. A key result is that the relative intensity in the blue wing of H-alpha spanning a range of velocities bear -150 km/s is modulated at a period of 2.98 +/- 0.4 days. We identify the 2.98 day period with the rotation period of the star. We also find that the occurrence of the periodic modulation of the mass outflow is episodic and most evident during a 2 week sequence of nightly observations. We find two other intervals where the periodic spectroscopic variability is likely present, although at a lower level of significance at a lower level of significance. The variability is otherwise stochastic in nature. The Mg II resonance lines exhibit clear variability that is most pronounced in the blue wing of the k line. A comparison of the Mg II k line profile with H-alpha profiles obtained nearly simultaneous yields no apparent correlation between the variable features in each line. The profile shapes of the Mg II h and k lines are generally indicative of formation in a wind. An analysis of the principal features that appear in the H-alpha profile set suggests that the line is composed of contributions from an enhanced chromosphere; a relatively slow moving, dense, optically thick component of a stellar wind formed relatively close to the star; and an optically thin, high-velocity, expanding stellar wind located further away from the star. An investigation of possible correlations among the principal features in the series of H-alpha profiles suggests that as the density in the wind increases, the wind may become more unstable to large turbulence. This may lead to a reduction in the wind bulk velocity, thus regulating the mass-loss rate. We also find that the position of the main absorption feature which is always present in the H-alpha profiles is not correlated with its depth, indicating that optical depth and wind velocity are not correlated in the denser portions of the wind.

Giampapa, Mark S.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗

Intercomparison of Atmospheric Radiation Schemes for the Lower Martian Atmosphere

Interpretation of the vast data flow from the recent and ongoing Mars missions requires extensive modeling activity to support the investigations and conclusions drawn from the observations. We have initiated a campaign to compare atmospheric radiation schemes for the lower Martian atmosphere with leastcompromise reference line-by-line calculations (Meadows and Crisp 1996, Crisp 1986). By using fixed reference conditions we are introducing an intercomparison framework to enable atmospheric modeling teams to compare their modeling results with the reference line-by-line calculations, and to improve their models. The atmosphere of Mars is thin and therefore its response to local radiative forcing is strong when compared to that of the Earth s troposphere. Airborne dust and also the water vapor give significant contributions to the radiative fluxes. Thus models of the Martian atmosphere should have fairly accurate radiative transfer (RT) algorithms. This is specifically important for the net radiative flux at the surface, which is the main driver for the surface temperature evolution, since the surface turbulent fluxes are small due to the thin atmosphere (e.g. Haberle et al. 1993, Savij arvi 1991a, Savij arvi 1991b, Savij arvi 1999). However, the accuracy of the RT schemes has not been rigorously tested in many cases as reference calculations have been rare. We introduce a set of reference cases for Martian atmospheric conditions with profiles for temperature, pressure, dust optical depth (Conrath 1975, Clancy et al. 1995, OckertBell et al. 1997, Smith et al. 2000), water, CO2,O3,CO and O2 (Table 1). The reference atmosphere uses a diurnally and globally averaged thermal structure derived from Mariner 9 IRIS observations acquired during late southern summer. Thermal structures for nighttime and daytime conditions will be included shortly. Computationally expensive reference simulations were performed by applying a spectrum resolving (line-by-line multiple scattering) model, SRM, to give the least-compromise base for comparisons. We call for other groups to join in future comparisons. This resembles the International Comparison of Radiation Codes for Climate Models (ICRCCM; e.g. Ellingson and Fouquart 1991, Ellingson et al. 1991). The ICRCCM led to many improvements in the Earth GCM radiation codes. The intercomparison framework presented in this paper is a step toward that direction for Mars. We have already five modeling teams participating in this intercomparison exercice. Additional teams investigating the Martian atmosphere are cordially invited to join this effort of comparing the radiation schemes in fixed conditions. The results achieved will be reported and some additional features to be included in the intercomparison will be discussed.

A-M Harri↗

Cirrus Clouds Optical, Microphysical and Radiative Properties Observed During Crystal-Face Experiment: I. A Radar-Lidar Retrieval System

A method of retrieving cloud microphysical properties using combined observations from both cloud radar and lidar is introduced. This retrieval makes use of an improvement to the traditional optimal estimation retrieval method, whereby a series of corrections are applied to the state vector during the search for an iterative solution. This allows faster convergence to a solution and is less processor intensive. The method is first applied to a synthetic cloud t o demonstrate its validity, and it is shown that the retrieval reliably reproduces vertical profiles of ice water content. The retrieval method is then applied to radar and lidar observations from the CRYSTAL-FACE experiment, and vertical profiles of ice crystal diameter, number concentration, and ice water content are retrieved for a cirrus cloud layers observed one day of that experiment. The validity of the relationship between visible extinction coefficient and radar reflectivity was examined. While synthetic tests showed such a functional relationship, the measured data only partially supported such a conclusion. This is due to errors in the forward model (as explained above) as well as errors in the data sets, including possible mismatch between lidar and radar profiles or errors in the optical depth. Empirical relationships between number concentrations and mean particle diameter were also examined. The results indicate that a distinct and robust relationship exists between these retrieved quantities and it is argued that such a relationship is more than an artifact of the retrieval process offering insight into the nature of the microphysical processes taking place in cirrus.

Mitrescu, C.↗

Global Estimates of PBL Depth from Space-Borne LIDAR

The planetary boundary layer (PBL) is responsible for communicating the exchange of energy, moisture, momentum, pollutants, and aerosols between the surface and the free atmosphere and is therefore crucial to many studies of the atmosphere. Unfortunately, there have historically been few observations of this important layer due to the complexity involved in its measurement. However, with the advent of more advanced satellites, global measurements of the PBL are now becoming possible. The PBL is often characterized by a high concentration of aerosols within the layer and low level clouds capping it and these are observable from space. The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite was launched in 2006 with the intention of observing aerosols and clouds and was the first space-based LIDAR optimized for this purpose. CALIPSO observations are therefore well suited to observing the depth of the PBL. Since it was launched, CALIPSO has been making nearly continuous measurements enabling a global picture of PBL depth. We plan to present a global PBL depth product and how it evolves throughout the year. The product is able to identify deeper PBL depths in the summer hemisphere over land and deeper depths along the northern hemisphere oceanic storm tracks in winter associated with cold air traveling over warm water. Large seasonal cycles are also evident in the subtropical desert locations among other features. In addition, comparisons will be made between several estimates of PBL depth based on turbulent intensity, meteorology profiles, and aerosol profiles from the GEOS5 model.

McGrath-Spangler, Erica lynn↗

Ice-Over-Water Cloud Identification in an Artificial Neural Network Approach

An artificial neural network (ANN) algorithm, employing several Aqua MODIS channels, the retrieved cloud phase and total cloud visible optical depth, and temperature and humidity vertical profiles is trained to detect multilayer (ML) ice-over-water cloud systems identified by matched 2008 CloudSat and CALIPSO (CC) data. The trained MLANN was applied to 2009 MODIS data resulting in combined ML and single layer detection accuracies of 87% (89%) and 86% (89%) for snow-free (snow-covered) regions during the day and night, respectively. When corrected for the viewing-zenith-angle dependence of each parameter, the ML fraction detected is relatively invariant across the swath. Compared to the CC ML variability, the MLANN is robust seasonally and interannually, and produces similar distribution patterns over the globe, except in the polar regions. Additional research is needed to conclusively evaluate the VZA dependence and further improve the MLANN accuracy. This approach should greatly improve the monitoring of cloud vertical structure using operational passive sensors.

MODIS↗

The Representation of Aerosols in GMAO’s Newest Reanalyses

Over the past few years, NASA’s Global Modelling and Assimilation Office has been working on the configuration and production of three new reanalysis products, GEOS-IT, GiOcean, and MERRA-21C. GEOS-IT, or the Goddard Earth Observing System for Instrument Teams, is a 3D variational data assimilation system that runs in a near real time framework however retrospectively provides data back through 1998 to deliver a consistent view of the Earth-atmosphere system for the production of observational NASA products. Retrospective production for GEOS-IT is complete and the meteorology has since been used to produce the one way weakly coupled GiOcean reanalysis. Due to differences in the atmospheric model, particularly related to scavenging, aerosols are not identical in GEOS-IT and GiOcean. MERRA-21C, or the Modern Era Restrospective analysis for Research and Applications in the 21st century, is a hybrid 4D ensemble variational system at a finer horizontal resolution of 0.25 degrees. Although different in their intended use, and therefore configuration, these systems prominently feature coupling between meteorology and aerosols. The differences and similarities in the aerosol configuration between the three systems will be discussed, covering biomass burning and anthropogenic emissions as well as observations used for the assimilation of aerosol optical depth. A large emphasis will be placed on the version of the underlying aerosol module, GOCART, which underwent a complete refactoring and the addition of radiatively active brown carbon between GEOS-IT and MERRA-21C. Independent observations will be used to evaluate the performance of aerosols in both reanalyses, focusing on aerosol optical depth, surface particulate matter, and vertical profiles of aerosol backscatter.

Allison Collow↗

Comparison of satellite-inferred (SAGE II) aerosol optical depths with corresponding airborne sunphotometer optical depths

Solar radiation measurements were made with an airborne multiwavelength autotracking sunphotometer detector off the northeastern coast of Brazil on April 6, 1985. These measurements are used to validate solar radiation measurements from the Stratospheric Aerosol and Gas Experiment (SAGE II). The SAGE II data below 17.5 km are dominated by attenuation due to thunderstorm-generated cirrus clouds. When integrated above 17.5 km, however, the aerosol extinction profiles from SAGE II yield aerosol optical depths that are consistent with the corresponding sunphotometer values above 10.1 km.

Livingston, John M.↗

Apollo 16 core 60013/14 as a product of path 1 and path 2 regolith evolution processes

We have made a petrographic study of 12 samples along the 62 cm profile of Apollo 16 core 60013/14 taken at Station 10 near the lunar module. The basic core description and maturity variations have been discussed earlier. This core is important because it enables a stratigraphic comparison with nearby cores 60001-7 and 60009/10. These cores form a triangle with sides 35-40 m. FMR data show a maturity-match between 60009/10 and 60013/14. Unfortunately, segments of the core 60001-7 are too disturbed for any cm scale comparison. We present new data on the grain size distribution and the modal petrography for six samples from six levels of 60013 using the same approach and rationale used for 60014. We have combined the new data with our previous data on 60014 to plot the variations of grain size distributions and grain type abundances, in the 500-1000 m and 90-150 micrometers size fractions, with depth. Agglutinate abundances match the I sub s/FeO profile and both show a significant break at a depth of about 45 cm, with lower maturity material below that point. Grain size distributions show that for all analyzed samples, the mean grain size is below 90 micrometers. The two lowest samples are somewhat more coarse-grained than the rest of the core. The break at about 45 cm is the most prominent feature of the core and was described during the initial dissection as a pronounced change in color from dark to light. It now seems clear that this break is primarily one of maturity rather than of bulk chemical composition.

Basu, A.↗

Handheld White Light Interferometer for Measuring Defect Depth in Windows

Accurate quantification of defects (scratches and impacts) is vital to the certification of flight hardware and other critical components. The amount of damage to a particular component contributes to the performance, reliability, and safety of a system, which ultimately affects the success or failure of a mission or test. The launch-commit criteria on a Space Shuttle Orbiter window are governed by the depth of the defects that are identified by a visual inspection. This measurement of a defect is not easy to obtain given the environment, size of the defect, and location of the window(s). The determination of depth has typically been performed by taking a mold impression and measuring the impression with an optical profiling instrument. Another method of obtaining an estimate of the depth is by using a refocus microscope. To use a refocus microscope, the surface of the glass and bottom of the defect are, in turn, brought into focus by the operator. The amount of movement between the two points corresponds to the depth of the defect. The refocus microscope requires a skilled operator and has been proven to be unreliable when used on Orbiter windows. White light interferometry was chosen as a candidate to replace the refocus microscope. The White Light Interferometer (WLI) was developed to replace the refocus microscope as the instrument used for measuring the depth of defects in Orbiter windows. The WLI consists of a broadband illumination source, interferometer, detector, motion control, displacement sensor, mechanical housing, and support electronics. The illumination source for the WLI is typically a visible light emitting diode (LED) or a near-infrared superluminescent diode (SLD) with power levels of less than a milliwatt. The interferometer is a Michelson configuration consisting of a 1-in. (2.5-cm) cube beam splitter, a 0.5-in. (1.3-cm) optical window as a movable leg (used to closely match the return intensity of the fixed leg from the window), and a mirrored prism to fold the optics into the mechanical housing. The detector may be one of many C-mount CCD (charge-coupled device) cameras. Motion is provided by a commercial nanostepping motor with a serial interface. The displacement sensor is a custom device specifically designed for this application. The mechanical housing and support electronics were designed to integrate the various components into an instrument that could be physically handled by a technician and easily transported.

Youngquist, Robert↗

Scattering of emitted radiation from inhomogeneous and nonisothermal layers

A parametric study is performed for the exiting monochromatic intensities scattered from finite plane-parallel inhomogeneous layers that are driven solely by a distribution of thermal sources. Intensities are obtained by invariantly imbedding the standard and thermal scattering functions. The single-scattering albedo and the Henyey-Greenstein phase-function parameter are varied independently, and both linear and exponential profiles are considered. Linear temperature profiles are used, including temperature inversions. The resulting intensities, as a function of the direction cosine of propagation, are discussed from a remote-sensing point of view. For an isothermal and homogeneous medium, the gross characteristics of the exiting intensity, represented by its overall slope, mean value (magnitude), and an interior maximum value, can be related to the total optical depth, single-scattering albedo, and phase function, respectively. For a homogeneous medium, linearly decreasing (in the line of sight) temperature profiles tend to obscure the phase-function information and decrease the apparent optical depth. On the other hand, linearly increasing temperature profiles tend to retain phase-function information and increase the apparent optical depth. Temperature inversion profiles give intensities very similar to those for purely linear profiles.

Bergstrom, R. W.↗

The Analysis and Modeling of the Deployment of NASA's X-38 Parafoil

NASA is currently working on developing the X-38 Crew Return Vehicle for the Internati6nal Space Station (ISS). The aerodynamics of the X-38 and its enormous 7500 square ft parafoil are extremeJy complicated and must be understood and modeled accurately before the vehicle is actually used for human space flight. A large amount of analysis has been performed on the steady state regions of the flight profile when the parafoil is already open. However, the deployment region of the flight profile has not been analyzed in much depth and is a particularly complicated dynamic flight period. During deployment, the X-38 system has its parafoil open up to its full area in five stages over a period of 20-30 seconds. The vehicle transitions from angles of attack around 90 deg to 0 deg in about two seconds and is extremely dynamic. This paper will detail the analysis work done to generate accurate flight profile characteristics of the deployment from flight data from vehicle test drops. Additionally, a computer simulation tool for steady state flight was modified to analyze and model the parafoil deployment. Eventually, a computer model was produced that generated accurate representations of the parafoil deployment dynamics for the vehicle test drops to date.

Ranier, Chris↗

Surface profile changes of scuffed bearing surfaces

A phase locked interference microscope capable of resolving depth differences to 30 A and planar displacements of 6000 A was constructed for the examination of the profiles of bearing surfaces without physical contact. This instrument was used to determine surface chemical reactivity by applying a drop of dilute alcoholic hydrochloric acid and measuring the profile of the solid surface before and after application of this probe. Scuffed bearing surfaces reacted much faster than unscuffed ones, but bearing surfaces which had been previously exposed to lubricants containing an organic chloride reacted much more slowly. In a separate series of experiments, a number of stainless steel plates were heated in a nitrogen atmosphere to different temperatures and their reactivity examined later at room temperature. The change of surface contour as a result of the probe reaction followed an Arrhenius type relation with respect to heat treatment temperature. This result could have implications on the scuffing mechanism.

Lauer, J. L.↗

Terrestrial geophysics in the SeaRISE project

Some areas of research in the SeaRISE project are briefly discussed. They are as follows: (1) Radar Sounding serves multiple purposes. The most general and obvious is mapping ice thickness and the surface and bedrock topography of the ice sheet. (2) The purpose of Seismic Shooting, in addition to water depth measurements on floating ice, is to provide information about the internal physical characteristics of the ice sheet, the rock beneath it, and the interface between the two. (3) Passive Seismic monitoring of microearthquakes can be used to study brittle fracture within the ice or the rock beneath it. Common parameters available from these studies are fault location, orientation, and displacement, as well as the size of the rupture area, stress drop, and energy released. (4) There is a large contrast in Electrical Resistivity between ice or permafrost on the one hand and liquid water or wet rock on the other hand. Thus, electrical resistivity profiles have the ability of revealing the depth to the melting point, whether it is found at the base of the ice or in the subglacial rock. (5) Gravity anomalies, especially combined with seismic measurements, are an effective tool for determining deeper crustal structure. Anomalies averaged over extensive areas are useful also for their potential to reveal isostatic imbalance, which is a measure of average glacial change over the last several hundred years.

Bentley, C. R.↗