The pre-drilling prediction of (paleo)-stress helps defining fracture and fault patterns which can either help or hinder hydrocarbon production, and helps to decipher the tectonic history of the crust. At the moment, paleostress studies are mostly based on the detailed mapping of faults and movement directions in the field. This limits the applicability of this approach to areas where the rocks of interest are outcropping and excludes sedimentary basins. In this work, a method is developed and tested to determine the paleostress stratigraphy solely based on 3D reflection data, and the initial work for a similar approach for ductile evaporites is presented. For the determination of paleostress, knowledge of the orientation and slip direction of the individual faults in the set is required. Both in field and seismic surveys it is relatively easy to determine the fault orientation. Determining the slip direction however remains a challenge in seismic data. In field studies, measurement of the slickenside orientation provides a quick and easy slip direction determination, but these structures are too small to be resolved in seismic data. Due to the converging pattern of slip lines, the measured slickenline does not need to represent the tectonic movement direction. Furthermore, since faults generally rupture in a number of seismic events distributed over the fault plane, contradicting slip directions formed during the same tectonic phase can overprint each other, complicating the interpretation of field based slip determination. The easiest way to determine the slip direction on faults is by connecting two points that were on direct opposite sides of the fault, prior to deformation. Unfortunately, no faulted channels or lineaments were present in the data set. By careful mapping of the Allan lines however, the slip directions of a number of faults, over a number of time periods were constrained. This data allowed for the determination of the paleostress stratigraphy of the NW corner of the Groningen High in the ...