Applications of wavelet-based compression to multidimensional Earth science data

A data compression algorithm involving vector quantization (VQ) and the discrete wavelet transform (DWT) is applied to two different types of multidimensional digital earth-science data. The algorithms (WVQ) is optimized for each particular application through an optimization procedure that assigns VQ parameters to the wavelet transform subbands subject to constraints on compression ratio and encoding complexity. Preliminary results of compressing global ocean model data generated on a Thinking Machines CM-200 supercomputer are presented. The WVQ scheme is used in both a predictive and nonpredictive mode. Parameters generated by the optimization algorithm are reported, as are signal-to-noise (SNR) measurements of actual quantized data. The problem of extrapolating hydrodynamic variables across the continental landmasses in order to compute the DWT on a rectangular grid is discussed. Results are also presented for compressing Landsat TM 7-band data using the WVQ scheme. The formulation of the optimization problem is presented along with SNR measurements of actual quantized data. Postprocessing applications are considered in which the seven spectral bands are clustered into 256 clusters using a k-means algorithm and analyzed using the Los Alamos multispectral data analysis program, SPECTRUM, both before and after being compressed using the WVQ program

Simple and accurate modelling of the gravitational potential produced by thick and thin exponential discs

We present accurate models of the gravitational potential produced by a radially exponential disc mass distribution. The models are produced by combining three separate Miyamoto–Nagai discs. Such models have been used previously to model the disc of the Milky Way, but here we extend this framework to allow its application to discs of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disc treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disc by <0.4 per cent out to 4 disc scalelengths, and <1.9 per cent out to 10 disc scalelengths. We tabulate fitting parameters which facilitate construction of exponential discs for any scalelength, and a wide range of disc thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disc galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disc, and a discy dwarf galaxy

Can Light Echoes Account for the Slow Decay of Type IIn Supernovae?

The spectra of type IIn supernovae indicate the presence of apre-existing slow, dense circumstellar wind (CSW). If the CSW extends sufficiently far from the progenitor star, then dust formation should occur in the wind. The light from the supernova explosion will scatter off this dust and produce a light echo. Continuum emission seen after the peak will have contributions from both this echo as well as from the shock of the ejecta colliding with the CSW, with a fundamental question of which source dominates the continuum. We calculate the brightness of the light echo as a function of time for a range of dust shell geometries, and use our calculations to fit to the light curves of SN 1988Z and SN 1997ab, the two slowest declining IIn supernovae on record. We find that the light curves of both objects can be reproduced by the echo model. However, their rate of decay from peak, color at peak and their observed peak absolute magnitudes when considered together are inconsistent with the echo model. Furthermore, when the observed values of M$_{B}$ are corrected for the effects of dust scattering, the values obtained imply that these supernovae have unrealistically high luminosities. We conclude that light echoes cannot properly account for the slow decline seen in some IIn's, and that the shock interaction is likely to dominate the continuum emission.Comment: 15 pages, 9 figure

Congruent families and invariant tensors

Classical results of Chentsov and Campbell state that -- up to constant multiples -- the only $2$-tensor field of a statistical model which is invariant under congruent Markov morphisms is the Fisher metric and the only invariant $3$-tensor field is the Amari-Chentsov tensor. We generalize this result for arbitrary degree $n$, showing that any family of $n$-tensors which is invariant under congruent Markov morphisms is algebraically generated by the canonical tensor fields defined in an earlier paper

s-Process Nucleosynthesis in Advanced Burning Phases of Massive Stars

We present a detailed study of s-process nucleosynthesis in massive stars of solar-like initial composition and masses 15, 20,25, and 30 Msun. We update our previous results of s-process nucleosynthesis during the core He-burning of these stars and then focus on an analysis of the s-process under the physical conditions encountered during the shell-carbon burning. We show that the recent compilation of the Ne22(alpha,n)Mg25 rate leads to a remarkable reduction of the efficiency of the s-process during core He-burning. In particular, this rate leads to the lowest overproduction factor of Kr80 found to date during core He-burning in massive stars. The s-process yields resulting from shell carbon burning turn out to be very sensitive to the structural evolution of the carbon shell. This structure is influenced by the mass fraction of C12 attained at the end of core helium burning, which in turn is mainly determined by the C12(alpha,gamma)O16 reaction. The still present uncertainty in the rate for this reaction implies that the s-process in massive stars is also subject to this uncertainty. We identify some isotopes like Zn70 and Rb87 as the signatures of the s-process during shell carbon burning in massive stars. In determining the relative contribution of our s-only stellar yields to the solar abundances, we find it is important to take into account the neutron exposure of shell carbon burning. When we analyze our yields with a Salpeter Initial Mass Function, we find that massive stars contribute at least 40% to s-only nuclei with mass A 90, massive stars contribute on average ~7%, except for Gd152, Os187, and Hg198 which are ~14%, \~13%, and ~11%, respectively.Comment: 52 pages, 16 figures, accepted for publication in Ap

Integral Battery Power Limiting Circuit for Intrinsically Safe Applications

A circuit topology has been designed to guarantee the output of intrinsically safe power for the operation of electrical devices in a hazardous environment. This design uses a MOSFET (metal oxide semiconductor field-effect transistor) as a switch to connect and disconnect power to a load. A test current is provided through a separate path to the load for monitoring by a comparator against a preset threshold level. The circuit is configured so that the test current will detect a fault in the load and open the switch before the main current can respond. The main current passes through the switch and then an inductor. When a fault occurs in the load, the current through the inductor cannot change immediately, but the voltage drops immediately to safe levels. The comparator detects this drop and opens the switch before the current in the inductor has a chance to respond. This circuit protects both the current and voltage from exceeding safe levels. Typically, this type of protection is accomplished by a fuse or a circuit breaker, but in order for a fuse or a circuit breaker to blow or trip, the current must exceed the safe levels momentarily, which may be just enough time to ignite anything in a hazardous environment. To prevent this from happening, a fuse is typically current-limited by the addition of the resistor to keep the current within safe levels while the fuse reacts. The use of a resistor is acceptable for non-battery applications where the wasted energy and voltage drop across the resistor can be tolerated. The use of the switch and inductor minimizes the wasted energy. For example, a circuit runs from a 3.6-V battery that must be current-limited to 200 mA. If the circuit normally draws 10 mA, then an 18-ohm resistor would drop 180 mV during normal operation, while a typical switch (0.02 ohm) and inductor (0.97 ohm) would only drop 9.9 mV. From a power standpoint, the current-limiting resistor protection circuit wastes about 18 times more power than the switch and the inductor configuration. In the fault condition, both the resistor and the inductor react immediately. The resistor reacts by allowing more current to flow and dropping the voltage. Initially, the inductor reacts by dropping the voltage, and then by not allowing the current to change. When the comparator detects the drop in voltage, it opens the switch, thus preventing any further current flow. The inductor alone is not sufficient protection, because after the voltage drop has settled, the inductor would then allow the current to change, in this example, the current would be 3.7 A. In the fault condition, the resistor is flowing 200 mA until the fuse blows (anywhere from 1 ms to 100 s), while the switch and inductor combination is flowing about 2 A test current while monitoring for the fault to be corrected. Finally, as an additional safety feature, the circuit can be configured to hold the switch opened until both the load and source are disconnected

Making a Place for the Next Generation of Geoscientists

Early-career scientists from the AGU Mentoring Network discuss how the global pandemic has exacerbated long-standing issues with the availability of positions in and the diversity of the geosciences.</jats:p

The New Pretender: A Large UK Case Series of Retinal Injuries in Children Secondary to Handheld Lasers

Purpose: To characterize a large single-center series of retinal injuries in children secondary to handheld laser devices, with emphasis on potential prognostic factors. / Design: Retrospective case series. / Methods: Sixteen children (24 eyes) with retinal injuries secondary to handheld lasers were identified from our electronic patient record system. Case notes, digital fundus photography, and spectral-domain optical coherence tomography images were reviewed. / Results: The mean age of affected children was 12.7 years (range 9–16 years), with 12 male and 4 female subjects. Mean follow up was 5.4 months (range 1–23 months). Five children (31%) were referred as suspected retinal dystrophies. The mean logMAR visual acuity at presentation was 0.30 (20/40) (range −0.20 [20/12.5] to 1.6 [20/800]). Eleven children (69%; 15 eyes) had “mild” injuries with focal retinal disruption confined to the photoreceptor and ellipsoid layers; such injuries were associated with a better prognosis, the mean visual acuity at presentation being 0.10 (20/25). “Moderate” injuries were seen in 3 eyes of 2 children, with retinal disruption confined to the outer retinal layer but diffuse rather than focal in nature. Three patients (4 eyes) had “severe” injuries, with subfoveal outer retinal architecture loss and overlying hyperreflective material in inner retinal layers. / Conclusion: Retinal injuries secondary to handheld laser devices may be difficult to diagnose and are likely underreported. It is important that such data are in the public domain, so regulatory authorities recognize the importance of laser retinopathy as an avoidable cause of childhood visual impairment and take steps to minimize the incidence and impact of laser injuries

Singularity Free Inhomogeneous Models with Heat Flow

We present a class of singularity free exact cosmological solutions of Einstein's equations describing a perfect fluid with heat flow. It is obtained as generalization of the Senovilla class [1] corresponding to incoherent radiation field. The spacetime is cylindrically symmetric and globally regular.Comment: 6 pages, TeX, to appear in Class.Quant.Gra