31:T9aa,viewof nuclide = Inputs.select(sampleType == 1 ? new Map([ // select nuclides for seawater ["Cs-137", 33], ["C-14", 3], ["Tc-99", 15], ["Po-210", 47], ["Pu-239,240", 77], ["Pu-241", 70], ["Cs-134", 31], ['Pu-238', 67] ]) : sampleType == 2 ? biotaType == 1 ? new Map([ //select nuclides for fish ["Cs-137", 33], ["C-14", 3], ["Tc-99", 15], ["Po-210", 47], ["Pu-239,240", 77], ["Pu-241", 70], ["Cs-134", 31], ['Pu-238', 67]] ) : biotaType == 2 ? new Map([ // select nuclides for seaweed ["Cs-137", 33], ["C-14", 3], ["Tc-99", 15] ["Po-210", 47], ["Pu-239,240", 77], ["Pu-241", 70], ["Cs-134", 31], ['Pu-238', 67]] ) : new Map([ // select nuclides for shellfish ["Cs-137", 33], ["C-14", 3], ["Tc-99", 15], ["Po-210", 47], ["Pu-239,240", 77], ["Pu-241", 70], ["Cs-134", 31], ['Pu-238', 67]] ) : new Map([ // select nuclides for sediment ["Cs-137", 33], ["C-14", 3], ["Tc-99", 15], ["Po-210", 47], ["Pu-239,240", 77], ["Pu-241", 70], ["Cs-134", 31], ['Pu-238', 67]]), {label: html`

Nuclide Type:

`, sort: true, unique: true })32:T5f3,globalMap = Plot.plot({ width, color: { type: "threshold", scheme: "YlGnBu", domain: [0.25, 0.5, 0.75, 0.95, 0.99].map((d) => d3.quantile(valuesForLegend, d), ) }, projection: { type: "mercator", domain: circle }, marks: [ Plot.geo(land, { fill: "lightgrey", stroke: "Black", strokeWidth: 0.5 }), Plot.graticule(), Plot.geo(circle, { stroke: "#bbb", strokeWidth: 1, strokeDasharray: "5,3" }), Plot.dot( dataWithinRadius, Plot.hexbin( { fill: (g) => d3[aggType](g, (d) => d.activity), title: (g) => `${nuclidesLut[nuclide]} act. conc. (${aggType}): ${d3[aggType](g, (d) => d.activity)} ${whatUnit(sampleType)}` // `${nuclidesLut[nuclide]} activity in ${whatUnit(sampleType)}` }, { x: "longitude", y: "latitude", fillOpacity: 0.9, symbol: "hexagon", stroke: "black", strokeWidth: 0.5, strokeOpacity: 5, binWidth: 6, channels: { activity: ((d) => d.activity) }, tip: { fomart: { activity: ((d) => d.activity) } } } ) ), Plot.geo(EISPolygon), Plot.geo(WISPolygon), Plot.geo(WOIPolygon), Plot.dot(nuclearSitesdata, {x: "Longitude decimal", y: "Latitude decimal", symbol: "square", title: "Name", tip:true, r: 5, stroke: "red", fill: "orange"}), ] })33:T77b,timePlot1 = Plot.plot({ width, height: 500, // color: { scheme: "OrRd" }, // color: { // scheme: "Category10", // type: "categorical" // }, color: { domain: [ "Eastern region of the Irish Sea", "Western region of the Irish Sea", "West of Ireland" ], range: ["#1b9e77", "#7570b3", "#d95f02"] // Colorbind-safe }, y: {type: "log", grid: true, tickFormat: ".2" } , x: {label: "Sampling date"}, marks: [ Plot.dot( dataInBoxes, { x: "begPeriod", y: "activity", fill: "roi", stroke: "roi", r: 3, fillOpacity: 0.4, strokeWidth: 0.5, strokeOpacity: 0.4 } // Plot.hexbin( // { fill: "count" }, // { // x: "begPeriod", // y: "activity", // symbol: "hexagon", // fillOpacity: 0.9, // stroke: "black", // strokeWidth: 0.5, // strokeOpacity: 0.5, // binWidth: 7, // tip: true // } // ) ), Plot.lineY( dataInBoxes, Plot.windowY( { k: dataInBoxes.length > 0 ? dataInBoxes.length / ((dataInBoxes[dataInBoxes.length - 1].begPeriod - dataInBoxes[0].begPeriod) / (3600 * 24 * 365.25 * 1000)) : 1, reduce: "median" }, { x: "begPeriod", y: "activity", z: "roi", stroke: "roi", strokeWidth: 2 } ) ) // Plot.text( // dataInBoxes, // Plot.selectLast({ // x: "begPeriod", // y: "activity", // z: "roi", // fill: "roi", // textAnchor: "end", // text: "roi", // lineAnchor: "top", // lineWidth: 10 // // dx: 3, // // dy: -8 // }) // ) ] })34:T478,summaryTableData = { let dataToSummarise = []; let summaryTableData = []; let average, median, standardDev, maximum, minimum; let regions = [EISPolygon.properties.name,WISPolygon.properties.name,WOIPolygon.properties.name]; for (let i = 0; i < regions.length; i++) { dataToSummarise = dataInBoxes.filter((d) => d.roi == regions[i]) dataToSummarise = dataToSummarise.map((d) => d.activity) average = dataToSummarise.reduce((a, b) => a + b, 0) / dataToSummarise.length median = dataToSummarise.sort((a,b) => a.activity - b.activity)[Math.floor((dataToSummarise.length)/2)] maximum = Math.max(...dataToSummarise) minimum = Math.min(...dataToSummarise) standardDev = Math.sqrt(dataToSummarise.reduce((a, b) => a + (b - average)**2, 0) / dataToSummarise.length) ; summaryTableData[i] = {name:regions[i], N: dataToSummarise.length, average: isNaN(average) ? "-" : f(average), stddev: isNaN(standardDev) ? "-" : f(standardDev), min: dataToSummarise.length == 0 ? "" : f(minimum), max: dataToSummarise.length == 0 ? "" : f(maximum), median: isNaN(median) ? "-" : f(median) }; } return summaryTableData }35:[{"id":"35af7a55a7fbd0e15e45a26c53ab610343bab6242bed430ecffa5e8f4efcdc7a4ac835985a4984ff2b23a5fb267ec90c291dd48c49f92bcb4784ff2e479688ef","href":"https://static.observableusercontent.com/files/35af7a55a7fbd0e15e45a26c53ab610343bab6242bed430ecffa5e8f4efcdc7a4ac835985a4984ff2b23a5fb267ec90c291dd48c49f92bcb4784ff2e479688ef","name":"Nuclear Power Station Releases.xlsx","type":"application/vnd.openxmlformats-officedocument.spreadsheetml.sheet","lastModified":1723627204725.421,"size":9156},{"id":"9f08a592efbd0188c4cc2b30b0b917b4b672a3c8648ce55e280332db57576ca3115fe716dc6d9ae582f5b1fb4cc79f2abd1007a11320ec5ea4ebd067dfecfa3a","href":"https://static.observableusercontent.com/files/9f08a592efbd0188c4cc2b30b0b917b4b672a3c8648ce55e280332db57576ca3115fe716dc6d9ae582f5b1fb4cc79f2abd1007a11320ec5ea4ebd067dfecfa3a","name":"Operating Nuclear Sites@2.xlsx","type":"application/vnd.openxmlformats-officedocument.spreadsheetml.sheet","lastModified":1723627204725.421,"size":11415}] 2f:{"title":"Establishing a Baseline of Radioactivity in the Irish Marine Environment - Graphs Only","cells":[{"id":1649,"value":"# Establishing a Baseline of Radioactivity in the Irish Marine Environment - Graphs Only\n\n\n### By Aoife Roche","pinned":false,"hidden":false,"mode":"md"},{"id":1751,"value":"*Note*\n\nThis Notebook is linked to the MARIS database through the MARIS API.\n\n*What is an API?*\n\n- Application Programming Interface\n- The MARIS API was created by the IAEA to allow easier access to the database. The API acts as a link between the MARIS website and this notebook, allowing interaction between the two.","pinned":false,"hidden":false,"mode":"md"},{"id":1334,"value":"### User Interface\nClick on the dropdown arrows and radio buttons to view the data of you preference.","pinned":false,"hidden":false,"mode":"md"},{"id":136,"value":"viewof sampleType = Inputs.select(new Map([[\"Seawater\", 1], [\"Biota\", 2], [\"Sediment\", 3]]), {\n label: html`

Sample Type:

`,\n sort: true, \n unique: true\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":602,"value":"viewof biotaType = Inputs.select(sampleType != 2 ? new Map([[\"N/A\", -1]]) : new Map([[\"Fish\", 1], [\"Seaweed\", 2] , [\"Shellfish\", 3]]), {\n label: html`

Biota Type:

`,\n sort: true, \n unique: true\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":124,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":53,"value":"viewof aggType = Inputs.radio([\"max\", \"mean\", \"median\"], {\n value: \"max\",\n label: html`

Aggregation Type:

`,\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":98,"value":" // D3 and Plot expect coordinates to be specified as [longitude, latitude], so if you're setting your center point by grabbing coordinates from Google Maps, you'll need to reverse them.\nviewof kms = Inputs.range([50, 1000], {\n value: 500,\n step: 50,\n label: html`

Radius around the center (km)

`\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1340,"value":"### ${nuclidesLut[nuclide]} Activity Concentrations in ${whatSampleType(sampleType, biotaType)} (${whatUnit(sampleType)}) in the Irish Marine Environment","pinned":false,"hidden":false,"mode":"md"},{"id":57,"value":"globalMap.legend(\"color\", {\n label: `${nuclidesLut[nuclide]} Activity Concentrations in ${whatSampleType(sampleType, biotaType)} (${whatUnit(sampleType)})`, \n width: 320,\n tickFormat: \",.3~s\",\n\n \n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":38,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":1360,"value":"### Summary Table: ${nuclidesLut[nuclide]} Activity Concentrations in ${whatSampleType(sampleType, biotaType)} (${whatUnit(sampleType)})","pinned":false,"hidden":false,"mode":"md"},{"id":1663,"value":"\n \n \n \n \n \n \n \n \n \n \n \n ${summaryTableData.map(region => htl.html`\n \n \n \n \n \n \n `)}\n \n

Region of Interest	N	Arithmetic mean	Standard deviation	Median	Minimum - Maximum
${region[\"name\"]}	${region[\"N\"]}	${region[\"average\"]}	${region[\"stddev\"]}	${region[\"median\"]}	${region[\"min\"]} - ${region[\"max\"]}

","pinned":false,"hidden":false,"mode":"html"},{"id":1352,"value":"### Moving Average of ${nuclidesLut[nuclide]} Activity Concentrations in ${whatSampleType(sampleType, biotaType)} (${whatUnit(sampleType)})","pinned":false,"hidden":false,"mode":"md"},{"id":739,"value":"timePlot1.legend(\"color\", {\n label: html`

Regions of Interest

`,\n sort: true,\n unique: true\n // width: 320\n // tickFormat: \",.3~s\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":738,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":1356,"value":"### Boxplot of ${nuclidesLut[nuclide]} Activity Concentrations in ${whatSampleType(sampleType, biotaType)} (${whatUnit(sampleType)})","pinned":false,"hidden":false,"mode":"md"},{"id":1640,"value":"boxPlot.legend(\"color\", {\n label: html`

Regions of Interest

`,\n sort: true,\n unique: true})","pinned":false,"hidden":false,"mode":"ojs"},{"id":330,"value":"boxPlot = Plot.plot({\n // Colorbind-safe\n height: 300 ,\n marginLeft: 200,\n width,\n color: {\n domain: [\n \"Eastern region of the Irish Sea\",\n \"Western region of the Irish Sea\",\n \"West of Ireland\"\n ],\n range: [\"#1b9e77\", \"#7570b3\", \"#d95f02\"],\n },\n // x: { type: logscale ? \"log\": \"linear\"},\n x: {type: \"log\", grid: true, tickFormat: \".2\" } , \n y: {\n label: \"Region Of Interest\",\n domain: [\n \"Eastern region of the Irish Sea\",\n \"Western region of the Irish Sea\",\n \"West of Ireland\"\n ]\n },\n marks: [\n Plot.boxX(\n dataInBoxes,\n { x: \"activity\", y: \"roi\", tip: true, fill: \"roi\" }\n )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1530,"value":"## Radionuclide Releases from Nearby Nuclear Power Stations ","pinned":false,"hidden":false,"mode":"md"},{"id":1506,"value":"nuclearStationReleases = nuclearReleases","pinned":false,"hidden":false,"mode":"ojs"},{"id":1365,"value":"# *-Code Below Only-*","pinned":false,"hidden":false,"mode":"md"},{"id":505,"value":"### Required for Map","pinned":false,"hidden":false,"mode":"md"},{"id":44,"value":"land50m = d3.json(\"https://unpkg.com/world-atlas@2.0.2/land-50m.json \")","pinned":true,"hidden":false,"mode":"ojs"},{"id":42,"value":"land = topojson.feature(land50m, land50m.objects.land)","pinned":true,"hidden":false,"mode":"ojs"},{"id":511,"value":"### Get all Activity values for legend","pinned":false,"hidden":false,"mode":"md"},{"id":66,"value":"valuesForLegend = dataWithinRadius.map((d) => d.activity)","pinned":true,"hidden":false,"mode":"ojs"},{"id":521,"value":"### MARIS API","pinned":false,"hidden":false,"mode":"md"},{"id":523,"value":"MARIS API query parameters: \n- Select data by specific IHO (shapeTypeId == 1) area_id (hardcoded) in order to speed up the query. This includes areaId = 9999 (unattributed). Is there are better way?\n- sampleTypeId and nuclideId set by input.select control\n- Sampling date >= 1 January 2010\n- unitId Bq/m3 (seawater), Bq/kg, Bq/g C\n \nNotes:\n- Cannot select by detection, use post-processing filter to remove, e.g. \"ND\" values","pinned":false,"hidden":false,"mode":"md"},{"id":4,"value":"params = ({\n areaId: \"1912,2350,2351,2357,2389,3141,4283,9999\",\n sampleTypeId: sampleType,\n nuclideId: nuclide, // 137Cs\n samplingDateFrom: \"2010-01-01\",\n shapeTypeId: 1,\n // bioGroup: sampleType != 2 ? 0 : (biotaType == 1 ? 4 : biotaType == 2 ? 11 : \"14,12,2,13,6\") // not working - accepts integers, not strings (e.g. 4 instead of \"fish\" but query does not see to be correct. In the meantime filter in dataInterim\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":600,"value":"params","pinned":true,"hidden":false,"mode":"ojs"},{"id":6,"value":"dataRaw = getMeasurementData(params)","pinned":true,"hidden":false,"mode":"ojs"},{"id":528,"value":"Additional filtering not supported by the MARIS API: \n- Only include detected values for now (detection == \"=\")","pinned":false,"hidden":false,"mode":"md"},{"id":479,"value":"dataInterim = dataRaw.filter(\n (d) => d.detection != \"ND\" && (sampleType != 2 ? d.bioGroupId == 0 : (biotaType == 1 ? d.bioGroupId == 4 : biotaType == 2 ? d.bioGroupId == 11 : (d.bioGroupId == 14 || d.bioGroupId == 12 || d.bioGroupId == 2 || d.bioGroupId == 13 || d.bioGroupId == 6))))\n\n// ","pinned":true,"hidden":false,"mode":"ojs"},{"id":604,"value":"dataInterim","pinned":true,"hidden":false,"mode":"ojs"},{"id":481,"value":"dataInterim2 =dataInterim.map((d) => ({ ...d, activity: (d[\"detection\"] == \"<\") ? d[\"activity\"]*0.5 : d[\"activity\"]}))","pinned":true,"hidden":false,"mode":"ojs"},{"id":532,"value":"Rules for conversion of units: \n\nSeawater \n- Standardise on Bq/m3\n- \nBiota: \n- Standardise on Bq/kg f.w.\n- Assume that Bq/kg (undefined fresh or dry weight) is **fresh**\n- Convert dry weight values using LaMer average values for each biogroup\n\nSediment: \n- Standardise on Bq/kg d.w.\n- Assume that Bq/kg (undefined fresh or dry weight) is **dry**\n- Convert fresh weight values values using factor of 5\n\nTo do:\n- check if we can access percentwt\n- check if unitId = 13 Bq/kgC is used ","pinned":false,"hidden":false,"mode":"md"},{"id":420,"value":"### Data loading","pinned":false,"hidden":false,"mode":"md"},{"id":410,"value":"data = dataInterim2.map((d) => ({\n ...d,\n begPeriod: new Date(d.begPeriod),\n roi: getROI([d.longitude, d.latitude], [EISPolygon, WOIPolygon, WISPolygon])\n}))","pinned":true,"hidden":false,"mode":"ojs"},{"id":88,"value":"dataWithinRadius = data.filter(\n (d) => d3.geoDistance([d.longitude, d.latitude], center) < kms / 6371\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":603,"value":"dataWithinRadius","pinned":true,"hidden":false,"mode":"ojs"},{"id":489,"value":"dataInBoxes = dataWithinRadius.filter((d) => d.roi != \"Elsewhere\").sort((a,b) => a.begPeriod - b.begPeriod)","pinned":true,"hidden":false,"mode":"ojs"},{"id":419,"value":"### Filter data using circle","pinned":false,"hidden":false,"mode":"md"},{"id":77,"value":"center = [-5,53.5]","pinned":true,"hidden":false,"mode":"ojs"},{"id":71,"value":"circle = d3\n .geoCircle()\n .center(center)\n .radius((kms / (6371 * Math.PI * 2)) * 360)\n .precision(10)()","pinned":true,"hidden":false,"mode":"ojs"},{"id":415,"value":"### User-defined Regions of Interest (ROI)","pinned":false,"hidden":false,"mode":"md"},{"id":414,"value":"EISPolygon = turf.polygon(\n [\n [\n [-4.5, 52],\n [-2.5, 52],\n [-2.5, 55.5],\n [-4.5, 55.5],\n [-4.5, 52]\n ]\n ],\n { name: \"Eastern region of the Irish Sea\" }\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":412,"value":"WISPolygon = turf.polygon(\n [\n [\n [-6.5, 52],\n [-4.5, 52],\n [-4.5, 55.5],\n [-6.5, 55.5],\n [-6.5, 52]\n ]\n ],\n { name: \"Western region of the Irish Sea\" }\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":413,"value":"WOIPolygon = turf.polygon(\n [\n [\n [-12, 51],\n [-6.5, 51],\n [-6.5, 56],\n [-12, 56],\n [-12, 51]\n ]\n ],\n { name: \"West of Ireland\" }\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":411,"value":"// // Get Region Of Interest\nfunction getROI(point, regions) {\n const region = regions.find((region) => {\n return turf.booleanPointInPolygon(turf.point(point), region);\n });\n return region ? region.properties.name : \"Elsewhere\";\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":382,"value":"## Utils and Imports","pinned":false,"hidden":false,"mode":"md"},{"id":381,"value":"import { getMeasurementData, luts } from \"@maris-iaea/api\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":380,"value":"turf = require(\"@turf/turf@7\")","pinned":true,"hidden":false,"mode":"ojs"},{"id":435,"value":"function whatUnit(samptypeId) {\n switch(samptypeId) {\n case 1:\n return \"Bq/m3\"\n case 2:\n return \"Bq/kg f.w.\"\n case 3:\n return \"Bq/kg d.w.\"\n default:\n return \"Unknown\"\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1481,"value":"function whatSampleType(samptypeId, biotaTypeID){\n switch(samptypeId) {\n case 1:\n return \"Seawater\"\n case 3: \n return \"Sediment\"\n case 2:\n switch(biotaTypeID) {\n case 1:\n return \"Fish\"\n case 2:\n return \"Seaweed\"\n case 3:\n return \"Shellfish\"\n default:\n return \"Unknown\"\n }\n default:\n return \"Unknown\" \n } \n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":437,"value":"function getDryWetFactor(biogroupId) {\n switch(biogroupId) {\n case 11:\n return 0.167278970397589\n case 4:\n return 0.232254394570392\n case 2:\n return 0.226485185185185\n case 14:\n return 0.169374196277088\n case 13:\n return 0.232484162739303\n case 12:\n return 0.181976062091503 \n case 6:\n return 0.194611473702632\n case 5:\n return 0.460238713453779\n default:\n return 1\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":734,"value":"nuclidesLut = _.fromPairs(dataWithinRadius.map((d) => [d.nuclideId, d.nuSymbol]))","pinned":true,"hidden":false,"mode":"ojs"},{"id":733,"value":"## Nuclear sites","pinned":false,"hidden":false,"mode":"md"},{"id":732,"value":"nuclearSites = FileAttachment(\"Operating Nuclear Sites@2.xlsx\").xlsx()","pinned":true,"hidden":false,"mode":"ojs"},{"id":731,"value":"nuclearSitesdata = nuclearSites.sheet(0, { headers: true })","pinned":true,"hidden":false,"mode":"ojs"},{"id":1503,"value":"workbook = FileAttachment(\"Nuclear Power Station Releases.xlsx\").xlsx()","pinned":true,"hidden":false,"mode":"ojs"},{"id":1504,"value":"workbook.sheetNames","pinned":true,"hidden":false,"mode":"ojs"},{"id":1505,"value":"nuclearReleases = workbook.sheet(0, {\n headers: true,\n // range: \"A1:J10\"\n })","pinned":true,"hidden":false,"mode":"ojs"},{"id":1662,"value":"$34","pinned":true,"hidden":false,"mode":"ojs"},{"id":1661,"value":"regions = [EISPolygon.properties.name,WISPolygon.properties.name,WOIPolygon.properties.name];","pinned":true,"hidden":false,"mode":"ojs"},{"id":1660,"value":"f = d3.format(\".1f\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W35","stdlib":"1","titling":"auto"}