f EFFICIENT AND SCALABLE ALGORITHMS FOR THE INVERSION OF VERY LARGE AEM DATASETS

Recently the standard size of airborne electromagnetic (AEM) datasets has been growing rapidly, while many surveys at the same time focus on modern applications that require very high resolution. Often, the particular target structures of interest cannot be identified directly from the data and approximate inversions or data transform techniques also prove insufficient. This type of application includes environmental studies of aquifers, salinity and geology and also cases of relevance to the mining industry e.g. uranium exploration. In these cases highly accurate full non-linear inversion provides the only reliable solution, which is a computational challenge. In designing an inversion scheme for these types of problems it is further advantageous to utilize as many assumptions as possible. The target structures of interest can often be regarded locally one dimensional, quasi layered and spatially coherent which can be integrated in the inversion. If inverting multiple soundings at a time for a set of spatially constrained 1D models, information is allowed to propagate spatially and a vastly improved model result can be obtained. Ideally, entire surveys should be inverted at a time in huge constrained problems, but poor scaling properties of the underlying algorithms typically makes this impossible.