Fault damage zones, composed of sub-seismic deformation structures, are difficult to detect using seismic data. Still, they can be related to fault throw, which is widely measured in the subsurface. This research employs a multiscale approach that integrates outcrop studies with seismic reflection data to investigate the attributes of fault damage zones affecting porous sandstones. We provide an in-depth understanding of the 3D fault zone volume, investigating correlations between geometric attributes of faults (width of damage zone, frequency of subsidiaries structures, and fault termination) in the outcrop and the fault throw in the subsurface, with insights from deformation mechanisms within the fault zone. The methods encompass 1D scanlines to constrain the damage zone width on the outcrop. At the same time, the subsurface analysis uses 3D seismic data, seismic attributes, and deep learning neural network (DNN) fault volumes to interpret fault geometries and quantify fault throw at different depths. Results show that the area presents a complex fault zone with multiple fault sets, deformation bands, and fracture corridors trending mainly NE-SW, NW-SE, and E-W. The integration of surface and subsurface fault data enabled the identification of two portions in each fault set outcropping at the fault tip for E-W/ESE-WNW-striking faults and the central part of the fault for NE-SW-striking faults. Faults with the greatest length do not outcrop with the largest damage zone width since they are outcropping the tip of the fault. The parallel faults overlap their damage zones, increasing the deformation zones in the affected sandstones. Fault throw and damage zone width present a positive correlation. This relation is affected by fault segments and subsidiary faults.