Confluence Retirement

In an effort to consolidate USGS hosted Wikis, myUSGS’ Confluence service is scheduled for retirement on January 27th, 2023. The official USGS Wiki and collaboration space is now SharePoint. Please migrate existing spaces and content to the SharePoint platform and remove it from Confluence at your earliest convenience. If you need any additional information or have any concerns about this change, please contact Thank you for your prompt attention to this matter.
Skip to end of metadata
Go to start of metadata


There are probably lots of ways to this. Having had some recent (10/2015) conversations w/the fine folks at Natural Resources Canada, the idea of using the Horton stream order looks pretty interesting. If you have other ideas for ways to resolve this from fairly limited attribution (like from/to connectivity, or other geometric info), please make a section!
So, consider an example like that shown in the following image as what we're trying to identify:
Here are some examples of different kinds of ordering, pointed out by Denis Boutin (NRCAN) in the GRASS documentation:

Note that most of these schemes are built a non-divergent model of flow direction downslope. While this is fine for DEM-derived networks, it might require that the network be pre-processed (pruned, generally cleaned up) to enforce this assumption prior to applying the selected ordering scheme (from which we extract the "main stem").

Hack Ordering

The algorithm walks upstream from outlet segment, selecting to travel up the tributary with the longest accumulated flowlength. This is repeated until a headwater segment is reached. 

(from the GRASS documentation on stream ordering) 

"the main stream of every subcatchment is the longest stream or stream with highest accumulation rate."

In most cases, the main stream is the longest watercourse of the catchment, but in some cases, when a catchment consists of both rounded and elongated subcatchments these rules may not be maintained. The algorithm assigns 1 to every outlet stream and goes upstream according to the maximum flow accumulation of every branch. When it reaches an initial stream, it steps back to the first unassigned confluence. It assigns order 2 to unordered tributaries and again goes upstream to the next initial stream. The process repeats until all branches of all outlets are ordered.

Horton Ordering

This requires first calculating Strahler ordering and then making a second pass, wherein the algorithm walks upstream from outlet segment, selecting to travel up to the tributary with the highest stream order. In the case of a tie in the Strahler order of the tributaries, the tributary with the longest accumulated flowlength is selected (as in Hack ordering). The process repeats until all branches of all outlets are ordered.


Pfaffstetter Ordering

Kristine Verdin (USGS) is longtime champion of this kind of scheme. Probably would be a smart plan to talk w/her about applying this method.


Software Resources


  • No labels