Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heat Generation Rate”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Unsteady pressure loads in a generic high speed engine model

Unsteady pressure loads were measured along the top interior wall of a generic high-speed engine (GHSE) model undergoing performance tests in the combustion-Heated Scramjet Test Facility at the Langley Research Center. Flow to the model inlet was simulated at 72000 ft and a flight Mach number of 4. The inlet Mach number was 3.5 with a total temperature and pressure of 1640 R and 92 psia. The unsteady pressure loads were measured with 5 piezoresistive gages, recessed into the wall 4 to 12 gage diameters to reduce incident heat flux to the diaphragms, and distributed from the inlet to the combustor. Contributors to the unsteady pressure loads included boundary layer turbulence, combustion noise, and transients generated by unstart loads. Typical turbulent boundary layer rms pressures in the inlet ranged from 133 dB in the inlet to 181 dB in the combustor over the frequency range from 0 to 5 kHz. Downstream of the inlet exist, combustion noise was shown to dominate boundary layer turbulence noise at increased heat release rates. Noise levels in the isolator section increased by 15 dB when the fuel-air ratio was increased from 0.37 to 0.57 of the stoichiometric ratio. Transient pressure disturbances associated with engine unstarts were measured in the inlet and have an upstream propagation speed of about 7 ft/sec and pressure jumps of at least 3 psia.

Parrott, Tony L.↗
NASA NTRS2011

Online Oxide Contamination Measurement and Purification Demonstration

Liquid metal sodium-potassium (NaK) has advantageous thermodynamic properties indicating its use as a fission reactor coolant for a surface (lunar, martian) power system. A major area of concern for fission reactor cooling systems is system corrosion due to oxygen contaminants at the high operating temperatures experienced. A small-scale, approximately 4-L capacity, simulated fission reactor cooling system employing NaK as a coolant was fabricated and tested with the goal of demonstrating a noninvasive oxygen detection and purification system. In order to generate prototypical conditions in the simulated cooling system, several system components were designed, fabricated, and tested. These major components were a fully-sealed, magnetically-coupled mechanical NaK pump, a graphite element heated reservoir, a plugging indicator system, and a cold trap. All system components were successfully demonstrated at a maximum system flow rate of approximately 150 cc/s at temperatures up to 550 C. Coolant purification was accomplished using a cold trap before and after plugging operations which showed a relative reduction in oxygen content.

Bradley, D. E.↗
NASA NTRS1986

Effects of cloud/radiation interactions on the climatic temperature distribution

Results are presented from one- and two-dimensional modeling of climate perturbations induced by interactive cloud/radiation feedback phenomena. The one-dimensional simulation generated a vertical temperature profile caused by CO2 doubling in the atmosphere, with and without implementation of a cumulus parameterization scheme. The simulation showed that the surface temperature sensitivity in the tropics is higher in the presence of cumulus convection than in its absence. The two-dimensional model accounted for the latitudinal dependence of heating rates caused by cumulus convection and the interactions between the mean circulation and horizontal eddy transport of sensible and latent heat fluxes.

Ou, S.-C.↗
NASA NTRS2002

Water cycles in closed ecological systems: effects of atmospheric pressure

In bioregenerative life support systems that use plants to generate food and oxygen, the largest mass flux between the plants and their surrounding environment will be water. This water cycle is a consequence of the continuous change of state (evaporation-condensation) from liquid to gas through the process of transpiration and the need to transfer heat (cool) and dehumidify the plant growth chamber. Evapotranspiration rates for full plant canopies can range from ~1 to 10 L m-2 d-1 (~1 to 10 mm m-2 d-1), with the rates depending primarily on the vapor pressure deficit (VPD) between the leaves and the air inside the plant growth chamber. VPD in turn is dependent on the air temperature, leaf temperature, and current value of relative humidity (RH). Concepts for developing closed plant growth systems, such as greenhouses for Mars, have been discussed for many years and the feasibility of such systems will depend on the overall system costs and reliability. One approach for reducing system costs would be to reduce the operating pressure within the greenhouse to reduce structural mass and gas leakage. But managing plant growth environments at low pressures (e.g., controlling humidity and heat exchange) may be difficult, and the effects of low-pressure environments on plant growth and system water cycling need further study. We present experimental evidence to show that water saturation pressures in air under isothermal conditions are only slightly affected by total pressure, but the overall water flux from evaporating surfaces can increase as pressure decreases. Mathematical models describing these observations are presented, along with discussion of the importance for considering "water cycles" in closed bioregenerative life support systems.

NASA Center KSC↗
NASA NTRS2012

TRMM Latent Heating Retrieval and Comparisons with Field Campaigns and Large-Scale Analyses

Rainfall production is a fundamental process within the Earth's hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating (LH), is one of the principal sources of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The vertical distribution of LH has a strong influence on the atmosphere, controlling large-scale tropical circulations, exciting and modulating tropical waves, maintaining the intensities of tropical cyclones, and even providing the energetics of midlatitude cyclones and other mobile midlatitude weather systems. Moreover, the processes associated with LH result in significant non-linear changes in atmospheric radiation through the creation, dissipation and modulation of clouds and precipitation. Yanai et al. (1973) utilized the meteorological data collected from a sounding network to present a pioneering work on thermodynamic budgets, which are referred to as the apparent heat source (Q1) and apparent moisture sink (Q2). Yanai's paper motivated the development of satellite-based LH algorithms and provided a theoretical background for imposing large-scale advective forcing into cloud-resolving models (CRMs). These CRM-simulated LH and Q1 data have been used to generate the look-up tables used in LH algorithms. This paper examines the retrieval, validation, and application of LH estimates based on rain rate quantities acquired from the Tropical Rainfall Measuring Mission satellite (TRMM). TRMM was launched in November 1997 as a joint enterprise between the American and Japanese space agencies -- with overriding goals of providing accurate four-dimensional estimates of rainfall and LH over the global Tropics and subtropics equatorward of 35o. Other literature has acknowledged the achievement of the first goal of obtaining an accurate rainfall climatology. This paper describes the second major goal of obtaining credible LH estimates as well as their applications within TRMM's zone of coverage, the standard TRMM LH products, and areas for further improvement.

latent heating↗
NASA NTRS2003

Radiative Heat Loss Measurements During Microgravity Droplet Combustion in a Slow Convective Flow

Radiative heat loss from burning droplets in a slow convective flow under microgravity conditions is measured using a broad-band (0.6 to 40 microns) radiometer. In addition, backlit images of the droplet as well as color images of the flame were obtained using CCD cameras to estimate the burning rates and the flame dimensions, respectively. Tests were carried out in air at atmospheric pressure using n-heptane and methanol fuels with imposed forced flow velocities varied from 0 to 10 centimeters per second and initial droplet diameters varied from 1 to 3 millimeters. Slow convective flows were generated using three different experimental configurations in three different facilities in preparation for the proposed International Space Station droplet experiments. In the 2.2 Second Drop-Tower Facility a droplet supported on the leading edge of a quartz fiber is placed within a flow tunnel supplied by compressed air. In the Zero-Gravity Facility (five-second drop tower) a tethered droplet is translated in a quiescent ambient atmosphere to establish a uniform flow field around the droplet. In the KC 135 aircraft an electric fan was used to draw a uniform flow past a tethered droplet. Experimental results show that the burn rate increases and the overall flame size decreases with increases in forced-flow velocities over the range of flow velocities and droplet sizes tested. The total radiative heat loss rate, Q(sub r), decreases as the imposed flow velocity increases with the spherically symmetric combustion having the highest values. These observations are in contrast to the trends observed for gas-jet flames in microgravity, but consistent with the observations during flame spread over solid fuels where the burning rate is coupled to the forced flow as here.

Hicks, Michael C.↗
NASA NTRS1984

Energization of ions in the auroral plasma by broadband waves - Generation of ion conics

Nonresonant heating of plasma ions flowing along geomagnetic field lines by VLF waves is shown to generate ion conics. Ions with energies of less than 1 eV are proven capable of resonance with ionospheric VLF waves excited by electron beams with energies over 100 eV at lower hybrid frequencies (LHF) of 1-1.01. The variation of resonant ion energy with LHF wave frequency is formulated and an expression is derived for the heating rate. The heating rate is then applied to model ion acceleration along a dipole geomagnetic field in terms of the power spectral density of broadband waves, i.e., the VLF noise. Attention is given to ion energization at lower altitudes than the auroral acceleration region and examples of conic acceleration are cited.

Singh, N.↗
NASA NTRS1989

Experimental determination of stator endwall heat transfer

Local Stanton numbers were experimentally determined for the endwall surface of a turbine vane passage. A six vane linear cascade having vanes with an axial chord of 13.81 cm was used. Results were obtained for Reynolds numbers based on inlet velocity and axial chord between 73,000 and 495,000. The test section was connected to a low pressure exhaust system. Ambient air was drawn into the test section, inlet velocity was controlled up to a maximum of 59.4 m/sec. The effect of the inlet boundary layer thickness on the endwall heat transfer was determined for a range of test section flow rates. The liquid crystal measurement technique was used to measure heat transfer. Endwall heat transfer was determined by applying electrical power to a foil heater attached to the cascade endwall. The temperature at which the liquid crystal exhibited a specific color was known from a calibration test. Lines showing this specific color were isotherms, and because of uniform heat generation they were also lines of nearly constant heat transfer. Endwall static pressures were measured, along with surveys of total pressure and flow angles at the inlet and exit of the cascade.

Boyle, Robert J.↗
NASA NTRS1989

Experimental determination of stator endwall heat transfer

Local Stanton numbers were experimentally determined for the endwall surface of a turbine vane possage. A six vane linear cascade having vanes with an axial chord of 13.81 cm was used. Resutls were obtained for Reynolds numbers based on inlet velocity and axial chord between 75,000 and 495,000. The test section was connected to a low pressure exhaust system. Ambient air was drawn into the test section, inlet velocity was controlled up to a maximum of 59.4 m/sec. The effect of the inlet boundary layer thickness on the endwall heat transfer was determined for a range of test section flow rates. The liquid crystal measurement technique was used to measure heat transfer. Endwall heat transfer was determined by applying electrical power to a foil heater attached to the cascade endwall. The temperature at which the liquid crystal exhibited a specific color was known from a calibration test. Lines showing this specific color were isotherms, and because of uniform heat generation they were also lines of nearly constant heat transfer. Endwall static pressures were measured, along with surveys of total pressure and flow angles at the inlet and exit of the cascade.

Boyle, Robert J.↗
NASA NTRS1997

A Novel Silicon Micromachined Integrated MCM Thermal Management System

"Micromachining" is a chemical means of etching three-dimensional structures, typically in single- crystalline silicon. These techniques are leading toward what is coming to be referred to as MEMS (Micro Electro Mechanical Systems), where in addition to the ordinary two-dimensional (planar) microelectronics, it is possible to build three-dimensional n-ticromotors, electrically- actuated raicrovalves, hydraulic systems and much more on the same microchip. These techniques become possible because of differential etching rates of various crystallographic planes and materials used for semiconductor n-ticrofabfication. The University of Cincinnati group in collaboration with Karl Baker at NASA Lewis were the first to form micro heat pipes in silicon by the above techniques. Current work now in progress using MEMS technology is now directed towards the development of the next generation in MCM (Multi Chip Module) packaging. Here we propose to develop a complete electronic thermal management system which will allow densifica6on in chip stacking by perhaps two orders of magnitude. Furthermore the proposed technique will allow ordinary conu-nercial integrated chips to be utilized. Basically, the new technique involves etching square holes into a silicon substrate and then inserting and bonding commercially available integrated chips into these holes. For example, over a 100 1/4 in. by 1 /4 in. integrated chips can be placed on a 4 in. by 4 in. silicon substrate to form a Multi-Chip Module (MCM). Placing these MCM's in-line within an integrated rack then allows for three-diniensional stacking. Increased miniaturization of microelectronic circuits will lead to very high local heat fluxes. A high performance thermal management system will be specifically designed to remove the generated energy. More specifically, a compact heat exchanger with milli / microchannels will be developed and tested to remove the heat through the back side of this MCM assembly for moderate and high heat flux applications, respectively. The high heat load application of particular interest in mind is the motor controller developed by Martin Marietta for Nasa to control the thruster's directional actuators on space vechicles. Work is also proposed to develop highly advanced and improved porous wick structures for use in advanced heat loops. The porous wick will be micromachined from silicon using MEMS technology, thus permitting far superior control of pore size and pore distribution (over wicks made from sintered n-ietals), which in turn is expected to led to significantly improved heat loop performance.

Kazmierczak, M. J.↗
NASA NTRS2003

Spectral Retrieval of Latent Heating Profiles from TRMM PR Data: Comparison of Look-Up Tables

The primary goal of the Tropical Rainfall Measuring Mission (TRMM) is to use the information about distributions of precipitation to determine the four dimensional (i.e., temporal and spatial) patterns of latent heating over the whole tropical region. The Spectral Latent Heating (SLH) algorithm has been developed to estimate latent heating profiles for the TRMM Precipitation Radar (PR) with a cloud- resolving model (CRM). The method uses CRM- generated heating profile look-up tables for the three rain types; convective, shallow stratiform, and anvil rain (deep stratiform with a melting level). For convective and shallow stratiform regions, the look-up table refers to the precipitation top height (PTH). For anvil region, on the other hand, the look- up table refers to the precipitation rate at the melting level instead of PTH. For global applications, it is necessary to examine the universality of the look-up table. In this paper, we compare the look-up tables produced from the numerical simulations of cloud ensembles forced with the Tropical Ocean Global Atmosphere (TOGA) Coupled Atmosphere-Ocean Response Experiment (COARE) data and the GARP Atlantic Tropical Experiment (GATE) data. There are some notable differences between the TOGA-COARE table and the GATE table, especially for the convective heating. First, there is larger number of deepest convective profiles in the TOGA-COARE table than in the GATE table, mainly due to the differences in SST. Second, shallow convective heating is stronger in the TOGA COARE table than in the GATE table. This might be attributable to the difference in the strength of the low-level inversions. Third, altitudes of convective heating maxima are larger in the TOGA COARE table than in the GATE table. Levels of convective heating maxima are located just below the melting level, because warm-rain processes are prevalent in tropical oceanic convective systems. Differences in levels of convective heating maxima probably reflect differences in melting layer heights. We are now extending our study to simulations of other field experiments (e.g. SCSMEX and ARM) in order to examine the universality of the look-up table. The impact of look-up tables on the retrieved latent heating profiles will also be assessed.

Shige, Shoichi↗
NASA NTRSDate not supplied

Investigating Material Behavior in Atmospheric Entry Conditions: Arc-Jet Testing Insights from Meteorite Ablation to High-Temperature Coatings.

Arc-jet testing is an indispensable tool in elucidating the complex interactions materials undergo when subjected to the extreme thermal and mechanical stresses of atmospheric entry. By replicating these harsh conditions, arc-jets enable researchers to explore critical phenomena such as ablation, melting, and spallation, generating invaluable data that informs the development of sophisticated, physics-based models. These models are pivotal not only in predicting material performance for spacecraft re-entry systems but also in assessing the threat posed by celestial bodies as they encounter Earth’s atmosphere. One of the leading facilities for such investigations is the Hypersonic Materials Environmental Test System (HyMETS) at NASA’s Langley Research Center. A pioneering test campaign conducted at HyMETS focused on unraveling the ablation mechanisms of an ordinary chondrite meteorite (Tamdakht H5) and a terrestrial analog (basalt). These studies unveiled distinct material behaviors: Tamdakht demonstrated a remarkably stable melt flow, with mass loss primarily driven by the volatilization of elements such as iron and sodium, whereas basalt exhibited more aggressive surface degradation due to the rapid decomposition of hydrated minerals, culminating in significant spallation and surface material ejection. Beyond meteorite analysis, HyMETS has also been integral in the exploration of advanced Thermal Protection Systems (TPS) such as Phenolic Impregnated Carbon Ablator (PICA), a material with extensive flight heritage in missions like Stardust, Mars Science Laboratory, Mars 2020, and Osiris-REx. PICA, when coated with NuSil CV-1144-0—a polysiloxane resin designed to prevent particle shedding—undergoes a remarkable transformation upon heating. The resin pyrolyzes to form a thin, oxidation-resistant silicon oxycarbide layer that profoundly influences the material’s thermal response through a sophisticated four-stage process. Initially, the silicon oxycarbide acts as a formidable barrier, impeding reactive interactions with the PICA char and effectively suppressing surface temperature rise and material recession. However, as the heating persists, the protective layer decomposes via carbothermal reduction, exposing the underlying char, which in turn leads to a dramatic spike in surface temperature and accelerated material erosion. In the final stage, equilibrium is reached, with recession rates aligning with those of the virgin material. This intricate understanding of PICA-NuSil behavior under extreme conditions offers valuable insights, enhancing material response models and advancing the development of next-generation TPS for space exploration.

Chemistry↗
NASA NTRS1980

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗
NASA NTRS1987

A new approach for vibration control in large space structures

An approach for augmenting vibration damping characteristics in space structures with large panels is presented. It is based on generation of bending moments rather than forces. The moments are generated using bimetallic strips, suitably mounted at selected stations on both sides of the large panels, under the influence of differential solar heating, giving rise to thermal gradients and stresses. The collocated angular velocity sensors are utilized in conjunction with mini-servos to regulate the control moments by flipping the bimetallic strips. A simple computation of the rate of dissipation of vibrational energy is undertaken to assess the effectiveness of the proposed approach.

Kumar, K.↗
NASA NTRS1992

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗
NASA NTRS2006

Life Test Approach for Refractory Metal/Sodium Heat Pipes

Heat pipe life tests described in the literature have seldom been conducted on a systematic basis. Typically one or more heat pipes are built and tested for an extended period at a single temperature with simple condenser loading. This paper describes an approach to generate carefully controlled data that can conclusively establish heat pipe operating life with material-fluid combinations capable of extended operation. Approximately 10 years of operational life might be compressed into 3 years of laboratory testing through a combination of increased temperature and mass fluence. Two specific test series have been identified and include: investigation of long term corrosion rates based on the guidelines contained in ASTM G-68-80 (using 7 heat pipes); and investigation of corrosion trends in a cross correlation sequence at various temperatures and mass fluences based on a central composite test design (using 9 heat pipes). The heat pipes selected for demonstration purposes are fabricated from a Mo-44.5%Re alloy with a length of 0.3 meters and a diameter of 1.59 cm(to conserve material) with a condenser to evaporator length ratio of approximately 3. The wick is a crescent annular design formed from 400-mesh Mo-Re alloy material hot isostatically pressed to produce a final wick core of 20 microns or less.

Martin, James J.↗
NASA NTRS2010

Design of Refractory Metal Life Test Heat Pipe and Calorimeter

Heat pipe life tests have seldom been conducted on a systematic basis. Typically, one or more heat pipes are built and tested for an extended period at a single temperature with simple condenser loading. Results are often reported describing the wall material, working fluid, test temperature, test duration, and occasionally the nature of any failure. Important information such as design details, processing procedures, material assay, power throughput, and radial power density are usually not mentioned. We propose to develop methods to generate carefully controlled data that conclusively establish heat pipe operating life with material-fluid combinations capable of extended operation. The test approach detailed in this Technical Publication will use 16 Mo-44.5%Re alloy/sodium heat pipe units that have an approximate12-in length and 5/8-in diameter. Two specific test series have been identified: (1) Long-term corrosion rates based on ASTM-G-68-80 (G-series) and (2) corrosion trends in a cross-correlation sequence at various temperatures and mass fluences based on a Fisher multifactor design (F-series). Evaluation of the heat pipe hardware will be performed in test chambers purged with an inert purified gas (helium or helium/argon mixture) at low pressure (10-100 torr) to provide thermal coupling between the heat pipe condenser and calorimeter. The final pressure will be selected to minimize the potential for voltage breakdown between the heat pipe and radio frequency (RF) induction coil (RF heating is currently the planned method of powering the heat pipes). The proposed calorimeter is constructed from a copper alloy and relies on a laminar flow water-coolant channel design to absorb and transport energy

Martin, J. J.↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗