SEARCH · Engineering Papers

Results for “SLOPE”

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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Impact of Footprint Diameter and Off-Nadir Pointing on the Precision of Canopy Height Estimates from Spaceborne Lidar

A spaceborne lidar mission could serve multiple scientific purposes including remote sensing of ecosystem structure, carbon storage, terrestrial topography and ice sheet monitoring. The measurement requirements of these different goals will require compromises in sensor design. Footprint diameters that would be larger than optimal for vegetation studies have been proposed. Some spaceborne lidar mission designs include the possibility that a lidar sensor would share a platform with another sensor, which might require off-nadir pointing at angles of up to 16 . To resolve multiple mission goals and sensor requirements, detailed knowledge of the sensitivity of sensor performance to these aspects of mission design is required. This research used a radiative transfer model to investigate the sensitivity of forest height estimates to footprint diameter, off-nadir pointing and their interaction over a range of forest canopy properties. An individual-based forest model was used to simulate stands of mixed conifer forest in the Tahoe National Forest (Northern California, USA) and stands of deciduous forests in the Bartlett Experimental Forest (New Hampshire, USA). Waveforms were simulated for stands generated by a forest succession model using footprint diameters of 20 m to 70 m. Off-nadir angles of 0 to 16 were considered for a 25 m diameter footprint diameter. Footprint diameters in the range of 25 m to 30 m were optimal for estimates of maximum forest height (R(sup 2) of 0.95 and RMSE of 3 m). As expected, the contribution of vegetation height to the vertical extent of the waveform decreased with larger footprints, while the contribution of terrain slope increased. Precision of estimates decreased with an increasing off-nadir pointing angle, but off-nadir pointing had less impact on height estimates in deciduous forests than in coniferous forests. When pointing off-nadir, the decrease in precision was dependent on local incidence angle (the angle between the off-nadir beam and a line normal to the terrain surface) which is dependent on the off-nadir pointing angle, terrain slope, and the difference between the laser pointing azimuth and terrain aspect; the effect was larger when the sensor was aligned with the terrain azimuth but when aspect and azimuth are opposed, there was virtually no effect on R2 or RMSE. A second effect of off-nadir pointing is that the laser beam will intersect individual crowns and the canopy as a whole from a different angle which had a distinct effect on the precision of lidar estimates of height, decreasing R2 and increasing RMSE, although the effect was most pronounced for coniferous crowns.

Pang, Yong↗

NASA NTRS2012

Calibration of Two-dimensional Floodplain Modeling in the Atchafalaya River Basin Using SAR Interferometry

Two-dimensional (2D) satellite imagery has been increasingly employed to improve prediction of floodplain inundation models. However, most focus has been on validation of inundation extent, with little attention on the 2D spatial variations of water elevation and slope. The availability of high resolution Interferometric Synthetic Aperture Radar (InSAR) imagery offers unprecedented opportunity for quantitative validation of surface water heights and slopes derived from 2D hydrodynamic models. In this study, the LISFLOOD-ACC hydrodynamic model is applied to the central Atchafalaya River Basin, Louisiana, during high flows typical of spring floods in the Mississippi Delta region, for the purpose of demonstrating the utility of InSAR in coupled 1D/2D model calibration. Two calibration schemes focusing on Manning s roughness are compared. First, the model is calibrated in terms of water elevations at a single in situ gage during a 62 day simulation period from 1 April 2008 to 1 June 2008. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. The best-fit models show that the mean absolute errors are 3.8 cm for a single in situ gage calibration and 5.7 cm/46 days for InSAR water level calibration. The optimum values of Manning's roughness coefficients are 0.024/0.10 for the channel/floodplain, respectively, using a single in situ gage, and 0.028/0.10 for channel/floodplain the using SAR. Based on the calibrated water elevation changes, daily storage changes within the size of approx 230 sq km of the model area are also calculated to be of the order of 107 cubic m/day during high water of the modeled period. This study demonstrates the feasibility of SAR interferometry to support 2D hydrodynamic model calibration and as a tool for improved understanding of complex floodplain hydrodynamics

Hahn Chul Jung↗

NASA NTRS2012

Kepler Observations of Rapid Optical Variability in Active Galactic Nuclei

Over three quarters in 2010 - 2011, Kepler monitored optical emission from four active galactic nuclei (AGN) with approx 30 min sampling, > 90% duty cycle and approx < 0.1% repeatability. These data determined the AGN optical fluctuation power spectral density functions (PSDs) over a wide range in temporal frequency. Fits to these PSDs yielded power law slopes of -2.6 to -3.3, much steeper than typically seen in the X-rays. We find evidence that individual AGN exhibit intrinsically different PSD slopes. The steep PSD fits are a challenge to recent AGN variability models but seem consistent with first order MRI theoretical calculations of accretion disk fluctuations.

Mushotzky, R. F.↗

NASA NTRS2012

On Spectral Invariance of Single Scattering Albedo for Weakly Absorbing Wavelengths

The single scattering albedo omega (sub 0 lambda) in atmospheric radiative transfer is the ratio of the scattering coefficient to the total extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength A and droplet size r. In this presentation we will show that for water droplets at weakly absorbing wavelengths, the ratio omega (sub 0 lambda)(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo omega (sub 0 lambda) via one known spectrum omega (sub 0 lambda)(r(sub o)). We will provide a simple physical explanation of the discovered relationship. In addition to water droplets, similar linear relationships were found for the single scattering albedo of non-spherical ice crystals. The single scattering albedo $\omega _ {0\lambda }$ in atmospheric radiative transfer is the ratio of the scattering coefficient to the total extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, and thus the single scattering albedo, are functions of wavelength $\lambda $ and droplet size $r$. We show that for water droplets at weakly absorbing wavelengths, the ratio $\omega _ {0\lambda } (r)$/$\omega _ {0\lambda } (r_{0})$ of two single scattering albedo spectra for two different droplet sizes is a linear function of $\omega _{0\lambda }(r)$. The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo $\omega_{0\lambda }(r)$ via one known spectrum $\omega_{0\lambda }(r_{0})$. We provide a simple physical explanation of the discovered relationship. Similar linear relationships characterize the single scattering albedo of non-spherical ice crystals.

Marshak, Alexander↗

NASA NTRS2012

Progress of Multi-Beam Long Trace-Profiler Development

The multi-beam long trace profiler (LTP) under development at NASA s Marshall Space Flight Center[1] is designed to increase the efficiency of metrology of replicated X-ray optics. The traditional LTP operates on a single laser beam that scans along the test surface to detect the slope errors. While capable of exceptional surface slope accuracy, the LTP single beam scanning has slow measuring speed. As metrology constitutes a significant fraction of the time spent in optics production, an increase in the efficiency of metrology helps in decreasing the cost of fabrication of the x-ray optics and in improving their quality. Metrology efficiency can be increased by replacing the single laser beam with multiple beams that can scan a section of the test surface at a single instance. The increase in speed with such a system would be almost proportional to the number of laser beams. A collaborative feasibility study has been made and specifications were fixed for a multi-beam long trace profiler. The progress made in the development of this metrology system is presented.

Gubarev, Mikhail↗

NASA NTRS2012

Energy Spectra, Composition, and Other Properties of Ground-Level Events During Solar Cycle 23

We report spacecraft measurements of the energy spectra of solar protons and other solar energetic particle properties during the 16 Ground Level Events (GLEs) of Solar Cycle 23. The measurements were made by eight instruments on the ACE, GOES, SAMPBX, and STEREO spacecraft and extend from approximately 0.1 to approximately 500-700 MeV. All of the proton spectra exhibit spectral breaks at energies ranging from approximately 2 to approximately 46 MeV and all are well fit by a double power-law shape. A comparison of GLE events with a larger sample of other solar energetic particle (SEP) events shows that the typical spectral indices are harder in GLE events, with a mean slope of -3.18 at greater than 40 MeV/nuc. In the energy range 45 to 80 MeV/nucleon about approximately 50% of GLE events have properties in common with impulsive He-3-rich SEP events, including enrichments in Ne/O, Fe/O, Ne-22/Ne-20, and elevated mean charge states of Fe. These He-3 rich events contribute to the seed population accelerated by CME-driven shocks. An analysis is presented of whether highly-ionized Fe ions observed in five events could be due to electron stripping during shock acceleration in the low corona. Making use of stripping calculations by others and a coronal density model, we can account for events with mean Fe charge states of (Q(sub Fe) is approximately equal to +20 if the acceleration starts at approximately 1.24-1.6 solar radii, consistent with recent comparisons of CME trajectories and type-II radio bursts. In addition, we suggest that gradual stripping of remnant ions from earlier large SEP events may also contribute a highly-ionized suprathermal seed population. We also discuss how observed SEP spectral slopes relate to the energetics of particle acceleration in GLE and other large SEP events.

Mewaldt, R. A.↗

NASA NTRS2012

Characterization of the Morphometry of Impact Craters Hosting Polar Deposits in Mercury's North Polar Region

Earth-based radar images of Mercury show radar-bright material inside impact craters near the planet s poles. A previous study indicated that the polar-deposit-hosting craters (PDCs) at Mercury s north pole are shallower than craters that lack such deposits. We use data acquired by the Mercury Laser Altimeter on the MESSENGER spacecraft during 11 months of orbital observations to revisit the depths of craters at high northern latitudes on Mercury. We measured the depth and diameter of 537 craters located poleward of 45 N, evaluated the slopes of the northern and southern walls of 30 PDCs, and assessed the floor roughness of 94 craters, including nine PDCs. We find that the PDCs appear to have a fresher crater morphology than the non-PDCs and that the radar-bright material has no detectable influence on crater depths, wall slopes, or floor roughness. The statistical similarity of crater depth-diameter relations for the PDC and non-PDC populations places an upper limit on the thickness of the radar-bright material (< 170 m for a crater 11 km in diameter) that can be refined by future detailed analysis. Results of the current study are consistent with the view that the radar-bright material constitutes a relatively thin layer emplaced preferentially in comparatively young craters.

Talpe Matthieu↗

NASA NTRS2012

Hesperian-Amazonian Transition Mid-Latitude Valleys: Markers of a Late Martian Climate Optima?

Recently the inventory of fluvial features that have been dated to the late Hesperian to early Amazonian epoch has increased dramatically, including a reassessment of the ages of the large alluvial fans and deltas (e.g., Eberswalde) to this time period. Mid-latitude Valleys (MLVs) are distinct from the older, more integrated Noachian-Hesperian Valley Networks which are deeply dissected, are generally of much larger spatial extent, and are more degraded. Although some MLVs involve rejuvenation of older Valley Networks, many MLVs are carved into smooth or rolling slopes and intercrater terrain. The MLVs range from a few meters to < 300 m in width, with nearly parallel valley walls and planforms that are locally sinuous. Although the MLVs in Newton and Gorgonum basins extend from the basin rims up to 75 km into the basin interior, most MLVs are shorter and often discontinuous. The occurrence of widespread MLVs suggest the possibility of their formation during one or perhaps more regional to global climatic episodes, possibly due to melting of seasonal to long-term accumulations of snow and ice. Temperatures warm enough to cause extensive melting may have occurred during optimal orbital and obliquity configurations, perhaps in conjunction with intensive volcanism releasing moisture and greenhouse gasses, or as a result of a brief episode of warming from a large impact. The concentration of MLVs to the northern and western basin slopes of Newton and Gorgonum basins suggests a possible aspect control to ice accumulation or melting. MLV activity occurred about at the same time as formation of the major outflow channels. A possible scenario is that delivery of water to the northern lowlands provided, through evaporation and sublimation, water that temporarily accumulated in the mid-southern latitudes as widespread ice deposits whose partial melting formed the MLVs and small, dominantly ice-covered lakes.

Moore, Jeffrey↗

NASA NTRS2013

Constraints on the Compositions of Phobos and Deimos from Mineral Absorptions

The compositions of Phobos and Deimos have remained controversial despite multiple Earth- and space-based observations acquired during the last 40 years. Phobos is composed of at least two spectral units that are both dark yet distinct at visible to near infrared wavelenghts; a spectrally red-sloped "red" unit covers most of the moon and a less red-sloped "blue" unit is present in the ejecta of the approximately 9-km diameter impact crater Stickney [1,2]. Deimos is similar spectrally to Phobos' "red" unit [2]. Here we report results from mapping mineral absorptions on Phobos and Deimos using visible/near infrared observations from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM). We find evidence for an absorption feature at 0.65 m in the Phobos red unit and Deimos that is reproducible in observations from other instruments. The phase responsible is uncertain but may be a Fe-bearing phyllosilicate and/or graphite, consistent with the notion that Phobos and Deimos have compositions similar to CM carbonaceous chondrites [3].

Fraeman, A. A.↗

NASA NTRS2013

Rolling Bearing Life Prediction, Theory, and Application

A tutorial is presented outlining the evolution, theory, and application of rolling-element bearing life prediction from that of A. Palmgren, 1924; W. Weibull, 1939; G. Lundberg and A. Palmgren, 1947 and 1952; E. Ioannides and T. Harris, 1985; and E. Zaretsky, 1987. Comparisons are made between these life models. The Ioannides-Harris model without a fatigue limit is identical to the Lundberg-Palmgren model. The Weibull model is similar to that of Zaretsky if the exponents are chosen to be identical. Both the load-life and Hertz stress-life relations of Weibull, Lundberg and Palmgren, and Ioannides and Harris reflect a strong dependence on the Weibull slope. The Zaretsky model decouples the dependence of the critical shear stress-life relation from the Weibull slope. This results in a nominal variation of the Hertz stress-life exponent. For 9th- and 8th-power Hertz stress-life exponents for ball and roller bearings, respectively, the Lundberg- Palmgren model best predicts life. However, for 12th- and 10th-power relations reflected by modern bearing steels, the Zaretsky model based on the Weibull equation is superior. Under the range of stresses examined, the use of a fatigue limit would suggest that (for most operating conditions under which a rolling-element bearing will operate) the bearing will not fail from classical rolling-element fatigue. Realistically, this is not the case. The use of a fatigue limit will significantly overpredict life over a range of normal operating Hertz stresses. Since the predicted lives of rolling-element bearings are high, the problem can become one of undersizing a bearing for a particular application.

Zaretsky, Erwin V.↗

NASA NTRS2012

Retrieval of Vertical LAI Profiles Over Tropical Rain Forests using Waveform Lidar at La Selva, Costa Rica

This study explores the potential of waveform lidar in mapping the vertical and spatial distributions of leaf area index (LAI) over the tropical rain forest of La Selva Biological Station in Costa Rica. Vertical profiles of LAI were derived at 0.3 m height intervals from the Laser Vegetation Imaging Sensor (LVIS) data using the Geometric Optical and Radiative Transfer (GORT) model. Cumulative LAI profiles obtained from LVIS were validated with data from 55 ground to canopy vertical transects using a modular field tower to destructively sample all vegetation. Our results showed moderate agreement between lidar and field derived LAI (r2=0.42, RMSE=1.91, bias=−0.32), which further improved when differences between lidar and tower footprint scales (r2=0.50, RMSE=1.79, bias=0.27) and distance of field tower from lidar footprint center (r2=0.63, RMSE=1.36, bias=0.0) were accounted for. Next, we mapped the spatial distribution of total LAI across the landscape and analyzed LAI variations over different land cover types. Mean values of total LAI were 1.74, 5.20, 5.41 and 5.62 over open pasture, secondary forests, regeneration forests after selective-logging and old-growth forests respectively. Lastly, we evaluated the sensitivities of our LAI retrieval model to variations in canopy/ground reflectance ratio and to waveform noise such as induced by topographic slopes. We found for both, that the effects were not significant for moderate LAI values (about 4). However model derivations of LAI might be inaccurate in areas of high-slope and high LAI (about 8) if ground return energies are low. This research suggests that large footprint waveform lidar can provide accurate vertical LAI profile estimates that do not saturate even at the high LAI levels in tropical rain forests and may be a useful tool for understanding the light transmittance within these canopies.

Tang, Hao↗

NASA NTRS2013

The W-W02 Oxygen Fugacity Buffer at High Pressures and Temperatures: Implications for f02 Buffering and Metal-silicate Partitioning

Oxygen fugacity (fO2) controls multivalent phase equilibria and partitioning of redox-sensitive elements, and it is important to understand this thermodynamic parameter in experimental and natural systems. The coexistence of a metal and its oxide at equilibrium constitutes an oxygen buffer which can be used to control or calculate fO2 in high pressure experiments. Application of 1-bar buffers to high pressure conditions can lead to inaccuracies in fO2 calculations because of unconstrained pressure dependencies. Extending fO2 buffers to pressures and temperatures corresponding to the Earth's deep interior requires precise determinations of the difference in volume (Delta) V) between the buffer phases. Synchrotron x-ray diffraction data were obtained using diamond anvil cells (DAC) and a multi anvil press (MAP) to measure unit cell volumes of W and WO2 at pressures and temperatures up to 70 GPa and 2300 K. These data were fitted to Birch-Murnaghan 3rd-order thermal equations of state using a thermal pressure approach; parameters for W are KT = 306 GPa, KT' = 4.06, and αKT = 0.00417 GPa K-1. Two structural phase transitions were observed for WO2 at 4 and 32 GPa with structures in P21/c, Pnma and C2/c space groups. Equations of state were fitted for these phases over their respective pressure ranges yielding the parameters KT = 190, 213, 300 GPa, KT' = 4.24, 5.17, 4 (fixed), and αKT = 0.00506, 0.00419, 0.00467 GPa K-1 for the P21/c, Pnma and C2/c phases, respectively. The W-WO2 buffer (WWO) was extended to high pressure by inverting the W and WO2 equations of state to obtain phase volumes at discrete pressures (1-bar to 100 GPa, 1 GPa increments) along isotherms (300 to 3000K, 100 K increments). The slope of the absolute fO2 of the WWO buffer is positive with increasing temperature up to approximately 70 GPa and is negative above this pressure. The slope is positive along isotherms from 1000 to 3000K with increasing pressure up to at least 100 GPa. The WWO buffer is at a higher fO2 than the IW buffer at pressures lower than 40 GPa, and the magnitude of this difference decreases at higher pressures. This qualitatively indicates an increasingly lithophile character for W at higher pressures. The WWO buffer was quantitatively applied to W metal-silicate partitioning by using the WWO-IW buffer difference in combination with literature data on W metal-silicate partitioning to model the exchange coefficient (KD) for the Fe-W exchange reaction. This approach captures the pressure dependence of W metal-silicate partitioning using the WWO-IW buffer difference and models the activities of the components in the silicate and metallic phases using an expression of the Gibbs excess energy of mixing. Calculation of KD along a peridotite liquidus predicts a decrease in W siderophility at higher pressures that supports the qualitative behavior predicted by the WWO-IW buffer difference, and agrees with findings of others. Comparing the competing effects of temperature and pressure on W metal-silicate partitioning, our results indicate that pressure exerts a greater effect.

Shofner, G. A.↗

NASA NTRS2013

Status of Multi-beam Long Trace-profiler Development

The multi-beam long trace profiler (MB-LTP) is under development at NASA's Marshall Space Flight Center. The traditional LTPs scans the surface under the test by a single laser beam directly measuring the surface figure slope errors. While capable of exceptional surface slope accuracy, the LTP single beam scanning has slow measuring speed. Metrology efficiency can be increased by replacing the single laser beam with multiple beams that can scan a section of the test surface at a single instance. The increase in speed with such a system would be almost proportional to the number of laser beams. The progress for a multi-beam long trace profiler development will be presented.

Gubarev, Mikhail V.↗

NASA NTRS2013

Status of Multi-beam Long Trace-profiler Development

Gubarev, Mikhail V.↗

NASA NTRS2012

Size Distributions of Solar Flares and Solar Energetic Particle Events

We suggest that the flatter size distribution of solar energetic proton (SEP) events relative to that of flare soft X-ray (SXR) events is primarily due to the fact that SEP flares are an energetic subset of all flares. Flares associated with gradual SEP events are characteristically accompanied by fast (much > 1000 km/s) coronal mass ejections (CMEs) that drive coronal/interplanetary shock waves. For the 1996-2005 interval, the slopes (alpha values) of power-law size distributions of the peak 1-8 Angs fluxes of SXR flares associated with (a) >10 MeV SEP events (with peak fluxes much > 1 pr/sq cm/s/sr) and (b) fast CMEs were approx 1.3-1.4 compared to approx 1.2 for the peak proton fluxes of >10 MeV SEP events and approx 2 for the peak 1-8 Angs fluxes of all SXR flares. The difference of approx 0.15 between the slopes of the distributions of SEP events and SEP SXR flares is consistent with the observed variation of SEP event peak flux with SXR peak flux.

Solar Flares↗

NASA NTRS2014

Smooth Pursuit of Flicker-Defined Motion

We examined the pursuit response to stimuli defined by space-variant flicker of a dense random dot carrier pattern. On each frame, every element of the pattern could change polarity, with a probability given by a two-dimensional Gaussian distribution. A normal distribution produces a circular region of twinkle, while inverting the distribution results in a spot of static texture in a twinkling surround. In this latter case, the carrier texture could be stationary, or could move with the twinkle modulator, thereby producing first-order motion in the region of the spot. While the twinkle-defined spot produces a strong sensation of motion, the complementary stimulus defined by the absence of twinkle does not, when viewed peripherally, it appears to move in steps even when the generating distribution moves smoothly. We examined pursuit responses to these stimuli using two techniques: 1) the eye movement correlogram, obtained by cross-correlating eye velocity with the velocity of a randomly-moving stimulus; and 2) delayed visual feedback, where transient stabilization of a target can produce spontaneous oscillations of the eye, with a period empirically observed to vary linearly with the applied delay. Both techniques provide an estimate of the internal processing time, which can be as short as 100 milliseconds for a first-order target. Assessed by the correlogram method, the response to flicker-defined motion is delayed by more than 100 milliseconds, and significantly weaker (especially in the vertical dimension). When initially presented in the delayed feedback condition, purely saccadic oscillation is observed. One subject eventually developed smooth oscillations (albeit with significant saccadic intrusions), showing a period-versus-delay slope similar to that observed for first-order targets. This result is somewhat surprising, given that we interpret the slope of the period-versus-delay-function as reflecting the balance between position- and velocity-sensitive inputs to pursuit.

eye movements↗

NASA NTRS2013

Efficient Flowline Simulations of Ice Shelf-Ocean Interactions: Sensitivity Studies with a Fully Coupled Model

Thermodynamic flowline and plume models for the ice shelf-ocean system simplify the ice and ocean dynamics sufficiently to allow extensive exploration of parameters affecting ice-sheet stability while including key physical processes. Comparison between geophysically and laboratory-based treatments of ice-ocean interface thermodynamics shows reasonable agreement between calculated melt rates, except where steep basal slopes and relatively high ocean temperatures are present. Results are especially sensitive to the poorly known drag coefficient, highlighting the need for additional field experiments to constrain its value. These experiments also suggest that if the ice-ocean interface near the grounding line is steeper than some threshold, further steepening of the slope may drive higher entrainment that limits buoyancy, slowing the plume and reducing melting; if confirmed, this will provide a stabilizing feedback on ice sheets under some circumstances.

polar oceanography↗

NASA NTRS2013

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

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

Hydrogen Concentrations↗