A method for determining spatial distribution of properties of the Earth's subsurface includes obtaining seismic data over a survey area of the Earth's subsurface. Controlled source electromagnetic survey data are obtained over substantially the same survey area. An initial model of the Earth's subsurface for each of the seismic data and the electromagnetic data is generated. Further data may include gravity, magnetics, seismics any type and borehole data. Each model is optimized on at least one model parameter. Consistency is determined between the models; and the at least one model parameter is adjusted and the optimizing and determining consistency are repeated until the models are consistent. Constraints are successively derived from the data sets and also cross checked against reservoir data where available.