The initial stress field is determined as follows: A tensional stress acting N103E (nearly parallel to the SISZ; cf. DeMets et al. 1990) is assumed, due to ridge push or basal drag of the adjacent plates. This rifting induces shear stresses in the region of the transform fault. The stress magnitude, which is unknown, is set to a value that produces left-lateral shear stresses in E-W direction as large as the stress drop determined for the largest event (M=7.1) in the studied earthquake sequence.
Tensional stresses at both ridges are modelled as constantly being released to have zero values at the rifts. These are the major disturbances of the unknown background stresses.
On this initial field, the stress changes due to earthquakes are iteratively superposed as well as the stress changes due to further spreading at the ridge segments based on an opening of 2 cm/year. This value is taken from DeMets et al. (1990). As the simplest assumption, lacking other data, the spreading rate is taken to be constant during the modelled time period, even though this can be questioned as for instance the present debate on the stress increase in the New Madrid seismic zone shows.
The stress field before every event is thus the sum of the initial field, the stress drop of all preceding events, and the plate tectonic stress build-up since the starting time of the model, which is set to 1706, when the first event in the series occurred.
Results were calculated for 56x44 test-points covering 280 km in E-W direction and 220 km in N-S direction. Stresses were computed for a homogeneous half-space, as a starting model. Although surface stress changes are calculated, these should be representative for crustal stresses using values for the moduli, that are typical for oceanic crust  and not for sedimentary layers at the surface.