We report preliminary results of a global 3‐D ionospheric electron density reanalysis demonstration study during 2002–2011 based on multisource data assimilation. The monthly global ionospheric electron density reanalysis has been done by assimilating the quiet days ionospheric data into a data assimilation model constructed using the International Reference Ionosphere (IRI) 2007 model and a Kalman filter technique. These data include global navigation satellite system (GNSS) observations of ionospheric total electron content (TEC) from ground‐based stations, ionospheric radio occultations by CHAMP, GRACE, COSMIC, SAC‐C, Metop‐A, and the TerraSAR‐X satellites, and Jason‐1 and 2 altimeter TEC measurements. The output of the reanalysis are 3‐D gridded ionospheric electron densities with temporal and spatial resolutions of 1 h in universal time, 5° in latitude, 10° in longitude, and ∼30 km in altitude. The climatological features of the reanalysis results, such as solar activity dependence, seasonal variations, and the global morphology of the ionosphere, agree well with those in the empirical models and observations. The global electron content derived from the international GNSS service global ionospheric maps, the observed electron density profiles from the Poker Flat Incoherent Scatter Radar during 2007–2010, and f_oF₂ observed by the global ionosonde network during 2002–2011 are used to validate the reanalysis method. All comparisons show that the reanalysis have smaller deviations and biases than the IRI‐2007 predictions. Especially after April 2006 when the six COSMIC satellites were launched, the reanalysis shows significant improvement over the IRI predictions. The obvious overestimation of the low‐latitude ionospheric F region densities by the IRI model during the 23/24 solar minimum is corrected well by the reanalysis. The potential application and improvements of the reanalysis are also discussed.