Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical homogeneity”

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 127 records · Page 7

The effect of repeating gamma-ray bursts on V/V(sub max)

I consider the effect of repeating burst sources on the the mean value of V/V(sub max) statistic. I find that if the burst sources are distributed homogeneously in a d-dimensional space, the mean value of V/V(sub max) converges to d/(d+3) as long as the luminosity function is dependent of position. Choosing the brightest event from a cluster of k events shifts the luminosity function to higher luminosities, but if the original luminosity function is position independent, it remains so. Therefore the treatment of repeating events, if applied consistently, will not affect the effectiveness the mean value of V/V(sub max) as a test of burst homogeneity. The calculation of the mean value of V/V(sub max) for apparent repeating and non-repeating source populations will be raised by the incorrect classification of faint bursts. In conclusion, the current practice of calculating the mean value of V/V(sub max) using all bursts, even apparent repeaters, and treating multispike bursts as single bursts, is valid and consistent.

Band, David L.↗

Estimations of ABL fluxes and other turbulence parameters from Doppler lidar data

Techniques for extraction boundary layer parameters from measurements of a short-pulse CO2 Doppler lidar are described. The measurements are those collected during the First International Satellites Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). By continuously operating the lidar for about an hour, stable statistics of the radial velocities can be extracted. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation. Combining many measurements would normally reduce the error provided that, it is unbiased and uncorrelated. The nature of some of the algorithms however, is such that, biased and correlated errors may be generated even though the raw measurements are not. Data processing procedures were developed that eliminate bias and minimize error correlation. Once bias and error correlations are accounted for, the large sample size is shown to reduce the errors substantially. The principal features of the derived turbulence statistics for two case studied are presented.

Gal-Chen, Tzvi↗

The edge-on Galaxies in the Pan-STARRS survey (EGIPS)

ABSTRACT We present a catalogue of 16 551 edge-on galaxies created using the public DR2 data of the Pan-STARRS survey. The catalogue covers the three quarters of the sky above Dec. = −30°. The galaxies were selected using a convolutional neural network, trained on a sample of edge-on galaxies identified earlier in the SDSS survey. This approach allows us to dramatically improve the quality of the candidate selection and perform a thorough visual inspection in a reasonable amount of time. The catalogue provides homogeneous information on astrometry, SExtractor photometry, and non-parametric morphological statistics of the galaxies. The photometry is reliably for objects in the 13.8–17.4 r-band magnitude range. According to the HyperLeda data base, redshifts are known for about 63 per cent of the galaxies in the catalogue. Our sample is well separated into the red sequence and blue cloud galaxy populations. The edge-on galaxies of the red sequence are systematically Δ(g − i) ≈ 0.1 mag redder than galaxies oriented at an arbitrary angle to the observer. We found a variation of the galaxy thickness with the galaxy colour. The red sequence galaxies are thicker than the galaxies of the blue cloud. In the blue cloud, on average, thinner galaxies turn out to be bluer. In the future, based on this catalogue it is intended to explore the three-dimensional structure of galaxies of different morphologies, as well as to study the scaling relations for discs and bulges.

Makarov, D. (ORCID:0000000191103221)↗

Trace Element Analyses of Micron-Size Particles and Statistical Determination of Minimum Detection Limits

Savannah River National Laboratory (SRNL) has developed expertise in producing homogeneous, ca. 1 m-diameter spherical particles of mixed-element components, wherein dopants can be varied from a trace constituent (ppm) to wt.% concentrations. The samples used for this work are nickel-doped cerium oxide microspheres produced by SRNL. They were initially selected as analogs for plutonium-doped uranium oxide particles and analyzed as part of a larger study to evaluate whether electron probe microanalyzers (EPMA) can be used to characterize nuclear materials as an alternative or complementary method to mass spectrometers. The five samples used in this study contained nominal compositions of 0, 0.004, 0.04, 0.4 and 4 wt.% Ni. They were analyzed by both an Agilent 7900 Q-ICP-MS at SRNL and the JEOL JXA8530F Plus EPMA at the University of Minnesota. In addition to EPMA results (calibrated with high-precision Q-ICP-MS analyses) suggesting that the EPMA could address outstanding nuclear material characterization needs, these samples 1) showcase the ability of the EPMA to quantify not just trace concentrations, but trace concentrations in microparticles (1 m diameter, Fig. 1), and 2) offer a unique opportunity to evaluate the methodology for assessing the minimum detection limits of EPMA analyses.

McSwiggen, Peter [JEOL USA, 11 Dearborn Road, Peab↗

A Combined Al-Mg/Pb-Pb Age of the Solar System

Astrophysical models of planet formation and protoplanetary disk evolution demand precise and accurate timing of the sequence of events in the solar nebula, relative to a time t=0, usually taken to be during the short epoch of CAI (Ca-rich, Al-rich inclusion) formation. Most CAIs formed withlive26Al (mean-life τ26= 1.034 Myr [1]), with an abundance 26Al/27Al ≈ (26Al/27Al)SS= 5.23 × 10-5[2]. We adopt this as the widespread level of 26Al in the solar nebula at t=0. Assuming spatial homogeneity of 26Al, an inclusion that had less 26Al, (26Al/27Al)0, formed a time Δt26= τ26ln[(26Al/27Al)SS/ (26Al/27Al)0] after t=0.These ages are typical precise to within ±0.1 Myr. Igneous bulk meteorites and inclusions can be relatively dated by the Al-Mg chronometer, but only ifΔt26<6 Myr. The Pb-Pb system is useful as a longer relative chronometer. It yields absolute ages tPb using 207Pb/206Pb, 206Pb/204Pb, and 238U/235U ratios measured indifferent portions of a sample, assuming certain half-lives [4]. These absolute ages are uncertain to within ±9 Myr due to uncertainties in the 235U half-life[3], but times of formation ΔtPb= tCAI–tPb relative to t=0, are more precise(±0.5Myr),iftCAI can be found. Here, tCAI means the Pb-Pb age that would be measured in CAIs using the half-lives the community typically uses, if they achieved isotopic closure at t=0. Unfortunately, direct Pb-Pb dating of CAIs has not definitively determined tCAI. Based on four CAIs with canonical (26Al/27Al)0,[5,6] found tPb= 4567.30 ± 0.16 Myr. No other CAI ages with measured 238U/235U have been reported in the refereed literature, but there are hints of other CAIs with ages tPb= 4568.0 ± 0.2 Myr [7] and tPb= 4568.3 ± 0.2 Myr [8].It is unclear whether anyof these igneous type B CAIs isotopically closed at t=0 or represents tCAI. Instead of measurements, we advocate finding tCAI by minimizing the discrepancies between the Al-Mg and Pb-Pb chronometers. Assuming Δt26=ΔtPb, we find the implied t’CAI= tPb+Δt26, then define t*CAIas the weighted mean of the t’CAI. t*CAIis the best guess for the Pb-Pb age of t=0; the assumption of homogeneity is justified if the t’CAI cluster within errors around t*CAI. This statistical approach is similar to, but improves on, that of[9]. We find t*CAI= 4568.73 ± 0.16 Myr. Below we discuss our methodology and the implications of this age for CAIs, 1.4 Myr older than the reported and typically used age 4567.30±0.16 Myr. Methods: We base our estimate of t*CAIon five achondrites for which published (26Al/27Al)0and Pb-Pb ages exist: the quenched angrites D’Orbigny, Sahara 99555 (SAH 99555), and Northwest Africa (NWA) 1670; the pseudo-eucrite Asuka 881394; and the inner disk achondrite. All are “NC” (non-carbonaceous) achondrites that likely cooled quickly enough that the Al-Mg and Pb-Pb systems achieved isotopic closure simultaneously. We also considered the “CC” (carbonaceous chondrite-like) achondrites NWA 2796 and NWA 6704, butdo not include them in our fit. Al-Mg and Pb-Pb seem not to have closed simultaneously, possibly because formation in the outer disk from volatile-rich composition led to slower cooling. Of the 8 chondrules from NWA 5697 measured by [20,21], we also consider the 4 for which 238U/235U was measured: 2-C1, 5-C2, 3-C5, 11-C1.Depending on their post-formation thermal histories, the Al-Mg and Pb-Pb systems in chondrules may or may not have closed simultaneously. Table 1: (26Al/27Al)0, Pb-Pb ages of selected samplesSample(26Al/27Al)0/ 10-6RefPb-PbRefD’Orbigny3.98±0.15104563.43±0.19♮10-12SAH 995553.64±0.18104563.88±0.2712NWA16705.92±0.59104564.39±0.24*10Asuka 88139413.1±0.5613-154564.98±0.1715NWA 73253.03±0.14164563.4±2.616NWA 27963.94±0.16174562.89±0.5917NWA 67043.03±0.14184562.76±0.26192-C17.56±1.53204567.57±0.56*215-C27.04±1.51204567.54±0.52*213-C58.85±1.83204566.20±0.63*2111-C15.55±1.84204565.84±0.72*21*regression based on one subset of data points ♮weighted mean of two datasets Pb-Pb ages are proportional to the intercept of the line formed by linear regression of 207Pb/206Pb vs. 206Pb/204Pb data from various washes, leachates and residues of acid dissolution of a sample. Because contamination by terrestrial or primordial Pb is pervasive, some fractions must be excluded from regressions to ensure a fit with acceptable mean squares weighted deviation (MSWD). Usually points are excluded based on low [Pb], or low 206Pb/204Pb ratio(low radiogenic component), with single outliers identified [11,12,15,16,17]. In the starred examples (Table 1)and the case of 3 CAI Pb-Pb ages [5], up to half the points were excluded solely because did not fit a pre-determined line. This approach is vulnerable to confirmation bias and produces fits with low MSWD and too-low Pb-Pb age uncertainty. Regressing the same data points as [10], were produce the Pb-Pb age of NWA 1670 of 4564.39±0.24 Myr. But selecting other combinations of data points, other, equally valid, isochrons yield ages from 4563.77±0.21 Myr to 4564.64±0.23 Myr. Similar arguments apply to the Pb-Pb isochrons built by [21] for chondrules 2-C1 (we find 4567.33±0.44 to 4567.85±0.46 Myr), 5-C2 (4566.84±0.53 to 4567.70±0.44 Myr), 3-C5 (4565.84±0.54to 4567.04±0.54)and 11-C1 (4565.36±0.51 to 4565.74±0.45 Myr). Our adopted ages for these and NWA 1670 are listed in Table 2.Table 2. tCAI estimated from various components, using our regressions for the chondrules & NWA 1670.SampleΔt26(Myr)tPb(Myr)t’CAI(Myr)D’Orbigny5.05±0.044563.43±0.194568.48±0.19SAH 995555.14±0.054563.88±0.274569.02±0.27NWA16704.64±0.104564.21±0.634568.85±0.67Asuka 8813943.81±0.044564.98±0.174568.79±0.17NWA 73255.33±0.054563.4±2.64568.7±2.6NWA 27965.06±0.044562.89±0.594567.95±0.59NWA 67045.29±0.134562.76±0.264568.05±0.292-C12.00±0.214567.59±0.704569.59±0.725-C22.07±0.224567.23±0.914569.30±0.933-C51.84±0.214566.44±1.124568.28±1.1411-C12.32±0.344565.52±0.664567.84±0.73achondrite4568.72±0.16chondrules4568.76±0.58combined4568.73±0.16A weighted average of the five NC achondrites(or just D’Orbigny, SAH 99555 and Asuka 881394)yields t*CAI= 4568.72 ± 0.16Myr. All are consistent with this value to within 1.8σ, and MSWD=1.5. Including the 4 U-corrected chondrules, t*CAI= 4568.73± 0.16Myrwith MSWD=1.66, which is statistically significant. All chondrules and NC achondrites are consistent with this to within 1.8σ, (Figure 1).Figure 1. Al-Mg formation times after t=0 vs. Pb-Pb ages. The five NC achondrites and four chondrules are consistent with a Pb-Pb age of t=0 of 4568.7 Myr. Discussion: The data from achondrites and chondrules are consistent with a single Pb-Pb age at t=0, justifying the assumption of 26Al homogeneity. The age, 4568.7 Myr, is ≈1.4 Myr older than the commonly accepted Pb-Pb age of CAIs that formed with canonical 26Al/27Al at t=0 [3]. Others have interpreted the discrepancy to signify 26Al heterogeneity in the CAI-forming region[5,21]. We suggest instead that CAIs were exposed to transient heating events that reset the Pb-Pb system without disturbing the Al-Mg system. Notably, chondrules typically experienced transient heating at these times in the nebula [22]. If so, direct measurements of CAIs will not yield as reliable a Pb-Pb age of t=0 as statistical approaches like this and that of [9].References:[1] Auer et al. 2009. [2] Jacobsen, B et al. 2008, EPSL 272, 353-364. [3] Tissot, Fet al. 2017, GCA 213, 593-617.[4] Villa, I et al. 2016, GCA 172, 387-392.[5] Amelin, Y et al. 2010, EPSL 300, 343-350.[6] Connelly, Jet al. 2012, Science 338, 651.[7] Bouvier, Aet al. 2011, LPICo 1639, 9054. [8] Bouvier, Aand Wadhwa, M2010, Nat Geosci 3, 637-641.[9] Nyquist, Let al. 2009, GCA 73, 5115-5136.[10] Schiller et al. 2015. [11] Wadhwa & Brennecka 2012. [12] Tissot et al. 2017. [13] Nyquist et al. 2003. [14] Wadhwa et al. 2009. [15] Wimpenny et al. 2019, GCA 244, 478-501.[16] Koefoed et al. 2016, GCA 183, 31-45. [17] Bouvier, A et al. 2011, GCA 75, 5310-5323. [18] Sanborn, Met al. 2019, GCA 245, 577-596. [19] Amelin, Yet al. 2019, GCA 245, 628-642.[20] Bollard, Jet al. 2017, Sci Adv 3 ,e1700407. [21] Bollard, Jet al. 2019, GCA 260, 62-83.[22] Villeneuve, J et al. 2009, Science 325, 985

S. J. Desch↗

Modeling of the reactant conversion rate in a turbulent shear flow

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Frankel, S. H.↗

Ectomycorrhizal Community Structure and Soil Characteristics of Mature Lodgepole Pine (Pinus Contorta) and Adjacent Stands of Old Growth Mixed Conifer in Yellowstone National Park, Wyoming USA

Forest development patterns following disturbance are known to influence the physical and chemical attributes of soils at different points in time. Changes in soil resources are thought to have a corresponding effect on ectomycorrhizal (ECM) community structure. We used molecular methods to compare below-ground ECM species richness, composition, and abundance between adjacent stands of homogenous lodgepole pine and old growth mixed conifer in Yellowstone National Park (YNP). In each stand-type we collected soil cores to both identify mycorrhizae and assess soil chemistry. Although no statistical difference was observed in the mean number of ECM root tips per core between stand types, the total number of species identified (85 versus 35) and the mean number of species per core (8.8 +/- 0.6 versus 2.5 +/- 0.3) were significantly higher in lodgepole pine. Differences between the actual and estimated species richness levels indicated that these forest types support a high number of ECM species and that undersampling was severe. Species compositions were widely disparate between stands where only four species were shared out of a total of 116. Soil analysis also revealed that mixed conifer was significantly lower in pH, but higher in organic matter, potassium, phosphorus, and ammonium when compared to lodgepole pine stands. Species richness per core was correlated with these chemical data, however, analysis of covariance indicated that stand type was the only statistically significant factor in the observed difference in species richness. Our data suggest that ECM fungal richness increases as homogenous lodgepole pine stands grow and mature, but declines after Engelmann spruce and subalpine fir colonize. Despite difficulties linking species composition with soil chemistry, there are a variety of physical and chemical factors that could be influencing ECM community structure. Future field experiments are necessary to test some of the mechanisms potentially operating within this system.

Douglas, Robert B.↗

Labile Fe(III) supersaturation controls nucleation and properties of product phases from Fe(II)-catalyzed ferrihydrite transformation

We report Fe(II)-catalyzed ferrihydrite (Fh) transformation to more crystalline iron (oxyhydr)oxide phases is a widely occurring geochemical process which has been extensively studied as a function of Fe(II)/Fh ratios at fixed Fh loadings. However, recent isolation of an intermediate Fe(III) species resulting from Fe(II)-Fh contact that facilitates transformation by dissolution/reprecipitation suggests that the kinetics and properties of product phases will instead depend mostly on its rate of accumulation to a critical concentration, consistent with principles in the classical nucleation theory (CNT). This suggests a dependence both on the loading of Fe(II) on the surface, which controls the rate of labile Fe(III) formation, as well as the available volume of solution, which also impacts how fast it can achieve its critical concentration to nucleate product phases. To specifically examine the latter effect, here we studied transformation of 15 mg Fh in 1 mM FeSO 4 solutions at pH 7.2 in batch suspensions of 30 mL, 150 mL, and 450 mL volumes. Time-dependent concentrations of aqueous Fe(II), surface-associated Fe(II), and resulting labile Fe(III) were monitored along with bulk solids characterization as a function of time. Transmission electron microscopy (TEM) was used to visualize the evolution of phases at identical locations on TEM grids. The collective results show that the rates of Fh loss and emergence of product lepidocrocite (Lp) and goethite (Gt) as well as their phase proportions, nucleation mode and morphological properties depend directly on the rate of accumulation of the labile Fe(III) precursor to its critical concentration, which in our experiments was controlled simply by varying the available volume of solution into which it enters. Statistical analyses of TEM image data suggest that while both heterogeneous and homogeneous nucleation occurred in all experiments, the former was increasingly favored at lower Fh/solution ratio due to its lower nucleation barrier being more favorable at attendant lower supersaturations of Fe(III).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Digital preprocessing of SEASAT imagery

A model for radar image data is derived for use in developing enhancement techniques. This model describes the image data as the result of a multiplicative-convolved noise process. This information is then used to design a minimum mean square error (MMSE) filter. The resulting filter is implemented adaptively to change with local statistics. A radar image processing technique which provides the MMSE estimate inside homogeneous areas and tends to preserve edge structure is thus developed. This technique has been implemented and tested using digitally correlated SEASAT-A synthetic aperture radar (SAR) imagery.

Frost, V. S.↗

Evaluation of the radar response to land surfaces and volumes: Examination of theoretical models, target statistics, and applications

Four areas of L-band radar remote sensing of terrain were examined: (1) the behavior of the radar backscatter coefficient of distributed surface and volumes as a function of the targets' dielectric and geometric parameters and as a fnction of their physical parameters; (2) the correspondence of the angular behavior of the relative backscatter coefficient as extracted from SIR-B digital imagery and truck mounted L-band scatterometer measurements for about 100 fields; (3) the statistical behavior of SIR-B image density for targets that appear homogeneous on Thematic Mapper (TM) optical imagery and/or color IR photography; and (4) the applicability of SIR-B imagery both alone and in conjunction with TM imagery for the classification and monitoring of land cover and renewable resources.

Ulaby, F. T.↗

Climate Impacts of Ice Nucleation

Several different ice nucleation parameterizations in two different General Circulation Models (GCMs) are used to understand the effects of ice nucleation on the mean climate state, and the Aerosol Indirect Effects (AIE) of cirrus clouds on climate. Simulations have a range of ice microphysical states that are consistent with the spread of observations, but many simulations have higher present-day ice crystal number concentrations than in-situ observations. These different states result from different parameterizations of ice cloud nucleation processes, and feature different balances of homogeneous and heterogeneous nucleation. Black carbon aerosols have a small (0.06 Wm(exp-2) and not statistically significant AIE when included as ice nuclei, for nucleation efficiencies within the range of laboratory measurements. Indirect effects of anthropogenic aerosols on cirrus clouds occur as a consequence of increasing anthropogenic sulfur emissions with different mechanisms important in different models. In one model this is due to increases in homogeneous nucleation fraction, and in the other due to increases in heterogeneous nucleation with coated dust. The magnitude of the effect is the same however. The resulting ice AIE does not seem strongly dependent on the balance between homogeneous and heterogeneous ice nucleation. Regional effects can reach several Wm2. Indirect effects are slightly larger for those states with less homogeneous nucleation and lower ice number concentration in the base state. The total ice AIE is estimated at 0.27 +/- 0.10 Wm(exp-2) (1 sigma uncertainty). This represents a 20% offset of the simulated total shortwave AIE for ice and liquid clouds of 1.6 Wm(sup-2).

Gettelman, Andrew↗

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

BATSE software for the analysis of the gamma ray burst spatial distribution

The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.

Hakkila, Jon↗

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗

Mirroring in the Fokker-Planck coefficient for cosmic-ray pitch-angle scattering in homogeneous magnetic turbulence

The Fokker-Planck coefficient for pitch-angle scattering, appropriate for cosmic rays in homogeneous stationary magnetic turbulence is computed without making any specific assumptions concerning the statistical symmetries of the random field. The Fokker-Planck coefficient obtained can be used to compute the parallel diffusion coefficient for high-energy cosmic rays propagating in the presence of strong turbulence, or for low-energy cosmic rays in the presence of weak turbulence. Because of the generality of magnetic turbulence allowed for in the analysis, special interplanetary magnetic field features, such as discontinuities or particular wave modes, can be included rigorously.

Goldstein, M. L.↗

Comparison of wheat classification accuracy using different classifiers of the image-100 system

Classification results using single-cell and multi-cell signature acquisition options, a point-by-point Gaussian maximum-likelihood classifier, and K-means clustering of the Image-100 system are presented. Conclusions reached are that: a better indication of correct classification can be provided by using a test area which contains various cover types of the study area; classification accuracy should be evaluated considering both the percentages of correct classification and error of commission; supervised classification approaches are better than K-means clustering; Gaussian distribution maximum likelihood classifier is better than Single-cell and Multi-cell Signature Acquisition Options of the Image-100 system; and in order to obtain a high classification accuracy in a large and heterogeneous crop area, using Gaussian maximum-likelihood classifier, homogeneous spectral subclasses of the study crop should be created to derive training statistics.

Dejesusparada, N.↗

An interactive system for analysis of global cloud imagery

Synoptic images of the global cloud pattern composited from six contemporaneous satellites provide an unprecedented view of the global cloud field. Having horizontal resolution of about 0.5 deg and temporal resolution of 3 h, the global cloud imagery (GCI) resolves most of the variability of organized convection, including several harmonics of the diurnal cycle. Although the GCI has these attractive features, the dense and 3D nature of that data make it a formidable volume of information to treat in a practical and efficient manner. An interactive image-analysis system (IAS) has been developed to investigate the space-time variability of global cloud behavior. In the IAS, data, hardware, and software are integrated into a single system providing a variety of space-time covariance analyses in a menu-driven format. Owing to its customized architecture and certain homogeneous properties of the GCI, the IAS calculates such quantities effectively. Many covariance statistics are derived from 3D data with interactive speed, allowing the user to interrogate the archive iteratively in a single session. The 3D nature of those analyses and the speed with which they are performed distinguish the IAS from conventional image processing of 2D data.

Woodberry, Karen↗

Discussion on selected symbiotic stars

Because of its large variety of aspects, the symbiotic phenomenon is not very suitable for a statistical treatment. It is also not clear whether symbiotic stars really represent a homogeneous group of astrophysical objects or a collection of objects of different natures but showing similar phenomena. However we are especially interested in the symbiotic phenomenon, i.e., in those physical processes occurring in the atmosphere of each individual object and in their time dependence. Such a research can be performed through the detailed analysis of individual objects. This study should be done for a time long enough to cover all the different phases of their activity, in all the spectral ranges. Since the typical time scale of the symbiotic phenomena is up to several years and decades, this represents a problem since, for instance, making astronomy outside the visual region is a quite new field of research. It was a fortunate case that a few symbiotic stars (Z And, AG Dra, CH Cyg, AX Per, and PU Vul) had undergone remarkable light variations (or 'outbursts') in recent years, which could have been followed in the space ultraviolet with IUE, and simultaneously in the optical and IR with ground-based telescopes. But, in general, the time coverage of most of the symbiotic objects is too short to have a complete picture of their behavior. In this regard, one should recall Mayall's remark about the light curve of Z And: 'Z Andromedae is another variable that shows it will require several hundred years of observations before a good analysis can be made of its variations'. This pessimistic remark should be considered as a note of caution for those involved in the interpretation of the observations. We shall discuss a number of individual symbiotic stars for which the amount of observational data is large enough to draw a rather complete picture of their general behavior and to make consistent models. We shall especially illustrate the necessary steps toward an empirical model and take the discussion of the individual objects as a useful occasion to describe different techniques of diagnosis.

Viotti, Roberto↗