We present a momentum-resolved study of strong field multiple ionization of ionic targets. Using a deconvolution method we are able to reconstruct the electron momenta from the ion momentum distributions after multiple ionization up to four sequential ionization steps. This technique allows an accurate determination of the saturation intensity as well as of the electron release times during the laser pulse. The measured results are discussed in comparison to typically used models of over-the-barrier ionization and tunnel ionization.Atoms exposed to super-intense laser pulses can be ionized to high charge states. In the optical regime, the ionization probability depends highly nonlinear on the field strength. Therefore, for a pulsed field, ionization is concentrated in a narrow intensity and a correspondingly narrow time interval for each ionization step: There is virtually no ionization at lower intensity because of negligible ionization rate and no ionization at higher intensity because all atoms or ions of the respective charge state are already ionized. Accordingly, the intensity where ionization peaks is an important quantity for many research areas of intense laser matter interaction. In the literature, different definitions and names are used for the characterization for this or related phenomena (e.g. appearance intensity) [1,2]. We will use the term saturation intensity for this quantity. This is of particular interest, e.g. in laser plasma physics, where correct modeling of ionization over a large number of charge states, which determines the spatio-temporal plasma density, is essen-The accurate determination of the saturation intensity is an important problem for modeling strong-field laser matter interaction. Difficulties exist on both on the experimental and the theoretical side. For the latter, a frequently used approach for the estimation of the saturation intensity is the over-the-barrier ionization (OBI) [1], i.e. the intensity where the Coulomb barrier is lowered below the ground state. Another approach is based on the tunneling rate, i.e. the theories of Perelomov, Popov and Terentev (PPT) [4] or Ammosov, Delone and Krainov (ADK) [5]. Both approaches have inherent problems: For OBI, the entire concept is controversial [6,7], for tunneling it is known that the approximations required to derive the respective formulas are not valid at the intensities typical for saturation with femtosecond pulses [8][9][10]. Experimentally the problem has been investigated by measuring the intensity dependence of the charge resolved ion yield [11,12], thereby the reliable determination of the intensity is the major problem. A meaningful verification of the predictive power of one or the other theoretical model is non-trivial, since the actual field strength at which a specific ionization event occurred is, in general not accessible experimentally.In this article, we communicate an experiment and method of data evaluation that enables an in situ intensity measurement of the saturation intensity for ions of various...