Australia's lithospheric density field, and its isostatic equilibration

Alan Aitken, C. Altinay, L. Gross

    Research output: Contribution to journalArticle

    6 Citations (Scopus)

    Abstract

    © The Authors 2015. Density is a key driver of tectonic processes, but it is a difficult property to define well in the lithosphere because the gravity method is non-unique, and because converting to density from seismic velocity models, themselves non-unique, is also highly uncertain. Here we use a new approach to define the lithospheric density field of Australia, covering from 100°E to 165°E, from 5°N to 55?S and from the crust surface to 300 km depth. A reference model was derived primarily from the recently released Australian Seismological Reference Model, and refined further using additional models of sedimentary basin thickness and crustal thickness. A novel form of finite-element method based deterministic gravity inversion was applied in geodetic coordinates, implemented within the open-source escript modelling environment. Three spatial resolutions were modelled: half-, quarter- and eighth-degree in latitude and longitude, with vertical resolutions of 5, 2.5 and 1.25 km, respectively. Parameter sweeps for the key inversion regularization parameters show that parameter selection is not scale dependent. The sweep results also show that finer resolutions are more sensitive to the uppermost crust, but less sensitive to the mid- to lower-crust and uppermost mantle than lower resolutions. All resolutions show similar sensitivity below about 100 km depth. The final density model shows that Australia's lithospheric density field is strongly layered but also has large lateral density contrasts at all depths.Within the continental crust, the structure of the middle and lower crust differs significantly from the crystalline upper crust, suggesting that the tectonic processes or events preserved in the deep crust differ from those preserved in the shallower crust. The lithospheric mantle structure is not extensively modified from the reference model, but the results reinforce the systematic difference between the density of the oceanic and continental domains, and help identify subdivisions within each. The lithospheric static pressure field was resolved in 3D from the gravity and density fields. The pressure field model also highlights the fundamental difference between the oceanic and continental domains, with the former possessing lower pressure through most of the model. Overall pressure variability is large in the upper crust (60 MPa) but reduces significantly by -30 km elevation (20-30MPa). By -50 km elevation, thick lower-crust generates further disequilibria (25-35MPa) that are not compensated until -125 km elevation (10-20MPa). Beneath -125 km elevation higher pressure is observed in the continental domain, extending to the base of the model. This indicates a lithosphere that is to a large degree isostatically compensated near the base of the felsic-intermediate continental crust, and again near the theoretical base of mature oceanic lithosphere.
    Original languageEnglish
    Pages (from-to)1961-1976
    JournalGeophysical Journal International
    Volume203
    Issue number3
    DOIs
    Publication statusPublished - 2015

    Fingerprint

    crusts
    lower crust
    crust
    Gravitation
    lithosphere
    pressure field
    Tectonics
    gravity
    upper crust
    continental crust
    gravitation
    pressure distribution
    tectonics
    Earth mantle
    geodetic coordinates
    mantle structure
    inversions
    crustal thickness
    oceanic lithosphere
    lower mantle

    Cite this

    @article{d0bf13a4c560403794b00adc469378d1,
    title = "Australia's lithospheric density field, and its isostatic equilibration",
    abstract = "{\circledC} The Authors 2015. Density is a key driver of tectonic processes, but it is a difficult property to define well in the lithosphere because the gravity method is non-unique, and because converting to density from seismic velocity models, themselves non-unique, is also highly uncertain. Here we use a new approach to define the lithospheric density field of Australia, covering from 100°E to 165°E, from 5°N to 55?S and from the crust surface to 300 km depth. A reference model was derived primarily from the recently released Australian Seismological Reference Model, and refined further using additional models of sedimentary basin thickness and crustal thickness. A novel form of finite-element method based deterministic gravity inversion was applied in geodetic coordinates, implemented within the open-source escript modelling environment. Three spatial resolutions were modelled: half-, quarter- and eighth-degree in latitude and longitude, with vertical resolutions of 5, 2.5 and 1.25 km, respectively. Parameter sweeps for the key inversion regularization parameters show that parameter selection is not scale dependent. The sweep results also show that finer resolutions are more sensitive to the uppermost crust, but less sensitive to the mid- to lower-crust and uppermost mantle than lower resolutions. All resolutions show similar sensitivity below about 100 km depth. The final density model shows that Australia's lithospheric density field is strongly layered but also has large lateral density contrasts at all depths.Within the continental crust, the structure of the middle and lower crust differs significantly from the crystalline upper crust, suggesting that the tectonic processes or events preserved in the deep crust differ from those preserved in the shallower crust. The lithospheric mantle structure is not extensively modified from the reference model, but the results reinforce the systematic difference between the density of the oceanic and continental domains, and help identify subdivisions within each. The lithospheric static pressure field was resolved in 3D from the gravity and density fields. The pressure field model also highlights the fundamental difference between the oceanic and continental domains, with the former possessing lower pressure through most of the model. Overall pressure variability is large in the upper crust (60 MPa) but reduces significantly by -30 km elevation (20-30MPa). By -50 km elevation, thick lower-crust generates further disequilibria (25-35MPa) that are not compensated until -125 km elevation (10-20MPa). Beneath -125 km elevation higher pressure is observed in the continental domain, extending to the base of the model. This indicates a lithosphere that is to a large degree isostatically compensated near the base of the felsic-intermediate continental crust, and again near the theoretical base of mature oceanic lithosphere.",
    author = "Alan Aitken and C. Altinay and L. Gross",
    year = "2015",
    doi = "10.1093/gji/ggv396",
    language = "English",
    volume = "203",
    pages = "1961--1976",
    journal = "Geophysical Journal International",
    issn = "0956-540X",
    publisher = "Oxford University Press",
    number = "3",

    }

    Australia's lithospheric density field, and its isostatic equilibration. / Aitken, Alan; Altinay, C.; Gross, L.

    In: Geophysical Journal International, Vol. 203, No. 3, 2015, p. 1961-1976.

    Research output: Contribution to journalArticle

    TY - JOUR

    T1 - Australia's lithospheric density field, and its isostatic equilibration

    AU - Aitken, Alan

    AU - Altinay, C.

    AU - Gross, L.

    PY - 2015

    Y1 - 2015

    N2 - © The Authors 2015. Density is a key driver of tectonic processes, but it is a difficult property to define well in the lithosphere because the gravity method is non-unique, and because converting to density from seismic velocity models, themselves non-unique, is also highly uncertain. Here we use a new approach to define the lithospheric density field of Australia, covering from 100°E to 165°E, from 5°N to 55?S and from the crust surface to 300 km depth. A reference model was derived primarily from the recently released Australian Seismological Reference Model, and refined further using additional models of sedimentary basin thickness and crustal thickness. A novel form of finite-element method based deterministic gravity inversion was applied in geodetic coordinates, implemented within the open-source escript modelling environment. Three spatial resolutions were modelled: half-, quarter- and eighth-degree in latitude and longitude, with vertical resolutions of 5, 2.5 and 1.25 km, respectively. Parameter sweeps for the key inversion regularization parameters show that parameter selection is not scale dependent. The sweep results also show that finer resolutions are more sensitive to the uppermost crust, but less sensitive to the mid- to lower-crust and uppermost mantle than lower resolutions. All resolutions show similar sensitivity below about 100 km depth. The final density model shows that Australia's lithospheric density field is strongly layered but also has large lateral density contrasts at all depths.Within the continental crust, the structure of the middle and lower crust differs significantly from the crystalline upper crust, suggesting that the tectonic processes or events preserved in the deep crust differ from those preserved in the shallower crust. The lithospheric mantle structure is not extensively modified from the reference model, but the results reinforce the systematic difference between the density of the oceanic and continental domains, and help identify subdivisions within each. The lithospheric static pressure field was resolved in 3D from the gravity and density fields. The pressure field model also highlights the fundamental difference between the oceanic and continental domains, with the former possessing lower pressure through most of the model. Overall pressure variability is large in the upper crust (60 MPa) but reduces significantly by -30 km elevation (20-30MPa). By -50 km elevation, thick lower-crust generates further disequilibria (25-35MPa) that are not compensated until -125 km elevation (10-20MPa). Beneath -125 km elevation higher pressure is observed in the continental domain, extending to the base of the model. This indicates a lithosphere that is to a large degree isostatically compensated near the base of the felsic-intermediate continental crust, and again near the theoretical base of mature oceanic lithosphere.

    AB - © The Authors 2015. Density is a key driver of tectonic processes, but it is a difficult property to define well in the lithosphere because the gravity method is non-unique, and because converting to density from seismic velocity models, themselves non-unique, is also highly uncertain. Here we use a new approach to define the lithospheric density field of Australia, covering from 100°E to 165°E, from 5°N to 55?S and from the crust surface to 300 km depth. A reference model was derived primarily from the recently released Australian Seismological Reference Model, and refined further using additional models of sedimentary basin thickness and crustal thickness. A novel form of finite-element method based deterministic gravity inversion was applied in geodetic coordinates, implemented within the open-source escript modelling environment. Three spatial resolutions were modelled: half-, quarter- and eighth-degree in latitude and longitude, with vertical resolutions of 5, 2.5 and 1.25 km, respectively. Parameter sweeps for the key inversion regularization parameters show that parameter selection is not scale dependent. The sweep results also show that finer resolutions are more sensitive to the uppermost crust, but less sensitive to the mid- to lower-crust and uppermost mantle than lower resolutions. All resolutions show similar sensitivity below about 100 km depth. The final density model shows that Australia's lithospheric density field is strongly layered but also has large lateral density contrasts at all depths.Within the continental crust, the structure of the middle and lower crust differs significantly from the crystalline upper crust, suggesting that the tectonic processes or events preserved in the deep crust differ from those preserved in the shallower crust. The lithospheric mantle structure is not extensively modified from the reference model, but the results reinforce the systematic difference between the density of the oceanic and continental domains, and help identify subdivisions within each. The lithospheric static pressure field was resolved in 3D from the gravity and density fields. The pressure field model also highlights the fundamental difference between the oceanic and continental domains, with the former possessing lower pressure through most of the model. Overall pressure variability is large in the upper crust (60 MPa) but reduces significantly by -30 km elevation (20-30MPa). By -50 km elevation, thick lower-crust generates further disequilibria (25-35MPa) that are not compensated until -125 km elevation (10-20MPa). Beneath -125 km elevation higher pressure is observed in the continental domain, extending to the base of the model. This indicates a lithosphere that is to a large degree isostatically compensated near the base of the felsic-intermediate continental crust, and again near the theoretical base of mature oceanic lithosphere.

    U2 - 10.1093/gji/ggv396

    DO - 10.1093/gji/ggv396

    M3 - Article

    VL - 203

    SP - 1961

    EP - 1976

    JO - Geophysical Journal International

    JF - Geophysical Journal International

    SN - 0956-540X

    IS - 3

    ER -