Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a St\"ackel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1 and Pal 5 streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 3.47^{+0.95}_{-1.44} \times 10^{11} \ M_{\odot}$. When including the Orphan and Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+5.69}_{-1.07} \times 10^{11} \ M_{\odot}$, indicative of the hidden systematic errors in fits to individual streams. These systematics are likely due to insufficiencies in our St\"ackel model of the potential and to the limited phase space explored by the data set of 4 streams, so the larger uncertainty of the combined fit is the more robust measure of the actual uncertainty in the Milky Way potential.

We present a study of the membership of the Hyades open cluster, derive kinematically-modelled parallaxes of its members, and study the colour-absolute magnitude diagram of the cluster. We use Gaia DR1 Tycho-Gaia Astrometric Solution (TGAS) data complemented by Hipparcos-2 data for bright stars not contained in TGAS. We supplement the astrometric data with radial velocities collected from a dozen literature sources. By assuming that all cluster members move with the mean cluster velocity to within the velocity dispersion, we use the observed and the expected motions of the stars to determine individual cluster membership probabilities. We subsequently derive improved parallaxes through maximum-likelihood kinematic modelling of the cluster. This method has an iterative component to deal with 'outliers', caused for instance by double stars or escaping members. Our method extends an existing method and supports the mixed presence of stars with and without radial velocities. We find 251 candidate members, 200 of which have a literature radial velocity, and 70 of which are new candidate members with TGAS astrometry. The cluster is roughly spherical in its centre but significantly flattened at larger radii. The observed colour-absolute magnitude diagram shows a clear binary sequence. The kinematically-modelled parallaxes that we derive are a factor ~1.7 / 2.9 more precise than the TGAS / Hipparcos-2 values and allow to derive an extremely sharp main sequence. This sequence shows evidence for fine-detailed structure which is elegantly explained by the full spectrum turbulence model of convection.