Abstract The Nordland area in northern Norway is the seismically most active area on mainland Fennoscandia. It exhibits patterns of coastal extension, which contrasts with the first-order regional stress pattern that reflects compression aligned with the North Atlantic ridge push. The regional stress field has been considered to emanate from the interaction of ridge push and glacial isostatic adjustment; while the local stress pattern can be additionally influenced by gravitational, topographic stresses, as well as the flexural effects of erosion and sediment deposition. We employ finite element numerical models at a crustal scale to study the 3D stress field, using existing geometric constraints from previous geophysical studies. Internal body forces, induced by variations in density, topography or Moho depth, already yield significant deviatoric stresses. In the models tested, these can strongly influence the near-surface stress regime, in particular for the continental margin setting we are investigating. In addition, redistribution of rock mass, which occurred mainly under Pleistocene glaciation, can modify the stress field significantly on a semi-regional scale. We consider this process to be the main driver for the coastal extension, in particular in areas where erosion has been high.