Refine
Has Fulltext
- no (14)
Document Type
- Article (14) (remove)
Language
- English (14) (remove)
Is part of the Bibliography
- yes (14)
Keywords
- methods: data analysis (14) (remove)
Institute
Thermophysical modelling and parameter estimation of small solar system bodies via data assimilation
(2020)
Deriving thermophysical properties such as thermal inertia from thermal infrared observations provides useful insights into the structure of the surface material on planetary bodies. The estimation of these properties is usually done by fitting temperature variations calculated by thermophysical models to infrared observations. For multiple free model parameters, traditional methods such as least-squares fitting or Markov chain Monte Carlo methods become computationally too expensive. Consequently, the simultaneous estimation of several thermophysical parameters, together with their corresponding uncertainties and correlations, is often not computationally feasible and the analysis is usually reduced to fitting one or two parameters. Data assimilation (DA) methods have been shown to be robust while sufficiently accurate and computationally affordable even for a large number of parameters. This paper will introduce a standard sequential DA method, the ensemble square root filter, for thermophysical modelling of asteroid surfaces. This method is used to re-analyse infrared observations of the MARA instrument, which measured the diurnal temperature variation of a single boulder on the surface of near-Earth asteroid (162173) Ryugu. The thermal inertia is estimated to be 295 +/- 18 Jm(-2) K-1 s(-1/2), while all five free parameters of the initial analysis are varied and estimated simultaneously. Based on this thermal inertia estimate the thermal conductivity of the boulder is estimated to be between 0.07 and 0.12,Wm(-1) K-1 and the porosity to be between 0.30 and 0.52. For the first time in thermophysical parameter derivation, correlations and uncertainties of all free model parameters are incorporated in the estimation procedure that is more than 5000 times more efficient than a comparable parameter sweep.
Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature. Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays). Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
In high-resolution solar physics, the volume and complexity of photometric, spectroscopic, and polarimetric ground-based data significantly increased in the last decade, reaching data acquisition rates of terabytes per hour. This is driven by the desire to capture fast processes on the Sun and the necessity for short exposure times "freezing" the atmospheric seeing, thus enabling ex post facto image restoration. Consequently, large-format and high-cadence detectors are nowadays used in solar observations to facilitate image restoration. Based on our experience during the "early science" phase with the 1.5 m GREGOR solar telescope (2014–2015) and the subsequent transition to routine observations in 2016, we describe data collection and data management tailored toward image restoration and imaging spectroscopy. We outline our approaches regarding data processing, analysis, and archiving for two of GREGOR's post-focus instruments (see http://gregor.aip.de), i.e., the GREGOR Fabry–Pérot Interferometer (GFPI) and the newly installed High-Resolution Fast Imager (HiFI). The heterogeneous and complex nature of multidimensional data arising from high-resolution solar observations provides an intriguing but also a challenging example for "big data" in astronomy. The big data challenge has two aspects: (1) establishing a workflow for publishing the data for the whole community and beyond and (2) creating a collaborative research environment (CRE), where computationally intense data and postprocessing tools are colocated and collaborative work is enabled for scientists of multiple institutes. This requires either collaboration with a data center or frameworks and databases capable of dealing with huge data sets based on virtual observatory (VO) and other community standards and procedures.
Aims. We study the evolution of an arch filament system (AFS) and of its individual arch filaments to learn about the processes occurring in them. Methods. We observed the AFS at the GREGOR solar telescope on Tenerife at high cadence with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) in the He I 10 830 angstrom spectral range. The He I triplet profiles were fitted with analytic functions to infer line-of-sight (LOS) velocities to follow plasma motions within the AFS. Results. We tracked the temporal evolution of an individual arch filament over its entire lifetime, as seen in the He I 10 830 angstrom triplet. The arch filament expanded in height and extended in length from 13 ' to 21 '. The lifetime of this arch filament is about 30 min. About 11 min after the arch filament is seen in He I, the loop top starts to rise with an average Doppler velocity of 6 km s(-1). Only two minutes later, plasma drains down with supersonic velocities towards the footpoints reaching a peak velocity of up to 40 km s(-1) in the chromosphere. The temporal evolution of He I 10 830 angstrom profiles near the leading pore showed almost ubiquitous dual red components of the He I triplet, indicating strong downflows, along with material nearly at rest within the same resolution element during the whole observing time.
The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.
Context. On 27 April 2015, when comet 67P/Churyumov-Gerasimenko was at 1.76 au from the Sun and moving toward perihelion, the OSIRIS and VIRTIS-M instruments on board the Rosetta spacecraft simultaneously observed the evolving dust and gas coma during a complete rotation of the comet. Aims. We aim to characterize the spatial distribution of dust, H2O, and CO2 gas in the inner coma. To do this, we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the rate we calculated using a thermophysical model. Methods. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O emission band), and 4200 nm (CO2 emission band). We measured the average signal in a circular annulus to study the spatial variation around the comet, and in a sector of the annulus to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet center. Results. The spatial correlation between dust and water, both coming from the sunlit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0 degrees subsolar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production with respect to the value calculated using a simple thermophysical model is observed when the head lobe and regions of the southern hemisphere with strong seasonal variations are illuminated (subsolar longitude 270 degrees-50 degrees). A drastic decrease in dust production when the water production (both measured and from the model) displays a maximum occurs when typical northern consolidated regions are illuminated and the southern hemisphere regions with strong seasonal variations are instead in shadow (subsolar longitude 50 degrees-90 degrees). Possible explanations of these observations are presented and discussed.
Aims. We demonstrate the high multiplex advantage of crowded field 3D spectroscopy with the new integral field spectrograph MUSE by means of a spectroscopic analysis of more than 12 000 individual stars in the globular cluster NGC 6397. Methods. The stars are deblended with a point spread function fitting technique, using a photometric reference catalogue from HST as prior, including relative positions and brightnesses. This catalogue is also used for a first analysis of the extracted spectra, followed by an automatic in-depth analysis via a full-spectrum fitting method based on a large grid of PHOENIX spectra. Results. We analysed the largest sample so far available for a single globular cluster of 18 932 spectra from 12 307 stars in NGC 6397. We derived a mean radial velocity of v(rad) = 17.84 +/- 0.07 km s(-1) and a mean metallicity of [Fe/H] = -2.120 +/- 0.002, with the latter seemingly varying with temperature for stars on the red giant branch (RGB). We determine Teff and [Fe/H] from the spectra, and log g from HST photometry. This is the first very comprehensive Hertzsprung-Russell diagram (HRD) for a globular cluster based on the analysis of several thousands of stellar spectra, ranging from the main sequence to the tip of the RGB. Furthermore, two interesting objects were identified; one is a post-AGB star and the other is a possible millisecond-pulsar companion.
We present K-s-band light curves for 299 Cepheids in the Small Magellanic Cloud (SMC) of which 288 are new discoveries that we have identified using multi-epoch near-infrared photometry obtained by the VISTA survey of the Magellanic Clouds system (VMC). The new Cepheids have periods in the range from 0.34 to 9.1 d and cover the magnitude interval 12.9 <= currency sign < K-s > <= currency sign 17.6 mag. Our method was developed using variable stars previously identified by the optical microlensing survey OGLE. We focus on searching new Cepheids in external regions of the SMC for which complete VMC K-s-band observations are available and no comprehensive identification of different types of variable stars from other surveys exists yet.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
Improved measurements of the photospheric and chromospheric three-dimensional magnetic and flow fields are crucial for a precise determination of the origin and evolution of active regions. We present an illustrative sample of multi-instrument data acquired during a two-week coordinated observing campaign in August 2015 involving, among others, the GREGOR solar telescope (imaging and near-infrared spectroscopy) and the space missions Solar Dynamics Observatory (SDO) and Interface Region Imaging Spectrograph (IRIS). The observations focused on the trailing part of active region NOAA 12396 with complex polarity inversion lines and strong intrusions of opposite polarity flux. The GREGOR Infrared Spectrograph (GRIS) provided Stokes IQUV spectral profiles in the photospheric Si i.1082.7 nm line, the chromospheric He I lambda 1083.0 nm triplet, and the photospheric Ca I lambda 1083.9 nm line. Carefully calibrated GRIS scans of the active region provided maps of Doppler velocity and magnetic field at different atmospheric heights. We compare quick-look maps with those obtained with the " Stokes Inversions based on Response functions" (SIR) code, which furnishes deeper insight into the magnetic properties of the region. We find supporting evidence that newly emerging flux and intruding opposite polarity flux are hampering the formation of penumbrae, i.e., a penumbra fully surrounding a sunspot is only expected after cessation of flux emergence in proximity to the sunspots. (C) 2016 WILEY-VCH Verlag GmbH& Co.KGaA, Weinheim
The new generation of solar instruments provides better spectral, spatial, and temporal resolution for a better understanding of the physical processes that take place on the Sun. Multiple-component profiles are more commonly observed with these instruments. Particularly, the He i 10830 triplet presents such peculiar spectral profiles, which give information on the velocity and magnetic fine structure of the upper chromosphere. The purpose of this investigation is to describe a technique to efficiently fit the two blended components of the He i 10830 triplet, which are commonly observed when two atmospheric components are located within the same resolution element. The observations used in this study were taken on 2015 April 17 with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) attached to the 1.5-m GREGOR solar telescope, located at the Observatorio del Teide, Tenerife, Spain. We apply a double-Lorentzian fitting technique using Levenberg-Marquardt least-squares minimization. This technique is very simple and much faster than inversion codes. Line-of-sight Doppler velocities can be inferred for a whole map of pixels within just a few minutes. Our results show sub-and supersonic downflow velocities of up to 32 km s(-1) for the fast component in the vicinity of footpoints of filamentary structures. The slow component presents velocities close to rest. (C) 2016 WILEY-VCH Verlag GmbH& Co. KGaA, Weinheim
Context. The solar magnetic field is responsible for all aspects of solar activity. Thus, emergence of magnetic flux at the surface is the first manifestation of the ensuing solar activity. Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. The small active region NOAA 12118 emerged on 2014 July 17 and was observed one day later with the 1.5-m GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-Perot Interferometer (GFPI) were complemented by synoptic line-of-sight magnetograms and continuum images obtained with the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator (DAVE), respectively. Morphological image processing was employed to measure the photometric and magnetic area, magnetic flux, and the separation profile of the emerging flux region during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days provides a comprehensive view of growth and decay. It traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km s(-1) is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of dark pores. We detected many elongated rapidly expanding granules between the two major polarities, with dimensions twice as large as the normal granules.
We present results of the analysis of 70 RR Lyrae stars located in the bar of the Large Magellanic Cloud (LMC). Combining the spectroscopically determined metallicity of these stars from the literature with precise periods from the OGLE III catalog and multi-epoch K-s photometry from the VISTA survey of the Magellanic Clouds system, we derive a new near-infrared period-luminosity-metallicity (PLKsZ) relation for RR Lyrae variables. In order to fit the relation we use a fitting method developed specifically for this study. The zero-point of the relation is estimated two different ways: by assuming the value of the distance to the LMC and by using Hubble Space Telescope parallaxes of five RR Lyrae stars in the Milky Way (MW). The difference in distance moduli derived by applying these two approaches is similar to 0.2 mag. To investigate this point further we derive the PL(Ks)Z relation based on 23 MW RR Lyrae stars that had been analyzed in Baade-Wesselink studies. We compared the derived PL(Ks)Z relations for RR Lyrae stars in the MW and LMC. Slopes and zero-points are different, but still consistent within the errors. The shallow slope of the metallicity term is confirmed by both LMC and MW variables. The astrometric space mission Gaia is expected to provide a huge contribution to the determination of the RR Lyrae PL(Ks)Z relation; however, calculating an absolute magnitude from the trigonometric parallax of each star and fitting a PL(Ks)Z relation directly to period and absolute magnitude leads to biased results. We present a tool to achieve an unbiased solution by modeling the data and inferring the slope and zero-point of the relation via statistical methods.
The Galactic center is an interesting region for high-energy (0.1-100 GeV) and very-high-energy (E > 100 GeV) gamma-ray observations. Potential sources of GeV/TeV gamma-ray emission have been suggested, e.g., the accretion of matter onto the supermassive black hole, cosmic rays from a nearby supernova remnant (e.g., Sgr A East), particle acceleration in a plerion, or the annihilation of dark matter particles. The Galactic center has been detected by EGRET and by Fermi/LAT in the MeV/GeV energy band. At TeV energies, the Galactic center was detected with moderate significance by the CANGAROO and Whipple 10 m telescopes and with high significance by H.E.S.S., MAGIC, and VERITAS. We present the results from three years of VERITAS observations conducted at large zenith angles resulting in a detection of the Galactic center on the level of 18 standard deviations at energies above similar to 2.5 TeV. The energy spectrum is derived and is found to be compatible with hadronic, leptonic, and hybrid emission models discussed in the literature. Future, more detailed measurements of the high-energy cutoff and better constraints on the high-energy flux variability will help to refine and/or disentangle the individual models.