Phaseless VLBI mapping of compact extragalactic radio sources

The problem of phaseless aperture synthesis is of current interest in phase-unstable VLBI with a small number of elements when either the use of closure phases is not possible (a two-element interferometer) or their quality and number are not enough for acceptable image reconstruction by standard adaptive calibration methods. Therefore, we discuss the problem of unique image reconstruction only from the spectrum magnitude of a source. We suggest an efficient method for phaseless VLBI mapping of compact extragalactic radio sources. This method is based on the reconstruction of the spectrum magnitude for a source on the entire UV plane from the measured visibility magnitude on a limited set of points and the reconstruction of the sought-for image of the source by Fienup's method from the spectrum magnitude reconstructed at the first stage. We present the results of our mapping of the extragalactic radio source 2200+420 using astrometric and geodetic observations on a global VLBI array. Particular attention is given to studying the capabilities of a two-element interferometer in connection with the putting into operation of a Russian-made radio interferometer based on Quasar RT-32 radio telescopes.     

On the Tomographic Reconstruction of the 3D Electron Density for the Solar Corona from STEREO COR1 Data

We present for the first time a three-dimensional reconstruction of the electron density in the corona at distances from 1.5R _⊙ to 4R _⊙ using COR1 STEREO observations. The reconstruction is performed using a regularized tomography inversion method for two biweekly periods corresponding to Carrington Rotations 2058 and 2066. Images from the two STEREO spacecraft are used to compare the reconstructed density structures with coronal features located by triangulation. We find that the location of a bright tip of a helmet streamer obtained from the tomographic reconstruction is in good agreement with the location obtained by triangulation. The reconstructed density structure of the equatorial streamer belt is largely consistent with the variation of the current sheet derived from a potential magnetic field extrapolation for most of the equatorial region and for an MHD model of the corona. A zero-value density region in the reconstruction is identified with a low-density region seen in an EUVI image below the reconstruction domain.     

TheU,B, V photoelectric brightness profiles of NGC 362

The photoelectric brightness profiles inU, B, andV colors of the globular cluster NGC 362 are reported down to about 9 magnitudes below the central value, covering a radial distance of about 10 arc min. They were derived fromE-W scans, with a large rectangular diaphragm, giving the integrated profile of the cluster in about 6 min. of time. The inverse problem of determining the brightness profiles from integrated luminosities is solved by means of the regularization method. The agreement with the data available in literature (photoelectric brightness measurements in the inner 80 arc sec, bright stars distribution in the range 80-330 arc sec) is generally good. TheB-V andU-B radial color distributions are presented. It may seem that there is no color gradient, outside the central region dominated by red giants.     

Saturn's rings. I. Optical thickness of rings A,B, D and structure of ring B

A photometric study of high-resolution (～0″.3) plates of Saturn taken at the Lowell Observatory in 1943 and 1945 is presented. N-S scans were taken over both the planet and rings. The excess brightness due to the planet seen through the rings is found by taking the difference between the central meridian (CM) scans and scans displaced by 5″.7. Adopting a value for the albedo of the planet, it is possible to obtain the optical thickness, τ_CM(r). In particular, for the regions of maximum brightness in rings A and B, we find τ_CM(I_{A max}) = 0.38 ± 0.11 and τ_CM(I_{B max}) = 0.61 ± 0.11. Observations by Barnard made in 1890 show evidence of ring D, recently discovered by Guerin (1969). The value for the optical thickness of this ring is τ_D(I_{D max}) = 0.03 ± 0.01. Ring B exhibits a pronounced (7–10%) decrease in brightness from the extremity of the major axis to the CM. After considering several possible explanations, we conclude that the ring particles are nonspherical and are in synchronous rotation around the planet with their long axis toward it. The mean value for the ratio of major to minor axis for the particles at 15″ is $\text{(a/b)}\text{? 1.08}$. Because of the shape and orientation of the particles, the optical thickness at the extremity of the major axis and at the CM are different for any saturnicentric latitude B ≠ 90°. Under these circumstances, only a minimum value for τ at the extremity can be derived.     

Correlation of near-infrared and optical variability of NGC 4151 in 2008–2013

We investigate the correlation between infrared (JHKL) and optical (B) fluxes of the variable nucleus of the Seyfert galaxy NGC4151 using partially published data for the last 6 years (2008–2013). We find that the lag of flux in all the infrared bands is the same, 40±6 days, to within the measurement accuracy. Variability in the J and K bands is not quite simultaneous, perhaps due to the differing contributions of the accretion disk in these bands. The lag found for the K band compared with the B band is not significantly different from earlier values obtained for the period 2000–2007. However, finding approximately the same lags in all IR bands for 2008–2013 differs from previous results at earlier epochs when the lag increased with increasing wavelength. Examples of almost the same lag in different IR bands are known for some other active nuclei. In the case of NGC 4151 it appears that the relative lags between the IR bands may be different in different years. The available data, unfortunately, do not allow us to investigate a possible change in the lags during the test interval. We discuss our results in the framework of the standard model where the variable infrared radiation is mainly due to thermal re-emission from the part of the dusty torus closest to the central source. There is also a contribution of some IR emission from the accretion disk, and this contribution increases with decreasing wavelength.     

Evolution of the dusty rings of Uranus

We present high quality images of the uranian ring system, obtained in August 2002, October 2003, and July 2004 at 2.2 μm with the adaptive optics camera NIRC2 on the Keck II telescope. Using these data, we report the first detection in backscattered light of a ring (which we refer to as the ζ ring) interior to Uranus' known rings. This ring consists of a generally uniform sheet of dust between 37,850 and 41,350 km with an equivalent width (in 2004; or ), and extends inward to 32,600 km at a gradually decreasing brightness. This ring might be related to the Voyager ring R/1986 U 2, although both its location and extent differ. This could be attributed to a difference in observing wavelength and/or solar phase angle, or perhaps to temporal variations in the ring. Through careful modeling of the I/F of the individual rings at each ansa, we reveal the presence of narrow (few 100 km wide) sheets of dust between the δ and ε rings, and between rings 4 and α. We derived a typical anisotropy factor g≈0.7 in the scattering behavior of these particles. The spatial distribution and relative intensity of these dust sheets is different than that seen in Voyager images taken in forward scattered light, due either to a difference in observing wavelength, and/or solar phase angle or to changes over time. We may have detected the λ ring in one scan at , but other scans provided upper limits below this value. A single detection, however, would be consistent with azimuthal asymmetries known to exist in this ring. We further demonstrate the presence of azimuthal asymmetries in all rings. We confirm the eccentricity of ∼0.001 in rings 4, 5, 6, which in 2004 are ∼70 km closer to Uranus in the north (near periapse; lower I/F) than in the south. We find a global optical depth of τ∼0.3 in the main rings, and of τ=0.25±0.05 in the ε ring.     

New near-aphelion light curves of Comet 2P/Encke

We present new, near-aphelion, time series of photometry of Comet 2P/Encke in Cousins-R band. With these light curves we find that the dominant, synodic rotational periodicity is either P₀=11.079±0.009 h or 2P₀=22.158±0.012 h. This is in contrast to data from the 1980s published by others that are consistent with 15.08- and 22.6-h periods. Those periods do not satisfy our phased light curves, and also the 1980s data are not easily reconciled with our periods. This could be due to P/Encke having non-principal axis rotation or due to a drift in the rotation period caused by outgassing torques. We observed the comet at five epochs: July, August, September, and October 2001, and September 2002, and the comet was at times intrinsically brighter than expected for a bare nucleus, due to an apparent contribution from an unresolved coma. Three-quarters of the data were obtained in the second and fifth epochs, and we analyzed these two time series using both the phase-dispersion minimization and “WindowCLEAN” techniques. At both epochs and with both techniques strong periodicities were found near frequencies and . By then using visual inspection of the phased light curves to corroborate these frequencies, and by using the data from the other three epochs to properly align light curve features, we were able to derive P₀ and 2P₀ as the only solutions that satisfy all our observations. The periodicity due to f₁ is clearly seen in our data, but we cannot tell from our data alone whether it is a manifestation of the nucleus's shape, non-principal axis rotation, or both.     

Fast computation of 2D transfer functions from adaptive optics data

The use of atmospheric transfer functions is common in image reconstruction techniques such as speckle interferometry to calibrate the Fourier amplitudes of the reconstructed images. Thus, an accurate model is needed to ensure proper photometry in the reconstruction. The situation complicates when adaptive optics (AO) are used during data acquisition. I propose a novel technique to derive two‐dimensional transfer functions from data collected using AO simultaneously with the observations. The technique is capable to compute the relevant transfer functions within a short time for the prevailing atmospheric conditions and AO performance during data acquisition (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Red Wing Creek structure,North Dakota: Petrographical and geochemical studies,and confirmation of impact origin

Abstract The 9 km diameter Red Wing Creek structure, North Dakota, is located within the oil-rich Williston Basin at 47°36′N and 103°33′W. Earlier geophysical studies indicated that this subsurface structure has a central uplift, surrounded by an annular crater moat, and a raised rim. Breccias were encountered during drilling between ～2000 and 2800 m depth in the central uplift area, and the presence of shatter cone fragments in drill core samples was suggested to indicate an impact origin of the Red Wing Creek structure. We studied the petrographic and geochemical characteristics of samples of well cuttings from two boreholes at the center of the structure: the True Oil 22–27 Burlington Northern and True Oil 11–27 Burlington Northern wells. We found planar deformation features (PDFs) in quartz with up to three sets of different crystallographic orientations in sandstone- and siltstone-dominated samples from the True Oil 11–27 borehole. U-stage measurements of the crystallographic orientations of the PDFs showed the occurrence of the shock-characteristic (0001), and orientations, with a dominance of (0001) and orientations. The relative frequencies of the orientations indicate a shock pressure of at least 12–20 GPa. These results provide unambiguous evidence for shock metamorphism at Red Wing Creek and confirm that the structure was formed by impact.     

Flare SOL2012-07-06: On the Origin of the Circular Polarization Reversal Between 17 GHz and 34 GHz

The new generation of multiwavelength radioheliographs with high spatial resolution will employ microwave imaging spectropolarimetry to recover flare topology and plasma parameters in the flare sources and along the wave propagation paths. The recorded polarization depends on the emission mechanism and emission regime (optically thick or thin), the emitting particle properties, and propagation effects. Here, we report an unusual flare, SOL2012-07-06T01:37, whose optically thin gyrosynchrotron emission of the main source displays an apparently ordinary mode sense of polarization in contrast to the classical theory that favors the extraordinary mode. This flare produced copious nonthermal emission in hard X-rays and in high-frequency microwaves up to 80 GHz. It is found that the main flare source corresponds to an interaction site of two loops with greatly different sizes. The flare occurred in the central part of the solar disk, which allows reconstructing the magnetic field in the flare region using vector magnetogram data. We have investigated the three possible known reasons of the circular polarization sense reversal – mode coupling, positron contribution, and the effect of beamed angular distribution. We excluded polarization reversal due to contribution of positrons because there was no relevant response in the X-ray emission. We find that a beam-like electron distribution can produce the observed polarization behavior, but the source thermal density must be much higher than the estimate from to the X-ray data. We conclude that the apparent ordinary wave emission in the optically thin mode is due to radio wave propagation across the quasi-transverse (QT) layer. The abnormally high transition frequency (above 35 GHz) can be achieved reasonably low in the corona where the magnetic field value is high and transverse to the line of sight. This places the microwave source below this QT layer, i.e. very low in the corona.     

Detection of C2HD and the D/H ratio on Titan

We report here the first detection of mono-deuterated acetylene (acetylene-d1, C₂HD) in Titan's atmosphere from the presence of two of its emission bands at 678 and 519 cm⁻¹ as observed in CIRS spectral averages of nadir and limb observations taken between July 2004 and mid-2007. By using new laboratory spectra for this molecule, we were able to derive its abundance at different locations over Titan's disk. We find the C₂HD value () to be roughly constant with latitude from the South to about 45° N and then to increase slightly in the North, as is the case for C₂H₂. Fitting the 678 cm⁻¹ν₅ band simultaneously with the nearby C₂H₂ 729 cm⁻¹ν₅ band, allows us to infer a D/H ratio in acetylene on Titan with an average of the modal values of 2.09±0.45×¹⁰⁻⁴ from the nadir observations, the uncertainties being mainly due to the vertical profile used for the fit of the acetylene band. Although still subject to significant uncertainty, this D/H ratio appears to be significantly larger than the one derived in methane from the CH₃D band (upper limit of 1.5×¹⁰⁻⁴; Bézard, B., Nixon, C.A., Kleiner, I., Jennings, D.E., 2007. Icarus, 191, 397-400; Coustenis, A., Achterberg, R., Conrath, B., Jennings, D., Marten, A., Gautier, D., Bjoraker, G., Nixon, C., Romani, P., Carlson, R., Flasar, M., Samuelson, R.E., Teanby, N., Irwin, P., Bézard, B., Orton, G., Kunde, V., Abbas, M., Courtin, R., Fouchet, Th., Hubert, A., Lellouch, E., Mondellini, J., Taylor, F.W., Vinatier, S., 2007. Icarus 189, 35-62). From the analysis of limb data we infer D/H values of (at 54° S), (at 15° S), (at 54° N) and (at 80° N), which average to a mean value of 1.63±0.27×¹⁰⁻⁴.     

Reconstruction of the surface topography from single images with the photometric method

The possibility of reconstructing the surface topography from single images with the photometric method in the linear approximation is analyzed. The photometric method or surface topography Reconstruction employs a statistical approach to the problem formulation and is the most mathematically correct. This method allows determination of the most probable surface topography given specific observational data. When only one image is available, the photometric method is superior in comparison with the currently available photoclinometry. The processing of test surface topography with the photometric method shows that, under typical conditions, the error of surface relief reconstruction is of higher than 40% in terms of the standard deviation of the surface height. The surface relief of some Martian areas are reconstructed from HRSC images obtained by the Mars Express spacecraft. It is shown that the image-reconstructed surface topography is in good agreement with the topographic information for the same Martian areas obtained by the MOLA altimeter.     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.     

Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Saturn A ring surface mass densities from spiral density wave dispersion behavior

We have undertaken an analysis of the Voyager photopolarimeter (PPS) stellar occultation data of Saturn's A ring. The Voyager PPS observed the bright star δ Scorpii as it was occulted by Saturn's main rings during the spacecraft flyby of the Saturn system in 1981. The occultation measurement produced a ring profile with radial resolution of approximately 100 m, and radial structure is evident in the profile down to the resolution limit. We have applied an autoregressive technique to the data for estimating the power spectrum as a function of radius, which has allowed us to identify 40 spiral density waves in Saturn's A ring, associated with the strongest torques due to forcing from the moons. The majority of the detected waves are observed to disperse linearly in regions beginning a few kilometers from the resonance location. We have used the dispersion behavior for those waves to calculate local surface mass densities in the vicinity of each wave. We find that the inner three-quarters of the A ring (up to the beginning of the Encke gap) has an average surface mass density of , while the outer region has an average surface mass density of . The two regions have different mean surface mass densities with a significance of approximately 0.999993, as estimated with a T-statistic, which corresponds to about 4.5σ. While the mean optical depth of the A ring increases slightly with increasing distance from Saturn, we find that it is not significantly correlated with the surface mass density; the two quantities having a linear Pearson's correlation coefficient of r_corr≈−0.03. The variation of mass density, independent of PPS optical depth, is consistent with previous conjectures that the particle size distribution and composition are not constant across the entire A ring, particularly in the very outer portion. We estimate the mass of Saturn's A ring from our surface mass density estimates as 4.9×10²¹ gm, or 8.61×10⁻⁹ of the mass of Saturn, roughly equivalent to the mass of a 110-km diameter icy satellite. This mass is almost 25% smaller than estimates from previous studies, but is well within the expected errors of the derived mass densities. We also identified three previously unstudied features which exhibit linear dispersion. The strongest of these features is tentatively identified as the Janus 13:11 density wave. The other two features do not fall near any known satellite resonances and may represent density waves created by previously undetected satellites.     

Model-independent restoration of the brightness distribution across the disks of stars from high precision eclipsing light curves

Estimates are made of the accuracy with which the brightness distributions across the disks of stars can be reconstructed through analysis of data from high precision space-based photometry of classical eclipsing systems and observations of the transit of planets across stellar disks. The ill-posed reconstruction problem was solved on a compact set of monotonically nonincreasing, upwardly convex, non-negative functions. One of the difficulties with this method in the case of stars with thin photospheres is the poor convergence of the solution at the point where the brightness distribution has a discontinuity at the edge of the star's disk. Nevertheless, the use of this method for analysis of high precision observational data is justified, since it can be used to obtain an estimate of the limb darkening that is independent of any model assumptions. The reconstructed brightness distribution for the star HD 209458, for which the transit of a planet over its disk was observed with the HST space telescope, is in good agreement with the results of a nonlinear model fit. Translated from Astrofizika, Vol. 51, No. 4, pp. 595–606 (November 2008).     

Generalized projection theorem for mapping with radio interferometers

A theorem, which provides a relationship between the one-dimensional Fourier transform of a line section across two-dimensional (spatial) spectrum and that across the correspoding (brightness) distribution function, is proved. The theorem is then shown to be relevant in some problems in radio astronomy and possibly in other fields connected with image reconstructions from one-dimensional scans through objects or their spectra.     

Combining Gaia Windows: The Method and Simulations

During the 5-year-observational phase of Gaia several one-dimensional observations will be obtained from the same objects at different epochs and transit angles. Combining the information from all the transits to produce a single two-dimensional image is scientifically important as this will allow to detect much fainter stars in the vicinity of the main source than otherwise possible from single observations on-board the satellite or from individual transmitted patches. Most of the emerging objects will be physical companions, but some background/foreground objects will also be found. The combination procedure further increases the image resolution in the across scan-direction. The procedure consists of combining patches that are stretched back to 2-D windows, only possible since the centroid positions of the windows are known. The initial image window from which the patches originate is not square, but elongated along the scan-direction if the Astro AF11 field of view is used in the current calculations. Transit angles can take any value between 0^∘ and 360^∘. The window combining method is a simple mapping procedure that is fast and should handle all realistic complexities occurring during the Gaia data reduction. The method is not intended to be a perfect image restoration technique, but its purpose is to fulfill the scientific requirement of detecting faint neighboring sources existing in the close vicinity of the main targets. In this study we show examples of the combined windows and illustrate how faint companions around primary targets could be detected. We investigate the performances and limitations of the method and discuss possible complexities that might occur during the final data analysis.