Integrated simulation results of femtosecond laser ablation of copper were compared with new experimental data. The numerical analysis was performed using our newly developed FEMTO-2D computer package based on the solution of the two-temperature model. Thermal dependence of target optical and thermodynamic processes was carefully considered. The experimental work was conducted with our 40 fs 800 nm Ti:sapphire laser in the energy range from 0.14 mJ to 0.77 mJ. Comparison of measured ablation profiles with simulation predictions based on phase explosion criterion has demonstrated that more than one ablation mechanisms contribute to the total material removal even in the laser intensity range where explosive boiling is dominating. Good correlation between experimental and simulation results was observed for skin depth and hot electron diffusion depth – two parameters commonly considered to identify two ablation regimes in metal. Analysis of the development dynamics for electron–lattice coupling and electron thermal conduction allowed explaining different ablation regimes because of the interplay of the two parameters.