ABSTRACT We develop, validate and apply a forward model to estimate stellar atmospheric parameters (Teff, log g, and [Fe/H]), revised distances and extinctions for 220 million stars with XP spectra from Gaia DR3. Instead of using ab initio stellar models, we develop a data-driven model of Gaia XP spectra as a function of the stellar parameters, with a few straightforward built-in physical assumptions. We train our model on stellar atmospheric parameters from the LAMOST survey, which provides broad coverage of different spectral types. We model the Gaia XP spectra with all of their covariances, augmented by 2MASS and WISE photometry that greatly reduces degeneracies between stellar parameters, yielding more precise determinations of temperature and dust reddening. Taken together, our approach overcomes a number of important limitations that the astrophysical parameters released in Gaia DR3 faced, and exploits the full information content of the data. We provide the resulting catalogue of stellar atmospheric parameters, revised parallaxes, and extinction estimates, with all their uncertainties. The modelling procedure also produces an estimate of the optical extinction curve at the spectral resolution of the XP spectra (R ∼ 20–100), which agrees reasonably well with the R(V) = 3.1 CCM model. Remaining limitations that will be addressed in future work are that the model assumes a universal extinction law, ignores binary stars and does not cover all parts of the Hertzsprung–Russell Diagram (e.g. white dwarfs).