SEARCH · Engineering Papers

Results for “land surface temperature”

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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Hummingbird Comet Nucleus Analysis Mission

Hummingbird is a highly focused scientific mission, proposed to NASA s Discovery Program, designed to address the highest priority questions in cometary science-that of the chemical composition of the cometary nucleus. After rendezvous with the comet, Hummingbird would first methodically image and map the comet, then collect and analyze dust, ice and gases from the cometary atmosphere to enrich characterization of the comet and support landing site selection. Then, like its namesake, Hummingbird would carefully descend to a pre-selected surface site obtaining a high-resolution image, gather a surface material sample, acquire surface temperature and then immediately return to orbit for detailed chemical and elemental analyses followed by a high resolution post-sampling image of the site. Hummingbird s analytical laboratory contains instrumentation for a comprehensive molecular and elemental analysis of the cometary nucleus as well as an innovative surface sample acquisition device.

Kojiro, Daniel↗

NASA NTRS2010

The Contribution of Soil Moisture Information to Forecast Skill: Two Studies

This talk briefly describes two recent studies on the impact of soil moisture information on hydrological and meteorological prediction. While the studies utilize soil moisture derived from the integration of large-scale land surface models with observations-based meteorological data, the results directly illustrate the potential usefulness of satellite-derived soil moisture information (e.g., from SMOS and SMAP) for applications in prediction. The first study, the GEWEX- and ClIVAR-sponsored GLACE-2 project, quantifies the contribution of realistic soil moisture initialization to skill in subseasonal forecasts of precipitation and air temperature (out to two months). The multi-model study shows that soil moisture information does indeed contribute skill to the forecasts, particularly for air temperature, and particularly when the initial local soil moisture anomaly is large. Furthermore, the skill contributions tend to be larger where the soil moisture initialization is more accurate, as measured by the density of the observational network contributing to the initialization. The second study focuses on streamflow prediction. The relative contributions of snow and soil moisture initialization to skill in streamflow prediction at seasonal lead, in the absence of knowledge of meteorological anomalies during the forecast period, were quantified with several land surface models using uniquely designed numerical experiments and naturalized streamflow data covering mUltiple decades over the western United States. In several basins, accurate soil moisture initialization is found to contribute significant levels of predictive skill. Depending on the date of forecast issue, the contributions can be significant out to leads of six months. Both studies suggest that improvements in soil moisture initialization would lead to increases in predictive skill. The relevance of SMOS and SMAP satellite-based soil moisture information to prediction are discussed in the context of these studies.

Koster, Randal↗

NASA NTRS2015

Variability of Phenology and Fluxes of Water and Carbon with Observed and Simulated Soil Moisture in the Ent Terrestrial Biosphere Model (Ent TBM Version 1.0.1.0.0)

The Ent Terrestrial Biosphere Model (Ent TBM) is a mixed-canopy dynamic global vegetation model developed specifically for coupling with land surface hydrology and general circulation models (GCMs). This study describes the leaf phenology submodel implemented in the Ent TBM version 1.0.1.0.0 coupled to the carbon allocation scheme of the Ecosystem Demography (ED) model. The phenology submodel adopts a combination of responses to temperature (growing degree days and frost hardening), soil moisture (linearity of stress with relative saturation) and radiation (light length). Growth of leaves, sapwood, fine roots, stem wood and coarse roots is updated on a daily basis. We evaluate the performance in reproducing observed leaf seasonal growth as well as water and carbon fluxes for four plant functional types at five Fluxnet sites, with both observed and prognostic hydrology, and observed and prognostic seasonal leaf area index. The phenology submodel is able to capture the timing and magnitude of leaf-out and senescence for temperate broadleaf deciduous forest (Harvard Forest and Morgan- Monroe State Forest, US), C3 annual grassland (Vaira Ranch, US) and California oak savanna (Tonzi Ranch, US). For evergreen needleleaf forest (Hyytiäla, Finland), the phenology submodel captures the effect of frost hardening of photosynthetic capacity on seasonal fluxes and leaf area. We address the importance of customizing parameter sets of vegetation soil moisture stress response to the particular land surface hydrology scheme. We identify model deficiencies that reveal important dynamics and parameter needs.

Soil moisture↗

NASA NTRS1975

Continued support in the study of lunar and planetary surfaces

Radar observations of various planetary surfaces are discussed. A radar investigation of Mars was conducted in conjunction with the Viking landing site selection process. Quasi-specular scattering from the lunar surface was interpreted in terms of horizonal scale dependence upon observing wavelengths. Furthermore, the effect of the extremely high temperatures encountered on the surface of Venus upon the dielectric constant of geophysical materials, and hence on the interpretation of radar results, was considered. The use of radio and radar techniques for the study of Saturns rings was also investigated.

Howard, H. T.↗

NASA NTRS2012

Recent Upgrades to NASA SPoRT Initialization Datasets for the Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed several products for its NOAA/National Weather Service (NWS) partners that can initialize specific fields for local model runs within the NOAA/NWS Science and Training Resource Center Environmental Modeling System (EMS). The suite of SPoRT products for use in the EMS consists of a Sea Surface Temperature (SST) composite that includes a Lake Surface Temperature (LST) analysis over the Great Lakes, a Great Lakes sea-ice extent within the SST composite, a real-time Green Vegetation Fraction (GVF) composite, and NASA Land Information System (LIS) gridded output. This paper and companion poster describe each dataset and provide recent upgrades made to the SST, Great Lakes LST, GVF composites, and the real-time LIS runs.

Case, Jonathan L.↗

NASA NTRS2017

Combustion Joining of Regolith Tiles for In-Situ Fabrication of Launch/Landing Pads on the Moon and Mars

To mitigate dust problems during launch-landing operations in lunar and Mars missions, it is desired to build solid pads on the surface. Recently, strong tiles have been fabricated from lunar regolith simulants using high-temperature sintering. The present work investigates combustion joining of these tiles through the use of exothermic intermetallic reactions. Specifically, nickel aluminum (1:1 mole ratio) mixture was placed in a gap between the tiles sintered from JSC-1A lunar regolith simulant. Upon ignition by a laser, a self-sustained propagation of the combustion front over the mixture occurred. Joining was improved with increasing the tile thickness from 6.3 mm to 12.7 mm. The temperatures sufficient for melting the glass phase of JSC-1A were recorded for 12.7-mm tiles, which explains the observed better joining.

Ferguson, Robert E.↗

NASA NTRS2012

Incorporating Plant Phenology Dynamics in a Biophysical Canopy Model

The Multi-Layer Canopy Model (MLCan) is a vegetation model created to capture plant responses to environmental change. Themodel vertically resolves carbon uptake, water vapor and energy exchange at each canopy level by coupling photosynthesis, stomatal conductance and leaf energy balance. The model is forced by incoming shortwave and longwave radiation, as well as near-surface meteorological conditions. The original formulation of MLCan utilized canopy structural traits derived from observations. This project aims to incorporate a plant phenology scheme within MLCan allowing these structural traits to vary dynamically. In the plant phenology scheme implemented here, plant growth is dependent on environmental conditions such as air temperature and soil moisture. The scheme includes functionality that models plant germination, growth, and senescence. These growth stages dictate the variation in six different vegetative carbon pools: storage, leaves, stem, coarse roots, fine roots, and reproductive. The magnitudes of these carbon pools determine land surface parameters such as leaf area index, canopy height, rooting depth and root water uptake capacity. Coupling this phenology scheme with MLCan allows for a more flexible representation of the structure and function of vegetation as it responds to changing environmental conditions.

environmental changes↗

NASA NTRS1989

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

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

Jedlovec, Gary J.↗

NASA NTRS1984

Origin of the Greenland Ice Sheet

A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.

North, G. R.↗

NASA NTRSDate not supplied

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗

NASA NTRS2017

The PMIP4 Contribution to CMIP6-Part 4: Scientific Objectives and Experimental Design of the PMIP4-CMIP6 Last Glacial Maximum Experiments and PMIP4 Sensitivity Experiments

The Last Glacial Maximum (LGM, 21,000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Kageyama, Masa↗

NASA NTRS2009

The GLACE-2 Experiment

A major motivation for the study of the coupled land-atmosphere system is the idea that soil moisture anomalies may affect future meteorological variables through their effects on future surface energy and water budgets. If true, the accurate initialization of soil moisture in a subseasonal or seasonal forecast system may improve forecast skill, making the forecast products more valuable to society. The Global Land-Atmosphere Coupling Experiment (GLACE-2) project is examining, with a wide variety of models, the degree to which subseasonal (out to two months) precipitation and air temperature forecasts improve through the realistic initialization of soil moisture. For the first time ever, a global consensus should emerge regarding the value of land initialization for forecasts, perhaps motivating national forecast centers to make full use of land moisture initialization in their operations. Participants in GLACE-2 perform two series of forecasts, each consisting of 100 2-month forecast ensembles (10 members per ensemble) covering ten boreal spring and summer start-dates in each of the years 1986-1995. Series 1 utilizes realistic land surface state initialization, provided through a decadal offline simulation using realistic meteorological forcing, as provided by the Global Soil Wetness Project - Phase 2 (GSWP-2), a research activity of the Global Land-Atmosphere System Study (GLASS) of GEWEX. Series 2 is identical to Series 1 in every way except for the fact that it does not benefit from realistic land state initialization. Through the comparison of Series 1 and 2, we isolate the impact of land initialization on the forecasts. Optional extensions to these base runs include forecasts covering additional years, using alternative meteorological forcing for the land initialization. To date, GLACE-2 has garnered participation from eleven modeling groups, covering 13 atmospheric models. Analysis of available results is already well underway, both at the individual modeling institutions and at the NASA Goddard Space Flight Center, which is coordinating the project. Analysis focuses on two elements of the forecast problem: (i) the quantification of model-specific "predictability" (i.e., the degree to which simulated atmospheric chaos will foil a forecast, even under the assumption of"perfect" model physics, initialization data, and validation data) and its decay with time; and (ii) the quantification of forecast skill, determined through a comparison of predicted precipitation and air temperature against observations. Indeed, the specific contribution of land initialization to both these elements is isolated through a comparison of the Series 1 and 2 forecasts. We examine the two elements at four different forecast leads: 1-15, 16-30, 31-45, and 46-60 days. Statistics-based approaches for enhancing skill (essentially using observational statistics to reduce the impact of the models' climatic biases) are also tested. In the present talk, we provide an update of progress in GLACE-2, featuring quantifications of predictability and forecast skill for a number of the participating models and providing a "first look" at the desired consensus view of land impacts on subseasonal forecasts.

Koster, Randal↗

NASA NTRS1974

Development of the Viking Mars lander thermal control subsystem design

Two Viking spacecraft, each consisting of a lander capsule coupled to an orbiter, will be launched toward Mars during the summer of 1975. About a year later, the orbiters will go into orbit around Mars and the landers will descend to the surface for 90-day landed missions. The lander must withstand a wide variety of environmental and operational conditions during all phases of the mission, including prelaunch sterilization. On the surface of Mars, the lander internal temperatures must be controlled under widely varying thermal environments and atmospheric conditions. The lander thermal design is based on the maximum use of passive techniques and is integrated into the overall vehicle design and operation. The solutions to the unusual combinations of design problems and a summary of the results of full scale model testing under simulated mission conditions are presented.

Morey, T. F.↗

NASA NTRS1976

Applications of Skylab data to land use and climatological analysis

The author has identified the following significant results. Skylab study in Central Atlantic Regional Ecological Test Site encompassed two separate but related tasks: (1) evaluation of photographic sensors S190A and B as sources of land use data for planning and managing land resources in major metropolitan regions, and (2) evaluation of the multispectral scanner S192 used in conjunction with associated data and analytical techniques as a data source on urban climates and the surface energy balance. Photographs from the Skylab S190B earth terrain camera were of greatest interest in the land use analysis task; they were of sufficiently high resolution to identify and map many level 2 and 3 land use categories. After being corrected to allow for atmosphere effects, output from thermal and visible bands of the S192 was employed in constructing computer map plots of albedo and surface temperature.

Alexander, R. H.↗

NASA NTRS2010

Design and Preliminary Thermal Performance of the Mars Science Laboratory Rover Heat Exchangers

The challenging range of proposed landing sites for the Mars Science Laboratory Rover requires a rover thermal management system that is capable of keeping temperatures controlled across a wide variety of environmental conditions. On the Martian surface where temperatures can be as cold as -123 degrees Centigrade and as warm as 38 degrees Centigrade, the Rover relies upon a Mechanically Pumped Fluid Loop (MPFL) and external radiators to maintain the temperature of sensitive electronics and science instruments within a -40 degrees Centigrade to 50 degrees Centigrade range. The MPFL also manages significant waste heat generated from the Rover power source, known as the Multi Mission Radioisotope Thermoelectric Generator (MMRTG). The MMRTG produces 110 Watts of electrical power while generating waste heat equivalent to approximately 2000 Watts. Two similar Heat Exchanger (HX) assemblies were designed to both acquire the heat from the MMRTG and radiate waste heat from the onboard electronics to the surrounding Martian environment. Heat acquisition is accomplished on the interior surface of each HX while heat rejection is accomplished on the exterior surface of each HX. Since these two surfaces need to be at very different temperatures in order for the MPFL to perform efficiently, they need to be thermally isolated from one another. The HXs were therefore designed for high in-plane thermal conductivity and extremely low through-thickness thermal conductivity by using aerogel as an insulator inside composite honeycomb sandwich panels. A complex assembly of hand welded and uniquely bent aluminum tubes are bonded onto the HX panels and were specifically designed to be easily mated and demated to the rest of the Rover Heat Recovery and Rejection System (RHRS) in order to ease the integration effort. During the cruise phase to Mars, the HX assemblies serve the additional function of transferring heat from the Rover MPFL to the separate Cruise Stage MPFL so that heat generated deep inside the Rover can be dissipated via the Cruise Stage radiators. Significant fabrication challenges had to be overcome in order to make the HX design a reality. The cruise phase thermal performance of the Rover HXs was verified in the cruise phase system level thermal vacuum test that was performed at JPL in January of 2009. The Rover HXs were modeled in I-DEAS TMG and predictions are compared to actual data from the test.

Heat Rejection System↗

NASA NTRSDate not supplied

Land-Atmosphere Coupling Simulation and Its Role in Subseasonal-to-Seasonal Prediction

Land-atmosphere (L-A) coupling can significantly influence subseasonal-to-seasonal (S2S) prediction. During periods of strong L-A coupling, land-atmosphere feedbacks are expected to enhance the memory of the system and therefore also the predictability and prediction skill. This study aims to evaluate S2S prediction of ambient surface air temperature under conditions of strong versus weak L-A coupling in forecasts.

Yuna Lim↗

NASA NTRSDate not supplied

Engineering the Interface: Advanced Surface Technologies for Lunar Dust Management and Equipment Longevity

Through the Artemis program, NASA intends to develop a sustainable human foothold on the Moon, ultimately paving the way for crewed exploration of Mars. The Moon's hostile environment poses numerous obstacles, including exposure to radiation, temperature extremes, micrometeoroid threats, and particularly the persistent problem of lunar dust. Lunar dust impacts nearly every aspect of surface operations through adhesion and abrasion mechanisms, with contamination from anthropogenic activities (landing, rovers) far outweighing natural phenomena. Multiple adhesion pathways contribute to surface contamination in the lunar environment, including van der Waals forces, electrostatic forces, chemical reaction, and magnetic forces from elemental iron deposits. Sharp asperities from micrometeoroid bombardment and atmospheric absence increase interaction potential and enable mechanical interlocking. Low cohesion between dust particles exacerbates these challenges, as minimal interaction potential between dust and nearby surfaces overcomes particle cohesion, causing contamination. Lunar dust adhesion mitigation technologies can be categorized as either active, requiring external energy, or passive, relying on intrinsic material properties. Ultrasonic and electrodynamic technologies have been developed to the highest technology readiness level for active approaches. Passive strategies primarily focus on surface chemistry and topography modifications. At NASA Langley Research Center, approaches include surface migration agents to reduce surface energy, topographical modification using laser ablation patterning, and tailored surface conductivity to reduce intrinsic adhesion force. Performance has been evaluated using custom-built ultrasonic and centrifuge instruments. Plume-surface interactions from lunar landers can propel micrometer-sized particles at velocities up to 1000 m s-1.8 These particles pose risks to landers, habitats and infrastructure, leading to erosion, degradation, and reduced component lifespan. A panel recovered from Surveyor III was determined to have been severely abraded because of lunar dust displaced from the Apollo 12 lunar module that landed 160 m away. The performance of metallic surfaces has been evaluated via high velocity single particle impact using the laser-induced project impact test (LIPIT) facility at the University of Utah. Peridynamics modeling, a form of continuum mechanics that uses a nonlocal approach enabling greater simulation capabilities of crack initiation and fracture, has also been utilized to gain greater insight into material response during impact events. Lunar dust contamination challenges extend to power generation systems and moving equipment. Cables, rotation stages, and other mechanisms may experience limited range of motion and reduced lifetime due to dust infiltration. NASA Langley Research Center has evaluated traditional aerospace alloys, softgoods, wear resistant ceramics, and several polymer and polymer composite materials. Test methods have included traditional techniques like Taber abrasion testing, as well as designed test configurations developed in the DUSTE (dust, ultraviolet radiation, and space thermal environmental) chamber that reproduce mechanism functions in operational environment. Beyond laboratory experiments, several flight experiments have been conducted. Materials were exposed to the low Earth orbit environment on the Materials International Space Station Experiment (MISSE) and to the lunar surface environment through the Aegis Aerospace Regolith Adherence Characterization (RAC) payload and the Honeybee Robotics PlanetVac payload. Determining lunar dust's impact on surface exploration and habitation requires comprehensive experimental and computational capabilities combined with lessons learned from initial lunar activities. Identifying the greatest environmental challenges and developing mitigation technologies provides the clearest path toward successfully, expeditiously, and efficaciously completing NASA's mission. This presentation will discuss ongoing efforts at NASA Langley Research Center and collaborator contributions to these critical objectives.

Surface Engineering↗

NASA NTRS1998

Progress in AMSR Snow Algorithm Development

Advanced Microwave Scanning Radiometer (AMSR) will be flown on-board of the Japanese Advanced Earth Observing Satellite-II (ADEOS-II) and United States Earth Observation System (EOS) PM-1 satellite. AMSR is a passive microwave radiometer with frequency ranges from 6.9 GHz to 89 GHz. It scans conically with a constant incidence angle of 55 deg at the Earth's surface. The swath width is about 1600 km. With a large antenna, AMSR will provide the best spatial resolution of multi-frequency radiometer from space. This provides us an opportunity to improve the snow parameter retrieval. Accurate determination of snow parameters from space is a challenging effort. Over the years, many different techniques have been used to account for the complicated snow parameters such as the density, stratigraphy, snow grain size, temperature variation of the snow-pack. Forest type, fractional forest cover and land use type also need to be considered in developing an improved retrieval algorithm. However, snow is such a dynamic variable, snow-pack parameter keeps changing once the snow is deposited on the earth surface. Currently, NASDA and NASA are developing AMSR snow retrieval algorithms. These algorithms are now being carefully tested and evaluated using the SSM/I data. Due to limited snow-pack data available for comparison, this activity is progressing slowly. However, it is clear that in order to improve the snow retrieval algorithm, it is necessary to model the metamorphism history of the snow-pack.

Chang, Alfred↗