diff --git a/_freeze/index/execute-results/docx.json b/_freeze/index/execute-results/docx.json index f41d186..61fb894 100644 --- a/_freeze/index/execute-results/docx.json +++ b/_freeze/index/execute-results/docx.json @@ -1,7 +1,7 @@ { "hash": "8ccdb86df4a0f84495987258f19b87bf", "result": { - "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#f6229d6a .cell .markdown}\n## Introduction\n:::\n\n::: {.cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {.cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-docx/fig-timeline-output-1.png){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {.cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=12}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#485c10f2 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", + "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#9f477374 .cell .markdown}\n## Introduction\n:::\n\n::: {.cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {.cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-docx/fig-timeline-output-1.png){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {.cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=6}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#f319e497 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", "supporting": [ "index_files/figure-docx" ], diff --git a/_freeze/index/execute-results/html.json b/_freeze/index/execute-results/html.json index e70d10f..3109db8 100644 --- a/_freeze/index/execute-results/html.json +++ b/_freeze/index/execute-results/html.json @@ -1,9 +1,9 @@ { "hash": "8ccdb86df4a0f84495987258f19b87bf", "result": { - "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#42487957 .cell .markdown}\n## Introduction\n:::\n\n::: {#0dd43b78 .cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {#cell-fig-timeline .cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-html/fig-timeline-output-1.png){#fig-timeline width=466 height=115 fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {#5018ff84 .cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=9}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#e96955a2 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", + "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#965dd61c .cell .markdown}\n## Introduction\n:::\n\n::: {#9dfe3472 .cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {#cell-fig-timeline .cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-html/fig-timeline-output-1.png){#fig-timeline width=466 height=115 fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {#be9fb910 .cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=3}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#2bf43e82 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", "supporting": [ - "index_files/figure-html" + "index_files" ], "filters": [], "includes": {} diff --git a/_freeze/index/execute-results/tex.json b/_freeze/index/execute-results/tex.json index 71f1c32..0c51bcb 100644 --- a/_freeze/index/execute-results/tex.json +++ b/_freeze/index/execute-results/tex.json @@ -1,7 +1,7 @@ { "hash": "8ccdb86df4a0f84495987258f19b87bf", "result": { - "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#7b0a4893 .cell .markdown}\n## Introduction\n:::\n\n::: {.cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {.cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-pdf/fig-timeline-output-1.pdf){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.' fig-pos='H'}\n:::\n:::\n\n\n::: {.cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=18}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#1eb4bd78 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", + "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#b61468ee .cell .markdown}\n## Introduction\n:::\n\n::: {.cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {.cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-pdf/fig-timeline-output-1.pdf){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.' fig-pos='H'}\n:::\n:::\n\n\n::: {.cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=12}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#7d4d67ef .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", "supporting": [ "index_files/figure-pdf" ], diff --git a/_freeze/index/execute-results/xml.json b/_freeze/index/execute-results/xml.json index a1445d2..ec89c5f 100644 --- a/_freeze/index/execute-results/xml.json +++ b/_freeze/index/execute-results/xml.json @@ -1,7 +1,7 @@ { "hash": "8ccdb86df4a0f84495987258f19b87bf", "result": { - "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#321d9992 .cell .markdown}\n## Introduction\n:::\n\n::: {#a5dd21fc .cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {#cell-fig-timeline .cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-jats/fig-timeline-output-1.png){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {#86046ea2 .cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=15}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#b471fd7f .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", + "markdown": "---\ntitle: La Palma Earthquakes\nauthor:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n corresponding: true\n email: steve@curvenote.com\n roles:\n - Investigation\n - Project administration\n - Software\n - Visualization\n affiliations:\n - Curvenote\n - name: Rowan Cockett\n orcid: 0000-0002-7859-8394\n corresponding: false\n roles: []\n affiliations:\n - Curvenote\nkeywords:\n - La Palma\n - Earthquakes\nabstract: |\n In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). ...\nplain-language-summary: |\n Earthquake data for the island of La Palma from the September 2021 eruption is found ...\nkey-points:\n - A web scraping script was developed to pull data from the Instituto Geogràphico Nacional into a machine-readable form for analysis\n - Earthquake events on La Palma are consistent with the presence of both mantle and crustal reservoirs.\ndate: last-modified\nbibliography: references.bib\ncitation:\n container-title: Earth and Space Science\nnumber-sections: true\n---\n\n:::{#67a3a4e8 .cell .markdown}\n## Introduction\n:::\n\n::: {#c69ee907 .cell execution_count=1}\n``` {.python .cell-code .hidden}\nimport matplotlib.pyplot as plt\nimport numpy as np\neruptions = [1492, 1585, 1646, 1677, 1712, 1949, 1971, 2021]\n```\n:::\n\n\n::: {#cell-fig-timeline .cell execution_count=2}\n``` {.python .cell-code .hidden}\nplt.figure(figsize=(6, 1))\nplt.eventplot(eruptions, lineoffsets=0, linelengths=0.1, color='black')\nplt.gca().axes.get_yaxis().set_visible(False)\nplt.ylabel('')\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Timeline of recent earthquakes on La Palma](index_files/figure-jats/fig-timeline-output-1.png){#fig-timeline fig-alt='An event plot of the years of the last 8 eruptions on La Palma.'}\n:::\n:::\n\n\n::: {#b397d1c8 .cell execution_count=3}\n``` {.python .cell-code .hidden}\navg_years_between_eruptions = np.mean(np.diff(eruptions[:-1]))\navg_years_between_eruptions\n```\n\n::: {.cell-output .cell-output-display .hidden execution_count=9}\n```\n79.83333333333333\n```\n:::\n:::\n\n\n:::{#1962cec2 .cell .markdown}\nBased on data up to and including 1971, eruptions on La Palma happen every 79.8 years on average.\n\nStudies of the magma systems feeding the volcano, such as @marrero2019, have proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).\n\nEight eruptions have been recorded since the late 1400s (@fig-timeline).\n\nData and methods are discussed in @sec-data-methods.\n\nLet $x$ denote the number of eruptions in a year. Then, $x$ can be modeled by a Poisson distribution\n\n$$\np(x) = \\frac{e^{-\\lambda} \\lambda^{x}}{x !}\n$$ {#eq-poisson}\n\nwhere $\\lambda$ is the rate of eruptions per year. Using @eq-poisson, the probability of an eruption in the next $t$ years can be calculated.\n\n| Name | Year |\n|---------------------|------|\n| Current | 2021 |\n| Teneguía | 1971 |\n| Nambroque | 1949 |\n| El Charco | 1712 |\n| Volcán San Antonio | 1677 |\n| Volcán San Martin | 1646 |\n| Tajuya near El Paso | 1585 |\n| Montaña Quemada | 1492 |\n\n: Recent historic eruptions on La Palma {#tbl-history}\n\n@tbl-history summarises the eruptions recorded since the colonization of the islands by Europeans in the late 1400s.\n\n![Map of La Palma](images/la-palma-map.png){#fig-map}\n\nLa Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands (@fig-map).\n\n\n\n\n{{< embed notebooks/data-screening.qmd#fig-spatial-plot >}}\n\n\n\n\n\n\n\n@fig-spatial-plot shows the location of recent Earthquakes on La Palma.\n\n## Data & Methods {#sec-data-methods}\n\n## Conclusion\n\n## References {.unnumbered}\n\n::: {#refs}\n:::\n:::\n\n", "supporting": [ "index_files/figure-jats" ], diff --git a/_freeze/index/figure-ipynb/cell-8-output-1.png b/_freeze/index/figure-ipynb/cell-8-output-1.png deleted file mode 100644 index 1494caa..0000000 Binary files a/_freeze/index/figure-ipynb/cell-8-output-1.png and /dev/null differ diff --git a/_freeze/index/figure-pdf/fig-timeline-output-1.pdf b/_freeze/index/figure-pdf/fig-timeline-output-1.pdf index 6aeea53..e239730 100644 Binary files a/_freeze/index/figure-pdf/fig-timeline-output-1.pdf and b/_freeze/index/figure-pdf/fig-timeline-output-1.pdf differ diff --git a/_freeze/notebooks/data-screening/execute-results/html.json b/_freeze/notebooks/data-screening/execute-results/html.json index 9950c05..2dc0881 100644 --- a/_freeze/notebooks/data-screening/execute-results/html.json +++ b/_freeze/notebooks/data-screening/execute-results/html.json @@ -1,7 +1,7 @@ { "hash": "d87b6d9ac2228b736962871bfd9167c0", "result": { - "markdown": "---\ntitle: Data Screening\ndescription: Screening and filtering data from the IGN catalog down to the relevant events for La Palma\nauthors:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n - name: Charlotte Wickham\n note: Edited original notebook by Steve Purves\n---\n\n::: {#6c6b8b20 .cell execution_count=1}\n``` {}\nimport pandas as pd\nimport matplotlib.pyplot as plt\n```\n:::\n\n\n:::{#a6b6d890 .cell .markdown}\nImport `la-palma.csv` - a clean subset of [`catalogoComunSV_1663233588717.csv`](https://github.com/Notebooks-Now/submission-quarto-full/blob/main/data/catalogoComunSV_1663233588717.csv):\n:::\n\n::: {#09c76fc6 .cell execution_count=2}\n``` {}\ndf = pd.read_csv('la-palma.csv')\ndf['DateTime'] = pd.to_datetime(df['DateTime'])\ndf = df.sort_values(by=['DateTime'], ascending=True)\ndf.head()\n```\n\n::: {.cell-output .cell-output-display execution_count=2}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)MagnitudeDateTime
11346-17.834928.534626.01.62017-03-09 23:44:06+00:00
11343-17.845928.549127.02.02017-03-10 00:16:10+00:00
11344-17.886328.500820.02.12017-03-10 00:16:11+00:00
11345-17.865728.520430.01.62017-03-10 03:20:26+00:00
11342-17.715628.59850.01.62017-08-21 02:06:55+00:00
\n
\n```\n:::\n:::\n\n\n::: {#e8ff1eca .cell execution_count=3}\n``` {}\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=3}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#64dcb61f .cell execution_count=4}\n``` {}\ndf.plot.scatter(x=\"Longitude\", y=\"Latitude\", figsize=(12,12), grid=\"on\");\n```\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-html/cell-5-output-1.png){width=978 height=947}\n:::\n:::\n\n\n:::{#ff18cec9 .cell .markdown}\n### Filter \n\nAlso some rows seem to be missing depth measurements, (see differences in the **count** for columns above) - take only rows with valid depths\n:::\n\n::: {#a3612d29 .cell execution_count=5}\n``` {}\ndf = df[df['Depth(km)'].notna()]\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=5}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#265563e8 .cell execution_count=6}\n``` {}\nplt.figure(figsize=(20,6))\nplt.subplot(1,3,1)\ndf[\"Latitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,2)\ndf[\"Longitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,3)\ndf[\"Depth(km)\"].plot.hist(bins=50, density=True)\n```\n\n::: {.cell-output .cell-output-display execution_count=6}\n```\n\n```\n:::\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-html/cell-7-output-2.png){width=1565 height=485}\n:::\n:::\n\n\n:::{#387bb661 .cell .markdown}\n## Spatial Plot\n\nScatter plot the spatial locations of events\n:::\n\n::: {#cell-fig-spatial-plot .cell execution_count=7}\n``` {}\n#| label: fig-spatial-plot\n#| fig-cap: Locations of earthquakes on La Palma since 2017.\n#| fig-alt: A scatterplot of earthquake locations plotting latitude against longitude.\nfrom matplotlib import colormaps\ncmap = colormaps['viridis_r']\nax = df.plot.scatter(x=\"Longitude\", y=\"Latitude\", \n s=47-df[\"Depth(km)\"], c=df[\"Magnitude\"], \n figsize=(12,10), grid=\"on\", cmap=cmap)\ncolorbar = ax.collections[0].colorbar\ncolorbar.set_label(\"Magnitude\")\n\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Locations of earthquakes on La Palma since 2017.](data-screening_files/figure-html/fig-spatial-plot-output-1.png){#fig-spatial-plot width=917 height=799 fig-alt='A scatterplot of earthquake locations plotting latitude against longitude.'}\n:::\n:::\n\n\n:::{#d3cde14f .cell .markdown}\n### Timeline Plot\n\nScatter plot the event time series and look for any quantization issues. Have times & dates been loaded correctly?\n\n\nax = df.plot.scatter(x='DateTime', y='Depth(km)', figsize=(20,8))\nax.set_ylim(50,0);\n:::\n\n", + "markdown": "---\ntitle: Data Screening\ndescription: Screening and filtering data from the IGN catalog down to the relevant events for La Palma\nauthors:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n - name: Charlotte Wickham\n note: Edited original notebook by Steve Purves\n---\n\n::: {#791d860d .cell execution_count=1}\n``` {}\nimport pandas as pd\nimport matplotlib.pyplot as plt\n```\n:::\n\n\n:::{#93584cbe .cell .markdown}\nImport `la-palma.csv` - a clean subset of [`catalogoComunSV_1663233588717.csv`](https://github.com/Notebooks-Now/submission-quarto-full/blob/main/data/catalogoComunSV_1663233588717.csv):\n:::\n\n::: {#895af1da .cell execution_count=2}\n``` {}\ndf = pd.read_csv('la-palma.csv')\ndf['DateTime'] = pd.to_datetime(df['DateTime'])\ndf = df.sort_values(by=['DateTime'], ascending=True)\ndf.head()\n```\n\n::: {.cell-output .cell-output-display execution_count=16}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)MagnitudeDateTime
11346-17.834928.534626.01.62017-03-09 23:44:06+00:00
11343-17.845928.549127.02.02017-03-10 00:16:10+00:00
11344-17.886328.500820.02.12017-03-10 00:16:11+00:00
11345-17.865728.520430.01.62017-03-10 03:20:26+00:00
11342-17.715628.59850.01.62017-08-21 02:06:55+00:00
\n
\n```\n:::\n:::\n\n\n::: {#ab73c8ba .cell execution_count=3}\n``` {}\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=17}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#217b3449 .cell execution_count=4}\n``` {}\ndf.plot.scatter(x=\"Longitude\", y=\"Latitude\", figsize=(12,12), grid=\"on\");\n```\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-html/cell-5-output-1.png){width=978 height=947}\n:::\n:::\n\n\n:::{#5829bcef .cell .markdown}\n### Filter \n\nAlso some rows seem to be missing depth measurements, (see differences in the **count** for columns above) - take only rows with valid depths\n:::\n\n::: {#4246bd08 .cell execution_count=5}\n``` {}\ndf = df[df['Depth(km)'].notna()]\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=19}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#d0f22c6c .cell execution_count=6}\n``` {}\nplt.figure(figsize=(20,6))\nplt.subplot(1,3,1)\ndf[\"Latitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,2)\ndf[\"Longitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,3)\ndf[\"Depth(km)\"].plot.hist(bins=50, density=True)\n```\n\n::: {.cell-output .cell-output-display execution_count=20}\n```\n\n```\n:::\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-html/cell-7-output-2.png){width=1565 height=485}\n:::\n:::\n\n\n:::{#e1a22c58 .cell .markdown}\n## Spatial Plot\n\nScatter plot the spatial locations of events\n:::\n\n::: {#cell-fig-spatial-plot .cell execution_count=7}\n``` {}\n#| label: fig-spatial-plot\n#| fig-cap: Locations of earthquakes on La Palma since 2017.\n#| fig-alt: A scatterplot of earthquake locations plotting latitude against longitude.\nfrom matplotlib import colormaps\ncmap = colormaps['viridis_r']\nax = df.plot.scatter(x=\"Longitude\", y=\"Latitude\", \n s=47-df[\"Depth(km)\"], c=df[\"Magnitude\"], \n figsize=(12,10), grid=\"on\", cmap=cmap)\ncolorbar = ax.collections[0].colorbar\ncolorbar.set_label(\"Magnitude\")\n\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Locations of earthquakes on La Palma since 2017.](data-screening_files/figure-html/fig-spatial-plot-output-1.png){#fig-spatial-plot width=917 height=799 fig-alt='A scatterplot of earthquake locations plotting latitude against longitude.'}\n:::\n:::\n\n\n:::{#70dd56c7 .cell .markdown}\n### Timeline Plot\n\nScatter plot the event time series and look for any quantization issues. Have times & dates been loaded correctly?\n\n\nax = df.plot.scatter(x='DateTime', y='Depth(km)', figsize=(20,8))\nax.set_ylim(50,0);\n:::\n\n", "supporting": [ "data-screening_files" ], diff --git a/_freeze/notebooks/data-screening/execute-results/xml.json b/_freeze/notebooks/data-screening/execute-results/xml.json index 6669413..6d4c893 100644 --- a/_freeze/notebooks/data-screening/execute-results/xml.json +++ b/_freeze/notebooks/data-screening/execute-results/xml.json @@ -1,7 +1,7 @@ { "hash": "d87b6d9ac2228b736962871bfd9167c0", "result": { - "markdown": "---\ntitle: Data Screening\ndescription: Screening and filtering data from the IGN catalog down to the relevant events for La Palma\nauthors:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n - name: Charlotte Wickham\n note: Edited original notebook by Steve Purves\n---\n\n::: {#e3a1133d .cell execution_count=1}\n``` {}\nimport pandas as pd\nimport matplotlib.pyplot as plt\n```\n:::\n\n\n:::{#21324404 .cell .markdown}\nImport `la-palma.csv` - a clean subset of [`catalogoComunSV_1663233588717.csv`](https://github.com/Notebooks-Now/submission-quarto-full/blob/main/data/catalogoComunSV_1663233588717.csv):\n:::\n\n::: {#3588ccfa .cell execution_count=2}\n``` {}\ndf = pd.read_csv('la-palma.csv')\ndf['DateTime'] = pd.to_datetime(df['DateTime'])\ndf = df.sort_values(by=['DateTime'], ascending=True)\ndf.head()\n```\n\n::: {.cell-output .cell-output-display execution_count=9}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)MagnitudeDateTime
11346-17.834928.534626.01.62017-03-09 23:44:06+00:00
11343-17.845928.549127.02.02017-03-10 00:16:10+00:00
11344-17.886328.500820.02.12017-03-10 00:16:11+00:00
11345-17.865728.520430.01.62017-03-10 03:20:26+00:00
11342-17.715628.59850.01.62017-08-21 02:06:55+00:00
\n
\n```\n:::\n:::\n\n\n::: {#cf8656b6 .cell execution_count=3}\n``` {}\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=10}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#204e1533 .cell execution_count=4}\n``` {}\ndf.plot.scatter(x=\"Longitude\", y=\"Latitude\", figsize=(12,12), grid=\"on\");\n```\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-jats/cell-5-output-1.png){}\n:::\n:::\n\n\n:::{#ebdb3a1c .cell .markdown}\n### Filter \n\nAlso some rows seem to be missing depth measurements, (see differences in the **count** for columns above) - take only rows with valid depths\n:::\n\n::: {#e3111732 .cell execution_count=5}\n``` {}\ndf = df[df['Depth(km)'].notna()]\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=12}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#528f6a0d .cell execution_count=6}\n``` {}\nplt.figure(figsize=(20,6))\nplt.subplot(1,3,1)\ndf[\"Latitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,2)\ndf[\"Longitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,3)\ndf[\"Depth(km)\"].plot.hist(bins=50, density=True)\n```\n\n::: {.cell-output .cell-output-display execution_count=13}\n```\n\n```\n:::\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-jats/cell-7-output-2.png){}\n:::\n:::\n\n\n:::{#6bcba60c .cell .markdown}\n## Spatial Plot\n\nScatter plot the spatial locations of events\n:::\n\n::: {#cell-fig-spatial-plot .cell execution_count=7}\n``` {}\n#| label: fig-spatial-plot\n#| fig-cap: Locations of earthquakes on La Palma since 2017.\n#| fig-alt: A scatterplot of earthquake locations plotting latitude against longitude.\nfrom matplotlib import colormaps\ncmap = colormaps['viridis_r']\nax = df.plot.scatter(x=\"Longitude\", y=\"Latitude\", \n s=47-df[\"Depth(km)\"], c=df[\"Magnitude\"], \n figsize=(12,10), grid=\"on\", cmap=cmap)\ncolorbar = ax.collections[0].colorbar\ncolorbar.set_label(\"Magnitude\")\n\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Locations of earthquakes on La Palma since 2017.](data-screening_files/figure-jats/fig-spatial-plot-output-1.png){#fig-spatial-plot fig-alt='A scatterplot of earthquake locations plotting latitude against longitude.'}\n:::\n:::\n\n\n:::{#b9129c3f .cell .markdown}\n### Timeline Plot\n\nScatter plot the event time series and look for any quantization issues. Have times & dates been loaded correctly?\n\n\nax = df.plot.scatter(x='DateTime', y='Depth(km)', figsize=(20,8))\nax.set_ylim(50,0);\n:::\n\n", + "markdown": "---\ntitle: Data Screening\ndescription: Screening and filtering data from the IGN catalog down to the relevant events for La Palma\nauthors:\n - name: Steve Purves\n orcid: 0000-0002-0760-5497\n - name: Charlotte Wickham\n note: Edited original notebook by Steve Purves\n---\n\n::: {#43daa9a5 .cell execution_count=1}\n``` {}\nimport pandas as pd\nimport matplotlib.pyplot as plt\n```\n:::\n\n\n:::{#25be6d6d .cell .markdown}\nImport `la-palma.csv` - a clean subset of [`catalogoComunSV_1663233588717.csv`](https://github.com/Notebooks-Now/submission-quarto-full/blob/main/data/catalogoComunSV_1663233588717.csv):\n:::\n\n::: {#7e7b2b9e .cell execution_count=2}\n``` {}\ndf = pd.read_csv('la-palma.csv')\ndf['DateTime'] = pd.to_datetime(df['DateTime'])\ndf = df.sort_values(by=['DateTime'], ascending=True)\ndf.head()\n```\n\n::: {.cell-output .cell-output-display execution_count=23}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)MagnitudeDateTime
11346-17.834928.534626.01.62017-03-09 23:44:06+00:00
11343-17.845928.549127.02.02017-03-10 00:16:10+00:00
11344-17.886328.500820.02.12017-03-10 00:16:11+00:00
11345-17.865728.520430.01.62017-03-10 03:20:26+00:00
11342-17.715628.59850.01.62017-08-21 02:06:55+00:00
\n
\n```\n:::\n:::\n\n\n::: {#cca50ca5 .cell execution_count=3}\n``` {}\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=24}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#3f3b8176 .cell execution_count=4}\n``` {}\ndf.plot.scatter(x=\"Longitude\", y=\"Latitude\", figsize=(12,12), grid=\"on\");\n```\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-jats/cell-5-output-1.png){}\n:::\n:::\n\n\n:::{#319030b4 .cell .markdown}\n### Filter \n\nAlso some rows seem to be missing depth measurements, (see differences in the **count** for columns above) - take only rows with valid depths\n:::\n\n::: {#f02e8517 .cell execution_count=5}\n``` {}\ndf = df[df['Depth(km)'].notna()]\ndf.describe()\n```\n\n::: {.cell-output .cell-output-display execution_count=26}\n```{=html}\n
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
LongitudeLatitudeDepth(km)Magnitude
count11347.00000011347.00000011347.00000011347.000000
mean-17.84406228.56889314.0893452.432423
std0.0273450.0218998.6247870.705282
min-17.99580028.3102000.0000000.200000
25%-17.85410028.55800010.1000001.900000
50%-17.83870028.56400011.3000002.600000
75%-17.82970028.57430013.9000002.900000
max-17.70570028.69780046.9000005.100000
\n
\n```\n:::\n:::\n\n\n::: {#f31daa87 .cell execution_count=6}\n``` {}\nplt.figure(figsize=(20,6))\nplt.subplot(1,3,1)\ndf[\"Latitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,2)\ndf[\"Longitude\"].plot.hist(bins=20, density=True)\nplt.subplot(1,3,3)\ndf[\"Depth(km)\"].plot.hist(bins=50, density=True)\n```\n\n::: {.cell-output .cell-output-display execution_count=27}\n```\n\n```\n:::\n\n::: {.cell-output .cell-output-display}\n![](data-screening_files/figure-jats/cell-7-output-2.png){}\n:::\n:::\n\n\n:::{#2e9530fb .cell .markdown}\n## Spatial Plot\n\nScatter plot the spatial locations of events\n:::\n\n::: {#cell-fig-spatial-plot .cell execution_count=7}\n``` {}\n#| label: fig-spatial-plot\n#| fig-cap: Locations of earthquakes on La Palma since 2017.\n#| fig-alt: A scatterplot of earthquake locations plotting latitude against longitude.\nfrom matplotlib import colormaps\ncmap = colormaps['viridis_r']\nax = df.plot.scatter(x=\"Longitude\", y=\"Latitude\", \n s=47-df[\"Depth(km)\"], c=df[\"Magnitude\"], \n figsize=(12,10), grid=\"on\", cmap=cmap)\ncolorbar = ax.collections[0].colorbar\ncolorbar.set_label(\"Magnitude\")\n\nplt.show()\n```\n\n::: {.cell-output .cell-output-display}\n![Locations of earthquakes on La Palma since 2017.](data-screening_files/figure-jats/fig-spatial-plot-output-1.png){#fig-spatial-plot fig-alt='A scatterplot of earthquake locations plotting latitude against longitude.'}\n:::\n:::\n\n\n:::{#f0aaa739 .cell .markdown}\n### Timeline Plot\n\nScatter plot the event time series and look for any quantization issues. Have times & dates been loaded correctly?\n\n\nax = df.plot.scatter(x='DateTime', y='Depth(km)', figsize=(20,8))\nax.set_ylim(50,0);\n:::\n\n", "supporting": [ "data-screening_files/figure-jats" ], diff --git a/_freeze/site_libs/bootstrap/bootstrap-icons.css b/_freeze/site_libs/bootstrap/bootstrap-icons.css index c3bfba6..285e444 100644 --- a/_freeze/site_libs/bootstrap/bootstrap-icons.css +++ b/_freeze/site_libs/bootstrap/bootstrap-icons.css @@ -1,5 +1,5 @@ /*! - * Bootstrap Icons v1.10.5 (https://icons.getbootstrap.com/) + * Bootstrap Icons v1.11.1 (https://icons.getbootstrap.com/) * Copyright 2019-2023 The Bootstrap Authors * Licensed under MIT (https://github.com/twbs/icons/blob/main/LICENSE) */ @@ -8,7 +8,7 @@ font-display: block; font-family: "bootstrap-icons"; src: -url("./bootstrap-icons.woff?1fa40e8900654d2863d011707b9fb6f2") format("woff"); +url("./bootstrap-icons.woff?2820a3852bdb9a5832199cc61cec4e65") format("woff"); } .bi::before, @@ -1979,3 +1979,100 @@ url("./bootstrap-icons.woff?1fa40e8900654d2863d011707b9fb6f2") format("woff"); .bi-sina-weibo::before { content: "\f8ca"; } .bi-tencent-qq::before { content: "\f8cb"; } .bi-wikipedia::before { content: "\f8cc"; } +.bi-alphabet-uppercase::before { content: "\f2a5"; } +.bi-alphabet::before { content: "\f68a"; } +.bi-amazon::before { content: "\f68d"; } +.bi-arrows-collapse-vertical::before { content: "\f690"; } +.bi-arrows-expand-vertical::before { content: "\f695"; } +.bi-arrows-vertical::before { content: "\f698"; } +.bi-arrows::before { content: "\f6a2"; } +.bi-ban-fill::before { content: "\f6a3"; } +.bi-ban::before { content: "\f6b6"; } +.bi-bing::before { content: "\f6c2"; } +.bi-cake::before { content: "\f6e0"; } +.bi-cake2::before { content: "\f6ed"; } +.bi-cookie::before { content: "\f6ee"; } +.bi-copy::before { content: "\f759"; } +.bi-crosshair::before { content: "\f769"; } +.bi-crosshair2::before { content: "\f794"; } +.bi-emoji-astonished-fill::before { content: "\f795"; } +.bi-emoji-astonished::before { content: "\f79a"; } +.bi-emoji-grimace-fill::before { content: "\f79b"; } +.bi-emoji-grimace::before { content: "\f7a0"; } +.bi-emoji-grin-fill::before { content: "\f7a1"; } +.bi-emoji-grin::before { content: "\f7a6"; } +.bi-emoji-surprise-fill::before { content: "\f7a7"; } +.bi-emoji-surprise::before { content: "\f7ac"; } +.bi-emoji-tear-fill::before { content: "\f7ad"; } +.bi-emoji-tear::before { content: "\f7b2"; } +.bi-envelope-arrow-down-fill::before { content: "\f7b3"; } +.bi-envelope-arrow-down::before { content: "\f7b8"; } +.bi-envelope-arrow-up-fill::before { content: "\f7b9"; } +.bi-envelope-arrow-up::before { content: "\f7be"; } +.bi-feather::before { content: "\f7bf"; } +.bi-feather2::before { content: "\f7c4"; } +.bi-floppy-fill::before { content: "\f7c5"; } +.bi-floppy::before { content: "\f7d8"; } +.bi-floppy2-fill::before { content: "\f7d9"; } +.bi-floppy2::before { content: "\f7e4"; } +.bi-gitlab::before { content: "\f7e5"; } +.bi-highlighter::before { content: "\f7f8"; } +.bi-marker-tip::before { content: "\f802"; } +.bi-nvme-fill::before { content: "\f803"; } +.bi-nvme::before { content: "\f80c"; } +.bi-opencollective::before { content: "\f80d"; } +.bi-pci-card-network::before { content: "\f8cd"; } +.bi-pci-card-sound::before { content: "\f8ce"; } +.bi-radar::before { content: "\f8cf"; } +.bi-send-arrow-down-fill::before { content: "\f8d0"; } +.bi-send-arrow-down::before { content: "\f8d1"; } +.bi-send-arrow-up-fill::before { content: "\f8d2"; } +.bi-send-arrow-up::before { content: "\f8d3"; } +.bi-sim-slash-fill::before { content: "\f8d4"; } +.bi-sim-slash::before { content: "\f8d5"; } +.bi-sourceforge::before { content: "\f8d6"; } +.bi-substack::before { content: "\f8d7"; } +.bi-threads-fill::before { content: "\f8d8"; } +.bi-threads::before { content: "\f8d9"; } +.bi-transparency::before { content: "\f8da"; } +.bi-twitter-x::before { content: "\f8db"; } +.bi-type-h4::before { content: "\f8dc"; } +.bi-type-h5::before { content: "\f8dd"; } +.bi-type-h6::before { content: "\f8de"; } +.bi-backpack-fill::before { content: "\f8df"; } +.bi-backpack::before { content: "\f8e0"; } +.bi-backpack2-fill::before { content: "\f8e1"; } +.bi-backpack2::before { content: "\f8e2"; } +.bi-backpack3-fill::before { content: "\f8e3"; } +.bi-backpack3::before { content: "\f8e4"; } +.bi-backpack4-fill::before { content: "\f8e5"; } +.bi-backpack4::before { content: "\f8e6"; } +.bi-brilliance::before { content: "\f8e7"; } +.bi-cake-fill::before { content: "\f8e8"; } +.bi-cake2-fill::before { content: "\f8e9"; } +.bi-duffle-fill::before { content: "\f8ea"; } +.bi-duffle::before { content: "\f8eb"; } +.bi-exposure::before { content: "\f8ec"; } +.bi-gender-neuter::before { content: "\f8ed"; } +.bi-highlights::before { content: "\f8ee"; } +.bi-luggage-fill::before { content: "\f8ef"; } +.bi-luggage::before { content: "\f8f0"; } +.bi-mailbox-flag::before { content: "\f8f1"; } +.bi-mailbox2-flag::before { content: "\f8f2"; } +.bi-noise-reduction::before { content: "\f8f3"; } +.bi-passport-fill::before { content: "\f8f4"; } +.bi-passport::before { content: "\f8f5"; } +.bi-person-arms-up::before { content: "\f8f6"; } +.bi-person-raised-hand::before { content: "\f8f7"; } +.bi-person-standing-dress::before { content: "\f8f8"; } +.bi-person-standing::before { content: "\f8f9"; } +.bi-person-walking::before { content: "\f8fa"; } +.bi-person-wheelchair::before { content: "\f8fb"; } +.bi-shadows::before { content: "\f8fc"; } +.bi-suitcase-fill::before { content: "\f8fd"; } +.bi-suitcase-lg-fill::before { content: "\f8fe"; } +.bi-suitcase-lg::before { content: "\f8ff"; } +.bi-suitcase::before { content: "\f900"; } +.bi-suitcase2-fill::before { content: "\f901"; } +.bi-suitcase2::before { content: "\f902"; } +.bi-vignette::before { content: "\f903"; } diff --git a/_freeze/site_libs/bootstrap/bootstrap-icons.woff b/_freeze/site_libs/bootstrap/bootstrap-icons.woff index 6e72a59..dbeeb05 100644 Binary files a/_freeze/site_libs/bootstrap/bootstrap-icons.woff and b/_freeze/site_libs/bootstrap/bootstrap-icons.woff differ diff --git a/_freeze/site_libs/quarto-html/quarto.js b/_freeze/site_libs/quarto-html/quarto.js index c3935c7..91be522 100644 --- a/_freeze/site_libs/quarto-html/quarto.js +++ b/_freeze/site_libs/quarto-html/quarto.js @@ -18,25 +18,12 @@ const layoutMarginEls = () => { // clear the top margin so we recompute it marginChild.style.marginTop = null; const top = marginChild.getBoundingClientRect().top + window.scrollY; - console.log({ - childtop: marginChild.getBoundingClientRect().top, - scroll: window.scrollY, - top, - lastBottom, - }); if (top < lastBottom) { const margin = lastBottom - top; marginChild.style.marginTop = `${margin}px`; } const styles = window.getComputedStyle(marginChild); const marginTop = parseFloat(styles["marginTop"]); - - console.log({ - top, - height: marginChild.getBoundingClientRect().height, - marginTop, - total: top + marginChild.getBoundingClientRect().height + marginTop, - }); lastBottom = top + marginChild.getBoundingClientRect().height + marginTop; } }