In this paper we describe the first combined observations executed by the Mars Express and Trace Gas Orbiter missions since the start of the TGO operational phase in April 2018 until June 2019. Also included is the science opportunity analysis that has been performed by the Science Operation Centres and instrument teams to identify the observation opportunities until the end of 2020. These results provide a valuable contribution to the scientific community by enabling collaborations within the instrument teams and enhance the scientific outcome of both missions. This information is also valuable to other Mars missions, which may be interested in observing these locations for wider scientific collaboration.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

Doppler tomography is a well-known method in astrophysics to image the accretion flow, often in the shape of thin discs, in compact binary stars. As accretion discs rotate, all emitted line radiation is Doppler-shifted. In fast-ion D-alpha (FIDA) spectroscopy measurements in magnetically confined plasma, the D-alpha-photons are likewise Doppler-shifted ultimately due to gyration of the fast ions. In either case, spectra of Doppler-shifted line emission are sensitive to the velocity distribution of the emitters. Astrophysical Doppler tomography has lead to images of accretion discs of binaries revealing bright spots, spiral structures, and flow patterns. Fusion plasma Doppler tomography has lead to an image of the fast-ion velocity distribution function in the tokamak ASDEX Upgrade. This image matched numerical simulations very well. Here we discuss achievements of the Doppler tomography approach, its promise and limits, analogies and differences in astrophysical and fusion plasma Doppler tomography, and what can be learned by comparison of these applications.

In a previous paper we investigated means for constraining the mass distribution of cluster galaxies by weak lensing. We concluded that a comprehensive method should treat the lensing effects of individual cluster galaxies and those resulting from a general cluster component simultaneously. To this end we now develop a non-parametric maximum likelihood cluster reconstruction algorithm that can implicitly take into account the presence of cluster galaxies. The method includes an entropy-like regularization prescription and directly uses the ellipticities of individual source galaxy images as observables rather than relying on an averaged ellipticity field. The mass distribution of cluster galaxies is described by parametrized models. For each set of galaxy parameters the cluster reconstruction algorithm allows to determine the best representation of the global underlying cluster component that is consistent with the presence of the cluster galaxies and the observed image ellipticities of background galaxies. Tests with simulations yielded convincing and robust results. We applied the method to a WFPC2 image of the cluster Cl0939+4713 and obtained a detection of the lensing effects of luminous elliptical cluster galaxies. We consider this application as a successful test of our technique. However, the small size of the image we analysed does not yet allow to draw strong conclusions regarding the mass distribution of cluster galaxies.

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the average mass distribution of cluster galaxies and potentially to test the extent of their dark matter haloes as a function of the density of their environment. The observable image distortions are a consequence of the interplay between the effects of a global cluster mass distribution and the perturbations due to individual cluster galaxies. Starting from a reconstruction of the cluster mass distribution with conventional techniques, we apply a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies. From simulations this approach is found to be reliable as long as the galaxies including their dark matter haloes only contribute a small fraction to the total mass of the system. If their haloes are extended, the galaxies contain a substantial mass fraction. In this case our method is still applicable in the outer regions of clusters, where the surface mass density is low, but yields biased estimates of the parameters describing the mass profiles of the cluster galaxies in the central part of the cluster. In that case it will be necessary to resort to more sophisticated strategies by modelling cluster galaxies and an underlying global mass distribution simultaneously. We conclude that galaxy-galaxy lensing in clusters provides a unique means to probe the presence and extent of dark haloes of cluster galaxies.

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the mass distribution of cluster galaxies and potentially to test the extent of their dark matter halos as a function of the density of the environment. Here I describe simulations of observational data and present a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies.

We propose a new technique to measure the time delay of radio-loud gravitational lens systems, which does not rely on the excessive use of interferometric observations. Instead, the method is based on single-dish flux density monitoring of the (unresolved) lens system's total lightcurve, combined with additional interferometric measurements of the flux density ratio at a few epochs during that monitoring period. The basic idea of the method is to reconstruct the individual image lightcurves from the observed total lightcurve by assuming a range of potential values for the time delay and the magnification ratio of the images. It is then possible to single out the correct reconstruction, and therefore determine the time delay, by checking the consistency of the reconstructed individual lightcurves with the additional interferometric observations. We performed extensive numerical simulations of synthetic lightcurves to investigate the dependence of the performance of this method on various parameters which are involved in the problem. Probably the most promising candidates for applying the method (and also for determining the Hubble constant) are lens systems consisting of multiply imaged compact sources and an Einstein ring, such as B0218+357 from which some of the parameters used for our simulations were adopted.