Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cubic”

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 181 records · Page 10

Inverse solution of Kepler's equation for hyperbolic orbits

An algorithm is presented for efficient inverse solution of Kepler's equation for hyperbolic orbits. It is shown that an expansion of Barker's equation into a bicubic polynomial provides a good approximation to obtain accurate starting values for rapid numerical solution of Kepler's equation. In the approximate equation a cubic in normalized elapsed flight time from pericenter is set equal to a cubic in a function S of eccentricity and true anomaly. The initial estimate of S to use in an iteration formula is obtained by evaluating the cubic in normalized flight time and finding in most cases the single real root of the other cubic. This initial estimate has an accuracy corresponding to values of true anomaly in error by less than 0.5 degrees generally.

Boltz, Frederick W.↗

The role of grain boundaries in hydrogen diffusion in metals at 25 C

The effect of grain size on hydrogen diffusion at 25 C was examined for 4340 steel (body-centered cubic) and for Inconel 718 (face-centered cubic). It was found that the effect of grain size is important for body-centered cubic structures, but plays a much less important role in face centered cubic structures. Accurate measurements of hydrogen desorption coefficients during hydrogen desorption show that these are not greatly different for both types of structures.

Danford, M. D.↗

Variability of Fram Strait Ice Flux and North Atlantic Oscillation

An important term in the mass balance of the Arctic Ocean sea ice is the ice export. We estimated the winter sea ice export through the Fram Strait using ice motion from satellite passive microwave data and ice thickness data from moored upward looking sonars. The average winter area flux over the 18-year record (1978-1996) is 670,000 square km, approximately 7% of the area of the Arctic Ocean. The winter area flux ranges from a minimum of 450,000 sq. km in 1984 to a maximum of 906,000 sq km in 1995. The daily, monthly and interannual variabilities of the ice area flux are high. There is an upward trend in the ice area flux over the 18-year record. The average winter volume flux over the winters of October 1990 through May 1995 is 1745 cubic km ranging from a low of 1375 cubic km in 1990 to a high of 2791 cubic km in 1994. The sea-level pressure gradient across the Fram Strait explains more than 80% of the variance in the ice flux over the 18-year record. We use the coefficients from the regression of the time-series of area flux versus pressure gradient across the Fram Strait and ice thickness data to estimate the summer area and volume flux. The average 12-month area flux and volume flux are 919,000 sq km and 2366 cubic km. We find a significant correlation (R =0.86) between the area flux and positive phases of the North Atlantic Oscillation (NAO) index over the months of December through March. Correlation between our six years of volume flux estimates and the NAO index gives R =0.56. During the high NAO years, a more intense Icelandic low increases the gradient in the sea-level pressure by almost 1 mbar across the Fram Strait thus increasing the atmospheric forcing on ice transport. Correlation is reduced during the negative NAO years because of decreased dominance of this large-scale atmospheric pattern on the sea-level pressure gradient across the Fram Strait. Additional information is contained in the original.

Kwok, Ron↗

A Side by Side Comparison of Filter-Based PM(sub 2.5) Measurements at a Suburban Site: A Closure Study

Reliable determination of the effects of air quality on public health and the environment requires accurate measurement of PM(sub 2.5) mass and the individual chemical components of fine aerosols. This study seeks to evaluate PM(sub 2.5) measurements that are part of a newly established national network by comparing them with a more conventional sampling system. Experiments were carried out during 2002 at a suburban site in Maryland, United States, where two samplers from the U.S. Environmental Protection Agency (USEPA) Speciation Trends Network: Met One Speciation Air Sampling System STNS and Thermo Scientific Reference Ambient Air Sampler STNR, two Desert Research Institute Sequential Filter Samplers DRIF, and a continuous TEOM monitor (Thermo Scientific Tapered Element Oscillating Microbalance) were sampling air in parallel. These monitors differ not only in sampling configuration but also in protocol-specific sample analysis procedures. Measurements of PM(sub 2.5) mass and major contributing species were well correlated among the different methods with r-values > 0.8. Despite the good correlations, daily concentrations of PM(sub 2.5) mass and major contributing species were significantly different at the 95% confidence level from 5 to 100% of the time. Larger values of PM(sub 2.5) mass and individual species were generally reported from STNR and STNS. The January STNR average PM(sub 2.5) mass (8.8 (micro)g/per cubic meter) was 1.5 (micro)g/per cubic meter larger than the DRIF average mass. The July STNS average PM(sub 2.5) mass (27.8 (micro)g/per cubic meter) was 3.8 (micro)g/per cubic meter larger than the DRIF average mass. These differences can only be partially accounted for by known random errors. Variations in flow control, face velocity, and sampling artifacts likely influence the measurement of PM(sub 2.5) speciation and mass closure. Simple statistical tests indicate that the current uncertainty estimates used in the STN network may underestimate the actual uncertainty.

Haines, Jennifer C.↗

A Vegetation Correction Methodology for Time Series Based Soil Moisture Retrieval From C-band Radar Observations

A methodology is presented to correct backscatter (sigma(sup 0)) observations for the effect of vegetation. The proposed methodology is based on the concept that the ratio of the surface scattering over the total amount of scattering (sigma(sup 0)(sub soil)/sigma(sup 0)) is only affected by the vegetation and can be described as a function of the vegetation water content. Backscatter observations sigma(sup 0) from the soil are not influenced by vegetation. Under bare soil conditions (sigma(sup 0)(sub soil)/sigma(sup 0)) equals 1. Under low to moderate biomass and soil moisture conditions, vegetation affects the observed sigma(sup 0) through absorption of the surface scattering and contribution of direct scattering by the vegetation itself. Therefore, the contribution of the surface scattering is smaller than the observed total amount of scattering and decreases as the biomass increases. For dense canopies scattering interactions between the soil surface and vegetation elements (e.g. leaves and stems) also become significant. Because these higher order scattering mechanisms are influenced by the soil surface, an increase in (sigma(sup 0)(sub soil)/sigma(sup 0)) may be observed as the biomass increases under densely vegetated conditions. This methodology is applied within the framework of time series based approach for the retrieval of soil moisture. The data set used for this investigation has been collected during a campaign conducted at USDA's Optimizing Production Inputs for Economic and Environmental Enhancement OPE-3) experimental site in Beltsville, Maryland (USA). This campaign took place during the corn growth cycle from May 10th to 0ctober 2nd, 2002. In this period the corn crops reached a vegetation water content of 5.1 kg m(exp -2) at peak biomass and a soil moisture range varying between 0.00 to 0.26 cubic cm/cubic cm. One of the deployed microwave instruments operated was a multi-frequency (C-band (4.75 GHz) and L-band (1.6 GHz)) quad-polarized (HH, HV, VV, VH) radar which was mounted on a 20 meter long boom. In the OPE-3 field campaign, radar observations were collected once a week at nominal times of 8 am, 10 am, 12 noon and 2 pm. During each data run the radar acquired sixty independent measurements within an azimuth of 120 degrees from a boom height of 12.2 m and at three different incidence angles (15,35, and 55 degrees). The sixty observations were averaged to provide one backscatter value for the study area and its accuracy is estimated to be 51.0 dB. For this investigation the C-band observations have been used. Application of the proposed methodology to the selected data set showed a well-defined relationship between (sigma(sup 0)(sub soil)/sigma(sup 0)) and the vegetation water content. It is found that this relationship can be described with two experimentally determined parameters, which depend on the sensing configuration (e.g. incidence angle and polarization). Through application of the proposed vegetation correction methodology and the obtained parameterizations, the soil moisture retrieval accuracy within the framework of a time series based approach is improved from 0.033 to 0.032 cubic cm/cubic cm, from 0.049 to 0.033 cubic cm/cubic cm and from 0.079 to 0.047 cubic cm/cubic cm for incidence angles of 15,35 and 55 degrees, respectively. Improvement in soil moisture retrieval due to vegetation correction is greater at larger incidence angles (due to the increased path length and larger vegetation effects on the surface signal at the larger angles).

Joseph, Alicia T.↗

Graded Index Silicon Geranium on Lattice Matched Silicon Geranium Semiconductor Alloy

A lattice matched silicon germanium (SiGe) semiconductive alloy is formed when a {111} crystal plane of a cubic diamond structure SiGe is grown on the {0001} C-plane of a single crystalline Al2O3 substrate such that a <110> orientation of the cubic diamond structure SiGe is aligned with a <1,0,-1,0> orientation of the {0001} C-plane. A lattice match between the substrate and the SiGe is achieved by using a SiGe composition that is 0.7223 atomic percent silicon and 0.2777 atomic percent germanium. A layer of Si(1-x), ,Ge(x) is formed on the cubic diamond structure SiGe. The value of X (i) defines an atomic percent of germanium satisfying 0.2277<X<1.0,(ii) is approximately 0.2777 where the layer of Si(1-x)Ge(x)interfaces with the cubic diamond structure SiGe, and (iii) increases linearly with the thickness of the layer of Si(1-x)Ge(x).

Park, Yeonjoon↗

Altitude Variation of the CO2 (V2)-O Quenching Rate Coefficient in Mesosphere and Lower Thermosphere

Among the processes governing the energy balance in the mesosphere and lower thermosphere (mlt), the quenching of CO2(N2) vibrational levels by collisions with oxygen atoms plays an important role. However, the k(CO2-O) values measured in the lab and retrieved from atmospheric measurements vary from 1.5 x 10(exp -12) cubic centimeters per second through 9.0 x 10(exp -12) cubic centimeters per second that requires further studying. In this work we used synergistic data from a ground based lidar and a satellite infrared radiometer to estimate K(CO2-O). We used the night- and daytime temperatures between 80 and 110 km measured by the colorado state university narrow-band sodium (Na) lidar located at fort collins, colorado (41N, 255E) as ground truth of the saber/timed nearly simultaneous (plus or minus 10 minutes) and common volume (within plus or minus 1 degree in latitude, plus or minus 2 degrees in longitude) observations. For each altitude in 80-110 km interval we estimate an "optimal" value of K(CO2-O) needed to minimize the discrepancy between the simulated 15 mm CO2 radiance and that measured by the saber/timed instrument. The K(CO2-O) obtained in this way varies in altitude from 3.5 x 10(exp -12) cubic centimeters per second at 80 km to 5.2 x 10(exp -12) cubic centimeters pers second for altitudes above 95 km. We discuss this variation of the rate constant and its impact on temperature retrievals from 15 mm radiance measurements and on the energy budget of mlt.

Feofilovi, Artem↗

Assessment of the SMAP Passive Soil Moisture Product

The National Aeronautics and Space Administration (NASA) Soil Moisture Active Passive (SMAP) satellite mission was launched on January 31, 2015. The observatory was developed to provide global mapping of high-resolution soil moisture and freeze-thaw state every two to three days using an L-band (active) radar and an L-band (passive) radiometer. After an irrecoverable hardware failure of the radar on July 7, 2015, the radiometer-only soil moisture product became the only operational Level 2 soil moisture product for SMAP. The product provides soil moisture estimates posted on a 36 kilometer Earth-fixed grid produced using brightness temperature observations from descending passes. Within months after the commissioning of the SMAP radiometer, the product was assessed to have attained preliminary (beta) science quality, and data were released to the public for evaluation in September 2015. The product is available from the NASA Distributed Active Archive Center at the National Snow and Ice Data Center. This paper provides a summary of the Level 2 Passive Soil Moisture Product (L2_SM_P) and its validation against in situ ground measurements collected from different data sources. Initial in situ comparisons conducted between March 31, 2015 and October 26, 2015, at a limited number of core validation sites (CVSs) and several hundred sparse network points, indicate that the V-pol Single Channel Algorithm (SCA-V) currently delivers the best performance among algorithms considered for L2_SM_P, based on several metrics. The accuracy of the soil moisture retrievals averaged over the CVSs was 0.038 cubic meter per cubic meter unbiased root-mean-square difference (ubRMSD), which approaches the SMAP mission requirement of 0.040 cubic meter per cubic meter.

passive↗

Characterization of Particulate Matter Profiling and Alveolar Deposition from Biomass Burning in Northern Thailand: The 7-SEAS Study

Biomass burning (BB) frequently occurs in SouthEast Asia (SEA), which significantly affects the air quality and could consequently lead to adverse health effects. The aim of this study was to characterize particulate matter (PM) and black carbon (BC) emitted from BB source regions in SEA and their potential of deposition in the alveolar region of human lungs. A 31-day characterization of PM profiling was conducted at the Doi Ang Khang (DAK) meteorology station in northern Thailand in March 2013. Substantial numbers of PM (10147 +/- 5800 # per cubic centimeter) with a geometric mean diameter (GMD) of 114.4 +/- 9.2 nm were found at the study site. The PM of less than 2.5 micron in aerodynamic diameter (PM sub 2.5) hourly-average mass concentration was 78.0 +/- 34.5 per cubic microgram whereas the black carbon (BC) mass concentration was 4.4 +/- 2.6 micrograms per cubic meter. Notably, high concentrations of nanoparticle surface area (100.5 +/- 54.6 square micrometers per cubic centimeter) emitted from biomass burning can be inhaled into the human alveolar region. Significant correlations with fire counts within different ranges around DAK were found for particle number, the surface area concentration of alveolar deposition, and BC. In conclusion, biomass burning is an important PM source in SEA, particularly nanoparticles, which has high potency to be inhaled into the lung environment and interact with alveolar cells, leading to adverse respiratory effects. The fire counts within 100 to 150 km shows the highest Pearson's r for particle number and surface area concentration. It suggests 12 to 24 hr could be a fair time scale for initial aging process of BB aerosols. Importantly, the people lives in this region could have higher risk for PM exposure.

biomass burning↗

TEMPERATURE CONSTRAINT OF THE MISTASTIN LAKE IMPACT STRUCTURE IMPACT MELT ROCKS

Impact melt is a product of hypervelocity cratering events formed by the instantaneous melting of near-surface target rocks. The temperatures of impact melt upon formation are believed to be superheated, vastly exceeding the liquidus temperature of igneous rocks formed via endogenic processes. However, quantitative understanding of the formation and duration of superheated temperatures remains incomplete. To date, only one datum on the hottest post-shock temperature (>2370 °C) has been derived by Timms et al. (2017) [1]. These authors studied the microstructures and crystallographic orientations preserved within a rim of vermicular ZrO2 of a zircon grain entrained within impact glass at the Mistastin Lake impact structure, Canada. Diagnostic crystallographic orientation indicated the presence of cubic zirconia, requiring a minimum temperature for the impact melt of >2370°C from the dissociation of zircon to cubic zirconia and liquid SiO2. With only one temperature datum, it cannot be stated whether the temperature of the impact melt sheet was homogenous across the entire structure or if the thermal history of the melt varied with spatial distribution. In this study, we present a preliminary investigation of 70 zircon grains, 22 of which exhibit rims of vermicular ZrO2 dissociation textures, similar to those studied by Timms et al. (2017). To acquire a more accurate temperature profile representative of impact melt, zircon crystals were collected from different types of impact melt bearing units, including additional samples of the impact glass studied by Timms et al. (2017). Of the 22 zircon grains with dissociation textures, at least 4 grains in the impact glass show diagnostic evidence of pre-existing cubic zirconia, and 1 grain from a glass-bearing breccia exhibits granular zircon textures, possible evidence of a transition from high-P to high-T conditions. These results show the heterogeneous post-shock temperature distribution of superheated impact melt within the Mistastin impact structure. [1] Timms, N.E., Erickson, T.M., Zanetti, M.R., Pearce, M.A., Cayron, C., Cavosie, A.J., Reddy, S.M., Wittmann, A. and Carpenter, P.K., 2017. Cubic zirconia in> 2370° C impact melt records Earth's hottest crust. Earth and Planetary Science Letters, 477, pp.52-58.

rock↗

The thermal expansion and high temperature transformation of GeSe

The thermal expansion of GeSe has been studied above room temperature up to the melting point of 670 plus or minus 5 C by X-ray diffraction techniques using a 190 mm Unicam high temperature camera. The thermal expansion of the crystallographic axes is linear with a distinct change in the expansion coefficients for all axes above 400 C. The relative changes in the axes indicate a rearrangement of the structure towards cubic symmetry with increasing temperature. The transformation of GeSe from the orthorhombic to a normal NaCl-type structure is observed at 651 plus or minus 5 C. The lattice parameter of the cubic form of GeSe is a 0 ? 5.730 plus or minus 0.003 A at 656 C. The GeSe lattice remains cubic up to the melting point.

Wiedemeier, H.↗

Infrared radiative transfer through a regular array of cuboidal clouds

Infrared radiative transfer through a regular array of cuboidal clouds is studied and the interaction of the sides of the clouds with each other and the ground is considered. The theory is developed for black clouds and is extended to scattering clouds using a variable azimuth two-stream approximation. It is shown that geometrical considerations often dominate over the microphysical aspects of radiative transfer through the clouds. For example, the difference in simulated 10 micron brightness temperature between black isothermal cubic clouds and cubic clouds of optical depth 10, is less than 2 deg for zenith angles less than 50 deg for all cloud fractions when viewed parallel to the array. The results show that serious errors are made in flux and cooling rate computations if broken clouds are modeled as planiform. Radiances computed by the usual practice of area-weighting cloudy and clear sky radiances are in error by 2 to 8 K in brightness temperature for cubic clouds over a wide range of cloud fractions and zenith angles. It is also shown that the lapse rate does not markedly affect the exiting radiances for cuboidal clouds of unit aspect ratio and optical depth 10.

HARSHVARDHAN↗

Performing spectroscopic and specific heat studies of improper ferroelectrics

The results of infrared measurements on Ni-Br, Cu-Cl, and Fe-I boracite improper ferroelectrics and far infrared measurements of Ni-Br boracite are presented. The boracites have the general formula X3B7O3Y, where X = divalent metal and Y = halogen. They undergo a first order phase transition from a high temperature paraelectric phase with cubic symmetry to a ferroelectric phase with orthorhombic symmetry. The boracites are "improper ferroelectrics" since the spontaneous polarization is not the primary order parameter in the cubic-orthorhombic phase transition. Current understanding of these materials is that the primary order parameter is associated with a doubly degenerate zone-boundary phonon in the cubic phase. The degenerate critical modes become homogeneous and split into the A sub 1 and A sub 2 modes in the orthorhombic phase, doubling the volume of the primitive cell. An harmonic coupling between the softing A sub 1 and a low frequency A sub 1 optic mode induces a spontaneous polarization as a secondary effect in the ferroelectric phase. This secondary non-critical nature of the ferroelectric mode earns these materials the "improper" title and is responsible for their unique properties and high figure of merit in detector use.

Coleman, L. B.↗

Non-noble catalysts and catalyst supports for phosphoric acid fuel cells

Tungsten carbide, which is active for hydrogen oxidation, is CO tolerant and has a hexagonal structure is discussed. Titanium carbide is inactive and has a cubic structure. Four different samples of the cubic alloys W sub x-1Ti sub XC sub 1-y were found to be active and CO tolerant. When the activities of these cubic alloys are weighted by the reciprocal of the square to those of highly forms of WC. They offer important insight into the nature of the active sites on W-C anode catalysts for use in phosphoric acid fuel cells.

Mcalister, A. J.↗

Rational-spline approximation with automatic tension adjustment

An algorithm for weighted least-squares approximation with rational splines is presented. A rational spline is a cubic function containing a distinct tension parameter for each interval defined by two consecutive knots. For zero tension, the rational spline is identical to a cubic spline; for very large tension, the rational spline is a linear function. The approximation algorithm incorporates an algorithm which automatically adjusts the tension on each interval to fulfill a user-specified criterion. Finally, an example is presented comparing results of the rational spline with those of the cubic spline.

Schiess, J. R.↗

Comparing the jerk with other global models of the geomagnetic field from 1960 to 1978

About 3300 satellite values of geomagnetic intensity and about 700 observatory values of annual mean magnetic vector components from 1960 to 1978 were fitted by three global models of the geomagnetic field B. Each model includes a spatially constant external field whose time dependence is a constant plus another constant times the Dst index, and each model accepts a time-independent station correction at each observatory. The time dependence of the internal Gauss coefficients is either cubic, quintic, or biquadratic (two independent quadratics, one before and one after January 1, 1970); and g1(0) also has an induced term proportional to the Dst index. The rms residual of the data fit is the same for the cubic and biquadratic models and insignificantly smaller for the quintic model. The quintic and biquadratic models have 1164 adjustable parameters, and the cubic has 1038. At a high level of significance the parameters of the best fitting biquadratic rule out a physical model for the magnetic impulse of 1969 in which the level surfaces of electrical conductivity in the lower mantle are approximately spherical, and the radial magnetic field at the core-mantle boundary goes from one quadratic time dependence to another in a year or less.

Backus, G. E.↗

The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.

Blackaby, Nicholas↗

Deep Basalt Aquifers in Orcus Patera, Elysium Basin Mars: Perspectives for Exobiology Exploration

Direct indicators of shorelines, spillways, and terraces allowed to determine the extent of the Elysium Paleolake between the contour-lines 1000 and 500 m below the Martian datum. The Elysium Paleolake is bordered north by Orcus Patera (14N/181W), which lies west of the Tartarus Montes and Tartarus Colles. The Orcus Patera displays an ellipse-shaped collapsed caldera of 360-km long and 100-km wide. Viking topographic data show that the bottom of the caldera is located at 2500 below the Martian datum, and surrounded by a steep-walled ram art which crest is located at about 0 m elevation. Considering the localization of Orcus Patera in the Elysium paleolake, its altimetry, and the magmatic origin of this caldera, we propose the existence of a paleolake in Orcus Patera generated (a) by juvenile water from magma during the Noachian period, and (b) by intermittent influx of the Elysium Basin from Hesperian to Amazonian. Results are encouraging to consider this site as a potential high-energy source environment for microbial communities. are circumscribed by a 50-km wide lava field mapped as Noachian material. The structure of Orcus Patera represents the record of material erupted from a magmatic reservoir. The caldera is enclosed by steep inner walls (25% measured from topographic data), values which could be in agreement with the presence of a deep magmatic reservoir, as suggested by the typology of Crumpler et.al. The depth of the caldera might be due to the collapse of the magma reservoir, and the release of gases accompanying the magma thermal evolution. Origins of water for the paleolake(s): The water that generated a paleolake in Orcus Patera may have come from two origins: (1) Juvenile water: Plescia and Crips estimated a magma H20 content by weight between 0.5% and 1.5% using for the first value a comparison with terrestrial basalt, and for the second values from a Martian meteorite. The amount of H20 can be estimated by the volume of erupted lava, and the lava content of the caldera. In this study, we adopt a water content of 1%. The total volume of magma that has been contained in the caldera, and the volume of lava contained in the observed lava field is about 110 x 10(exp 6) cubic km, that gives a total volume of 1.10 x 10(exp 6) cubic km of water. The juvenile water expelled by the overpressure within the magma chamber charged with desolved water-vapor may have moved into the crust. The decrease in overburden pressure led to bubble formation. The ascent of these bubbles generated a pressurization of the magma, which was sufficient to fracture the overlaying magma layer, (2) Water from Elysium paleolake. During the Amazonian, the rise of the Elysium paleolake level generated an overspilling that supplied the caldera with water. The southern portion of the crest shows a deep gap 12-km wide at -1500 m elevation, locating the gap between 500 to 1000 in below the assumed water of Elysium paleolake, thus facilitating the influx of Elysium paleolake water into Orcus Patera. Bathymetric calculations give a floor area of 25,500 sq km at -2000 m elevation, and a water volume of 42,000 cubic km, with a lake-level at -1500 m. A substantial amount of water may have percolated through the fractured lava, and part of the volume may have overspilled the northern crest of Orcus Patera to debouch in the Tartarus Montes region. We envision the formation of a subsurface aqueous environment in basaltic rocks at the contact of the two water-source origins, possibly the percolating surface lake water, and more likely the juvenile water. Similarly to terrestrial calderas, Orcus Patera might be surrounded by ring-fractures caused by the collapse of the magma chamber that followed the release of gases. These ring-fractures may have been covered later by sedimentation in the caldera (lacustrine, aeolian, and volcanic), and by mass wasting. The detumescence of the magma in the caldera, and the vesiculation of the juvenile water may have operated simultaneously. Comparatively to terrestrial melts, Martian iron-rich melts are denser. This greater density implies greater effusion rates (eight-times terrestrial values), and larger fissuration widths (two-times terrestrial ones). With increasing vesiculation of magma, the bubbles interact with one-another because there are of similar pressure. They make a magma froth at the contact with the caldera surface, and on the walls of the fractures. In the saturated magma, froth, where the volume ratio of gases-to-liquid is about 4:1, the bubbles form a huge surface area of interconnected spaces. Bubbles near the caldera surface disrupt the magma, and fragmentation takes place, which moves downward through the magma column. On Earth, the bubbles are likely to grow between 1 and 50 mm in diameter due to the difference between the magma surface tension, and the bubble supersaturation pressure. The Martian low-pressure at surface level is likely to accelerate the expansion of the bubbles, and increase their final diameter and number, creating more voids in the magma. The strong magma froth with enclosed juvenile water bubbles interconnected with exsolved gas bubbles constitute a potential geothermal environment for geochemical energy production from basalt and water that does not require excessive temperatures. This process can start at +20C. Similar types of environments have been shown on Earth as potential energy sources for microbial metabolism, and could have provided deep aqueous basaltic niches for possible Martian microorganisms, even geologically recently. During the Amazonian, combination of volcanism and water activity still existed on Mars. Moreover, this type of potential niches open ways for investigation of possible oases of extinct or extant life, not only on paleolakes, and surface hydrothermalism spring areas, but also all large systems of fossae, which combine hydrologic and volcanic activities, and which provide an energy source, and an underground shelter to prevent surface UV bombardment. Additional information contained in the original.

Grin, E. A.↗