Effect of non-uniform station coverage on the inversion for earthquake rupture history for a Haskell-type source model

In order to determine how reliably one can invert accelerograms to determin e the rupture process details, when the station configuration is less than optimal, we use the vertical component of synthetic accelerograms for a Has kell-type earthquake rupture model, at stations in the vicinity of a dip-sl ip fault and solve the inverse problem. Of the various station configuratio ns used, one is a uniform distribution and the others are very non-uniform. Faults of two different aspect ratios are considered. We mainly use much l arger spatial and temporal cell sizes in the inversion than we use to const ruct the artificial data. The fault mechanism and the fault area are taken as known in the inversions. To solve the inverse problem, we use the method of linear programming and stabilize the solution by the use of physical co nstraints. The constraints of positivity of the slip rates on the fault is used in all cases in this study. In some cases, additional physical constra ints such as preassigning the final moment, the rupture speed, and so on, a re also used. We find that using a cell size almost double the wavelength o f interest, we are able to reproduce the solution of the problem, even when we add a small amount of random noise to the artificial data, provided the source medium structure is known. We show that the best station configurat ion is when the stations are on the hanging wall, due to the fact that they provide the best illumination of the fault surface. This provides an incen tive to install permanent ocean bottom strong ground motion stations in sub duction zones. We also analyzed the effect of the rupture propagation direc tion on the results of the inversion showing that even four stations are su fficient to retrieve the rupture process if they are in the forward directi on of the rupture propagation; the results for this case are better than wh en the four stations are placed in the backward direction, even when their positions are such that they illuminate the fault in exactly the same way a s the four stations in the forward direction. Thus azimuthal distribution a nd the resulting illumination of the fault as well as the relation of the p osition of the stations to the direction of rupture propagation are more im portant than simply the number of stations. Finally, we find that propel kn owledge of source medium structure is essential to recover the source proce ss details reliably and that poor knowledge of crustal structure cannot be compensated by adding stations or by additional constraints.