Refine
Document Type
- Article (14)
- Doctoral Thesis (1)
Language
- English (15)
Is part of the Bibliography
- yes (15)
Keywords
- methods: data analysis (15) (remove)
Institute
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
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.
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.
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.
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.
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.
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.
Most of the baryonic matter in the Universe resides in a diffuse gaseous phase in-between galaxies consisting mostly of hydrogen and helium. This intergalactic medium (IGM) is distributed in large-scale filaments as part of the overall cosmic web. The luminous extragalactic objects that we can observe today, such as galaxies and quasars, are surrounded by the IGM in the most dense regions within the cosmic web. The radiation of these objects contributes to the so-called ultraviolet background (UVB) which keeps the IGM highly ionized ever since the epoch of reionization.
Measuring the amount of absorption due to intergalactic neutral hydrogen (HI) against extragalactic background sources is a very useful tool to constrain the energy input of ionizing sources into the IGM. Observations suggest that the HI Lyman-alpha effective optical depth, τ_eff, decreases with decreasing redshift, which is primarily due to the expansion of the Universe. However, some studies find a smaller value of the effective optical depth than expected at the specific redshift z~3.2, possibly related to the complete reionization of helium in the IGM and a hardening of the UVB. The detection and possible cause of a decrease in τ_eff at z~3.2 is controversially debated in the literature and the observed features need further explanation.
To better understand the properties of the mean absorption at high redshift and to provide an answer for whether the detection of a τ_eff feature is real we study 13 high-resolution, high signal-to-noise ratio quasar spectra observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) at the Very Large Telescope (VLT). The redshift evolution of the effective optical depth, τ_eff(z), is measured in the redshift range 2.7≤z≤3.6. The influence of metal absorption features is removed by performing a comprehensive absorption-line-fitting procedure.
In the first part of the thesis, a line-parameter analysis of the column density, N, and Doppler parameter, b, of ≈7500 individually fitted absorption lines is performed. The results are in good agreement with findings from previous surveys.
The second (main) part of this thesis deals with the analysis of the redshift evolution of the effective optical depth. The τ_eff measurements vary around the empirical power law τ_eff(z)~(1+z)^(γ+1) with γ=2.09±0.52. The same analysis as for the observed spectra is performed on synthetic absorption spectra. From a comparison between observed and synthetic spectral data it can be inferred that the uncertainties of the τ_eff values are likely underestimated and that the scatter is probably caused by high-column-density absorbers with column densities in the range 15≤logN≤17. In the real Universe, such absorbers are rarely observed, however. Hence, the difference in τ_eff from different observational data sets and absorption studies is most likely caused by cosmic variance. If, alternatively, the disagreement between such data is a result of an too optimistic estimate of the (systematic) errors, it is also possible that all τ_eff measurements agree with a smooth evolution within the investigated redshift range. To explore in detail the different analysis techniques of previous studies an extensive literature comparison to the results of this work is presented in this thesis.
Although a final explanation for the occurrence of the τ_eff deviation in different studies at z~3.2 cannot be given here, our study, which represents the most detailed line-fitting analysis of its kind performed at the investigated redshifts so far, represents another important benchmark for the characterization of the HI Ly-alpha effective optical depth at high redshift and its indicated unusual behavior at z~3.2.
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.
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