modeling marine magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back...

17
Group A Wei Wang & Yongfei Wang Oct. 28 2013 Modeling marine Magnetic anomalies

Upload: others

Post on 22-Mar-2021

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

Group A Wei Wang & Yongfei Wang

Oct. 28 2013

Modeling marine Magnetic anomalies

Page 2: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

o  Background o  Basic theory for forward and inverse

problem o  Numerical results and inversion o  Conclusion

Outline

Page 3: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

In the late 1950's, scientists mapped the present-day magnetic field generated by rocks on the floor of the Pacific Ocean. The volcanic rocks which make up the sea floor have magnetization because, as they cool, magnetic minerals within the rock align to the Earth's magnetic field. The intensity of the magnetic field they measured was very different from the intensity they had calculated.

Background

Figure-1 Oceanic structure Figure-2 Magnetic anomalies distribution and signals (http://volcano.oregonstate.edu/)

Sandwell,2001

Page 4: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

Basic theory to calculate magnetic anomalies To calculate the anomalous scalar field on the sea surface due to thin magnetic stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB = −∇U∇2U = 0∇2U = µ0∇•M

We define a scalar potential U and a magnetization vector M. The magnetic anomaly ΔB is the negative gradient of the potential. The potential satisfies Laplace’s equation above the source layer and is satisfies Poisson’s equation within the source layer.

Sandwell,2001

Page 5: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  Use the following equation relating the fourier transform of the magnetic anomaly to the fourier transform of the magnetic timescale.

•  This earth filter attenuates both long and short wavelengths so it acts like a band-pass filter. In the space domain it modifies the shape of the square-wave reversal pattern as sketched in the following diagram.

Sandwell,2001

Page 6: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

Numerical Results 1) Earth Filter

Sandwell, 2001

Page 7: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  2). Skewness function:

•  Periodic function •  Different skewness values

Page 8: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  3). Observed paleo-magnetic field

•  Symmetric pattern •  Skewness

Page 9: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  4). Grid Search Method: •  Method:

– Separate the parameter domain into to grids – Minimum value of misfit function

•  Assumption: – Speading rate kept the same – Skewness kept the same

•  Three parameters: – Mean crust thickness: 1-9 km

•  Not sensitive, set to 3 km – Skewness: -180 – 180 degree

•  Grid 10 degree – Spreading rate: 30 – 70 km/my

•  Grid 2 km/my

Page 10: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  5). Misfit function

•  a). Difference between observed and synthetic paleo-magnetic field.

•  b). Emphersize the phase difference to fit the polarity misfit.

L. Zhao & D. Helmberger 1993

!"" = (!!"# ! ∗ !!"#(!)!!!"!!! )

( !!"#(!) !)!!!"!!! ( !!"#(!)

!)!!!"!!!

!!!!!!!!!

Page 11: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  Misfit function value map

Page 12: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  Best Fit Solution: –  Misfit function: 0.326 –  Skewness: -40 degree –  Spreading rate: 46 km/my

Page 13: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  Repeat the experiment with half time range

•  Misfit function value map

Page 14: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  Best Fit Solution: –  Misfit function: 0.62 ( almost double ) –  Skewness: 0 degree ( vary with time ) –  Spreading rate: 46 km/my ( not vary with time )

Page 15: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

•  1. The synthetic and observation fit well within 10 my and misfit after 10 my

•  2. Mean crust thickness: 3 km (no constraint)

•  3. Spreading rate is 46 km/my •  4. Skewness is -40 degree,

– not located above the polarity

Conclusion

Page 16: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ

Discussion •  1. The two assumptions:

– The speading rate may be change a little ? – The skewness changes a lot?

•  2. Noise in data: – constrain the accuracy – constrain the meshes of grid search

•  3. Future-work – Change the skewness with time for better

accuracy

Page 17: Modeling marine Magnetic anomalies · 2013. 10. 30. · stripes on the seafloor, we go back Poisson’s equation relating magnetic field to magnetization. ΔB=−∇U ∇2U=0 ∇2U=µ