In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.

The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars and can trap or scatter stars close to orbital resonances. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy, ring-like structures in the disk and correlations between different components of the stellar velocity. These perturbations can also cause 'phase wrapping' signatures in the disk, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations, especially in the vertical direction, and signatures of incomplete phase mixing in the disk. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

The nature and origin of the Galactic warp represent one of the open questions posed by Galactic evolution. Thanks to Gaia high precision absolute astrometry, steps towards the understanding of the warp's dynamical nature can be made. Indeed, proper motions for long-lived stable warp are expected to show measurable trends in the component vertical to the galactic plane. Within this context, we search for the kinematic warp signal in the first Gaia data release (DR1). By analyzing distant spectroscopically-identified OB stars in the Hipparcos subset in Gaia DR1, we find that the kinematic trends cannot be explained by a simple model of a long-lived warp. We therefore discuss possible scenarios for the interpretation of the obtained results. We also present current work in progress to select a larger sample of OB star candidates from the Tycho-Gaia Astrometric Solution (TGAS) subsample in DR1, and delineate the points that we will be addressing in the near future.

Recent ALMA observations indicate that CW Leo, the closest carbon-rich AGB star to the Sun, might have a low-mass stellar companion. We present archival ground-based astrometric measurements of CW Leo obtained within the context of the Torino Parallax Program and with $>6$ yr ($1995 - 2001$) of time baseline. The residuals to a single-star solution show significant curvature, and they are strongly correlated with the well-known $I$-band photometric variations due to stellar pulsations. We describe successfully the astrometry of CW Leo with a variability-induced motion (VIM) + acceleration model. We obtain proper motion and parallax of the center-of-mass of the binary, the former in fair agreement with recent estimates, the latter at the near end of the range of inferred distances based on indirect methods. The VIM + acceleration model results allow us to derive a companion mass in agreement with that inferred by ALMA, they point towards a somewhat longer period than implied by ALMA, but are not compatible with much longer period estimates. These data will constitute a fundamental contribution towards the full understanding of the orbital architecture of the system when combined with Gaia astrometry, providing a $\sim 25$\,yr time baseline.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The mechanism responsible for the warp of our Galaxy, as well as its dynamical nature, continues to remain unknown. With the advent of high precision astrometry, new horizons have been opened for detecting the kinematics associated with the warp and constraining possible warp formation scenarios for the Milky Way. The aim of this contribution is to establish whether the first Gaia data release (DR1) shows significant evidence of the kinematic signature expected from a long-lived Galactic warp in the kinematics of distant OB stars. As the first paper in a series, we present our approach for analyzing the proper motions and apply it to the sub-sample of Hipparcos stars. We select a sample of 989 distant spectroscopically-identified OB stars from the New Reduction of Hipparcos, of which 758 are also in Gaia DR1, covering distances from 0.5 to 3 kpc from the Sun. We develop a model of the spatial distribution and kinematics of the OB stars from which we produce the probability distribution functions of the proper motions, with and without the systematic motions expected from a long-lived warp. A likelihood analysis is used to compare the expectations of the models with the observed proper motions from both Hipparcos and Gaia DR1. We find that the proper motions of the nearby OB stars are consistent with the signature of a warp, while those of the more distant stars (parallax<1 mas) are not. The kinematics of our sample of young OB stars suggests that systematic vertical motions in the disk cannot be explained by a simple model of a stable long-lived warp. The Galactic warp may either be a transient feature, or additional phenomena are acting on the gaseous component of the Milky Way, causing systematic vertical motions that are masking the expected warp signal. A larger and deeper sample of stars with Gaia astrometry will be needed to constrain the dynamical nature of the Galactic warp.

At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. We summarize Gaia DR1 and provide illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Gaia DR1 consists of: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the Hipparcos and Tycho-2 catalogues and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set,consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ~3000 Cepheid and RR Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas/yr for the proper motions. A systematic component of ~0.3 mas should be added to the parallax uncertainties. For the subset of ~94000 Hipparcos stars in the primary data set, the proper motions are much more precise at about 0.06 mas/yr. For the secondary astrometric data set, the typical uncertainty of the positions is ~10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ~0.03 mag over the magnitude range 5 to 20.7. Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

The GSC-II is an all-sky database of objects derived from the uncompressed DSS that the STScI has created from the Palomar and UK Schmidt survey plates and made available to the community. Like its predecessor (GSC-I), the GSC-II was primarily created to provide guide star information and observation planning support for HST. This version, however, is already employed at some of the ground-based new-technology telescopes such as GEMINI, VLT, and TNG, and will also be used to provide support for the JWST and Gaia space missions as well as LAMOST, one of the major ongoing scientific projects in China. Two catalogs have already been extracted from the GSC-II database and released to the astronomical community. A magnitude-limited (R=18.0) version, GSC2.2, was distributed soon after its production in 2001, while the GSC2.3 release has been available for general access since 2007. The GSC2.3 catalog described in this paper contains astrometry, photometry, and classification for 945,592,683 objects down to the magnitude limit of the plates. Positions are tied to the ICRS; for stellar sources, the all-sky average absolute error per coordinate ranges from 0.2" to 0.28" depending on magnitude. When dealing with extended objects, astrometric errors are 20% worse in the case of galaxies and approximately a factor of 2 worse for blended images. Stellar photometry is determined to 0.13-0.22 mag as a function of magnitude and photographic passbands (B,R,I). Outside of the galactic plane, stellar classification is reliable to at least 90% confidence for magnitudes brighter than R=19.5, and the catalog is complete to R=20.

Gaia is an ambitious space astrometry mission of ESA with a main objective to map the sky in astrometry and photometry down to a magnitude 20 by the end of the next decade. While the mission is built and operated by ESA and an industrial consortium, the data processing is entrusted to a consortium formed by the scientific community, which was formed in 2006 and formally selected by ESA one year later. The satellite will downlink around 100 TB of raw telemetry data over a mission duration of 5 years from which a very complex iterative processing will lead to the final science output: astrometry with a final accuracy of a few tens of microarcseconds, epoch photometry in wide and narrow bands, radial velocity and spectra for the stars brighter than 17 mag. We discuss the general principles and main difficulties of this very large data processing and present the organisation of the European Consortium responsible for its design and implementation.

Gaia is ESA's ambitious space astrometry mission the main objective of which is to astrometrically and spectro-photometrically map 1000 Million celestial objects (mostly in our galaxy) with unprecedented accuracy. The announcement of opportunity for the data processing will be issued by ESA late in 2006. The Gaia Data Processing and Analysis Consortium (DPAC) has been formed recently and is preparing an answer. The satellite will downlink close to 100 TB of raw telemetry data over 5 years. To achieve its required accuracy of a few 10s of Microarcsecond astrometry, a highly involved processing of this data is required. In addition to the main astrometric instrument Gaia will host a Radial Velocity instrument, two low-resolution dispersers for multi-color photometry and two Star Mappers. Gaia is a flying Giga Pixel camera. The various instruments each require relatively complex processing while at the same time being interdependent. We describe the overall composition of the DPAC and the envisaged overall architecture of the Gaia data processing system. We shall delve further into the core processing - one of the nine, so-called, coordination units comprising the Gaia processing system.

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined.

A large-scale three-dimensional model of Galactic extinction is presented based on the Galactic dust distribution model of Drimmel and Spergel (2001). The extinction A_V to any point within the Galactic disk can be quickly deduced using a set of three-dimensional cartesian grids. Extinctions from the model are compared to empirical extinction measures, including lines-of-sight in and near the Galactic plane using optical and NIR extinction measures; in particular we show how extinction can be derived from NIR color-magnitude diagrams in the Galactic plane to a distance of 8 kiloparsec.

We present parallax determinations for six white dwarf candidates in the Torino Observatory Parallax Program. The absolute parallaxes are found with precisions at the 2-3 milliarcsecond level. For WD 1126+185 we find a distance incompatible with being a white dwarf, implying an incorrect classification. For WD 2216+484 we find our distance is consistent with a simple DA white dwarf rather than a composite system as previously proposed in the literature. In general it is found that the published photometric distance is an overestimate of the distance found here.

We present a proper motion survey over about 200 square degrees towards the NGP, based on the material used for the construction of the GSC-II, that we are using to study the vertical structure and kinematics of the Galaxy. In particular, we measured the rotation velocity of the halo up to 10 kpc above the galactic plane traced by a sample of RR Lyr{\ae} and BHB giants for which radial velocities were used to recover the complete distribution of the spatial velocities. Finally, the impact of astrometric and spectroscopic GAIA observation are discussed.

A probable sky map of the emission from a short-lived isotrope produced by massive stars is presented. The model is based on the nonaxisymmetric component of a dust distribution model developed to reproduce Galactic FIR emission. Features seen in COMPTEL observations are qualitatively reproduced.

We present a three dimensional model for the Milky Way fit to the far-infrared (FIR) and near-infrared (NIR) data from the COBE/DIRBE instrument for Galactic latitudes |b| < 30 degrees and to within 20 degrees of the Galactic center. Because of the low optical depth at 240 microns, the FIR emission traces the distribution of galactic dust in the Galaxy. We model the dust distribution as due to three components: a warped exponential disk with scale length 0.28 R_sun and a flaring scale height, a spiral arm component with four arms as traced by Galactic HII regions, and the local (Orion) arm which produces prominent emission features at galactic longitude l approx. 80 and -100 degrees. A Cosmic Infrared Background of 1.07 MJy/sr is recovered, consistent with previous determinations. The dust distribution is then used to calculate absorption in J and K, and the stellar emission in these wavebands is modeled with two components: a warped exponential disk with a scale length of 0.28 R_sun and a spiral arm component dominated by two arms. This small scale length is consistent with a maximal disk model for our Galaxy, which is inconsistent with the cuspy dark matter halos predicted in CDM models. We find different amplitudes for the warp in the stars and dust, which starts within the Solar Circle.

Using far-infrared 240 micron and near-infrared K band data from the COBE/DIRBE instrument, we model the Galactic stellar and dust distribution. Making the assumption that the Galaxy is transparent in the 240 micron band, the dust emission is modeled using the following components: a warped exponential disk of scale length 0.26 R_sun, a spiral arm component as mapped by HII regions, and a feature coinciding with the local (Orion) arm. The dust distribution is used to calculate absorption in the K band, and the stellar emission is likewise modeled with a warped exponential disk, with a scale length of 0.29 R_sun, and a spiral arm component. Models of the K band emission in the Galactic plane indicate that in this waveband a two arm spiral dominates the nonaxisymmetric emission. The warp is evident in both the dust and stellar component, and is found to start within the Solar Circle.

Emission profiles of the Galactic plane in K and at $240\micron$ are presented, and features associated with the tangents of the spiral arms are identified. In the K band, which traces stellar emission and suffers little from absorbtion, features associated with the arm tangents indicate that a two-armed logarithmic spiral dominates the nonaxisymmetric structure of the Milky Way. In contrast, the $240\micron$ emission from dust entrained in the interstellar gas is consistent with a four-armed model, in concordance with radio data and optical spiral tracers. This suggests that the non-axisymmetric mass perturbation responsible for the four-armed spiral structure in the gas has a two rather than four-armed structure.

The kinematics of distant OB stars perpendicular to the Galactic plane, inferred from proper motions in the Hipparcos catalogue, are analysed and compared to the kinematic signature that would be induced by a long-lived Galactic warp. Previously we reported that the kinematics of the OB stars toward the anticenter were inconsistent with the assumption of a long-lived warp (Nature 392, 471), showing negative systematic motions as opposed to the expected positive motions. Using a larger sample of OB stars, improved distances, and a more accurate model for the uncertainties, we confirm our previous result for a larger range of galactocentric radii. However, we note that the new model for errors in the photometric distances reveal an important bias that causes the observed systematic vertical motions to be smaller than their true values. Using synthetic catalogues we investigate the effect of this bias on the systematic vertical motions in conjunction with the possibility of a smaller warp amplitude and a warp precession rate, both of which can also lead to smaller systematic motions. Taken together these three effects can work to produce negative observed systematic vertical motions, similar to those detected in the data, though only with both excessively high precession rates (-25 km/sec/kpc) and very large photometric errors (1 magnitude).