Isothermal Atmosphere Research Articles

Theoretical analysis of the atmospheric entry of sub-mm meteoroids of MgxCa[formula omitted]CO3 composition

Current models allow to reliably simulate mechanical and thermal phenomena associated with a micrometeor passage through the Earth’s atmosphere. However, these models have rarely been applied to materials other than those most common in meteorites, such as silicates and metals. A particular case that deserves attention is the one of micrograins made of minerals, in particular carbonates, which have been associated, in meteorites, with organic molecules. Carbonates are known for their decomposition in vacuum at moderate temperatures, and they might contribute to the thermal protection of organic matter. In this work, a model with non isothermal atmosphere, power balance, evaporation, ablation, radiation losses and stoichiometry, is proposed. This paper includes the very first calculations for meteoroids with a mixed carbonate composition. Results show that the carbonate fraction of these objects always go to zero at high altitudes except for grazing entries, where the reached temperature is lower and some carbonate remains unreacted. For all entry conditions, peculiar temperature curves are obtained due to the decomposition process. Furthermore, a significant impact of decomposition cooling on the temperature peak is observed for some grazing entry cases. Although specific solutions used in these calculations can be improved, this work sets a definite model and a basis for future research on sub-mm grains of relatively volatile minerals entering the Earth’s atmosphere.

High-precision multiwavelength eclipse photometry of the ultra-hot gas giant exoplanet WASP-103 b

We present 16 occultation and three transit light curves for the ultra-short period hot Jupiter WASP-103b, in addition to five new radial velocity measurements. We combine these observations with archival data and perform a global analysis of the resulting extensive data set, accounting for the contamination from a nearby star. We detect the thermal emission of the planet in both the z΄ and KS bands, the measured occultation depths being 699±110 ppm (6.4σ) and |$3567_{-350}^{+400}$| ppm (10.2σ), respectively. We use these two measurements, together with recently published HST/WFC3 data, to derive joint constraints on the properties of WASP-103b's dayside atmosphere. On one hand, we find that the z΄ band and WFC3 data are best fit by an isothermal atmosphere at 2900 K or an atmosphere with a low H2O abundance. On the other hand, we find an unexpected excess in the KS band measured flux compared to these models, which requires confirmation with additional observations before any interpretation can be given. From our global data analysis, we also derive a broad-band optical transmission spectrum that shows a minimum around 700 nm and increasing values towards both shorter and longer wavelengths. This is in agreement with a previous study based on a large fraction of the archival transit light curves used in our analysis. The unusual profile of this transmission spectrum is poorly matched by theoretical spectra and is not confirmed by more recent observations at higher spectral resolution. Additional data, in both emission and transmission, are required to better constrain the atmospheric properties of WASP-103b.

Hubble PanCET: an isothermal day-side atmosphere for the bloated gas-giant HAT-P-32Ab

We present a thermal emission spectrum of the bloated hot Jupiter HAT-P-32Ab from a single eclipse observation made in spatial scan mode with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST). The spectrum covers the wavelength regime from 1.123 to 1.644 microns which is binned into 14 eclipse depths measured to an averaged precision of 104 parts-per million. The spectrum is unaffected by a dilution from the close M-dwarf companion HAT-P-32B, which was fully resolved. We complemented our spectrum with literature results and performed a comparative forward and retrieval analysis with the 1D radiative-convective ATMO model. Assuming solar abundance of the planet atmosphere, we find that the measured spectrum can best be explained by the spectrum of a blackbody isothermal atmosphere with Tp = 1995 +/- 17K, but can equally-well be described by a spectrum with modest thermal inversion. The retrieved spectrum suggests emission from VO at the WFC3 wavelengths and no evidence of the 1.4 micron water feature. The emission models with temperature profiles decreasing with height are rejected at a high confidence. An isothermal or inverted spectrum can imply a clear atmosphere with an absorber, a dusty cloud deck or a combination of both. We find that the planet can have continuum of values for the albedo and recirculation, ranging from high albedo and poor recirculation to low albedo and efficient recirculation. Optical spectroscopy of the planet's day-side or thermal emission phase curves can potentially resolve the current albedo with recirculation degeneracy.

Haze heats Pluto's atmosphere yet explains its cold temperature.

Pluto's atmosphere is cold and hazy. Recent observations have shown it to be much colder than predicted theoretically, suggesting an unknown cooling mechanism. Atmospheric gas molecules, particularly water vapour, have been proposed as a coolant; however, because Pluto's thermal structure is expected to be in radiative-conductive equilibrium, the required water vapour would need to be supersaturated by many orders of magnitude under thermodynamic equilibrium conditions. Here we report that atmospheric hazes, rather than gases, can explain Pluto's temperature profile. We find that haze particles have substantially larger solar heating and thermal cooling rates than gas molecules, dominating the atmospheric radiative balance from the ground to an altitude of 700 kilometres, above which heat conduction maintains an isothermal atmosphere. We conclude that Pluto's atmosphere is unique among Solar System planetary atmospheres, as its radiative energy equilibrium is controlled primarily by haze particles instead of gas molecules. We predict that Pluto is therefore several orders of magnitude brighter at mid-infrared wavelengths than previously thought-a brightness that could be detected by future telescopes.

Well-Posedness of a Problem of Propagation of Nonlinear Acoustic-Gravity Waves in the Atmosphere Generated by Surface Pressure Variations

Currently there are many international microbarograph networks for high-resolution recording of wave pressure variations on the Earth’s surface. This arouses interest in wave propagation in the atmosphere generated by atmospheric pressure variations. A full system of nonlinear hydrodynamic equations for atmospheric gases with lower boundary conditions in the form of wavelike pressure variations on the Earth’s surface is considered. Since the wave amplitudes near the Earth’s surface are small, linearized equations are used in the analysis of well-posedness of the problem. With the help of a wave energy functional method, it is shown that in the non-dissipative case the solution to the boundary value problem is uniquely determined by the variable pressure field on the Earth’s surface. The corresponding dissipative problem is well-posed if, in addition to the pressure field, appropriate conditions on the velocity and temperature on the Earth’s surface are given. In the case of an isothermal atmosphere, the problem admits analytical solutions that are harmonic in the variables x and t. A good agreement between the numerical and analytical solutions is obtained. The study shows that the temperature and density can rapidly vary at the lower boundary of the boundary value problem. An example of solving the three-dimensional problem with variable pressure on the Earth’s surface taken from experimental observations is given. The developed algorithms and computer programs can be used to simulate atmospheric waves generated by pressure variations on the Earth’s surface.

Exoplanet phase curves at large phase angles. Diagnostics for extended hazy atmospheres

At optical wavelengths, Titan's brightness for large Sun-Titan-observer phase angles significantly exceeds its dayside brightness. The brightening that occurs near back-illumination is due to moderately large haze particles in the moon's extended atmosphere that forward-scatter the incident sunlight. Motivated by this phenomenon, here we investigate the forward scattering from currently known exoplanets, its diagnostics possibilities, the observational requirements to resolve it, and potential implications. An analytical expression is derived for the amount of starlight forward-scattered by an exponential atmosphere that takes into account the finite angular size of the star. We use this expression to tentatively estimate how prevalent this phenomenon may be. Based on numerical calculations that consider exoplanet visibility, we identify numerous planets with predicted out-of-transit forward scattering signals of up to tens of parts-per-million provided that aerosols of $>$1 micron size form over an extended vertical region near the optical radius level. We propose that the interpretation of available optical phase curves should be revised to constrain the strength of this phenomenon that might provide insight into aerosol scale heights and particle sizes. For the relatively general atmospheres considered here, forward scattering reduces the transmission-only transit depth by typically less than the equivalent to a scale height. For short-period exoplanets the finite angular size of the star severely affects the amount of radiation scattered towards the observer at mid-transit.

Nonlinear Effects in Three-minute Oscillations of the Solar Chromosphere. I. An Analytical Nonlinear Solution and Detection of the Second Harmonic

Abstract The vertical propagation of nonlinear acoustic waves in an isothermal atmosphere is considered. A new analytical solution that describes a finite-amplitude wave of an arbitrary wavelength is obtained. Although the short- and long-wavelength limits were previously considered separately, the new solution describes both limiting cases within a common framework and provides a straightforward way of interpolating between the two limits. Physical features of the nonlinear waves in the chromosphere are described, including the dispersive nature of low-frequency waves, the steepening of the wave profile, and the influence of the gravitational field on wavefront breaking and shock formation. The analytical results suggest that observations of three-minute oscillations in the solar chromosphere may reveal the basic nonlinear effect of oscillations with combination frequencies, superposed on the normal oscillations of the system. Explicit expressions for a second-harmonic signal and the ratio of its amplitude to the fundamental harmonic amplitude are derived. Observational evidence of the second harmonic, obtained with the Fast Imaging Solar Spectrograph, installed at the 1.6 m New Solar Telescope of the Big Bear Observatory, is presented. The presented data are based on the time variations of velocity determined from the Na i D2 and Hα lines.

The theory of transmission spectra revisited: a semi-analytical method for interpreting WFC3 data and an unresolved challenge

The computation of transmission spectra is a central ingredient in the study of exoplanetary atmospheres. First, we revisit the theory of transmission spectra, unifying ideas from several workers in the literature. Transmission spectra lack an absolute normalization due to the a priori unknown value of a reference transit radius, which is tied to an unknown reference pressure. We show that there is a degeneracy between the uncertainty in the transit radius, the assumed value of the reference pressure (typically set to 10 bar) and the inferred value of the water abundance when interpreting a WFC3 transmission spectrum. Second, we show that the transmission spectra of isothermal atmospheres are nearly isobaric. We validate the isothermal, isobaric analytical formula for the transmission spectrum against full numerical calculations and show that the typical errors are ~0.1% (~10 ppm) within the WFC3 range of wavelengths for temperatures of 1500 K (or higher). Third, we generalize the previous expression for the transit radius to include a small temperature gradient. Finally, we analyze the measured WFC3 transmission spectrum of WASP-12b and demonstrate that we obtain consistent results with the retrieval approach of Kreidberg et al. (2015) if the reference transit radius and reference pressure are fixed to assumed values. The unknown functional relationship between the reference transit radius and reference pressure implies that it is the product of the water abundance and reference pressure that is being retrieved from the data, and not just the water abundance alone. This degeneracy leads to a limitation on how accurately we may extract molecular abundances from transmission spectra using WFC3 data alone. Finally, we compare our study to that of Griffith (2014) and discuss why the degeneracy was missed in previous retrieval studies. [abridged]

Spectral line intensity and polarization in gas-dusty medium

We study the standard spectral line radiative transfer equation for media consisting of resonant atoms and non-resonant components (the dust grains and the atoms without considered spectral transition). Our goal is to study the intensity and polarization of the resonance radiation emerging from semi-infinite atmosphere. Using the known technique of resolvent matrices, we obtain the exact solution of vectorial radiative transfer equation for various sources of non-polarized radiation in semi-infinite atmosphere. Homogeneous, linear increasing and exponentially decreasing sources are considered. Recall, that homogeneous sources correspond to the isothermal atmosphere and the exponential ones correspond to the incident radiation with different angles of inclination.

Hybrid Galerkin numerical modelling of elastodynamics and compressible Navier–Stokes couplings: applications to seismo-gravito acoustic waves

We introduce a hybrid Galerkin modelling tool for the nonlinear acoustic and gravity wave propagation in planetary atmospheres coupled through topography to a solid medium. We rely on a 2-D spectral-element technique to model linear visco-elastic solid media and couple it to a discontinuous Galerkin method for the atmosphere modelled by the fully nonlinear Navier–Stokes equations. Significant benefits of such a method are, first, its versatility because it handles both acoustic and gravity waves in the same simulation, second, it enables one to observe nonlinear effects as convection or wave-breaking and, finally, it allows one to study the impact of ground-atmosphere coupling for waves propagating from seismic sources. Simulations are performed for 2-D isothermal atmosphere models with complex wind and viscosity profiles. We validate the computations by comparing them to finite-difference solutions, already validated in a previous paper. Specific benchmark validation cases are considered for both acoustic and gravity waves subject to viscosity variations, wind duct and nonlinear wave breaking. We apply this tool to study acoustic and gravity waves generated by a strong seismic source and its nonlinear breaking in the upper atmosphere.

Importance of Cross-redistribution in Scattering Polarization of Spectral Lines:The Cases of 3P−3S Triplets of Mg i and Ca i

Abstract Scattering on a multi-level atomic system has dominant contributions from resonance and Raman scattering. While initial and final levels are the same for resonance scattering, they are different for Raman scattering. The frequency redistribution for resonance scattering is described by the usual partial frequency redistribution functions of Hummer, while that for Raman scattering is described by cross-redistribution (XRD) function. In the present paper, we investigate the importance of XRD on linear polarization profiles of 3P−3S triplets of Mg i and Ca i formed in an isothermal one-dimensional atmosphere. We show that XRD produces significant effects on the linear polarization profiles when the wavelength separations between the line components of the multiplet are small, like in the cases of Mg i b and Ca i triplets.

A clear link connecting the troposphere and ionosphere: ionospheric reponses to the 2015 Typhoon Dujuan

The global navigation satellite system (GNSS) total electron content (TEC) sequences were used to capture the arrival time and location of the ionosphere disturbances in response to the 2015 Typhoon Dujuan. After removing the de-trended TEC variation, the clear ionosphere disturbances on the typhoon landing day could be distinguished, and these disturbances disappeared from the TEC sequences before and after the typhoon landing day. The foF2 data observed by Xiamen ionosonde station also show ionosphere disturbances. Based on the advantages of GNSS multi-point observations, the disturbances horizontal velocity in the ionosphere were estimated according to the linear theory for a dispersion relation of acoustic gravity waves (AGWs) in an isothermal atmosphere. The average horizontal velocity (\(\sim \)240 m/s) and the radial velocity (\(\sim \)287 m/s) were used in the two-dimensional grid search for the origin point on the Earth’s surface. The origin area was determined to be on the eastern side of Taiwan. Lastly, a possible physical mechanism is discussed in this study. When typhoons land on Taiwan, the severe convective storms and the drag effect from the Central Mountains create an ideal location for development of AGWs. Topographic conditions, like the high lapse rate, contribute to the formation of AGWs, which then propagates into the ionosphere altitude.

Coupling a Gas Chromatography Unit to an Inductively Coupled Plasma Mass Spectrometer

Although the eminent threat of a terrorist group detonating an improvised nuclear device (IND) in an urban environment is low, it is crucial that countries develop modern nuclear forensic capabilities to expedite response in a post-detonation scenario. In particular, new instruments need to be created to shorten dissolution time, expedite chemical separation, and improve forensic analysis of the nuclear melt glass that is created during the detonation of the device. To expedite this process, an instrument was designed to thermally couple a gas chromatograph (GC) to a time-of-flight inductively coupled plasma time-of-flight mass spectrometer (ICPTOFMS) In order to couple these two instruments, another instrument was designed to provide an isothermal atmosphere between the GC and TOFICPMS to expedite rapid gas separations processes. By using gas separations instead of the current wet chemistry processes, the required separation and analysis time of the melt glass significantly decreases. The new instrument would also provide a more detailed analysis of the elemental and isotopic composition of the melt glass. By completing these tasks simultaneously, this significantly decreases the required time to conduct these separations and improves the elemental and isotopic analysis.

EXPLORING BIASES OF ATMOSPHERIC RETRIEVALS IN SIMULATED JWST TRANSMISSION SPECTRA OF HOT JUPITERS

ABSTRACT With a scheduled launch in 2018 October, the James Webb Space Telescope (JWST) is expected to revolutionize the field of atmospheric characterization of exoplanets. The broad wavelength coverage and high sensitivity of its instruments will allow us to extract far more information from exoplanet spectra than what has been possible with current observations. In this paper, we investigate whether current retrieval methods will still be valid in the era of JWST, exploring common approximations used when retrieving transmission spectra of hot Jupiters. To assess biases, we use 1D photochemical models to simulate typical hot Jupiter cloud-free atmospheres and generate synthetic observations for a range of carbon-to-oxygen ratios. Then, we retrieve these spectra using TauREx, a Bayesian retrieval tool, using two methodologies: one assuming an isothermal atmosphere, and one assuming a parameterized temperature profile. Both methods assume constant-with-altitude abundances. We found that the isothermal approximation biases the retrieved parameters considerably, overestimating the abundances by about one order of magnitude. The retrieved abundances using the parameterized profile are usually within 1σ of the true state, and we found the retrieved uncertainties to be generally larger compared to the isothermal approximation. Interestingly, we found that by using the parameterized temperature profile we could place tight constraints on the temperature structure. This opens the possibility of characterizing the temperature profile of the terminator region of hot Jupiters. Lastly, we found that assuming a constant-with-altitude mixing ratio profile is a good approximation for most of the atmospheres under study.

POLARIZED LINE FORMATION IN NON-MONOTONIC VELOCITY FIELDS

ABSTRACT For a correct interpretation of the observed spectro-polarimetric data from astrophysical objects such as the Sun, it is necessary to solve the polarized line transfer problems taking into account a realistic temperature structure, the dynamical state of the atmosphere, a realistic scattering mechanism (namely, the partial frequency redistribution—PRD), and the magnetic fields. In a recent paper, we studied the effects of monotonic vertical velocity fields on linearly polarized line profiles formed in isothermal atmospheres with and without magnetic fields. However, in general the velocity fields that prevail in dynamical atmospheres of astrophysical objects are non-monotonic. Stellar atmospheres with shocks, multi-component supernova atmospheres, and various kinds of wave motions in solar and stellar atmospheres are examples of non-monotonic velocity fields. Here we present studies on the effect of non-relativistic non-monotonic vertical velocity fields on the linearly polarized line profiles formed in semi-empirical atmospheres. We consider a two-level atom model and PRD scattering mechanism. We solve the polarized transfer equation in the comoving frame (CMF) of the fluid using a polarized accelerated lambda iteration method that has been appropriately modified for the problem at hand. We present numerical tests to validate the CMF method and also discuss the accuracy and numerical instabilities associated with it.

First-order analytic propagation of satellites in the exponential atmosphere of an oblate planet

The paper offers the fully analytic solution to the motion of a satellite orbiting under the influence of the two major perturbations, due to the oblateness and the atmospheric drag. The solution is presented in a time-explicit form, and takes into account an exponential distribution of the atmospheric density, an assumption that is reasonably close to reality. The approach involves two essential steps. The first one concerns a new approximate mathematical model that admits a closed-form solution with respect to a set of new variables. The second step is the determination of an infinitesimal contact transformation that allows to navigate between the new and the original variables. This contact transformation is obtained in exact form, and afterwards a Taylor series approximation is proposed in order to make all the computations explicit. The aforementioned transformation accommodates both perturbations, improving the accuracy of the orbit predictions by one order of magnitude with respect to the case when the atmospheric drag is absent from the transformation. Numerical simulations are performed for a low Earth orbit starting at an altitude of 350 km, and they show that the incorporation of drag terms into the contact transformation generates an error reduction by a factor of 7 in the position vector. The proposed method aims at improving the accuracy of analytic orbit propagation and transforming it into a viable alternative to the computationally intensive numerical methods.

Similarity analysis of the streamer zone of Blue jets

Multiple observations of Blue Jets (BJ) and Gigantic Blue Jets (GBJ) show that BJ and GBJ emits a fan of streamers similar to a laboratory leader. Moreover, in the exponential atmosphere those long streamers grow preferentially upward, producing a narrow cone confined by the aperture angle. It was also noticed that BJ are similar to the streamer zone of a leader (streamer corona) and the modeling studies based on the streamers fractal structure were conducted. Objective of this paper is to study the fractal dimension of the bunch of streamer channels emitted by BJ and GBJ, at different altitude and under the varying electric field. This similarity analysis has been done in three steps: first we described the dendritic structure of streamers in corona discharge applying the fractal theory. Then using this model and the data from existing laboratory experiments we obtained the fractal dimension of the branching streamer channels. Finally the model was validated by the observations of BJ available from the literature.

Exploring the atmosphere using smartphones

The characteristics of the inner layer of the atmosphere, the troposphere, are determinant for Earth's life. In this experience we explore the first hundreds of meters using a smartphone mounted on a quadcopter. Both the altitude and the pressure are obtained using the smartphone's sensors. We complement these measures with data collected from the flight information system of an aircraft. The experimental results are compared with the International Standard Atmosphere and other simple approximations: isothermal and constant density atmospheres.

Oscillations of a vertically stratified dissipative atmosphere. I. Solution above source

A method of construction of solution for acoustic-gravity waves (AGW) above a wave source, taking dissipation throughout the atmosphere into account (dissipative solution above source, DSAS), is proposed. The method is to combine three solutions for three parts of the atmosphere: an analytical solution for the upper isothermal part and numerical solutions for the real non-isothermal dissipative atmosphere in the middle part and for the real non-isothermal small dissipation atmosphere in the lower one. In this paper the method has been carried out for the atmosphere with thermal conductivity but without viscosity. The heights of strong dissipation and the total absorption index in the regions of weak and average dissipation are found. For internal gravity waves the results of test calculations for an isothermal atmosphere and calculations for a real non-isothermal atmosphere are shown in graphical form. An algorithm and appropriate code to calculate DSAS, taking dissipation due to finite thermal conductivity into account throughout the atmosphere, are developed. The results of test DSAS calculations for an everywhere isothermal atmosphere are given. The calculation results for DSAS for the real non-isothermal atmosphere are also presented. A method for construction of the 2×2 Green's matrix fully taking dissipation into account and allowing us to find disturbance from some source of AGW in the atmosphere is proposed.

Approximate Method for Estimating Gaseous Loss at Very Low Angles Along Paths of Finite Length and Earth-Space Paths

This paper presents an approximate method for estimating gaseous absorption at very low elevation angles along paths of finite length or Earth-space paths, taking into account both the effects of the Earth curvature and refraction. An exponential atmosphere is assumed. Gaseous absorption is evaluated by breaking its integral representation into a number of smaller pieces that are evaluated in closed form by the leading term of their asymptotic expansions; there is no need for numerical integration. The new method extends the validity of Annex 2 of Recommendation ITU-R P.676-10 to path elevation angles below 5°. Comparisons with the line-by-line method of the same Recommendation, using radiosonde data as input, show good to very good agreement at 20, 30, and 350 GHz for an Arctic location, a tropical area in Australia, and two midlatitude locations-one in the Rocky Mountains and the other closer to sea level. Some limitations of the approximate method over very high mountains, where surface pressure is lower than the 780-hPa reference of the water vapor attenuation model of Annex 2, are also discussed.