Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADVECTION”

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 775 records · Page 43

On Animating 2D Velocity Fields

A velocity field. even one that represents a steady state flow implies a dynamical system. Animated velocity fields is an important tool in understanding such complex phenomena. This paper looks at a number of techniques that animate velocity fields and propose two new alternatives, These are texture advection and streamline cycling. The common theme among these techniques is the use of advection on some texture to generate a realistic animation of the velocity field. Texture synthesis and selection for these methods are presented. Strengths and weaknesses of the techniques are also discussed in conjunction with several examples.

Kao, David↗

On Animating 2D Velocity Fields

A velocity field, even one that represents a steady state flow, implies a dynamical system. Animated velocity fields is an important tool in understanding such complex phenomena. This paper looks at a number of techniques that animate velocity fields and propose two new alternatives. These are texture advection and streamline cycling. The common theme among these techniques is the use of advection on some texture to generate a realistic animation of the velocity field. Texture synthesis and selection for these methods are presented. Strengths and weaknesses of the techniques are also discussed in conjunctions with several examples.

Kao, David↗

Algebraic Nonoverlapping Domain Decomposition Methods for Stabilized FEM and FV Discretizations

We consider preconditioning methods for convection dominated fluid flow problems based on a nonoverlapping Schur complement domain decomposition procedure for arbitrary triangulated domains. The triangulation is first partitioned into a number of subdomains and interfaces which induce a natural 2 x 2 partitioning of the p.d.e. discretization matrix. We view the Schur complement induced by this partitioning as an algebraically derived coarse space approximation. This avoids the known difficulties associated with the direct formation of an effective coarse discretization for advection dominated equations. By considering various approximations of the block factorization of the 2 x 2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 using MPI message passing protocol. A number of 2-D CFD calculations will be presented for both scalar advection-diffusion equations and the Euler equations discretized using stabilized finite element and finite volume methods. These results show very good scalability of the preconditioner for various discretizations as the number of processors is increased while the number of degrees of freedom per processor is fixed.

Barth, Timothy J.↗

Are Stronger North-Atlantic Southwesterlies the Forcing to the Late-Winter Warming in Europe?

We examine a possible mechanism leading to late-winter warming and thus to an early spring in Europe. From the NCEP Reanalysis, we extract for the years 1948-1999 ocean-surface winds over the eastern North Atlantic, and air temperatures at the surface, T(sub s), and at the 500 mb level, T(sub 500), in late-winter and spring. T(sub s) is extracted at six European locations, all at 50.5 N, ranging in longitude from 1.9 E (northeastern France) to 26.2 E (Ukraine). To quantify the advection of maritime air into Europe, we evaluate for 3-pentad groups the Index I(sub na) of the southwesterlies at 45 N; 20 W: I(sub na) is the average wind speed at this point if the direction is from the quadrant 180-270 deg (when the direction is different, the contribution counts as zero). In late winter correlations C(sub it) between the Index I(sub na) and the temperature T(sub s) are substantial, up to the 0.6 level, in western Europe (but weaker correlations for Poland and Ukraine). C(sub it) drops sharply by mid-March, taking occasionally negative values subsequently. This drop in C(sub it) indicates that maritime air advection is no longer associated closely with the surface-air warming, the role of immolation becomes important, and thus the drop in C(sub it) marks the arrival of spring. Correlations C(sub i delta) between I(sub na) and our lapse-rate parameter delta, the difference between T(sub s) and T(sub 500), indicate that the flow of warm maritime-air from the North Atlantic into this 'corridor' at 50.5 N is predominantly at low tropospheric level. By computing the best linear fit to I(sub na) and T(sub s), the trends for the period 1948-1999 are evaluated. The trends are appreciable in the second half of February and the first half of March. Our 3-pentad analysis points to the interval from mid-February to mid-March as the end-of-winter period in which the southwesterlies over the eastern North Atlantic become stronger and the surface-air temperatures in Europe rise markedly, the lapse rate becomes steeper, and concurrently the longitudinal temperature gradient between the Somme (France) and the Oder (Germany/Poland border) is reduced by 0.8 C, that is, by 20% of its 1948 value. Our thesis, that the observed late-winter warming and the corollary advancement of spring in Europe resulted at least in part from stronger southwesterlies over the North Atlantic, merits further investigations.

Otterman, J.↗

Lidar Observations of Tropospheric Aerosols Over Northeastern South Africa During the ARREX and SAFARI-2000 Dry Season Experiments

During the ARREX-1999 and SAFARI-2000 Dry Season experiments a micropulse lidar (523 nm) instrument was operated at the Skukuza Airport in northeastern South Africa. The Mar was collocated with a diverse array of passive radiometric equipment. For SAFARI-2000 the processed Mar data yields a daytime time-series of layer mean/derived aerosol optical properties, including extinction-to-backscatter ratios and vertical extinction cross-section profile. Combined with 523 run aerosol optical depth and spectral Angstrom exponent calculations from available CIMEL sun-photometer data and normalized broadband flux measurements the temporal evolution of the near surface aerosol layer optical properties is analyzed for climatological trends. For the densest smoke/haze events the extinction-to-backscatter ratio is found to be between 60-80/sr, and corresponding Angstrom exponent calculations near and above 1.75. The optical characteristics of an evolving smoke event from SAFARI-2000 are extensively detailed. The advecting smoke was embedded within two distinct stratified thermodynamic layers, causing the particulate mass to advect over the instrument array in an incoherent manner on the afternoon of its occurrence. Surface broadband flux forcing due to the smoke is calculated, as is the evolution in the vertical aerosol extinction profile as measured by the Han Finally, observations of persistent elevated aerosol during ARREX-1999 are presented and discussed. The lack of corroborating observations the following year makes these observation; both unique and noteworthy in the scope of regional aerosol transport over southern Africa.

Campbell, James R.↗

Anomalous Heat Budgets in the Interior Pacific Ocean on Seasonal- to -Timescales and Gyre Spacescales

This study quantifies uncertainties in closing the seasonal cycle of diabatic heat storage over the Pacific Ocean from 20 degrees S to 60 degrees N through the synthesis of World Ocean Circulation Experiment (WOCE) products over 7 years from 1993-1999. We utilize WOCE reanalysis products from the following sources: diabatic heat storage (DHS) from the Scripps Institution of Oceanography (SIO); near-surface geostrophic and Ekman currents from the Earth and Space Research (ESR); and air-sea heat fluxes from Comprehensive Ocean-Atmosphere Data Set (COADS), National Centers for Environmental Prediction (NCEP), and European Center for Mid-Range Weather Forecasts (ECMWF). We interpolate these products onto a common grid, allowing the seasonal cycle of DHS to be modeled for comparison with that observed. Everywhere latent heat flux residuals dominate sensible heat flux residuals and shortwave heat flux residuals dominate longwave heat flux residuals, both comparable in magnitude to the residual horizontal heat advection. We find the root-mean-square (RMS) of the differences between observed and model residual DHS tendencies to be less than 15 W per square meters everywhere except in the Kuroshio extension. Comparable COADS and NCEP products perform better than ECMWF products in the extra-tropics, while the NCEP product performs best in the tropics. Radiative and turbulent air-sea heat flux residuals computed from ship-born measurements perform better than those computed from satellite cloud and wind measurements. Since the RMS differences derive largely from biases in measured wind speed and cloud fraction, least-squares minimization is used to correct the residual Ekman heat advection and air-sea heat flux. Minimization reduces RMS differences less than 5 W per square meters except in the Kuroshio extension, suggesting how winds, clouds, and exchange coefficients in the NCEP, ECMWF, and ESR products can be improved.

White, Warren↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

CODA: Coordinated Observations of Dynamics in Aurora

The grant funds were provided to carry out chemical tracer wind and turbulence measurements in the mesosphere and lower thermosphere as part of the Coordinated Observations of Dynamics in Aurora (CODA) experiment for which Dr. Andy Christensen of the Aerospace Corp. was the Principal Investigator. Clemson University was a co-investigator institution on the project and was tasked with designing, building, and testing a total of four chemical tracer release payloads which were used to release trimethyl aluminum (TMA) to measure the winds and turbulent diffusion. All four payloads were flown from the rocket range at Poker Flat, Alaska, in support of instrumented payloads flown from the same location. The CODA experiment was designed to investigate the changes in the chemical composition in the lower thermosphere during an auroral substorm event. In particular, the objective was to determine the relative importance of two competing mechanisms. The first mechanism is associated with the advection of air from the polar cap into the auroral oval. A flow in that direction is typical in the postmidnight sector of the oval. As the air moves from the region of negligible forcing in the polar cap to the region of strong forcing and heating in the oval, changes in the composition are expected to occur. Alternatively, enhanced local mixing will also cause such changes. Since the auroral forcing accelerates the neutral winds in the lower E region, but not at lower altitudes, the forcing results in large shears that can be large enough to be unstable in the Richardson number sense. The unstable shears are expected to generate turbulence and therefore enhanced mixing. Such mixing will also create composition changes similar to what is observed. In the CODA experiment an instrumented payload designed and built by the Aerospace Corporation was used to measure the composition and plasma parameters needed to define the auroral forcing. In addition, three chemical tracer payloads were to be flown along three trajectories widely separated in azimuth to map out the winds and wind gradients. The measurements would thus provide information about the local composition changes above the launch site as well as detailed information about the advection and location of stable or unstable shears both locally and over a broader area in the vicinity of the launch site.

Larsen, M. F.↗

Various Numerical Applications on Tropical Convective Systems Using a Cloud Resolving Model

In recent years, increasing attention has been given to cloud resolving models (CRMs or cloud ensemble models-CEMs) for their ability to simulate the radiative-convective system, which plays a significant role in determining the regional heat and moisture budgets in the Tropics. The growing popularity of CRM usage can be credited to its inclusion of crucial and physically relatively realistic features such as explicit cloud-scale dynamics, sophisticated microphysical processes, and explicit cloud-radiation interaction. On the other hand, impacts of the environmental conditions (for example, the large-scale wind fields, heat and moisture advections as well as sea surface temperature) on the convective system can also be plausibly investigated using the CRMs with imposed explicit forcing. In this paper, by basically using a Goddard Cumulus Ensemble (GCE) model, three different studies on tropical convective systems are briefly presented. Each of these studies serves a different goal as well as uses a different approach. In the first study, which uses more of an idealized approach, the respective impacts of the large-scale horizontal wind shear and surface fluxes on the modeled tropical quasi-equilibrium states of temperature and water vapor are examined. In this 2-D study, the imposed large-scale horizontal wind shear is ideally either nudged (wind shear maintained strong) or mixed (wind shear weakened), while the minimum surface wind speed used for computing surface fluxes varies among various numerical experiments. For the second study, a handful of real tropical episodes (TRMM Kwajalein Experiment - KWAJEX, 1999; TRMM South China Sea Monsoon Experiment - SCSMEX, 1998) have been simulated such that several major atmospheric characteristics such as the rainfall amount and its associated stratiform contribution, the Qlheat and Q2/moisture budgets are investigated. In this study, the observed large-scale heat and moisture advections are continuously applied to the 2-D model. The modeled cloud generated from such an approach is termed continuously forced convection or continuous large-scale forced convection. A third study, which focuses on the respective impact of atmospheric components on upper Ocean heat and salt budgets, will be presented in the end. Unlike the two previous 2-D studies, this study employs the 3-D GCE-simulated diabatic source terms (using TOGA COARE observations) - radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the Ocean mixed-layer (OML) model.

Shie, C.-L.↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Nonlinear Gulf Stream Interaction with the Deep Western Boundary Current System: Observations and a Numerical Simulation

Gulf Stream (GS) separation near its observed Cape Hatteras (CH) separation location, and its ensuing path and dynamics, is a challenging ocean modeling problem. If a model GS separates much farther north than CH, then northward GS meanders, which pinch off warm core eddies (rings), are not possible or are strongly constrained by the Grand Banks shelfbreak. Cold core rings pinch off the southward GS meanders. The rings are often re-absorbed by the GS. The important warm core rings enhance heat exchange and, especially, affect the northern GS branch after GS bifurcation near the New England Seamount Chain. This northern branch gains heat by contact with the southern branch water upstream of bifurcation, and warms the Arctic Ocean and northern seas, thus playing a major role in ice dynamics, thermohaline circulation and possible global climate warming. These rings transport heat northward between the separated GS and shelf slope/Deep Western Boundary Current system (DWBC). This region has nearly level time mean isopycnals. The eddy heat transport convergence/divergence enhances the shelfbreak and GS front intensities and thus also increases watermass transformation. The fronts are maintained by warm advection by the Florida Current and cool advection by the DWBC. Thus, the GS interaction with the DWBC through the intermediate eddy field is climatologically important.

Dietrich, David E.↗

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

Use of GOES, SSM/I, TRMM Satellite Measurements Estimating Water Budget Variations in Gulf of Mexico - Caribbean Sea Basins

This study presents results from a multi-satellite/multi-sensor retrieval system designed to obtain the atmospheric water budget over the open ocean. A combination of 3ourly-sampled monthly datasets derived from the GOES-8 5-channel Imager, the TRMM TMI radiometer, and the DMSP 7-channel passive microwave radiometers (SSM/I) have been acquired for the combined Gulf of Mexico-Caribbean Sea basin. Whereas the methodology has been tested over this basin, the retrieval system is designed for portability to any open-ocean region. Algorithm modules using the different datasets to retrieve individual geophysical parameters needed in the water budget equation are designed in a manner that takes advantage of the high temporal resolution of the GOES-8 measurements, as well as the physical relationships inherent to the TRMM and SSM/I passive microwave measurements in conjunction with water vapor, cloud liquid water, and rainfall. The methodology consists of retrieving the precipitation, surface evaporation, and vapor-cloud water storage terms in the atmospheric water balance equation from satellite techniques, with the water vapor advection term being obtained as the residue needed for balance. Thus, the intent is to develop a purely satellite-based method for obtaining the full set of terms in the atmospheric water budget equation without requiring in situ sounding information on the wind profile. The algorithm is validated by cross-checking all the algorithm components through multiple- algorithm retrieval intercomparisons. A further check on the validation is obtained by directly comparing water vapor transports into the targeted basin diagnosed from the satellite algorithms to those obtained observationally from a network of land-based upper air stations that nearly uniformly surround the basin, although it is fair to say that these checks are more effective m identifying problems in estimating vapor transports from a leaky operational radiosonde network than in verifying the transport estimates determined from the satellite algorithm system Total columnar atmospheric water budget results are presented for an extended annual cycle consisting of the months of October-97, January-98, April-98, July-98,October-98, and January 1999. These results are used to emphasize the changing relationship in E-P, as well as in the varying roles of storage and advection in balancing E-P both on daily and monthly time scales and on localized and basin space scales. Results from the algorithm-to-algorithm intercomparisons are also presented in the context of sensitivity testing to help understand the intrinsic uncertainties in evaluating the water budget terms by an all-satellite algorithm approach.

Smith, Eric A.↗

Simulated convective systems using a cloud resolving model: Impact of large-scale temperature and moisture forcing using observations and GEOS-3 reanalysis

The GCE (Goddard Cumulus Ensemble) model, which has been developed and improved at NASA Goddard Space Flight Center over the past two decades, is considered as one of the finer and state-of-the-art CRMs (Cloud Resolving Models) in the research community. As the chosen CRM for a NASA Interdisciplinary Science (IDS) Project, GCE has recently been successfully upgraded into an MPI (Message Passing Interface) version with which great improvement has been achieved in computational efficiency, scalability, and portability. By basically using the large-scale temperature and moisture advective forcing, as well as the temperature, water vapor and wind fields obtained from TRMM (Tropical Rainfall Measuring Mission) field experiments such as SCSMEX (South China Sea Monsoon Experiment) and KWAJEX (Kwajalein Experiment), our recent 2-D and 3-D GCE simulations were able to capture detailed convective systems typical of the targeted (simulated) regions. The GEOS-3 [Goddard EOS (Earth Observing System) Version-3] reanalysis data have also been proposed and successfully implemented for usage in the proposed/performed GCE long-term simulations (i.e., aiming at producing massive simulated cloud data -- Cloud Library) in compensating the scarcity of real field experimental data in both time and space (location). Preliminary 2-D or 3-D pilot results using GEOS-3 data have generally showed good qualitative agreement (yet some quantitative difference) with the respective numerical results using the SCSMEX observations. The first objective of this paper is to ensure the GEOS-3 data quality by comparing the model results obtained from several pairs of simulations using the real observations and GEOS-3 reanalysis data. The different large-scale advective forcing obtained from these two kinds of resources (i.e., sounding observations and GEOS-3 reanalysis) has been considered as a major critical factor in producing various model results. The second objective of this paper is therefore to investigate and present such an impact of large-scale forcing on various modeled quantities (such as hydrometeors, rainfall, and etc.). A third objective is to validate the overall GCE 3-D model performance by comparing the numerical results with sounding observations, as well as available satellite retrievals.

Shie, C.-L.↗

Monthly-Diurnal Water Budget Variability Over Gulf of Mexico-Caribbean Sea Basin from Satellite Observations

This study presents results from a multi-satellite/multi-sensor retrieval system design d to obtain the atmospheric water budget over the open ocean. A combination of hourly-sampled monthly datasets derived from the GOES-8 5-channel Imager, the TRMM TMI radiometer, and the DMSP 7-channel passive microwave radiometers (SSM/I) have been acquired for the combined Gulf of Mexico-Caribbean Sea basin. Whereas the methodology has been tested over this basin, the retrieval system is designed for portability to any open-ocean region. Algorithm modules using the different datasets to retrieve individual geophysical parameters needed in the water budget equation are designed in a manner that takes advantage of the high temporal resolution of the GOES-8 measurements, as well as the physical relationships inherent to the TRMM and SSM/I passive microwave measurements in conjunction with water vapor, cloud liquid water, and rainfall. The methodology consists of retrieving the precipitation, surface evaporation, and vapor-cloud water storage terms in the atmospheric water balance equation from satellite techniques, with the water vapor advection term being obtained as the residue needed for balance. Thus, the intent is to develop a purely satellite-based method for obtaining the full set of terms in the atmospheric water budget equation without requiring in situ sounding information on the wind profile. The algorithm is validated by cross-checking all the algorithm components through multiple-algorithm retrieval intercomparisons. A further check on the validation is obtained by directly comparing water vapor transports into the targeted basin diagnosed from the satellite algorithms to those obtained observationally from a network of land-based upper air stations that nearly uniformly surround the basin, although it is fair to say that these checks are more effective in identifying problems in estimating vapor transports from a "leaky" operational radiosonde network than in verifying the transport estimates determined from the satellite algorithm system. Total columnar atmospheric water budget results are presented for an extended annual cycle consisting of the months of October-97, January-98, April-98, July-98,October-98, and January- 1999. These results are used to emphasize the changing relationship in E-P, as well as in the varying roles of storage and advection in balancing E-P both on daily and monthly time scales and on localized and basin space scales. Results from the algorithm-to-algorithm intercomparisons are also presented in the context of sensitivity testing to help understand the intrinsic uncertainties in evaluating the water budget terms by an all-satellite algorithm approach.

Smith, E. A.↗

Stable Low Cloud Phase II: Nocturnal Event Study

This report describes the work done by the Applied Meteorology Unit (AMU) in developing a database of nights that experienced rapid (< 90 minutes) low cloud formation in a stable atmosphere, resulting in ceilings at the Shuttle Landing Facility (TTS) that violated Space Shuttle Flight Rules (FR). This work is the second phase of a similar AMU task that examined the same phenomena during the day. In the first phase of this work, the meteorological conditions favoring the rapid formation of low ceilings include the presence of any inversion below 8000 ft, high relative humidity (RH) beneath the inversion and a clockwise turning of the winds from the surface to the middle troposphere (-15000 ft). The AMU compared and contrasted the atmospheric and thermodynamic conditions between nights with rapid low ceiling formation and nights with low ceilings resulting from other mechanisms. The AMU found that there was little to discern between the rapidly-forming ceiling nights and other low ceiling nights at TTS. When a rapid development occurred, the average RH below the inversions was 87% while non-events had an average RH of 79%. One key parameter appeared to be the vertical wind profile in the Cape Canaveral, FL radiosonde (XMR) sounding. Eighty-three percent of the rapid development events had veering winds with height from the surface to the middle troposphere (-15,000 ft) while 61% of the non-events had veering winds with height. Veering winds indicate a warm-advection regime, which supports large-scale rising motion and ultimately cloud formation in a moist environment. However, only six of the nights (out of 86 events examined) with low cloud ceilings had an occurrence of rapidly developing ceilings. Since only 7% rapid development events were observed in this dataset, it is likely that rapid low cloud development is not a common occurrence during the night, or at least not as common as during the day. In the AMU work on the daytime rapid low cloud development (Case and Wheeler 2005), nearly 30% of the low cloud ceiling cases investigated were identified as rapidly developing events. Forecasters at the Space Meteorology Group (SMG) issue 30 to 90 minute forecasts for low cloud ceilings at TTS to support Space Shuttle landings. Mission verification statistics have shown ceilings to be the number one forecast challenge. More specifically, forecasters at SMG are concerned with any rapidly developing clouds/ceilings below 8000 R in a stable, capped thermodynamic environment. Therefore, the AMU was tasked to examine archived events of rapid stable cloud formation resulting in ceilings below 8000 ft, and document the atmospheric regimes favoring this type of cloud development. The AMU examined the cool season months of November to March during the years of 1994-2005 for nights that had low-level inversions and rapid, stable low cloud formation that resulted in ceilings violating the Space Shuttle FR. The AMU wrote and modified existing code to identify inversions from the evening and morning XMR radiosonde during the cool season and output pertinent sounding information. They parsed all days with cloud ceilings below 8000 ft at TTS, forming a database of possible rapidly-developing low ceiling events. Nights with precipitation or noticeable fog burn-off situations were excluded from the database. Only the nighttime hours were examined for possible ceiling development events since the daytime events were examined in the first phase of this work. The report presents one sample case of rapidly-developing low cloud ceilings. The case depicts the representative meteorological and thermodynamic characteristics of such events. The case also illustrates how quickly the cloud decks can develop, sometimes forming in 30 minutes or less. The report also summarizes the composite meteorological conditions for 6 event nights with rapid low cloud ceiling formation and 80 non-events nights consisting of advection or widespread low cloud ceilings. The teorological conditions were quite similar for both the event and non-event nights, since both types of nights experienced low cloud ceilings. Both types of nights had a relatively moist environment beneath an inversion based below 8000 ft.

Bauman, William H., III↗

A Stellar-mass Black Hole in the Ultra-luminous X-ray Source M82 X-1

We have analyzed the archival XMM-Newton data of the archetypal Ultra-Luminous X-ray Source (ULX) M82 X-1 with an LO5 ksec exposure when the source was in the steady state. Thanks to the high photon statistics from the large effective area and long exposure, we were able to discriminate different X-ray continuum spectral models. Neither the standard accretion disk model (where the radial dependency of the disk effective temperature is T(r) proportional to r(sup -3/4)) nor a power-law model gives a satisfactory fit. In fact, observed curvature of the M82 X-1 spectrum was just between those of the two models. When the exponent of the radial dependence (p in T(r) proportional to r(sup -P)) of the disk temperature is allowed to be free, we obtained p = 0.61 (sup +0.03)(sub -0.02). Such a reduction of p from the standard value 3/4 under extremely high mass accretion rates is predicted from the accretion disk theory as a consequence of the radial energy advection. Thus, the accretion disk in M82 X-1 is considered to be in the Slim disk state, where an optically thick Advection Dominant Accretion Flow (ADAF) is taking place. We have applied a theoretical slim disk spectral model to M82 X-1, and estimated the black hole mass approximately equal to 19 - 32 solar mass. We conclude that M82 X-1 is a stellar black hole which has been produced through evolution of an extremely massive star, shining at a several times the super-Eddington luminosity.

Okajima, Takashi↗