{"version":"desci-nodes-0.2.0","components":[{"id":"root","name":"root","type":"data","payload":{"cid":"bafybeig5rcvd6qmpzjkpm3s7y7yubjhb7mspvvhtu5yfb6h7i5erburfea","path":"root"}},{"id":"9aedbec2-a910-4de9-9568-9dee8e354ea3","name":"cremons-et-al-2022","type":"code","payload":{"title":"Software for Simulating Lunar Surface Hydration Measurements for Multispectral Lidar at 3 µm","licenseType":"creative commons zero v1.0 universal","description":"Python reproduction of Cremons, Daniel R., “Software for Simulating Lunar Surface Hydration Measurements for Multispectral Lidar at 3 µm”. Zenodo, Jul. 01, 2022. doi: 10.5281/zenodo.6025377.","path":"root/cremons-et-al-2022","cid":"pending"}},{"id":"0b998a5a-abae-4e9e-97df-df66be47e98b","name":"cremons-et-al-2022","type":"code","payload":{"url":"bafybeie7fo4b7ofoqs27gnkd4ld7e2gdllqxn6wxxz7dusrixgl7sisrxm","path":"root/cremons-et-al-2022","externalUrl":"https://github.com/philiplinden/cremons-et-al-2022"},"starred":false},{"id":"85d347ca-cc3e-4d56-9051-ad186860199e","name":"Earth and Space Science - 2022 - Cremons.pdf","type":"pdf","payload":{"annotations":[{"startX":0.2780753902230222,"startY":0.49195640075382924,"endX":0.5088727563078891,"endY":0.509197780064174,"pageIndex":0,"id":"3b8a934f-5798-47eb-8e45-bc20f5f7dce4","text":"The DeSci node published as a reproduction of this work using Python instead of Matlab. The goal is to reproduce the figures and results from the data published by Cremons, et. al that accompanied their paper.\n\n\nThe goal of this project is to learn by doing, and to reproduce the figures in a manner than is more broadly accessible. Comments and opinions expressed in the code and supporting text in this notebook are my own.\n\n\n- Philip Linden, MoonDAO\n","rects":[{"startX":0.2780753902230222,"startY":0.49195640075382924,"endX":0.5088727563078891,"endY":0.509197780064174,"pageIndex":0}]},{"startX":0.05068948567423321,"startY":0.515983048330974,"endX":0.2454383963900261,"endY":0.5829584962628743,"pageIndex":0,"id":"05b393c0-162c-4a8d-852b-78f91087ffdd","text":"> Cremons, D. R., & Honniball, C. I.(2022). Simulated lunar surface hydrationmeasurements using multispectral lidarat 3 μm. Earth and Space Science,9, e2022EA002277. https://doi.org/10.1029/2022EA002277\n\n\n\nThis work was originally published in _Earth and Space Science in 2022_(doi: 10.5281/zenodo.6025377). The manuscript, data, and Matlab code is included unaltered from the originally published materials in the `src` directory.\n\n\nFigures and algorithms presented in this manuscript by Cremons, et. al. were reproduced with Python code written by Philip Linden. The reproduction is included in the `repro` directory.\n\n","rects":[{"startX":0.05069106275010837,"startY":0.515983048330974,"endX":0.2140751870279488,"endY":0.5258352650797424,"pageIndex":0},{"startX":0.05069421690185869,"startY":0.5274089498472917,"endX":0.2454383963900261,"endY":0.5372611665960602,"pageIndex":0},{"startX":0.05069096418286617,"startY":0.5388313561237503,"endX":0.23344966272153939,"endY":0.5486835728725188,"pageIndex":0},{"startX":0.05069549827600726,"startY":0.5502614669611887,"endX":0.24438875382786446,"endY":0.5616838732376472,"pageIndex":0},{"startX":0.05068948567423321,"startY":0.5616838732376472,"endX":0.2443757429518944,"endY":0.5715360899864158,"pageIndex":0},{"startX":0.05069648394842924,"startY":0.5731062795141059,"endX":0.18414977802780708,"endY":0.5829584962628743,"pageIndex":0}]},{"startX":0.05068715521377259,"startY":0.6419472975449843,"endX":0.1634057142949295,"endY":0.651537743719903,"pageIndex":0,"id":"f972499a-8ca5-4f0d-a75e-d0c31a11137f","text":"Philip Linden independently reproduced the results and adapted code to Python for MoonDAO as a Creative Commons Non Commercial work.\n","rects":[{"startX":0.05068715521377259,"startY":0.6419472975449843,"endX":0.1634057142949295,"endY":0.651537743719903,"pageIndex":0}]},{"startX":0.574184988063854,"startY":0.9010709189429073,"endX":0.8825266820941582,"endY":0.9119083799164872,"pageIndex":3,"id":"32c75b85-d89e-4b25-89d3-e8957f4845ee","text":"> Figure 1. (a) Regolith and pyroxene spectra as described in Table 1.\n\n\nThis figure was reproduced in Code Cell [4] of the IPy notebook.\n\n[results.ipynb](#/code/3e71b606-f76a-497e-90df-f3d6011ad09a?path=root%2Fcremons-et-al-2022%2Frepro%2Fresults.ipynb)\n\nSo far, so good. The imported datasets match what is shown in Figure 1a of the manuscript after limiting the x-axis to the spectrum of interest.\n\nIt's not clear to me what the units of reflectance are in this data, as it is unlabeled in the source data files. I looked around briefly in the source dataset ([NASA RELAB](https://sites.brown.edu/relab/lunar-rock-and-mineral-characterization-database/)) and related publications ([source](https://www.researchgate.net/publication/253034096_RELAB_Reflectance_Experiment_Laboratory_A_NASA_multiuser_spectroscopy_facility)) but didn't see anything obvious. This does not undermine the analysis, though.\n\nI observed the features described in the manuscript in the plotted data.\n","rects":[{"startX":0.574184988063854,"startY":0.9010709189429073,"endX":0.8825266820941582,"endY":0.9119083799164872,"pageIndex":3}]},{"startX":0.27784375720385796,"startY":0.9010709565261312,"endX":0.909098143544044,"endY":0.9372932170999458,"pageIndex":3,"id":"cb505f0f-2103-4b72-96b8-da2575ee34b7","text":"> Figure 1. (b) Laboratory MORB glass reflectance spectra used to simulate surface hydration in the reflectance simulations. The rawMORB spectra were normalized to each other at a wavelength 2.6 μm.\n\n\n\nThis figure was reproduced in Code Cell [10] of the IPy notebook.\n\n[results.ipynb](#/code/3e71b606-f76a-497e-90df-f3d6011ad09a?path=root%2Fcremons-et-al-2022%2Frepro%2Fresults.ipynb)\n\nThe MORB reflectance over the spectrum of interest (1-4 micron) looks very different from the \"processed\" traces in Figure 1b. In the authors' code the MORB spectra are normalized to the 22 ppm reflectance spectrum. The \"raw\" heated MORB spectra also contains strange (to me) features below 1.5 microns.\n\nThe manuscript attributes these features to the presence of unquantified hydration and organic doublets in this part of the spectrum, and substitutes parts of the reflectance spectra with linearly fitted data.\n\nThe first data processing operation that is performed is to replace the spectra below 1.6 micron with another reflectance spectrum. It's a bit strange to _replace_ a portion of the observations with other observations, but in this case it feels analagous to an HDR composite photo.\n\nApparent artifacts had to be removed below 2.6 micron for the figure to match the manuscript. This is done by simply replacing the data below 2.6 micron with the least hydrated reflectance spectrum, 22 ppm.\n\nEven the 22 ppm MORB spectrum has odd behavior between 1.0-1.5 microns, though. This region of the spectrum is replaced with the observations that focused on 1.0-1.6 microns (`Morb_D38A_Low_wavelength.txt`), which also has to be normalized to mesh with the 22 ppm reflectance spectrum at the boundary between those datasets.\n\nWhile there is some pretty heavy processing done on the spectra below 2.6 micron, the focus of this study is the 3-4 micron range. The processing is done to show relative characteristics between spectra without including other features that may distract from the analysis itself.","rects":[{"startX":0.8858470182148243,"startY":0.9010709565261312,"endX":0.901368500410913,"endY":0.9119083799164868,"pageIndex":3},{"startX":0.27784769989354585,"startY":0.9137633751178608,"endX":0.909098143544044,"endY":0.9246007985082163,"pageIndex":3},{"startX":0.27784375720385796,"startY":0.9264557937095903,"endX":0.6521845981010015,"endY":0.9372932170999458,"pageIndex":3}]},{"startX":0.05069323122943674,"startY":0.34495691008168494,"endX":0.4722735564730364,"endY":0.3938587733677456,"pageIndex":4,"id":"c2f8eee5-5138-4842-a645-ca74a9d1f772","text":"> Figure 2. A Subset of interpolated MORB single scattering albedo spectra.The interpolated spectra correspond to total water values of (from low tohigh): 0, 150, 300, 450, 600, 750, 900, 1,050, 1,200 ppm, and 1,350, and1,500 ppm.\n\n\n\nThis figure was reproduced in Code Cell [16] of the IPy notebook.\n\n[results.ipynb](#/code/3e71b606-f76a-497e-90df-f3d6011ad09a?path=root%2Fcremons-et-al-2022%2Frepro%2Fresults.ipynb)\n","rects":[{"startX":0.05069796245706222,"startY":0.34495691008168494,"endX":0.4722735564730364,"endY":0.3569180718783675,"pageIndex":4},{"startX":0.05069323122943674,"startY":0.3576441797716865,"endX":0.43733827025372035,"endY":0.3684816031620421,"pageIndex":4},{"startX":0.050695695410491674,"startY":0.3703366359466403,"endX":0.4374921337187907,"endY":0.3811740593369959,"pageIndex":4},{"startX":0.05069426618547981,"startY":0.38302134997739,"endX":0.11153844485263015,"endY":0.3938587733677456,"pageIndex":4}]},{"startX":0.2778382374382949,"startY":0.6073502601661117,"endX":0.9492457640980082,"endY":0.6506076493286732,"pageIndex":4,"id":"48666a4a-93d1-469e-b690-8998ca41f5b2","text":"> To create the spectral mixtures, each endmember reflectance spectrum was first converted to single-scattering albedo (SSA) using the methods of Hapke (2012) and least-squares minimization.\n\n\nThe Matlab implementation of this method had to be modified slightly in order to reproduce the author's results. The Matlab code that reproduced the author's findings is included here:\n\n[Hydrated_Regolith_Spectra_Generator_and_Retrieval.m](#/code/5acc09bf-f463-42f3-a378-feb197a39379?path=root%2Fcremons-et-al-2022%2Frepro%2FHydrated_Regolith_Spectra_Generator_and_Retrieval.m)\n","rects":[{"startX":0.6552378169953125,"startY":0.6073502601661117,"endX":0.9491012645209466,"endY":0.6201581344228658,"pageIndex":4},{"startX":0.2778382374382949,"startY":0.6225750364106275,"endX":0.9492457640980082,"endY":0.6353829106673816,"pageIndex":4},{"startX":0.27783892740899024,"startY":0.6377997750719191,"endX":0.4689681040059418,"endY":0.6506076493286732,"pageIndex":4}]},{"startX":0.2778438557711002,"startY":0.6530245395716776,"endX":0.9485816180200815,"endY":0.6810571524897233,"pageIndex":4,"id":"867d4161-1831-499b-b654-e31f2a10d470","text":"> the SSA in the Hapkemodel by (Lucey, 1998)\n\n\nThis model was successfully reproduced in Python.\n\n[hapke.py](#/code/0b998a5a-abae-4e9e-97df-df66be47e98b?path=repro%2Fhapke.py)\n","rects":[{"startX":0.8211477761380862,"startY":0.6530245395716776,"endX":0.9485816180200815,"endY":0.6658324138284317,"pageIndex":4},{"startX":0.2778438557711002,"startY":0.6682493158161935,"endX":0.42523184729804975,"endY":0.6810571524897233,"pageIndex":4}]},{"startX":0.05069559684324946,"startY":0.5875856841139017,"endX":0.9269636533268324,"endY":0.6111078967014555,"pageIndex":14,"id":"885436b1-8130-448c-931c-482505b421eb","text":"> (b) Comparison of interpolated MORB spectra (obtained from a linear ESPAT fit) and laboratory MORB spectra at the same total watercontent for each. The total water content for each laboratory spectrum is given in the legend.\n\n\n\nThis figure is reproduced in Code Cell [17] of the IPy notebook.\n\n[results.ipynb](#/code/3e71b606-f76a-497e-90df-f3d6011ad09a?path=root%2Fcremons-et-al-2022%2Frepro%2Fresults.ipynb)\n","rects":[{"startX":0.2026795324557647,"startY":0.5875856841139017,"endX":0.9269636533268324,"endY":0.5984231075042573,"pageIndex":14},{"startX":0.05069559684324946,"startY":0.6002703981446513,"endX":0.5425196701044411,"endY":0.6111078967014555,"pageIndex":14}]},{"startX":0.2778376460348418,"startY":0.8241605993562146,"endX":0.49293679885835984,"endY":0.839431565383385,"pageIndex":15,"id":"a3306f07-a1f7-4e07-ae3f-2c9a6c68af4c","text":"The complete dataset originally published with this manuscript in _Earth and Space Science 2022_ is included in this node in the `data` directory.\n","rects":[{"startX":0.2778376460348418,"startY":0.8241605993562146,"endX":0.49293679885835984,"endY":0.839431565383385,"pageIndex":15}]}],"path":"root/cremons-et-al-2022/Earth and Space Science - 2022 - Cremons.pdf","cid":"bafybeig42ber4ai7swxzmakvsh7kh2taukr32rohvb5xfo5z77htzmuduy","title":"Simulated Lunar Surface Hydration Measurements Using Multispectral Lidar at 3 μm (Manuscript)","licenseType":"creative commons zero v1.0 universal"},"starred":true},{"id":"3e71b606-f76a-497e-90df-f3d6011ad09a","name":"results.ipynb","type":"code","payload":{"path":"root/cremons-et-al-2022/repro/results.ipynb","cid":"bafybeihzbh2jdwzka5unfwvpdlf6fyfwp4imtemekqhalcyet7yhcsp7zm","title":"Reproduced Results","licenseType":"mit license"},"starred":true},{"id":"5acc09bf-f463-42f3-a378-feb197a39379","name":"Hydrated_Regolith_Spectra_Generator_and_Retrieval.m","type":"code","payload":{"path":"root/cremons-et-al-2022/repro/Hydrated_Regolith_Spectra_Generator_and_Retrieval.m","cid":"bafkreigv3cimir4wpo5lqubwrskd3hmeomn3wvn5vsi4qg5rdjhkoecgmi"},"starred":true}],"authors":[{"name":"Philip Linden","orcid":"","role":"Author"},{"name":"Casey Honniball","orcid":"0000-0001-8248-8991","role":"Author"},{"name":"Daniel R. Cremons","orcid":"0000-0003-2211-1329","role":"Author"}],"title":"Software for Simulating Lunar Surface Hydration Measurements for Multispectral Lidar at 3 µm","defaultLicense":"CC BY","researchFields":["geology"],"dpid":{"prefix":"beta","id":"137"}}