The laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) technique is considered versatile for multielemental analysis and imaging because it is easy to handle, compatible with different types of solid samples, requires minimal sample preparation, and provides high spatial resolution and sensitivity. One of the challenges of imaging analysis is to obtain accurate and precise spatial information of elements distribution in the sample, so the optimization of the laser ablation (LA) parameters is essential. In this context, this study aimed to optimize the LA parameters for direct analysis of speleothem samples. Laser intensity, frequency, and spot diameter were evaluated through multivariate experimental design and multi-response data, the influence of ablation scan speed and the use of Ca-44 as an internal standard (IS) for the qualitative image of Zn-66, Ba-137, 55Mn, Fe-57, Sr-88, 6(0)Ni, and Mg-26 distribution in the sample were also evaluated. The multivariate optimization revealed positive interactions between the parameters evaluated, i.e., the greater the laser intensity, LA frequency, and spot diameter, the greater the analyte signal and, thus the sensitivity. Therefore, 90% laser ablation intensity, 20 Hz repetition rate, and 100 mu m spot diameter were selected. In the scan speed evaluation, the images obtained with 40 and 20 mu m s(-1) were very similar for all isotopes. The use of 44Ca as IS did not impact the resolution of the images. The use of Ca-44 can provide important information about the speleothem formation and paleoclimate changes.