The Fourth IPCC Assessment Report concluded that ice sheet flow models, in their current state, were unable to provide accurate forecast for the increase of polar ice sheet discharge and the ...associated contribution to sea level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice flow models, and as a result a significant number of new ice sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger-scale problems, earning the status of an ice sheet model. Here, we summarise almost 10 yr of development performed by different groups. Elmer/Ice solves the full-Stokes equations, for isotropic but also anisotropic ice rheology, resolves the grounding line dynamics as a contact problem, and contains various basal friction laws. Derived fields, like the age of the ice, the strain rate or stress, can also be computed. Elmer/Ice includes two recently proposed inverse methods to infer badly known parameters. Elmer is a highly parallelised code thanks to recent developments and the implementation of a block preconditioned solver for the Stokes system. In this paper, all these components are presented in detail, as well as the numerical performance of the Stokes solver and developments planned for the future.
Aims.In this paper, we present an innovative data reduction method for single-mode interferometry. It has been specifically developed for the AMBER instrument, the three-beam combiner of the Very ...Large Telescope Interferometer, but it can be derived for any single-mode interferometer. Methods.The algorithm is based on a direct modelling of the fringes in the detector plane. As such, it requires a preliminary calibration of the instrument in order to obtain the calibration matrix that builds the linear relationship between the interferogram and the interferometric observable, which is the complex visibility. Once the calibration procedure has been performed, the signal processing appears to be a classical least-square determination of a linear inverse problem. From the estimated complex visibility, we derive the squared visibility, the closure phase, and the spectral differential phase. Results.The data reduction procedures have been gathered into the so-called amdlib software, now available for the community, and are presented in this paper. Furthermore, each step in this original algorithm is illustrated and discussed from various on-sky observations conducted with the VLTI, with a focus on the control of the data quality and the effective execution of the data reduction procedures. We point out the present limited performances of the instrument due to VLTI instrumental vibrations which are difficult to calibrate.
Aims. We present the first NIR spectro-interferometry of the LBV η Carinae. The observations were performed with the AMBER instrument of the ESO Very Large Telescope Interferometer (VLTI) using ...baselines from 42 to 89 m. The aim of this work is to study the wavelength dependence of η Car's optically thick wind region with a high spatial resolution of 5 mas (11 AU) and high spectral resolution. Methods. The observations were carried out with three 8.2 m Unit Telescopes in the K-band. The raw data are spectrally dispersed interferograms obtained with spectral resolutions of 1500 (MR-K mode) and 12 000 (HR-K mode). The MR-K observations were performed in the wavelength range around both the $\ion{He}{i}$ 2.059 μm and the Brγ 2.166 μm emission lines, the HR-K observations only in the Brγ line region. Results. The spectrally dispersed AMBER interferograms allow the investigation of the wavelength dependence of the visibility, differential phase, and closure phase of η Car. In the K-band continuum, a diameter of $4.0\pm0.2$ mas (Gaussian FWHM, fit range 28–89 m baseline length) was measured for η Car's optically thick wind region. If we fit Hillier et al. (2001, ApJ, 553, 837) model visibilities to the observed AMBER visibilities, we obtain 50% encircled-energy diameters of 4.2, 6.5 and 9.6 mas in the 2.17$\,\mu$m continuum, the $\ion{He}{i}$, and the Brγ emission lines, respectively. In the continuum near the Brγ line, an elongation along a position angle of $120\degr\pm15\degr$ was found, consistent with previous VINCI/VLTI measurements by van Boekel et al. (2003, A&A, 410, L37). We compare the measured visibilities with predictions of the radiative transfer model of Hillier et al. (2001), finding good agreement. Furthermore, we discuss the detectability of the hypothetical hot binary companion. For the interpretation of the non-zero differential and closure phases measured within the Brγ line, we present a simple geometric model of an inclined, latitude-dependent wind zone. Our observations support theoretical models of anisotropic winds from fast-rotating, luminous hot stars with enhanced high-velocity mass loss near the polar regions.
The young stellar object MWC 297 is an embedded B1.5Ve star exhibiting strong hydrogen emission lines and a strong near-infrared continuum excess. This object has been observed with the VLT ...interferometer equipped with the AMBER instrument during its first commissioning run. AMBER/VLTI is currently the only near infrared interferometer that can observe spectrally dispersed visibilities. MWC 297 has been spatially resolved in the continuum with a visibility of $0.50^{+0.08}_{-0.10}$ as well as in the Brγ emission line where the visibility decreases to $0.33\pm0.06$. This change in the visibility with wavelength can be interpreted by the presence of an optically thick disk responsible for the visibility in the continuum and of a stellar wind traced by the Brγ emission line and whose apparent size is 40% larger. We validate this interpretation by building a model of the stellar environment that combines a geometrically thin, optically thick accretion disk model consisting of gas and dust, and a latitude-dependent stellar wind outflowing above the disk surface. The continuum emission and visibilities obtained from this model are fully consistent with the interferometric AMBER data. They agree also with existing optical, near-infrared spectra and other broad-band near-infrared interferometric visibilities. We also reproduce the shape of the visibilities in the Brγ line as well as the profile of this line obtained at an higher spectral resolution with the VLT/ISAAC spectrograph, and those of the Hα and Hβ lines. The disk and wind models yield a consistent inclination of the system of approximately 20°. A picture emerges in which MWC 297 is surrounded by an equatorial flat disk that is possibly still accreting and an outflowing wind that has a much higher velocity in the polar region than at the equator. The AMBER/VLTI unique capability of measuring spectral visibilities therefore allows us for the first time to compare the apparent geometry of a wind with the disk structure in a young stellar system.
Main modes of variability of the Antarctic tropospheric circulation (500 hPa geopotential height) and precipitation are identified through their empirical orthogonal functions (EOF). This is done by ...combining various sources of information, including meteorological analyses and forecasts (NCEP and ECMWF), atmospheric general circulation model (LMDZ) simulations, and satellite data (GPCP). Unlike previous similar work on circulation variability, the mode analyses are restricted to the Antarctic region. The main modes that relate the Antarctic region to the mid and tropical latitudes, e.g. in association with ENSO, are nonetheless clearly identified and thus robust. The contribution of the sea-surface or of the circumpolar Antarctic atmospheric dynamics to the occurrence and to the chronology of these modes is evaluated through various atmospheric model simulations. EOF analyses results are somewhat less stable, across the various datasets, and more noisy for precipitation than for circulation. Yet, through moisture advection considerations, the two most significant precipitation modes can be well related to the three main modes of circulation variability. The signatures of both the Southern Oscillation Index (SOI) and the Antarctic Oscillation Index (AOI) are found in one same precipitation mode, suggesting that they have a substantially common spatial structure. In addition, the relative strength of the signature of the AOI and SOI appears to change in time. In particular, the signature of the SOI was weak in the 1980s precipitations, but turned very strong in the 1990s. Common spatial patterns and variable strength in time may explain why hints of an ENSO signature in Antarctic precipitation have been reported but not unequivocally demonstrated so far.PUBLICATION ABSTRACT
The grounding line is a key element of coastal outlet glaciers, acting on their dynamics. Accurately knowing its position is fundamental for both modelling the glacier dynamics and establishing a ...benchmark for later change detection. Here we map the grounding line of the Astrolabe Glacier in East Antarctica (66°41' S, 140°05' E), using both hydrostatic and tidal methods. The first method is based on new surface and ice thickness data from which the line of buoyant floatation is found. The second method uses kinematic GPS measurements of the tidal response of the ice surface. By detecting the transitions where the ice starts to move vertically in response to the tidal forcing we determine control points for the grounding line position along GPS profiles. Employing a two-dimensional elastic plate model, we compute the rigid short-term behaviour of the ice plate and estimate the correction required to compare the kinematic GPS control points with the previously determined line of floatation. These two approaches show consistency and lead us to propose a grounding line for the Astrolabe Glacier that significantly deviates from the lines obtained so far from satellite imagery.
Aims.We study the geometry and kinematics of the circumstellar environment of the Be star κ CMa in the Brγ emission line and its nearby continuum. Methods.We use the AMBER/VLTI instrument operating ...in the K band, which provides a spatial resolution of about 6 mas with a spectral resolution of 1500, to study the kinematics within the disk and to infer its rotation law. To obtain more kinematical constraints we also use a high spectral resolution Paβ line profile obtain in December 2005 at the Observatorio do Pico do Dios, Brazil and we compile $V/R$ line profile variations and spectral energy distribution data points from the literature. Results.Using differential visibilities and differential phases across the Brγ line we detect an asymmetry in the disk. Moreover, we found that κ CMa seems difficult to fit within the classical scenario for Be stars, illustrated recently by α Arae observations, i.e. a fast rotating B star close to its breakup velocity surrounded by a Keplerian circumstellar disk with an enhanced polar wind. We discuss the possibility that κ CMa is a critical rotator with a Keplerian rotating disk and examine whether if the detected asymmetry can be interpreted within the “one-armed” viscous disk framework.
Aims.This paper describes the design goals and engineering efforts that led to the realization of AMBER (Astronomical Multi BEam combineR) and to the achievement of its present performance. ...Methods.On the basis of the general instrumental concept, AMBER was decomposed into modules whose functions and detailed characteristics are given. Emphasis is put on the spatial filtering system, a key element of the instrument. We established a budget for transmission and contrast degradation through the different modules, and made the detailed optical design. The latter confirmed the overall performance of the instrument and defined the exact implementation of the AMBER optics. Results.The performance was assessed with laboratory measurements and commissionings at the VLTI, in terms of spectral coverage and resolution, instrumental contrast higher than 0.80, minimum magnitude of 11 in K, absolute visibility accuracy of 1%, and differential phase stability of 10-3 rad over one minute.
Context.Optical long-baseline interferometry is moving a crucial step forward with the advent of general-user scientific instruments that equip large aperture and hectometric baseline facilities, ...such as the Very Large Telescope Interferometer (VLTI). Aims.AMBER is one of the VLTI instruments that combines up to three beams with low, moderate and high spectral resolutions in order to provide milli-arcsecond spatial resolution for compact astrophysical sources in the near-infrared wavelength domain. Its main specifications are based on three key programs on young stellar objects, active galactic nuclei central regions, masses, and spectra of hot extra-solar planets. Methods.These key science goals led to scientific specifications, which were used to propose and then validate the instrument concept. AMBER uses single-mode fibers to filter the entrance signal and to reach highly accurate, multiaxial three-beam combination, yielding three baselines and a closure phase, three spectral dispersive elements, and specific self-calibration procedures. Results.The AMBER measurements yield spectrally dispersed calibrated visibilities, color-differential complex visibilities, and a closure phase allows astronomers to contemplate rudimentary imaging and highly accurate visibility and phase differential measurements. AMBER was installed in 2004 at the Paranal Observatory. We describe here the present implementation of the instrument in the configuration with which the astronomical community can access it. Conclusions.After two years of commissioning tests and preliminary observations, AMBER has produced its first refereed publications, allowing assessment of its scientific potential.
Aims.We investigate the origin of the ${\rm Br}\gamma$ emission of the Herbig Ae star HD 104237 on Astronomical Unit (AU) scales. Methods.Using AMBER/VLTI at a spectral resolution $\mathcal{R}=1500$ ...we spatially resolve the emission in both the ${\rm Br}\gamma$ line and the adjacent continuum. Results.The visibility does not vary between the continuum and the ${\rm Br}\gamma$ line, even though the line is strongly detected in the spectrum, with a peak intensity 35% above the continuum. This demonstrates that the line and continuum emission have similar size scales. We assume that the K-band continuum excess originates in a “puffed-up” inner rim of the circumstellar disk, and discuss the likely origin of ${\rm Br}\gamma$. Conclusions.We conclude that this emission most likely arises from a compact disk wind, launched from a region 0.2–0.5 AU from the star, with a spatial extent similar to that of the near infrared continuum emission region, i.e., very close to the inner rim location.
PDF
