45:["$","title",null,{"children":"Rio Grande near Del Norte | Peter Donovan | Observable"}] 4d:T5e2,brushFilterX = function ( extent, { defaultExtent, marginTop = 10, marginRight = 20, marginBottom = 20, marginLeft = 65, width = document.body.clientWidth, height = 80 } = {} ) { const x = d3 .scaleTime() .domain(extent) .range([marginLeft, width - marginRight]); const xAxis = (g) => g .attr("transform", `translate(0,${height - marginBottom + 1})`) .call(d3.axisBottom(x).ticks(width < 710 ? 6 : 10)); const svg = d3 .create("svg") .attr("width", width) .attr("height", height) .attr("viewBox", [0, 0, width, height]); svg .selectAll("text") .data(["adjust grey rectangle to select date range for ALL CHARTS"]) .join("text") .text((d) => d) .style("font-size", 8) .style("font-family", "sans-serif") .attr("dx", 10) .attr("dy", (d, i) => marginTop + 10 + i * 10); const brush = d3 .brushX() .extent([ [marginLeft, marginTop], [width - marginRight, height - marginBottom] ]) .on("brush end", brushed); svg.append("g").call(xAxis); let b = svg.append("g").call(brush); if (defaultExtent) b.call(brush.move, defaultExtent.map(x)); function brushed(event) { svg.node().value = event.selection === null ? null : event.selection.map(x.invert); svg.node().dispatchEvent(new Event("input", { bubbles: true })); const selection = event.selection; } svg.node().value = defaultExtent ? defaultExtent : null; return svg.node(); }4e:[{"id":"f381e09dc69d3e6faef65aef295c42db92c91b9ae38540ba9f60973198d757922ee28d36beb3272b5e8268ac862754e9e0e08f381e7ca99b595313422d850692","href":"https://static.observableusercontent.com/files/f381e09dc69d3e6faef65aef295c42db92c91b9ae38540ba9f60973198d757922ee28d36beb3272b5e8268ac862754e9e0e08f381e7ca99b595313422d850692","name":"delnorteDWR.csv","type":"text/csv","lastModified":1739471305079.58,"size":22659},{"id":"528183bb7ba50721795b164fb9c5e3df2cc7d28a55db1ad75d938eec2e4451f8b45a0e4eb1682fd9d7b9578dfe3c0320045989f826101fddab29cfde93429705","href":"https://static.observableusercontent.com/files/528183bb7ba50721795b164fb9c5e3df2cc7d28a55db1ad75d938eec2e4451f8b45a0e4eb1682fd9d7b9578dfe3c0320045989f826101fddab29cfde93429705","name":"delnorte.csv","type":"text/csv","lastModified":1739471313843.864,"size":1330267}] 49:{"id":"1b1fa9805c1cdd22","slug":"rio-grande-near-del-norte","title":"Rio Grande near Del Norte","thumbnail":"cfe2f2d54b7234fea58930a7cb8c962876ce563bd7cd10355bd3e2eb2eaa2f49","access_level":"public","create_time":"$D2025-02-13T13:48:43.347Z","update_time":"$D2025-04-04T21:03:23.333Z","view_time":null,"first_public_time":"$D2025-02-14T16:07:17.658Z","has_unpublished_edits":null,"fork":false,"stars":0,"starred":false,"owner":{"type":"team","id":"e2b593d20ef3388a","login":"managingwholes","name":"Peter Donovan","avatar_url":"https://avatars-next.observableusercontent.com/6ff849db058c39c81342ebf90aa84a9f6a536ac9ded98087e309e6a9ce577bf5","tier":"starter_2024","bio":"","followed_by_you":false,"follows_you":false},"version":354,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"0550a7b94f6b94ca","login":"managingwholes","name":"Peter Donovan","avatar_url":"https://avatars-next.observableusercontent.com/6ff849db058c39c81342ebf90aa84a9f6a536ac9ded98087e309e6a9ce577bf5","tier":null,"bio":"","followed_by_you":false,"follows_you":false}],"members":[{"type":"individual","id":"0550a7b94f6b94ca","login":"managingwholes","name":"Peter Donovan","avatar_url":"https://avatars-next.observableusercontent.com/6ff849db058c39c81342ebf90aa84a9f6a536ac9ded98087e309e6a9ce577bf5","tier":null,"role":"editor"},{"type":"team","id":"e2b593d20ef3388a","login":"managingwholes","name":"Peter Donovan","avatar_url":"https://avatars-next.observableusercontent.com/6ff849db058c39c81342ebf90aa84a9f6a536ac9ded98087e309e6a9ce577bf5","tier":"starter_2024","role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"titling":"auto","title":"Rio Grande near Del Norte","cells":[{"id":0,"value":"# Rio Grande near Del Norte\n\n

Data from USGS, DWR. Del Norte gauge (previously USGS 08220000, data from 1890-2023 with some gaps) now operated by State Dept. of Water Resources: https://dwr.state.co.us/Tools/Stations/RIODELCO

\n\n

To SHIFT DATE RANGE: Drag the gray rectangle atop the timeline, or its left or right edges, to shift date ranges on charts below.

","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":57,"value":"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":23,"value":"viewof limitsDelNorte = brushFilterX(\n d3.extent(delnorteUSGScleaned.concat(delnorteDWRCleaned), (d) => d.date),\n { defaultExtent: [d3.isoParse(\"1996-01-01\"), d3.isoParse(new Date())] }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":9,"value":"delNorteGauge = Plot.plot({\n x: {\n label: \"date\",\n domain: limitsDelNorte || undefined\n },\n y: {\n type: \"sqrt\",\n label: \"cfs\"\n },\n grid: true,\n title: \"Rio Grande near Del Norte: daily cfs\",\n caption: md`Daily mean flows of Rio Grande near Del Norte in cubic feet per second from from USGS and DWR. Mouse over line to read values.`,\n height: 250,\n width,\n marginLeft: 50,\n marks: [\n Plot.ruleY([0]),\n Plot.lineY(delnorteCombined, {\n x: \"date\",\n y: \"cfs\",\n stroke: \"black\",\n tip: true\n }),\n Plot.areaY(delnorteCombined, {\n x: \"date\",\n y: \"cfs\",\n fill: \"steelblue\",\n fillOpacity: 0.15\n })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":131,"value":"delNorteAF = Plot.plot({\n x: {\n label: \"year\",\n domain: limitsDelNorte || undefined\n },\n y: {\n // type: \"sqrt\",\n label: \"April-July discharge AF\"\n },\n grid: true,\n title: \"Volume of Apr-Jul discharge by year\",\n caption: md`Acre-feet of April-July discharge by year, charted by midpoints. 5-year rolling mean in red.`,\n height: 200,\n width,\n marginLeft: 50,\n marks: [\n Plot.ruleY([0]),\n Plot.lineY(midpoints, {\n x: \"date\",\n y: \"AprJulAF\",\n stroke: \"blue\",\n tip: true\n }),\n Plot.lineY(\n midpoints,\n Plot.windowY({\n k: 5,\n curve: \"basis\",\n x: \"date\",\n y: \"AprJulAF\",\n stroke: \"red\",\n clip: true\n })\n )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":297,"value":"midpointChart = Plot.plot({\n x: {\n label: \"midpoint date\",\n domain: limitsDelNorte || undefined\n },\n y: {\n // type: \"sqrt\",\n label: \"midpoint day of year\"\n },\n grid: true,\n title: \"Timing of midpoint of Apr-Jul discharge by year\",\n caption: md`Showing the date and day of year when half the volume of Apr-Jul discharge had occurred. 5-year rolling mean in red. One hypothesis is that dust on snow is contributing to earlier runoff.`,\n height: 200,\n width,\n marginLeft: 50,\n marks: [\n // Plot.ruleY([0]),\n Plot.lineY(midpoints, {\n x: \"date\",\n y: \"doy\",\n stroke: \"blue\",\n tip: true\n }),\n Plot.lineY(\n midpoints,\n Plot.windowY({\n k: 5,\n curve: \"basis\",\n x: \"date\",\n y: \"doy\",\n stroke: \"red\",\n clip: true\n })\n )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":324,"value":"## Functions and imports","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":128,"value":"midpoints = aq // using the arquero library, which has many data manipulation functions\n .from(delnorteCombined)\n .filter((d) => d.doy > 90 && d.doy < 213) // filters data from April 1 - July 31\n .groupby(\"year\")\n .derive({ cumsum: aq.rolling((d) => Math.floor(op.sum(d[\"af\"]))) }) // cumulative sum of acre feet for period\n .derive({ AprJulAF: (d) => Math.floor(op.sum(d[\"af\"])) }) // sum of acre feet for period\n .derive({ diff: (d) => op.abs(d[\"cumsum\"] - d[\"AprJulAF\"] / 2) }) // absolute value of difference from half of Apr-Jul discharge\n .orderby(\"diff\") // sorts by diff, ascending, within each year\n .filter((d) => op.rank() === 1) // selects the lowest diff for each year, closest to midpoint\n .orderby(\"year\")\n .objects() // yields array of objects, 1 for each year, including date and doy of midpoint","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":50,"value":"// combines cleaned versions of downloaded USGS and DWR data, adds acre feet, day of year, and full year\ndelnorteCombined = delnorteUSGScleaned.concat(delnorteDWRCleaned).map((d) => ({\n date: d[\"date\"],\n cfs: d[\"cfs\"],\n af: Math.floor(d[\"cfs\"] * 1.983),\n doy: moment(d.date).dayOfYear() + 1, // the plus 1 allows for time zone shift?\n year: d[\"date\"].getFullYear()\n}))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":41,"value":"delnorteUSGScleaned = delnorteUSGS.map((d) => ({\n date: d[\"date\"],\n cfs: d[\"cfs\"]\n}))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":7,"value":"delnorteDWRCleaned = delnorteDWR.map((d) => ({\n date: new Date(\n d[\"date\"].split(\"/\")[2] + // re-orders the date string\n \"-\" +\n d[\"date\"].split(\"/\")[0] +\n \"-\" +\n d[\"date\"].split(\"/\")[1]\n ),\n cfs: d[\"cfs\"]\n}))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":34,"value":"delnorteDWR = FileAttachment(\"delnorteDWR.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":32,"value":"delnorteUSGS = FileAttachment(\"delnorte.csv\").csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":61,"value":"import { aq, op } from \"@uwdata/arquero\" // good library for grouping, reducing, etc.\n// docs: https://observablehq.com/@uwdata/arquero-cookbook","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":67,"value":"moment = require(\"moment\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":12,"value":"$4d","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W4e","stdlib":"1"}} 44:null