Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed estimation”

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 541 records · Page 30

An Update on Oceanic Precipitation Rate and its Zonal Distribution in Light of Advanced Observations from Space

This study contributes to the estimation of the global mean and zonal distribution of oceanic precipitation rate using complementary information from advanced precipitation measuring sensors and provides an independent reference to assess current precipitation products. Precipitation estimates from the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) and CloudSat cloud profiling radar (CPR) were merged, as the two complementary sensors yield an unprecedented range of sensitivity to quantify rainfall from drizzle through the most intense rates. At higher latitudes, where TRMM PR does not exist, precipitation estimates from Aqua's Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) complemented CloudSat CPR to capture intense precipitation rates. The high sensitivity of CPR allows estimation of snow rate, an important type of precipitation at high latitudes, not directly observed in current merged precipitation products. Using the merged precipitation estimate from the CloudSat, TRMM, and Aqua platforms (this estimate is abbreviated to MCTA), the authors' estimate for 3-yr (2007-09) nearglobal (80degS-80degN) oceanic mean precipitation rate is approx. 2.94mm/day. This new estimate of mean global ocean precipitation is about 9% higher than that of the corresponding Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) value (2.68mm/day) and about 4% higher than that of the Global Precipitation Climatology Project (GPCP; 2.82mm/day). Furthermore, MCTA suggests distinct differences in the zonal distribution of precipitation rate from that depicted in GPCPand CMAP, especially in the Southern Hemisphere.

profiling radar↗

Intercomparison of Precipitation Estimates over the Arctic Ocean and Its Peripheral Seas from Reanalyses

Precipitation over the Arctic Ocean has a significant impact on the basin-scale freshwater and energy budgets but is one of the most poorly constrained variables in atmospheric reanalyses. Precipitation controls the snow cover on sea ice, which impedes the exchange of energy between the ocean and atmosphere, inhibiting sea ice growth. Thus, accurate precipitation amounts are needed to inform sea ice modeling, especially for the production of thickness estimates from satellite altimetry freeboard data. However, obtaining a quantitative estimate of the precipitation distribution in the Arctic is notoriously difficult because of a number of factors, including a lack of reliable, long-term in situ observations; difficulties in remote sensing over sea ice; and model biases in temperature and moisture fields and associated uncertainty of modeled cloud microphysical processes in the polar regions. Here, we compare precipitation estimates over the Arctic Ocean from eight widely used atmospheric reanalyses over the period 2000–16 (nominally the ‘‘new Arctic’’). We find that the magnitude, frequency, and phase of precipitation vary drastically, although interannual variability is similar. Reanalysis-derived precipitation does not increase with time as expected; however, an increasing trend of higher fractions of liquid precipitation (rainfall) is found. When compared with drifting ice mass balance buoys, three reanalyses (ERA-Interim, MERRA, and NCEP R2) produce realistic magnitudes and temporal agreement with observed precipitation events, while two products [MERRA, version 2 (MERRA-2), and CFSR] show large, implausible magnitudes in precipitation events. All the reanalyses tend to produce overly frequent Arctic precipitation. Future work needs to be undertaken to determine the specific factors in reanalyses that contribute to these discrepancies in the new Arctic.

Mcpartland, Linette N.↗

Estimation of Snow Parameters from Dual-Wavelength Airborne Radar

Estimation of snow characteristics from airborne radar measurements would complement In-situ measurements. While In-situ data provide more detailed information than radar, they are limited in their space-time sampling. In the absence of significant cloud water contents, dual-wavelength radar data can be used to estimate 2 parameters of a drop size distribution if the snow density is assumed. To estimate, rather than assume, a snow density is difficult, however, and represents a major limitation in the radar retrieval. There are a number of ways that this problem can be investigated: direct comparisons with in-situ measurements, examination of the large scale characteristics of the retrievals and their comparison to cloud model outputs, use of LDR measurements, and comparisons to the theoretical results of Passarelli(1978) and others. In this paper we address the first approach and, in part, the second.

Liao, Liang↗

Understanding Climate Policy Data Needs. NASA Carbon Monitoring System Briefing: Characterizing Flux Uncertainty, Washington D.C., 11 January 2012

Climate policy in the United States is currently guided by public-private partnerships and actions at the local and state levels. This mitigation strategy is made up of programs that focus on energy efficiency, renewable energy, agricultural practices and implementation of technologies to reduce greenhouse gases. How will policy makers know if these strategies are working, particularly at the scales at which they are being implemented? The NASA Carbon Monitoring System (CMS) will provide information on carbon dioxide fluxes derived from observations of earth's land, ocean and atmosphere used in state of the art models describing their interactions. This new modeling system could be used to assess the impact of specific policy interventions on CO2 reductions, enabling an iterative, results-oriented policy process. In January of 2012, the CMS team held a meeting with carbon policy and decision makers in Washington DC to describe the developing modeling system to policy makers. The NASA CMS will develop pilot studies to provide information across a range of spatial scales, consider carbon storage in biomass, and improve measures of the atmospheric distribution of carbon dioxide. The pilot involves multiple institutions (four NASA centers as well as several universities) and over 20 scientists in its work. This pilot study will generate CO2 flux maps for two years using observational constraints in NASA's state-of -the-art models. Bottom-up surface flux estimates will be computed using data-constrained land and ocean models; comparison of the different techniques will provide some knowledge of uncertainty in these estimates. Ensembles of atmospheric carbon distributions will be computed using an atmospheric general circulation model (GEOS-5), with perturbations to the surface fluxes and to transport. Top-down flux estimates will be computed from observed atmospheric CO2 distributions (ACOS/GOSAT retrievals) alongside the forward-model fields, in conjunction with an inverse approach based on the CO2 model of GEOS ]Chem. The forward model ensembles will be used to build understanding of relationships among surface flux perturbations, transport uncertainty and atmospheric carbon concentration. This will help construct uncertainty estimates and information on the true spatial resolution of the top-down flux calculations. The relationship between the top-down and bottom-up flux distributions will be documented. Because the goal of NASA CMS is to be policy relevant, the scientists involved in the flux modeling pilot need to understand and be focused on the needs of the climate policy and decision making community. If policy makers are to use CMS products, they must be aware of the modeling effort and begin to design policies that can be evaluated with information. Improving estimates of carbon sequestered in forests, for example, will require information on the spatial variability of forest biomass that is far more explicit than is presently possible using only ground observations. Carbon mitigation policies being implemented by cities around the United States could be designed with the CMS data in mind, enabling sequential evaluation and subsequent improvements in incentives, structures and programs. The success of climate mitigation programs being implemented in the United States today will hang on the depth of the relationship between scientists and their policy and decision making counterparts. Ensuring that there is two-way communication between data providers and users is important for the success both of the policies and the scientific products meant to support them..

Brown, Molly E.↗

Estimating risk for earth-satellite attenuation prediction

Annual cumulative distributions of attenuation measurements and of rain-rate measurements were obtained from several locations in Europe and the United States. They were analyzed to estimate the year-to-year variability to be associated with a prediction of the expected cumulative distribution. Two models were constructed to estimate variability, an ad hoc model that summarized a number of observations when compared to model predictions and a probabilistic model that applied the ideas of order statistics to the prediction problem when the number of independent attenuation or rain-rate events in a sample year could be estimated. Based on these models, the statistical uncertainty in a model prediction may be estimated. The estimation procedure also provides an answer to the question of the number of years of observation needed to provide an estimate of the empirical distribution with a specified statistical uncertainty.

Crane, Robert K.↗

Stochastic Filtering of Rain Profiles Using Radar, Surface-Referenced Radar, or Combined Radar-Radiometer Measurements

This paper describes a computationally efficient nearly optimal Bayesian algorithm to estimate rain (and drop size distribution) profiles, given a radar reflectivity profile at a single attenuating wavelength. In addition to estimating the averages of all the mutually ambiguous combinations of rain parameters that can produce the data observed, the approach also calculates the n-ns uncertainty in its estimates (this uncertainty thus quantifies "the amount of ambiguity" in the "solution"). The paper also describes a more general approach that can make estimates based on a radar reflectivity profile together with an approximate measurement of the path-integrated attenuation, or a radar reflectivity profile and a set of passive microwave brightness temperatures. This more general "combined" algorithm is currently being adapted for the Tropical Rainfall Measuring Mission.

Haddad, Ziad S.↗

Probabilistic Study Conducted on Sensor-Based Engine Life Calculation

Turbine engine life management is a very complicated process to ensure the safe operation of an engine subjected to complex usage. The challenge of life management is to find a reasonable compromise between the safe operation and the maximum usage of critical parts to reduce maintenance costs. The commonly used "cycle count" approach does not take the engine operation conditions into account, and it oversimplifies the calculation of the life usage. Because of the shortcomings, many engine components are regularly pulled for maintenance before their usable life is over. And, if an engine has been running regularly under more severe conditions, components might not be taken out of service before they exceed their designed risk of failure. The NASA Glenn Research Center and its industrial and academic partners have been using measurable parameters to improve engine life estimation. This study was based on the Monte Carlo simulation of 5000 typical flights under various operating conditions. First a closed-loop engine model was developed to simulate the engine operation across the mission profile and a thermomechanical fatigue (TMF) damage model was used to calculate the actual damage during takeoff, where the maximum TMF accumulates. Next, a Weibull distribution was used to estimate the implied probability of failure for a given accumulated cycle count. Monte Carlo simulations were then employed to find the profiles of the TMF damage under different operating assumptions including parameter uncertainties. Finally, probabilities of failure for different operating conditions were analyzed to demonstrate the importance of a sensor-based damage calculation in order to better manage the risk of failure and on-wing life.

Guo, Ten-Huei↗

Ionizing Radiation Environment on the International Space Station: Performance vs. Expectations for Avionics and Material

The role of structural shielding mass in the design, verification, and in-flight performance of International Space Station (ISS), in both the natural and induced orbital ionizing radiation (IR) environments, is reported. Detailed consideration of the effects of both the natural and induced ionizing radiation environment during ISS design, development, and flight operations has produced a safe, efficient manned space platform that is largely immune to deleterious effects of the LEO ionizing radiation environment. The assumption of a small shielding mass for purposes of design and verification has been shown to be a valid worst-case approximation approach to design for reliability, though predicted dependences of single event effect (SEE) effects on latitude, longitude, SEP events, and spacecraft structural shielding mass are not observed. The Figure of Merit (FOM) method over predicts the rate for median shielding masses of about 10g/cm(exp 2) by only a factor of 3, while the Scott Effective Flux Approach (SEFA) method overestimated by about one order of magnitude as expected. The Integral Rectangular Parallelepiped (IRPP), SEFA, and FOM methods for estimating on-orbit (Single Event Upsets) SEU rates all utilize some version of the CREME-96 treatment of energetic particle interaction with structural shielding, which has been shown to underestimate the production of secondary particles in heavily shielded manned spacecraft. The need for more work directed to development of a practical understanding of secondary particle production in massive structural shielding for SEE design and verification is indicated. In contrast, total dose estimates using CAD based shielding mass distributions functions and the Shieldose Code provided a reasonable accurate estimate of accumulated dose in Grays internal to the ISS pressurized elements, albeit as a result of using worst-on-worst case assumptions (500 km altitude x 2) that compensate for ignoring both GCR and secondary particle production in massive structural shielding.

Koontz, Steven L.↗

Proportion estimation using prior cluster purities

The prior distribution of CLASSY component purities is studied, and this information incorporated into maximum likelihood crop proportion estimators. The method is tested on Transition Year spring small grain segments.

Terrell, G. R.↗

Estimating Total Electron Content Using 1,000+ GPS Receivers

A computer program uses data from more than 1,000 Global Positioning System (GPS) receivers in an Internet-accessible global network to generate daily estimates of the global distribution of vertical total electron content (VTEC) of the ionosphere. This program supersedes an older program capable of processing readings from only about 200 GPS receivers. This program downloads the data via the Internet, then processes the data in three stages. In the first stage, raw data from a global subnetwork of about 200 receivers are preprocessed, station by station, in a Kalman-filter-based least-squares estimation scheme that estimates satellite and receiver differential biases for these receivers and for satellites. In the second stage, an observation equation that incorporates the results from the first stage and the raw data from the remaining 800 receivers is solved to obtain the differential biases for these receivers. The only remaining error sources for which an account cannot be given are multipath and receiver noise contributions. The third stage is a postprocessing stage in which all the processed data are combined and used to generate new data products, including receiver differential biases and global and regional VTEC maps and animations.

Komjathy, Attila↗

A Tracking Sun Photometer Without Moving Parts

This innovation is small, lightweight, and consumes very little electricity as it measures the solar energy attenuated by gases and aerosol particles in the atmosphere. A Sun photometer is commonly used on the Earth's surface, as well as on aircraft, to determine the solar energy attenuated by aerosol particles in the atmosphere and their distribution of sizes. This information is used to determine the spatial and temporal distribution of gases and aerosols in the atmosphere, as well as their distribution sizes. The design for this Sun photometer uses a combination of unique optics and a charge coupled device (CCD) array to eliminate moving parts and make the instrument more reliable. It could be selfcalibrating throughout the year. Data products would be down-welling flux, the direct-diffuse flux ratio, column abundance of gas phase constituents, aerosol optical depth at multiple-wavelengths, phase functions, cloud statistics, and an estimate of the representative size of atmospheric particles. These measurements can be used to obtain an estimate of aerosol size distribution, refractive index, and particle shape. Incident light is received at a light-reflecting (inner) surface, which is a truncated paraboloid. Light arriving from a hemispheric field of view (solid angle 2 steradians) enters the reflecting optic at an entrance aperture at, or adjacent to, the focus of the paraboloid, and is captured by the optic. Most of this light is reflected from an inner surface. The light proceeds substantially parallel to the paraboloid axis, and is detected by an array detector located near an exit aperture. Each of the entrance and exit apertures is formed by the intersection of the paraboloid with a plane substantially perpendicular to the paraboloid axis. Incident (non-reflected) light from a source of limited extent (the Sun) illuminates a limited area on the detector array. Both direct and diffuse illumination may be reflected, or not reflected, before being received on the detector array. As the Sun traverses a path in the sky over some time interval, the track of the Sun can be traced on the detector array. A suitably modified Sun photometer might be used to study the dynamics of an environment on another planet or satellite with an atmosphere.

Strawa, Anthony W.↗

Determination of Eros Physical Parameters for Near Earth Asteroid Rendezvous Orbit Phase Navigation

Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Miller, J. K.↗

Jovian Plasmas Torus Interaction with Europa. Plasma Wake Structure and Effect of Inductive Magnetic Field: 3D Hybrid Kinetic Simulation

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to a variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo Orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy et al., 2007; Shematovich et al., 2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider the models with Oþ þ and Sþ þ background plasma, and various betas for background ions and electrons, and pickup electrons. The majority of O2 atmosphere is thermal with an extended non-thermal population (Cassidy et al., 2007). In this paper, we discuss two tasks: (1) the plasma wake structure dependence on the parameters of the upstream plasma and Europa's atmosphere (model I, cases (a) and (b) with a homogeneous Jovian magnetosphere field, an inductive magnetic dipole and high oceanic shell conductivity); and (2) estimation of the possible effect of an induced magnetic field arising from oceanic shell conductivity. This effect was estimated based on the difference between the observed and modeled magnetic fields (model II, case (c) with an inhomogeneous Jovian magnetosphere field, an inductive magnetic dipole and low oceanic shell conductivity).

Jovian magnetosphere↗

Stratospheric aerosol increase after eruption of Pinatubo observed with lidar and aureolemeter

An increase in the amount of stratospheric aerosol due to the Pinatubo eruption (June 12-15, 1991, 15.14 deg N, 120.35 deg E) was observed from the end of June, 1991 by a lidar in NIES (National Institute for Environmental Studies), Tsukuba (36.0 deg N, 140.1 deg E). After large fluctuations in summer of 1991, the amount of the aerosols increased in mid-September as a result of enhanced transportation from the subtropical region. In autumn and winter of 1991, dense aerosol layers were continuously observed. Aureolemeter (scanning spectral radiometer) measurements were also carried out with lidar measurements and columnar size distribution of stratospheric aerosols was estimated for some cases. Collaborative measurements with the lidar and aureolemeter provided some information on height distribution of the surface area of aerosols in late 1991.

Hayashida, Sachiko↗

Temperature-time distribution and thermal stresses on the RTG fins and shell during water cooling

Radioisotope thermoelectric generator (RTG) packages designed for space missions generally do not require active cooling. However, the heat they generate cannot remain inside of the launch vehicle bay and requires active removal. Therefore, before the Shuttle bay door is closed, the RTG coolant tubes attached to the heat rejection fins must be filled with water, which will circulate and remove most of the heat from the cargo bay. There is concern that charging a system at initial temperature around 200 C with water at 24 C can cause unacceptable thermal stresses in the RTG shell and fins. A computer model is developed to estimate the transient temperature distribution resulting from such charging. The thermal stresses resulting from the temperature gradients do not exceed the elastic deformation limit for the material. Since the simplified mathematical model for thermal stresses tends to overestimate stresses, it is concluded that the RTG can be cooled by introducing water at 24 C to the initially hot fin coolant tubes while the RTG is in the Shuttle cargo bay.

Turner, R. H.↗

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗

A New Bond Albedo for Performing Orbital Debris Brightness to Size Transformations

We have developed a technique for estimating the intrinsic size distribution of orbital debris objects via optical measurements alone. The process is predicated on the empirically observed power-law size distribution of debris (as indicated by radar RCS measurements) and the log-normal probability distribution of optical albedos as ascertained from phase (Lambertian) and range-corrected telescopic brightness measurements. Since the observed distribution of optical brightness is the product integral of the size distribution of the parent [debris] population with the albedo probability distribution, it is a straightforward matter to transform a given distribution of optical brightness back to a size distribution by the appropriate choice of a single albedo value. This is true because the integration of a powerlaw with a log-normal distribution (Fredholm Integral of the First Kind) yields a Gaussian-blurred power-law distribution with identical power-law exponent. Application of a single albedo to this distribution recovers a simple power-law [in size] which is linearly offset from the original distribution by a constant whose value depends on the choice of the albedo. Significantly, there exists a unique Bond albedo which, when applied to an observed brightness distribution, yields zero offset and therefore recovers the original size distribution. For physically realistic powerlaws of negative slope, the proper choice of albedo recovers the parent size distribution by compensating for the observational bias caused by the large number of small objects that appear anomalously large (bright) - and thereby skew the small population upward by rising above the detection threshold - and the lower number of large objects that appear anomalously small (dim). Based on this comprehensive analysis, a value of 0.13 should be applied to all orbital debris albedo-based brightness-to-size transformations regardless of data source. Its prima fascia genesis, derived and constructed from the current RCS to size conversion methodology (SiBAM Size-Based Estimation Model) and optical data reduction standards, assures consistency in application with the prior canonical value of 0.1. Herein we present the empirical and mathematical arguments for this approach and by example apply it to a comprehensive set of photometric data acquired via NASA's Liquid Mirror Telescopes during the 2000-2001 observing season.

Mulrooney, Mark K.↗

Shallow and Deep Latent Heating Modes Over Tropical Oceans Observed with TRMM PR Spectral Latent Heating Data

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of its formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the Tropics. Associated with this rainfall is a vast amount of heat, which is known as latent heat. It arises mainly due to the phase change of water vapor condensing into liquid droplets; three-fourths of the total heat energy available to the Earth's atmosphere comes from tropical rainfall. In addition, fresh water provided by tropical rainfall and its variability exerts a large impact upon the structure and motions of the upper ocean layer. Three-dimensional distributions of latent heating estimated from Tropical Rainfall Measuring Mission Precipitation Radar (TRMM PR)utilizing the Spectral Latent Heating (SLH) algorithm are analyzed. Mass-weighted and vertically integrated latent heating averaged over the tropical oceans is estimated as approx.72.6 J/s (approx.2.51 mm/day), and that over tropical land is approx.73.7 J/s (approx.2.55 mm/day), for 30degN-30degS. It is shown that non-drizzle precipitation over tropical and subtropical oceans consists of two dominant modes of rainfall systems, deep systems and congestus. A rough estimate of shallow mode contribution against the total heating is about 46.7 % for the average tropical oceans, which is substantially larger than 23.7 % over tropical land. While cumulus congestus heating linearly correlates with the SST, deep mode is dynamically bounded by large-scale subsidence. It is notable that substantial amount of rain, as large as 2.38 mm day-1 in average, is brought from congestus clouds under the large-scale subsiding circulation. It is also notable that even in the region with SST warmer than 28 oC, large-scale subsidence effectively suppresses the deep convection, remaining the heating by congestus clouds. Our results support that the entrainment of mid-to-lower-tropospheric dry air, which accompanies the large-scale subsidence is the major factor suppressing the deep convection. Therefore, representation of the realistic entrainment is very important for proper reproduction of precipitation distribution and resultant large-scale circulation.

Takayabu, Yukari N.↗