Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional probability distribution”

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

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗

Shock recovery analogs and the origin of mesosiderites

The origin of mesosiderites, which consist of approximately equal-weight proportions of Fe-Ni metal and silicates (gabbros, basalts, orthopyroxenites, dunites), remains an interesting and complex problem in meteoritics. There is general agreement that multiple impact events were probably involved in the formation of these brecciated stony-iron meteorites, but given the heterogeneity of mesosiderites, additional processes have been invoked to explain the unique and intricate textural and compositional makeup of mesosiderites. We conducted a series of shock recovery experiments to test the impact event(s) scenario. The results indicated significant similarities between the shocked analogs and many mesosiderites. We have compared our analogs with a suite of thin sections of Barea mesosiderite. I have conducted a series of flash heating experiments in which equal-weight proportions of gabbro and stainless steel (SS304) powders were compressed into small charges and heated under reducing conditions for short times. These experiments were used to bracket localized, peak postshock temperatures in our analog shots and to compare the mixing relations between the silicate and metal. The shock recovery experiments used porous metal-silicate powder starting mixtures, therefore our experiments are most analogous to an impact scenario where the target is an asteroidal regolith surface composed of a loose mixture of Fe-Ni metal and heterogeneous silicates. Analog experiments may really describe a secondary impact process similar to the late-stage, localized impact melting event. This leaves one of the crucial questions about mesosiderite genesis unanswered, namely what is the source of the Fe-Ni metal that is so intimately distributed in these meteorites?

Rowan, L. R.↗

Theoretical Basis for the Surface Spectral Reflectance Relationships Used in the MODIS Aerosol Algorithm

The analysis of data from the MODIS instrument on the Terra platform to derive global distribution of aerosols assumes a set of relationships between the blue, rho (sub blue), the red, rho (sub red), and 2.1 micrometers, rho (sub 2.1), spectral channels. These relations have been established from a series of measurements indicating that rho (sub blue) approximately 0.5 rho (sub red) approximately 0.25 rho (sub 2.1). Here we use a model to describe the transfer of radiation through a vegetation canopy composed of randomly oriented leaves to assess the theoretical foundations for these relationships. The influence of varying fractional vegetation coverage is simulated simply as a linear combination of pure soil and pure vegetation conditions, also known as Independent Pixel Approximation (IPA). Calculations for a wide range of leaf area indices and vegetation fractions show that rho (sub blue) is consistently about 1/4 of rho (sub 2.1) as used by MODIS for the whole range of analyzed cases, except for very dark soils, such as those found in burn scars. For its part, the ratio rho (sub red)/rho (sub 2.1) varies from less than the empirically derived value of 1/2 for dense and dark vegetation (rho (sub 2.1) less than 0.1), to more than 1/2 for bright mixture of soil and vegetation. This is in agreement with measurements over uniform dense vegetation, but not with measurements over mixed dark scenes. In the later case, the discrepancy is probably mitigated by shadows due to uneven canopy and terrain on a large scale. It is concluded that the value of this ratio should ideally be made dependent on the land cover type in the operational processing of MODIS data, especially over dense forests.

Kaufman, Yoram J.↗

An X-ray and optical study of the supernova remnant W44

We report the results of a 8000 s observation of the supernova remnant W44 using the ROSAT Position Sensitive Proportional Counter (PSPC). The image shows the same centrally peaked morphology observed by the Einstein IPC and contrasts with the shell-like radio morphology. The eastern limb shows a lack of X-ray emission within the radio shell, probably due to the interaction between the Supernova Remnants (SNR) and a molecular cloud. No counterpart to the pulsar 1853 + 01 in W44 has been detected, with L(sub X) less than 1.3 x 10(exp 32) ergs/s in the 0.2 to 2.4 keV band. The spectral analysis of the central part of W44, combining EXOSAT ME and Einstein SSS data, shows that the shocked plasma has not reached ionization equilibrium. The best nonequilibrium fit to PSPC, ME, and SSS spectra gives Eta = 10(exp 51) ergs cm(exp -6), T(sub s) = 10(exp 7) K with T(sub e) = T(sub i), suggesting conditions are approaching ionization equilibrium. There is no evidence of enhanced abundances of Mg, Si, S, or Fe. The variation of temperature and column density was obtained region by region using the PSPC and Einstein IPC. The temperature is largely uniform over the remnant, but strong column density variations are found to be consistent with molecular clouds in the line of sight. An evaporation model with a two-phase interstellar medium structure of clumps and interclump gas (White & Long 1991) can explain the X-ray centrally peaked morphology of W44. The clumps remaining behind a SN shock provide a reservoir of material, and evaporat e to increase the density of X-ray emitting gas in the interior of a SNR. The uniform temperature distribution of W44 strongly supports the predictions of this model. In addition, mosaiced H alpha and (S II) images of W44, taken using the prime focus universal extragalactic instrument (PFUEI) camera on the Palomar 60 sec telescope, reveal the first discovery of optical filaments (both H alpha and (S II)) in the northwestern and southeastern portion of the remnant, within the X-ray emitting region. The optical filaments and the X-ray image showing locally brighter emission and clumps along the filaments suggest both are produced by the interaction between the supernova shock front and regions of enhanced ambient density.

Rho, Jeongee↗

Stratospheric Sulfuric Acid and Black Carbon Aerosol Measured During POLARIS and its Role in Ozone Chemistry

Stratospheric aerosol can affect the environment in three ways. Sulfuric acid aerosol have been shown to act as sites for the reduction of reactive nitrogen and chlorine and as condensation sites to form Polar Stratospheric Clouds, under very cold conditions, which facilitate ozone depletion. Recently, modeling studies have suggested a link between BCA (Black Carbon Aerosol) and ozone chemistry. These studies suggest that HNO3, NO2, and O3 may be reduced heterogeneously on BCA particles. The ozone reaction converts ozone to oxygen molecules, while HNO3 and NO2 react to form NOx. Finally, a buildup of BCA could reduce the single-scatter albedo of aerosol below a value of 0.98, a critical value that has been postulated to change the effect of stratospheric aerosol from cooling to warming. Correlations between measured BCA amounts and aircraft usage have been reported. Attempts to link BCA to ozone chemistry and other stratospheric processes have been hindered by questions concerning the amount of BCA that exists in the stratosphere, the magnitude of reaction probabilities, and the scarcity of BCA measurements. The Ames Wire Impactors (AWI) participated in POLARIS as part of the complement of experiments on the NASA ER-2. One of our main objectives was to determine the amount of aerosol surface area, particularly BCA, available for reaction with stratospheric constituents and assess if possible, the importance of these reactions. The AWI collects aerosol and BCA particles on thin Palladium wires that are exposed to the ambient air in a controlled manner. The samples are returned to the laboratory for subsequent analysis. The product of the AWI analysis is the size, surface area, and volume distributions, morphology and elemental composition of aerosol and BCA. This paper presents results from our experiments during POLARIS and puts these measurements in the context of POLARIS and other missions in which we have participated. It describes modifications to the AWI data analysis procedures in which the collection of BCA is modeled as a fractal aggregate. The new analysis results in an increase in BCA surface area of approximately 24 and an increase in mass of 7-10 from the previous method. For the current study, BCA surface area is used in computer models that attempt to predict measured NOx/NOy ratios and O3 depletion rates. Inclusion of the HNO3 reaction with BCA in one model tends to improve the agreement of calculated to measured NOx/NOy ratio. However, it was found that these trends are viable only if the reactions are catalytic.

Strawa, Anthony W.↗

Ventifact Formation in the Mojave Desert: Field Analogs for Martian Processes

Two field studies in the Mojave Desert, California, shed light on processes of ventifact formation. The field sites are located on a ridge at Little Cowhole Mountain, which lies approximately 12 km south of Baker, and on an unnamed ridge situated along the northern boundary of the Mojave River Sink (Rasor Road site). The rocks at Little Cowhole Mountain are a blue-grey marble/dolomite, whereas those at Rasor Road are Miocene volcanic rocks (basalt). At both sites the abrasive agent is a fine-grained aeolian sand which was probably derived largely from the Mojave River. There are minimal modem inputs of sand to either site: abrasion occurs as a result of unique climatic and topographic conditions which allow pre-existing sand to be recycled from one aspect of the ridge to the other. Climatic conditions are well suited for ventifact formation. Owing, to the dry climate (marked by low average relative humidity, infrequent dew, and low annual rainfall), rates of chemical weathering are low. Where resurfacing of the rocks by sand abrasion proceeds at a rate greater than weathering, the ventifacts are considered "active." Active ventifacts are found atop and straddling the ridge crests, in the zone of maximum wind velocity and sediment supply. Inactive ventifacts occur where modem weathering Processes exceed abrasion rates; principally on the basal two-thirds of the hillslope, where wind velocity and sediment supply are lower. At intermediate locations between the slope base and crest, ventifacts are either active or inactive, depending on local conditions. The presence of relict ventifacts at the study sites, as well as elsewhere in the eastern Mojave Desert, suggests that the conditions for venti fact formation must have been more intense and extensive in the past. Together, the studies illustrate that the processes that interact to form ventifacts are highly complex, and must be studied at many scales. Small-scale effects, such as local topography, plant cover, or even the spatial distribution of boulders, strongly influence the formation of each individual ventifact. Mesoscale effects (hillslope form, total availability of sediment, seasonality of winds, etc.) determine the distribution of overall erosional energy and the location of active and relict forms.

Laity, J. E.↗

The galactic globular cluster system

We explore correlations between various properties of Galactic globular clusters, using a database on 143 objects. Our goal is identify correlations and trends which can be used to test and constrain theoretical models of cluster formation and evolution. We use a set of 13 cluster parameters, 9 of which are independently measured. Several arguments suggest that the number of clusters still missing in the obscured regions of the Galaxy is of the order of 10, and thus the selection effects are probably not severe for our sample. Known clusters follow a power-law density distribution with a slope approximately -3.5 to -4, and an apparent core with a core radius approximately 1 kpc. Clusters show a large dynamical range in many of their properties, more so for the core parameters (which are presumably more affected by dynamical evolution) than for the half-light parameters. There are no good correlations with luminosity, although more luminous clusters tend to be more concentrated. When data are binned in luminosity, several trends emerge: more luminous clusters tend to have smaller and denser cores. We interpret this as a differential survival effect, with more massive clusters surviving longer and reaching more evolved dynamical states. Cluster core parameters and concentrations also correlate with the position in the Galaxy, with clusters closer to the Galactic center or plane being more concentrated and having smaller and denser cores. These trends are more pronounced for the fainter (less massive) clusters. This is in agreement with a picture where tidal shocks form disk or bulge passages accelerate dynamical evolution of clusters. Cluster metallicities do not correlate with any other parameter, including luminosity and velocity dispersion; the only detectable trend is with the position in the Galaxy, probably reflecting Zinn's disk-halo dichotomy. This suggests that globular clusters were not self-enriched systems. Velocity dispersions show excellent correlations with luminosity and surface brightness. Their origin is not well understood, but they may well reflect initial conditions of cluster formation, and perhaps even be used to probe the initial density perturbation spectrum on a approximately 10(exp 6) solar mass scale. Core radii and concentrations play a role of a 'second parameter' in these correlations. While a global manifold of cluster properties has a high statistical dimensionality (D greater than 4), a subset of structural, photometric, and dynamical parameters forms a statistically three-dimensional family, as expected from objects following King models; we propose to call this set of quantities the King Manifold. Some of the observed correlations may be usable as distance indicator relations for globular clusters.

Djorgovski, S.↗

Utilitarian models of the solar nebula

Models of the primitive solar nebula based on a combination of theory, observations of T Tauri stars, and global conservation laws are presented. The models describe the motions of nebular gas, mixing of interstellar material during the formation of the nebula, and evolution of thermal structure in terms of several characteristic parameters. The parameters describe key aspects of the protosolar cloud (its rotation rate and collapse rate) and the nebula (its mass relative to the Sun, decay time, and density distribution). For most applications, the models are heuristic rather than predicted. Their purpose is to provide a realistic context for the interpretation of solar system data, and to distinquish those nebular characteristics that can be specified with confidence, independently of the assumtions of particular models, form those that are poorly constrained. It is demonstrated that nebular gas typically experienced large radial excursions during the evolution of the nebula and that both inward and outward mean radial velocities on the order of meters per second occured in the terrestrial planet region, with inward velocities predominant for most ofthe evolution. However, the time history of disk size, surface density, and radial velocities are sensitive to the total angular momentun of the protosolar cloud, which cannot be constrained by purely theoretical considerations.It is shown that a certain amount of 'formational' mixing of interstellar material was an inevitable consequenc of nebular mass and angular momentum transport during protostellar collapse, regardless of the specific transport mechanisms invloved. Even if the protosolar cloud was initially homogeneous, this mixing was important because it had the effect of mingling presolar material that had experienced different degrees of thermal processing during collapse and passage through the accertion shock. Nebular thermal structure is less sensitive to poorly constrained parameters than is dynamical history. A simple criterion is derived for the condition that silicate grains are evaporated at midplane, and it is argued that this condition was probably fulfilled early in nebular history. Cooling of a hot nebula due tocoagulation of dust and consequent local reduction of optical depth is examined, and it is shown how such a process leads naturally to an enrichment of rock-forming elements in the gas phase.

Cassen, Patrick↗

Probability of Decompression Sickness and Venous Gas Emboli from 49 NASA Hypobaric Chamber Tests with Reference to Exploration Atmosphere

Introduction: Decompression sickness (DCS) is a complex biophysical event; it combines human perception of pain, for instance, and the presence of a gas phase in the tissues. Living tissues are complex and dynamic. Micronuclei and later bubbles may or may not form given what appears to be the same conditions. Even when bubbles grow, symptoms may or may not develop under what appears to be the same conditions. Therefore, at this time it is appropriate to consider DCS as a probabilistic rather than a deterministic event. Methods: Probabilistic models about hypobaric DCS and venous gas emboli (VGE) require a large amount of quality research data, a definition of decompression dose using physical and physiologic variables, and a flexible analytical approach that can quantify the association between each outcome and all covariates of interest (assuming independence between DCS and VGE) and then ultimately be extended to acknowledge dependencies between DCS and VGE. Our DCS and VGE data are from 1,031 hypobaric decompressions from 1983 to 2016. A total of 577 humans participated in 49 hypobaric chamber tests to evaluate denitrogenation procedures used by astronauts in the Space Shuttle and International Space Station programs. We defined decompression dose as the ratio of computed nitrogen tension in a theoretical 360-minute half-time compartment to ambient pressure, which accounts for denitrogenation and exposure pressure as well as explanatory variables such as age, sex, body mass index, and the presence or absence of ambulation as part of exercise at the exposure pressure. A parametric survival model, using a log-logistic distribution, was used to quantify the time to development of DCS, VGE, and Grade IV VGE. Results: Our survival estimates are applicable to simple hypobaric decompressions, such as depressurizations in 5 to 30 minutes to exposure pressures between 4 and 10 pounds per square inch absolute (psia) and after minutes to hours of denitrogenation, either under resting or exercise conditions to accelerated denitrogenation. The regressions are applicable to exposures between 2 to 6 hours and under conditions of ambulation or no ambulation as part of exercise at the test pressure. We estimate that an exposure to 4.3 psia with simulated extravehicular activity (EVA) that includes ambulation after equilibration to the exploration atmosphere at 8.2 psia with a 34% oxygen atmosphere will result in 3.1% DCS (1.8% to 5.2 95% confidence interval), 23.2% VGE (16.7 to 31.2%), and 8.5% Grade IV VGE (4.7 to 14.7%) in equal samples of men and women exposed for 6 hours. Discussion: Probabilistic models for DCS, VGE, and Grade IV VGE can be used to inform those that plan future EVAs. Their applications are useful to quantify the risk of DCS and VGE in astronauts that perform EVAs in low-pressure space suits while in space or while exploring the surfaces of the moon or Mars.

Johnny Conkin↗

Transient Reliability Analysis Capability Developed for CARES/Life

The CARES/Life software developed at the NASA Glenn Research Center provides a general-purpose design tool that predicts the probability of the failure of a ceramic component as a function of its time in service. This award-winning software has been widely used by U.S. industry to establish the reliability and life of a brittle material (e.g., ceramic, intermetallic, and graphite) structures in a wide variety of 21st century applications.Present capabilities of the NASA CARES/Life code include probabilistic life prediction of ceramic components subjected to fast fracture, slow crack growth (stress corrosion), and cyclic fatigue failure modes. Currently, this code can compute the time-dependent reliability of ceramic structures subjected to simple time-dependent loading. For example, in slow crack growth failure conditions CARES/Life can handle sustained and linearly increasing time-dependent loads, whereas in cyclic fatigue applications various types of repetitive constant-amplitude loads can be accounted for. However, in real applications applied loads are rarely that simple but vary with time in more complex ways such as engine startup, shutdown, and dynamic and vibrational loads. In addition, when a given component is subjected to transient environmental and or thermal conditions, the material properties also vary with time. A methodology has now been developed to allow the CARES/Life computer code to perform reliability analysis of ceramic components undergoing transient thermal and mechanical loading. This means that CARES/Life will be able to analyze finite element models of ceramic components that simulate dynamic engine operating conditions. The methodology developed is generalized to account for material property variation (on strength distribution and fatigue) as a function of temperature. This allows CARES/Life to analyze components undergoing rapid temperature change in other words, components undergoing thermal shock. In addition, the capability has been developed to perform reliability analysis for components that undergo proof testing involving transient loads. This methodology was developed for environmentally assisted crack growth (crack growth as a function of time and loading), but it will be extended to account for cyclic fatigue (crack growth as a function of load cycles) as well.

Nemeth, Noel N.↗

Critical review of electrical conductivity measurements and charge distribution analysis of magnesium oxide

The electrical conductivity sigma of MgO single crystals shows a sharp increase at 500-800 C, in particular of sigma surface, generally attributed to surface contamination. Charge Distribution Analysis (CDA), a new technique providing information on fundamental properties that was previously unavailable, allows for the determination of surface charges, their sign and associated internal electric field. Data on 99.99% purity, arc-fusion grown MgO crystals show that mobile charge carriers start to appear in the bulk of the MgO crystals between 200 and 400 C when sigma (measured by conventional techniques) is in t he 10(exp -14) to 10(exp -16) /omega/cm range. Above 500 C, as sigma increases to 10(exp -6) to 10(exp -7)/omega/cm, more charges appear giving rise to a strong positive surface charge supported by a strong internal field. This indicates that charges are generated in the bulk and diffuse to the surface by an internally controlled process. On the basis of their positive sign they are identified as holes (defect electrons). Because of the low cation content of these very pure MgO crystals, theses holes cannnot be associated with transition metal impurties. Instead, they are associated with the O(2-) sublattice, e.g. consist of O(-) states or positive holes. This conclusion is supported by magnetic susceptibility data showing the appearance of 1000 +/- 500 ppm paramagnetic species between 200-500 C. The magnetic data are consistent with strongly coupled, diamagnetic O(-) pairs below 200-500 C, chemically equivalent to peroxy anions, O2(2-), and probably associated with cation vacancies in the MgO matrix. The formation of O2(2-) in arc-fusion grown MgO crystals is very unexpected because of the highly reducing growth conditions. Their presence implies an internal redox reaction involving dissolved 'water' by which OH(-) pairs convert to O2(2-) plus H2 molecules. This redox conversion is supported by mass spectroscopic measurements of the H2 release from highly OH(-)-doped, finely divided MgO and by wet-chemical analysis of its oxidant concentration.

Freund, Friedemann↗

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗

HST eclipse mapping of dwarf nova OY Carinae in quiescence: An 'Fe II curtain' with Mach approx. = 6 velocity dispersion veils the white dwarf

Hubble Space Telescope (HST) observations of the eclipsing dwarf nova OY Car in its quiescent state are used to isolate the ultraviolet spectrum (1150-2500 A at 9.2 A Full Width at Half Maximum (FWHM) resolution) of the white dwarf, the accretion disk, and the bright spot. The white dwarf spectrum has a Stark-broadened photospheric L(alpha) absorption, but is veiled by a forest of blended Fe II features that we attribute to absorption by intervening disk material. A fit gives T(sub w) approx. = 16.5 x 10(exp 3) K for the white dwarf with a solar-abundance, log g = 8 model atmosphere, and T approx. = 10(exp 4) K, n(sub e) approx. = 10(exp 13)/cu cm, N(sub H) approx. = 10(exp 22) sq cm, and velocity dispersion delta V approx. = 60 km/s for the veil of homogeneous solar-abundance local thermodynamic equilibrium (LTE) gas. The veil parameters probably measure characteristic physical conditions in the quiescent accretion disk or its chromosphere. The large velocity dispersion is essential for a good fit; it lowers (chi square)/778 from 22 to 4. Keplerian shear can produce the velocity dispersion if the veiling gas is located at R approx. = 5 R(sub W) with (delta R)/R approx. = 0.3, but this model leaves an unobscured view to the upper hemisphere of the white dwarf, incompatible with absorptions that are up to 80% deep. The veiling gas may be in the upper atmosphere of the disk near its outer rim, but we then require supersonic (Mach approx. = 6) but sub-Keplerian (delta V/V(sub Kep) approx. = 0.07) velocity disturbances in this region to produce both the observed radial velocity dispersion and vertical motions sufficient to elevate the gas to z/R = cos i = 0.12. Such motions might be driven by the gas stream, since it may take several Kepler periods to reestablish the disk's vertical hydrostatic equilibrium. The temperature and column density of the gas we see as Fe II absorption in the ultraviolet are similar to what is required to produce the strong Balmer jump and line emissions seen in optical spectra of OY Car and similar quiescent dwarf novae. The outer accretion disk is detected at mid-eclipse with a spectrum that rises from 0.05 to 0.3 mJy between 2000 and 2500 A, consistent with combinations of cool blackbodies, blended Fe II emission lines, and Balmer continuum emission. The total disk flux density is 0.5 mJy at 2500 A, and this shallow disk eclipse implies a roughly flat surface brightness distribution. The bright spot, somewhat bluer than the disk, has a flux density rising from 0.05 to 0.15 mJy between 1600 and 2500 A. The C IV emission line has a broad shallow eclipse, but the radial velocity variations observed during the eclipse do not clearly distinguish between a disk or wind origin. The only possible indications of boundary layer emission are fast UV flares that appear to arise from near the central object -- not from the bright spot.

Horne, Keith↗

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗

A Report on Stochastic Fairness Queueing (SFQ) Experiments

SRI International (SRI) has developed an improved queueing algorithm, known as Stochastic Fairness Queueing (SFQ), for best-effort traffic (i.e., traffic that does not require any guaranteed service). SFQ is a probablistic variant of strict fair queueing where instead of a single queue being allocated per flow, a fixed number of queues are used and a hash function maps the IP source and destination to a particular queue. A seed to the hash function is also perturbed occasionally to help distribute the flows amongst different queues when more than one flow maps to the same queue during the lifetime of the flow. SFQ provides 'fair' access by trying to ensure that each flow from source to destination host obtains equal access to the available bandwidth. This report covers a series of experiments performed on DARTnet evaluating the behavior and performance of SFQ against a FIFO queueing discipline. These experiments were designed to show SFQ's advantages and performance, and include tests demonstrating: Fair utilization of available resources; Starvation prevention; Graceful degradation under overload conditions; and Resource usage. In general, the experiments do show that SFQ is better than FIFO queueing at allocating bandwidth equally among a set of flows. SFQ also prevents a stream from dominating the available bandwidth, which seems to be a tendency with FIFO queueing (i.e., if a flow demands more than its share of the available bandwidth, with FIFO queueing that stream receives a disproportionate amount when compared to flows demanding less than their share). Furthermore, SFQ seems to reward 'nice' users of the network by providing a lower variance in delay and more throughput when their resource demand is less than their available share. Both SFQ and FIFO queueing seem to degrade fairly well as the network becomes saturated and to recover well as the network becomes less congested. Not unexpectedly, FIFO queueing is a little more efficient than SFQ-the delays are less and the throughput slightly higher because SFQ requires more processing. However, the performance difference between the two queueing disciplines is relatively small. However, the experiments do point out some interesting behavior. FIFO queueing can behave better than SFQ with seed perturbation. We recommend further evaluation of the hash function and the seed perturbation technique. There are probably weaknesses in their current selection that cause this unexpected behavior. SFQ also seems to possess good scaling properties. To verify this, more experiments with a larger number of streams from more hosts need to be executed and examined, including the staggered introduction of streams. Staggering the streams may prove important, because graphs in the degradation experiment revealed some unexpected increases and decreases in throughput, which should be examined. This may again be due to the interaction of the hash function with the seed perturbation but it may also be related to some other unknown problem.

Denny, Barbara A.↗

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗

Convective instabilities of electromagnetic ion cyclotron waves in the outer magnetosphere

The path-integrated linear growth of electromagnetic ion cyclotron waves in the outer (L is greater than or equal to 7) magnetosphere is investigated using a realistic thermal plasma distribution with an additional anisotropic energetic ring current H(+) to provide free energy for instability. The results provide a realistic simulation of the recent Active Magneto- spheric Particle Tracer Explorers (AMPTE) observations. For conditions typical of the dayside magnetosphere, high plasma beta effects reduce the group velocity and significantly increase the spatial growth rates for left-handed polarized instabilities just below the helium gyrofrequency Omega(sub He(+)), and on the guided mode above Omega(sub He(+)) but below the cross over frequency omega(sub cr). Relatively high densities, typical of the afternoon local time sector, favor these low group velocity effects for predominantly field-aligned waves. Lower densities, typical of those found in the early morning local time sector, increase the group velocity but allow strong convective instabilities at high normalized frequencies well above Omega(sub He(+)). These waves are reflected in the magnetosphere and can exist for several equatorial transits without significant damping. They are left-handed polarized only on the first equatorial crossing and become linearly polarized for the remainder of the ray path. Consequently, these waves should be observed with basically linear polarization at all frequencies and all latitudes in the early morning local time sector. Wave growth below Omega(sub He(+)) is severely limited owing to the narrow bandwidth for instability and the small resonant path lengths. In the afternoon sector, where plasma densities can exceed 10(exp 7)/cu m, intense convective amplification is possible both above and below Omega(sub He(+)). Waves below Omega(sub He(+)) are not subject to reflection when the O(+) concentration is small and therefore should be observed with left-handed polarization near the equator and essentially linear polarization at higher latitudes. Since the He(+) concentration is usually large in the afternoon sector, guided mode waves above Omega(sub He(+)) reflect to form a background distribution with basically linear polarization. We suggest that the strong left-handed polarized emissions observed by AMPTE in the afternoon sector near the equator are probably due to strongly growing low group velocity waves at frequencies just below Omega(sub He(+)), and on the guided mode above Omega(sub He(+)).

Horne, Richard B.↗

Interactive Soil Dust Aerosol Model in the GISS GCM: Sensitivity of the Soil Dust Cycle to Radiative Properties of Soil Dust Aerosols - Part 1

The sensitivity of the soil dust aerosol cycle to the radiative forcing by soil dust aerosols is studied. Four experiments with the NASA/GISS atmospheric general circulation model, which includes a soil dust aerosol model, are compared, all using a prescribed climatological sea surface temperature as lower boundary condition. In one experiment, dust is included as dynamic tracer only (without interacting with radiation), whereas dust interacts with radiation in the other simulations. Although the single scattering albedo of dust particles is prescribed to be globally uniform in the experiments with radiatively active dust, a different single scattering albedo is used in those experiments to estimate whether regional variations in dust optical properties, corresponding to variations in mineralogical composition among different source regions, are important for the soil dust cycle and the climate state. On a global scale, the radiative forcing by dust generally causes a reduction in the atmospheric dust load corresponding to a decreased dust source flux. That is, there is a negative feedback in the climate system due to the radiative effect of dust. The dust source flux and its changes were analyzed in more detail for the main dust source regions. This analysis shows that the reduction varies both with the season and with the single scattering albedo of the dust particles. By examining the correlation with the surface wind, it was found that the dust emission from the Saharan/Sahelian source region and from the Arabian peninsula, along with the sensitivity of the emission to the single scattering albedo of dust particles, are related to large scale circulation patterns, in particular to the trade winds during Northern Hemisphere winter and to the Indian monsoon circulation during summer. In the other regions, such relations to the large scale circulation were not found. There, the dependence of dust deflation to radiative forcing by dust particles is probably dominated by physical processes with short time scales. The experiments show that dust radiative forcing can lead to significant changes both in the soil dust cycle and in the climate state. To estimate dust concentration and radiative forcing by dust more accurately, dust size distributions and dust single scattering albedo in the model should be a function of the source region, because dust concentration and climate response to dust radiative forcing are sensitive to dust radiative parameters.

Perlwitz, Jan↗