Greedy selection of optimal location of sensors for uncertainty reduction in seismic moment tensor inversion

We address an optimal sensor placement problem through Bayesian experimental design for seismic full waveform inversion for the recovery of the associated moment tensor. The objective is that of optimally choosing the location of the sensors (stations) from which to collect the observed data. The Shannon expected information gain is used as the objective function to search for the optimal network of sensors. A closed form for such objective is available due to the linear structure of the forward problem, as well as the Gaussian modeling of the observational errors and prior distribution. The resulting problem being inherently combinatorial, a greedy algorithm is deployed to sequentially select the sensor locations that form the best network for learning the moment tensor. Numerical results are presented to display the optimal network of sensors and how the uncertainty in the inferred seismic moment tensor contracts. The scenario of full three-dimensional velocity models or unknown earthquake source locations is treated as nuisance uncertainty, contributing to the overall uncertainty without being the focus of the inversion. This is addressed using a consensus approach over a set of realizations of the nuisance parameter. We analyzed the resulting network of stations for the moment tensor inversion under model misspecification, which reflects realistic data-generating processes.