A full vector paleomagnetic record, comprising directional data and relative paleointensity (rPI), were derived from 16 sediment cores recovered from the southeastern Black Sea. The obtained data were used to create a stack covering the time window between 68.9 and 14.5 ka. Age models are based on radiocarbon dating and correlations of warming/cooling cycles monitored by high-resolution X-ray fluorescence (XRF) elementary ratios as well as by ice-rafted debris (IRD) in Black Sea sediments to the sequence of 'Dansgaard-Oeschger' (D-O) events defined from the Greenland ice core oxygen isotope stratigraphy. The PSV stacks between 68.9 ka and 42.5 ka were calculated for 100-year time bins, as well as between 39.4 ka and 20 ka . The PSV stacks spanning from 42.5 to 39.4 ka, covering the Laschamps excursion and sampled with a special high-resolution sampling, were calculated for 40-year time bins. With plenty of samples, the PSV stacks between 20 ka and 14.5 ka were calculated for 50-year time bins. I chrm: characteristic inclination D chrm: characteristic declination alpha95: error range of inclination and declination by applying Fisher statistics ' Slope-rPI: relative paleointensity determined by slope NRM/ARM during alternating field demagnetization; NRM: natural remanent magnetization, ARM: anhysteretic remanent magnetization n: number of samples used in stacking for every time bins std: standar error range of relative paleointensity F(µT): total field intensity X': field component in North and South direction Y': field component in West and East direction Z': field componnet parallel to geocentric axial dipole Pla: latitude of virtual geomagnetic pole Plo: longitude of virtual geomagnetic pole Pi: paleosecular variation index