Spectrophotometric Variability of Mira

Data from the Pulkovo spectrophotometric data base on the absolute quasimonochromatic fluxes from oCet in the 320–1080 nm range are used to determine the physical parameters of this star in different phases of its light curve. The continuum emission layer is found to expand between the phases of the cycle corresponding to the rising and falling branches of the light curve. The average expansion velocity is 32 km/s. By the time the star’s brightness has fallen by roughly three magnitudes, its radius has increased by almost a factor of three. Over this same time the temperature of the layer has fallen from 3000 K to 2200 K. For this expansion velocity, the calculated mass rate loss is ⊙ M /year.__________Translated from Astrofizika, Vol. 48, No. 2, pp. 175–189 (May 2005).     

New photometric analysis of hot Jupiter: WASP-135b

We report photometric follow-up observations of WASP-135b obtained using the 1.23-m telescope at Calar Alto Observatory and 1.00-m telescope at TÜBİTAK National Observatory during the 2017 and 2018 observational seasons. Eight new transit light curves of WASP-135b were analyzed with jktebop code. The ratio of the planet radius to radius of host star, fractional radius of host star, and orbital inclination of WASP-135b were found to be 0.138 ± 0.002, 0.181 ± 0.008 and 82.44 ± 0.64 degrees, respectively. Planetary radius of WASP-135b was derived from transit parameters to be 1.075 ± 0.150 R_J. The transit ephemeris of WASP-135b was also updated using the maximum likelihood method (MLM). 165 well-known hot Jupiters (HJs) were selected from the Exoplanet Data Explorer database and the classification of these HJs together with WASP-135b, based on their equilibrium temperatures and Safronov numbers, is discussed.     

Spectroscopic and photometric study of the contact binary BO CVn

We present the results of the study of the contact binary system BO CVn. We have obtained physical parameters of the components based on combined analysis of new, multi-color light curves and spectroscopic mass ratio. This is the first time the latter has been determined for this object. We derived the contact configuration for the system with a very high filling factor of about 88%. We were able to reproduce the observed light curve, namely the flat bottom of the secondary minimum, only if a third light has been added into the list of free parameters. The resulting third light contribution is significant, about 20–24%, while the absolute parameters of components are: M₁ = 1.16, M₂ = 0.39, R₁ = 1.62 and R₂ = 1.00 (in solar units).The O-C diagram shows an upward parabola which, under the conservative mass transfer assumption, would correspond to a mass transfer rate of dM/dt = 6.3 × 10^?8M_⊙/yr, matter being transferred from the less massive component to the more massive one. No cyclic, short-period variations have been found in the O-C diagram (but longer-term variations remain a possibility).     

A photochemical model for the Venus atmosphere at 47–112 km

The model is intended to respond to the recent findings in the Venus atmosphere from the Venus Express and ground-based submillimeter and infrared observations. It extends down to 47 km for comparison with the kinetic model for the lower atmosphere (Krasnopolsky, V.A. [2007]. Icarus 191, 25–37) and to use its results as the boundary conditions. The model numerical accuracy is significantly improved by reduction of the altitude step from 2 km in the previous models to 0.5 km. Effects of the NUV absorber are approximated using the detailed photometric observations at 365 nm from Venera 14. The H₂O profile is not fixed but calculated in the model. The model involves odd nitrogen and OCS chemistries based on the detected NO and OCS abundances. The number of the reactions is significantly reduced by removing of unimportant processes. Column rates for all reactions are given, and balances of production and loss may be analyzed in detail for each species.The calculated vertical profiles of CO, H₂O, HCl, SO₂, SO, OCS and of the O₂ dayglow at 1.27 μm generally agree with the existing observational data; some differences are briefly discussed. The OH dayglow is ～30 kR, brighter than the OH nightglow by a factor of 4. The H + O₃ process dominates in the nightglow excitation and O + HO₂ in the dayglow, because of the reduction of ozone by photolysis. A key feature of Venus’ photochemistry is the formation of sulfuric acid in a narrow layer near the cloud tops that greatly reduces abundances of SO₂ and H₂O above the clouds. Delivery of SO₂ and H₂O through this bottleneck determines the chemistry and its variations above the clouds. Small variations of eddy diffusion near 60 km result in variations of SO₂, SO, and OCS at and above 70 km within a factor of ～30. Variations of the SO₂/H₂O ratio at the lower boundary have similar but weaker effect: the variations within a factor of ～4 are induced by changes of SO₂/H₂O by ±5%. Therefore the observed variations of the mesospheric composition originate from minor variations of the atmospheric dynamics near the cloud layer and do not require volcanism. NO cycles are responsible for production of a quarter of O₂, SO₂, and Cl₂ in the atmosphere. A net effect of photochemistry in the middle atmosphere is the consumption of CO₂, SO₂, and HCl from and return of CO, H₂SO₄, and SO₂Cl₂ to the lower atmosphere. These processes may be balanced by thermochemistry in the lower atmosphere even without outgassing from the interior, though the latter is not ruled out by our models. Some differences between the model and observations and the previous models are briefly discussed.     

Reading the climate record of the martian polar layered deposits

The martian polar regions have layered deposits of ice and dust. The stratigraphy of these deposits is exposed within scarps and trough walls and is thought to have formed due to climate variations in the past. Insolation has varied significantly over time and caused dramatic changes in climate, but it has remained unclear whether insolation variations could be linked to the stratigraphic record. We present a model of layer formation based on physical processes that expresses polar deposition rates of ice and dust in terms of insolation. In this model, layer formation is controlled by the insolation record, and dust-rich layers form by two mechanisms: (1) increased summer sublimation during high obliquity, and (2) variations in the polar deposition of dust modulated by obliquity variations. The model is simple, yet physically plausible, and allows for investigations of the climate control of the polar layered deposits (PLD). We compare the model to a stratigraphic column obtained from the north polar layered deposits (NPLD) (Fishbaugh, K.E., Hvidberg, C.S., Byrne, S., Russel, P.S., Herkenhoff, K.E., Winstrup, M., Kirk, R. [2010a]. Geophys. Res. Lett., 37, L07201) and show that the model can be tuned to reproduce complex layer sequences. The comparison with observations cannot uniquely constrain the PLD chronology, and it is limited by our interpretation of the observed stratigraphic column as a proxy for NPLD composition. We identified, however, a set of parameters that provides a chronology of the NPLD tied to the insolation record and consistently explains layer formation in accordance with observations of NPLD stratigraphy. This model dates the top 500 m of the NPLD back to ～1 million years with an average net deposition rate of ice and dust of 0.55 mm a^?1. The model stratigraphy contains a quasi-periodic ～30 m cycle, similar to a previously suggested cycle in brightness profiles from the NPLD (Laskar, J., Levrard, B., Mustard, F. [2002]. Nature, 419, 375–377; Milkovich, S., Head, J.W. [2005]. J. Geophys. Res. 110), but here related to half of the obliquity cycles of 120 and 99 kyr and resulting from a combination of the two layer formation mechanisms. Further investigations of the non-linear insolation control of PLD formation should consider data from other geographical locations and include radar data and other stratigraphic datasets that can constrain the composition and stratigraphy of the NPLD layers.     

A new class of charged analogues of Vaidya–Tikekar type super-dense star

First ever closed form solution for charged fluid sphere expressed by a space time with its hypersurfaces t= constant as spheroid is obtained for the case 0<K<1. The same is utilized to construct a superdense star with surface density 2×10¹⁴ gm/cm³. The star is seen to satisfy the reality and causality conditions for 0<K≤0.045 and possesses maximum mass and radius to be 0.065216M _Θ and 1.137496 km respectively. Moreover the interior of the star satisfy strong energy condition. However in the absence of the causality condition, the reality conditions are valid for a wider range 0<K≤0.13. The maximum mass and radius for the later case are 1.296798M _Θ and 2.6107 km respectively for the strong energy condition, while the said parameters for the weak energy condition read as 1.546269M _Θ and 2.590062 km respectively.     

Restitution coefficients and sticking velocities of a chondrule analogue colliding on a porous silica layer at impact velocities between 0.1 and 80 m s−1

To study the accretional growth of rimmed chondrules and their agglomerates in the solar nebula, we measured the restitution coefficients, ε, and the sticking velocities to a porous silica layer, v_c, by impacting the silica layer with a glass ball at velocities from 0.1 to 80 m s^?1. We used a porous silica layer covering a basalt block with thicknesses ranging from 1/5 of the glass ball radius to equal to the glass ball radius as a rimmed chondrule analogue, and the porosity of the silica layer was set to be 70%, 80%, 85%, and 90%. Collisional experiments were conducted by means of the free fall method or by the use of a spring gun or a gas gun, allowing us to vary the impact velocity. We used a laser displacement meter to estimate the impact and rebound velocities as well as the acceleration during the collision at impact velocities below 1 m s^?1. As a result, the sticking velocity, v_c, of 90%- and 85%- porosity layers with a thickness equal to 1/2 of the glass ball diameter was 0.44 and 2.4 m s^?1, respectively. On the other hand, we found a distinct barrier to sticking for smaller-porosity layers: the silicate layer with a porosity smaller than 80% never exhibited sticking at any impact velocity below 1 m s^?1. Instead, we observed a rebound effect with restitution coefficients larger than 0.2. In the case of a silica layer with a porosity smaller than 80%, we observed the sub-sticking condition defined by ε < 0.1 at velocities extending from 5 m s^?1 to 70 m s^?1.     

Spatial distributions and seasonal variations of features related to a venting process at high southern latitudes observed by the MOC camera

A distinctive terrain named cryptic region which is characterized by regions of low albedo and low temperature has been identified on the Martian south polar cap. In this zone, many fan- and spider-shaped features of km-scale appeared following the sublimation of the CO₂ frost layer. These peculiar features were apparently caused by a wind-blown system of dust-laden jets. During the warming period starting at L_s∼180°, the seasonal ice cap regresses and fans and spiders appear in sequence. These surface features are repeatable events that tend to occupy the same areas from year to year. In this study, we use the Mars Orbiter camera (MOC) narrow-angle images to produce a statistical study of the time distributions of the fans and spiders as functions of L_s and as functions of the topography. The time variations and spatial distributions of these features are further correlated with the CO₂ ice coverage measured by the Mars Orbiter laser altimeter (MOLA) instrument. We have documented that most of the fans are found in the early spring with L_s<230° and the fans and spiders coexist at L_s=250°±20°. It is also found that there is a strong dependence on latitude and altitude with fans and spiders most often observed at high latitude (>83°S) and high altitude (>2500 m). Our statistical result also indicates that the occurrence of fans is highly correlated with the thickness of the CO₂ frost thus providing support for the venting model.     

Quantitative estimation of helium-3 spatial distribution in the lunar regolith layer

³He (helium-3) in the lunar regolith implanted by the solar wind is one of the most valuable resources because of its potential as a fusion fuel. The abundance of ³He in the lunar regolith is related to solar wind flux, lunar surface maturity and TiO₂ content, etc. A model of solar wind flux, which takes account of variations due to shielding of the nearside when the Moon is in the Earth's magnetotail, is used to present a global distribution of relative solar wind flux over the lunar surface. Using Clementine UV/VIS multispectral data, the global distribution of lunar surface optical maturity (OMAT) and the TiO₂ content in the lunar regolith are calculated. Based on Apollo regolith samples, a linear relation between ³He abundance and normalized solar wind flux, optical maturity, and TiO₂ content is presented. To simulate the brightness temperature of the lunar surface, which is the mission of the Chinese Chang-E project's multichannel radiometers, a global distribution of regolith layer thickness is first empirically constructed from lunar digital elevation mapping (DEM). Then an inversion approach is presented to retrieve the global regolith layer thickness. It finally yields the total amount of ³He per unit area in the lunar regolith layer, which is related to the regolith layer thickness, solar wind flux, optical maturity and TiO₂ content, etc. The global inventory of ³He is estimated as 6.50×¹⁰⁸ kg, where 3.72×¹⁰⁸ kg is for the lunar nearside and 2.78×¹⁰⁸ kg is for the lunar farside.     

The refined physical properties of the transiting exoplanetary system WASP‐41

We present the first follow‐up study of the transiting system WASP‐41 after its discovery in 2011. The main goal of this paper is to refine the physical parameters of the system and search for possible signs of transit timing variations. The observations used for analysis were taken from the public archive of the Exoplanet Transit Database (ETD). Assuming three different limb darkening laws, we found the best‐fitting model and redetermined parameters of the system. Although system parameters obtained in this study were found to be in good agreement with the discovery study, the planetary radius determined here is notably smaller, namely 1.12^+0.06 _–0.07 R_Jup. The Safronov number Θ = 0.071 ± 0.002 and equilibrium temperature T_eq = 1271 ± 50 K were also determined. Both values indicate that the planet WASP‐41b belongs to Class I of transiting planets. No significant transit timing variations were detected. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Kinematics of Umbral Dots: Their Typical Area,Brightness and Lifetime

The statistical properties of area and brightness of umbral dots, formed in a sunspot umbra observed on 18 June 2004, are studied using an improved method of image segmentation and feature tracking algorithm. Central (peripheral) umbral dots have a typical size around 0.16 (0.17) arcsec. The brightness distribution of umbral dots shows a multi-population distribution. The brightness of the central umbral dots does not exceed 0.6 I _quiet. In most cases, the area of central (peripheral) umbral dots reaches its maximum 40 s (150 s) after their brightness maximum, respectively. The brightness of umbral dots does not show any directly proportional behavior with their lifetimes and areas. The temporal variations of some physical parameters of umbral dots are studied. It seems that these variations are not physical processes in an umbra and are correlated to the images quality (seeing variations). The lifetime distribution of umbral dots is an exponential function and central (peripheral) umbral dots have a typical half-life of about 300 s (180 s), respectively.     

Spatial variability,composition and thickness of the seasonal north polar cap of Mars in mid-spring

The planetary fourier spectrometer (PFS) for the Mars express mission (MEX) is an infrared spectrometer operating in the wavelength range from 1.2 to 45 μm by means of two spectral channels, called SWC (short wavelength channel) and LWC (long wavelength channel), covering, respectively, 1.2–5.5 and 5.5–45 μm.The middle-spring Martian north polar cap (Ls∼40°) has been observed by PFS/MEX in illuminated conditions during orbit 452. The SWC spectra are here used to study the cap composition in terms of CO₂ ice, H₂O ice and dust content. Significant spectral variation is noted in the cap interior, and regions of varying CO₂ ice grain sizes, water frost abundance, CO₂ ice cover and dust contamination can be distinguished. In addition, we correlate the infrared spectra with an image acquired during the same orbit by the OMEGA imaging spectrometer and with the altimetry from MOLA data. Many of the spectra variations correlate with heterogeneities noted in the image, although significant spectral variations are not discernible in the visible. The data have been divided into five regions with different latitude ranges and strong similarities in the spectra, and then averaged. Bi-directional reflectance models have been run with the appropriate lighting geometry and used to fit the observed data, allowing for CO₂ ice and H₂O ice grain sizes, dust and H₂O ice contaminations in the form of intimate granular mixtures and spatial mixtures.A wide annulus of dusty water ice surrounds the recessing CO₂ seasonal cap. The inner cap exhibits a layered structure with a thin CO₂ layer with varying concentrations of dark dust, on top of an H₂O ice underneath ground. In the best-fits, the ices beneath the top layer have been considered as spatial mixtures. The results are still very good everywhere in the spectral range, except where the CO₂ ice absorption coefficients are such that even a thin layer is enough to totally absorb the incoming radiation (i.e. the band is saturated). This only happens around 3800 cm⁻¹, inside the strong 2.7-μm CO₂ ice absorption band. The effect of finite snow depth has been investigated through a layered albedo model. The thickness of the CO₂ ice deposits increases with latitude, ranging from 0.5–1 g cm⁻² within region II to 60–80 g cm⁻² within the highest-latitude (up to 84°N) region V.Region I is at the cap edge and extends from 65°N to 72°N latitude. No CO₂ ice is present in this region, which consists of relatively large grains of water ice (20 μm), highly contaminated by dust (0.15 wt%). The adjacent region II is a narrow region [76–79°N] right at the edge of the north residual polar cap. This region is very distinct in the OMEGA image, where it appears to surround the whole residual cap. The CO₂ ice features are barely visible in these spectra, except for the strong saturated 2.7 μm band. It basically consists of a thin layer of 5-mm CO₂ ice on top of an H₂O ice layer with the same composition as region I. A third interesting region III is found all along the shoulder of the residual cap [79–81°N]. It extends over 1.5 km in altitude and over only 2° of latitude and consists of CO₂ ice with a large dust content. It is an admixture of CO₂ ice (3–4 mm), with several tens of ppm by mass of water ice and more than 2 ppt by mass of dust. The surface temperatures have been retrieved from the LWC spectra for each observation. We found an increase in the surface temperature in this region, indicating a spatial mixture of cold CO₂ ice and warmer dust/H₂O ice. Region IV is close to the top of the residual cap [81–84°N]; it is much brighter than region III, with a dust content 10 times lower than the latter. The CO₂ grain size is 3 mm and strong CO₂ ice features are present in the data, indicating a thicker CO₂ ice layer than in region II (1–2 g cm⁻²). The final region V is right at the top of the residual cap (⩾84°N). It is “pure” CO₂ ice (no dust) of 5 mm grain sizes, with 30 ppm by weight of water ice. The CO₂ ice features are very pronounced and the 2.7 μm band is saturated. The optical thickness is close to the semi-infinite limit (30–40 g cm⁻²). Assuming a snowpack density of 0.5 g cm⁻³, we get a minimum thickness of 1–2 cm for the top-layer of regions II and III, 4–10 cm for region IV, and ⩾60–80 cm thickness for region V. These values are in close agreement with several recent results for the south seasonal polar cap.These results should provide new, useful constraints in models of the Martian climate system and volatile cycles.     

Calculation and visualisation of future glacier extent in the Swiss Alps by means of hypsographic modelling

The observed rapid glacier wastage in the European Alps during the past 20 years already has strong impacts on the natural environment (rock fall, lake formation) as well as on human activities (tourism, hydro-power production, etc.) and poses several new challenges also for glacier monitoring. With a further increase of global mean temperature in the future, it is likely that Alpine glaciers and the high-mountain environment as an entire system will further develop into a state of imbalance. Hence, the assessment of future glacier geometries is a valuable prerequisite for various impact studies. In order to calculate and visualize in a consistent manner future glacier extent for a large number of individual glaciers (> 100) according to a given climate change scenario, we have developed an automated and simple but robust approach that is based on an empirical relationship between glacier size and the steady-state accumulation area ratio (AAR₀) in the Alps. The model requires digital glacier outlines and a digital elevation model (DEM) only and calculates new glacier geometries from a given shift of the steady-state equilibrium line altitude (ELA₀) by means of hypsographic modelling. We have calculated changes in number, area and volume for 3062 individual glacier units in Switzerland and applied six step changes in ELA₀ (from + 100 to + 600 m) combined with four different values of the AAR₀ (0.5, 0.6, 0.67, 0.75). For an AAR₀ of 0.6 and an ELA₀ rise of 200 m (400 m) we calculate a total area loss of − 54% (− 80%) and a corresponding volume loss of − 50% (− 78%) compared to the 1973 glacier extent. In combination with a geocoded satellite image, the future glacier outlines are also used for automated rendering of perspective visualisations. This is a very attractive tool for communicating research results to the general public. Our study is illustrated for a test site in the Upper Engadine (Switzerland), where landscape changes above timberline play an important role for the local economy. The model is seen as a first-step approach, where several parts can be (and should be) further developed.     

Photometry of the long period dwarf nova GY Hya

Although comparatively bright, the cataclysmic variable GY Hya has not attracted much attention in the past. As part of a project to better characterize such systems photometrically, we observed light curves in white light, each spanning several hours, at Bronberg Observatory, South Africa, in 2004 and 2005, and at the Observatório do Pico dos Dias, Brazil, in 2014 and 2016. These data permit to study orbital modulations and their variations from season to season. The orbital period, already known from spectroscopic observations of Peters and Thorstensen (2005), is confirmed through strong ellipsoidal variations of the mass donor star in the system and the presence of eclipses of both components. A refined period of 0.34723972 (6) days and revised ephemeries are derived. Seasonal changes in the average orbital light curve can qualitatively be explained by variations of the contribution of a hot spot to the system light together with changes of the disk radius. The amplitude of the ellipsoidal variations and the eclipse contact phases permit to put some constraints on the mass ratio, orbital inclination and the relative brightness of the primary and secondary components. There are some indications that the disk radius during quiescence, expressed in units of the component separation, is smaller than in other dwarf novae.     

The cyclic variation of the outburst recurrence time in GRS 1747–312

We present an analysis of the long-term evolution of outbursts in the neutron star soft X-ray transient GRS 1747–312. Observations taken from ASM/RXTE, in the 1.5–12 keV passband, are utilized. We reveal a cyclic behavior in the residuals of the outburst recurrence time with respect to the mean value of T_C = 136 ± 2 days. The profile of this cycle is approximately sinusoidal; the remaining cycle-to-cycle fluctuations possess a considerably smaller amplitude. We find that, although the peak flux of the outbursts displays a significant scatter at a given phase of the cycle, the most luminous outbursts occur after the longest T_C. The fluence displays a large scatter for the individual outbursts and tends to decrease with time. We argue that although the cycle-length of ～5.4 yr is compatible with that of the presumed magnetic activity of the late-type donor, it cannot be explained by variations of the mass outflow from the donor to the disk. In our interpretation, the stellar activity is translated to variations of T_C via interaction of the magnetic field of the spots on the donor with the magnetic field of the disk. This gives rise to a variable efficiency of the removal of the angular momentum from the quiescent disk during the activity cycle of the donor. This mechanism can be strengthened by accompanying variations of the radius of the optically thin advection-dominated accretion flow in quiescence. We show that the peak mass accretion rate onto the neutron star in the individual outbursts of GRS 1747–312 is considerably more stable than in two other similar systems with frequent outbursts, Aql X-1 and 4U 1608–52; this allows the cyclic modulation of T_C to show itself in GRS 1747–312.     

The Influence of Coronal Background Structure with Nested and Closed Magnetic Fields on the Triggering of CME

The influence of coronal streamer background with nested and closed magnetic fields on the of the triggering of coronal mass ejections (CMEs) is investigated in the meridian plane. In the coronal streamers’ background magnetic structure there are three small-scale closed magnetic fields, of which the middle one has a direction opposite to that of the global dipolar field of coronal streamers. The trigger model of CMEs emerges from beneath this small-scale closed magnetic field and possesses a concentric circular structure with radius of a = 0.1R_s (R_s being the solar radius). The direction of the magnetic field in the front half of the CME trigger model is opposite to that of the small-scale closed field and is the same as that of the streamers’ global dipolar field. As revealed by numerical simulation, when the ratio of the plasma pressure at the center of the CME trigger model to the boundary pressure is m ≤ 2, then the emerging model can trigger CMEs. When m < 2, then it cannot. The error in this critical value of 2 is less than 1%.     

An Investigation of the Physical Mechanisms of Orbital Period Variations of the Semidetached Binary UW Vir

The orbital period variations of the Algol-type semidetached binary UW Vir are analyzed. It is shown that in addition to a long-term rapid increase (dP/dt = + 1.37 × 10⁻⁶ day/year), its orbit period has a variation with the period of 62.3 years. Based on the basic physical parameters given by Brancewicz and Dworak in 1980, the physical mechanisms causing the orbital period variations are investigated. The analysis indicates that the periodical variation of orbital period can be interpreted by the light-travel time effect due to the presence of a third body with the mass of M₃ ≥ 0.94 M. As no observational information has been reported for this tertiary component, it might be a compact object (e.g., a white dwarf). The long-term increase of orbital period can be explained in terms of the mass transfer from the secondary to the primary component (dM₂/dt = 1.43 × 10⁻⁷ M/year). This is in agreement with the semidetached configuration of the system with a lobe-filling secondary component. But according to the evolution theory of binaries, the Algol-type semidetached binary UW Vir should be at the evolutionary stage of slow mass transfer on the nuclear-reaction timescale of the secondary component. However, the analysis shows that the timescale for the periodical variation of orbital period is much shorter than the nuclear-reaction timescale of the secondary component, but close to the thermodynamic timescale of the secondary. This reveals that: (1) This binary system is at the evolutionary stage of rapid mass transfer on the thermodynamic timescale of the secondary component; or (2) The circumstellar matter of the system makes a contribution to the rapid increase of orbital period via the angular momentum transfer.     

New light curve analysis and period changes of the overcontact binary EQ Tauri

A photometric CCD study of EQ Tauris was carried out in the R, V and B bands. These new data were analyzed by using the PHOEBE and Binary Marker 3 programs, which yielded the geometrical system parameters and mass and radius of both stars. We find these to confirm the results of Pribulla and Vanko (2002). Using these new data together with data from the literature, the orbital period variations were studied by using the Kalimeris method. We find the same timescale of variability of the orbital period as Pribulla and Vanko (2002), namely about 49 years.     

A multi-wavelength study of the 2009 impact on Jupiter: Comparison of high resolution images from Gemini,Keck and HST

Within several days of A. Wesley’s announcement that Jupiter was hit by an object on UT 19 July 2009, we observed the impact site with (1) the Hubble Space Telescope (HST) at UV through visible (225–924 nm) wavelengths, (2) the 10-m W.M. Keck II telescope in the near-infrared (1–5 μm), and (3) the 8-m Gemini-North telescope in the mid-infrared (7.7–18 μm). All observations reported here were obtained between 22 and 25 July 2009. Observations at visible and near-infrared wavelengths show that large (～0.75-μm radius) dark (imaginary index of refraction m_i ～ 0.01–0.1) particulates were deposited at atmospheric pressures between 10 and 200–300 mbar; analysis of HST-UV data reveals that in addition smaller-sized (～0.1 μm radius) material must have been deposited at the highest altitudes (～10 mbar). Differences in morphology between the UV and visible/near-IR images suggest three-dimensional variations in particle size and density across the impact site, which probably were induced during the explosion and associated events. At mid-infrared wavelengths the brightness temperature increased due to both an enhancement in the stratospheric NH₃ gas abundance and the physical temperature of the atmosphere. This high brightness temperature coincides with the center part of the impact site as seen with HST. This observation, combined with (published) numerical simulations of the Shoemaker-Levy 9 impacts on Jupiter and the Tunguska airburst on Earth, suggests that the downward jet from the terminal explosion probably penetrated down to the ～700-mbar level.     

Quartz–coesite–stishovite relations in shocked metaquartzites from the Vredefort impact structure,South Africa

Coesite and stishovite are developed in shock veins within metaquartzites beyond a radius of ~30 km from the center of the 2.02 Ga Vredefort impact structure. This work focuses on deploying analytical field emission scanning electron microscopy, electron backscattered diffraction, and Raman spectrometry to better understand the temporal and spatial relations of these silica polymorphs. α-Quartz in the host metaquartzites, away from shock veins, exhibits planar features, Brazil twins, and decorated planar deformation features, indicating a primary (bulk) shock loading of >5 < 35 GPa. Within the shock veins, coesite forms anhedral grains, ranging in size from 0.5 to 4 μm, with an average of 1.25 μm. It occurs in clasts, where it displays a distinct jigsaw texture, indicative of partial reversion to a less dense SiO₂ phase, now represented by microcrystalline quartz. It is also developed in the matrix of the shock veins, where it is typically of smaller size (<1 μm). Stishovite occurs as euhedral acicular crystals, typically <0.5 μm wide and up to 15 μm in length, associated with clast–matrix or shock vein margin–matrix interfaces. In this context, the needles occur as radiating or subparallel clusters, which grow into/over both coesite and what is now microcrystalline quartz. Stishovite also occurs as more blebby, subhedral to anhedral grains in the vein matrix (typically <1 μm). We propose a model for the evolution of the veins (1) precursory frictional melting in a microfault (~1 mm wide) generates a molten matrix containing quartz clasts. This is followed by (2) arrival of the main shock front, which shocks to 35 GPa. This generates coesite in the clasts and in the matrix. (3) On initial shock release, the coesite partly reverts to a less dense SiO₂ phase, which is now represented by microcrystalline quartz. (4) With continued release, stishovite forms euhedral needle clusters at solid–liquid interfaces and as anhedral crystals in the matrix. (5) With decreasing pressure–temperature, the matrix completes crystallization to yield a microcrystalline quasi-igneous texture comprising quartz–coesite–stishovite–kyanite–biotite–alkali feldspar and accessory phases. It is possible that the shock vein represents the locus of a thermal spike within the bulk shock, in which case there is no requirement for additional pressure (i.e., the bulk shock was ≃35 GPa). However, if that pressure was not realized from the main shock, then supplementary pressure excursions within the vein would have been required. These could have taken the form of localized reverberations from wave trapping, or implosion processes, including pore collapse, phase change–initiated volume reduction, and melt cavitation.