Model:Chi analysis tools: Difference between revisions
(Created page with "{{Model identity |Model also known as=Channel segment finder |Model type=Single }} {{Start models incorporated}} {{End a table}} {{Model identity2 |Categories=Terrestrial |Spa...") |
m (Text replacement - "{{End headertab}}" to "{{End headertab}} {{{{PAGENAME}}_autokeywords}}") |
||
(34 intermediate revisions by 2 users not shown) | |||
Line 6: | Line 6: | ||
{{End a table}} | {{End a table}} | ||
{{Model identity2 | {{Model identity2 | ||
| | |ModelDomain=Terrestrial | ||
|Spatial dimensions=1.5D | |Spatial dimensions=1.5D | ||
|Spatialscale=Regional-Scale, Landscape-Scale, Watershed-Scale | |Spatialscale=Regional-Scale, Landscape-Scale, Watershed-Scale | ||
|One-line model description=Tool for examining channel profiles in chi-elevation space using the integral method of channel analysis | |One-line model description=Tool for examining channel profiles in chi-elevation space using the integral method of channel analysis | ||
|Extended model description=This tool is used for examining bedrock channels. The tool is based on the assumption that the stream power incision model (SPIM) adequately describes channel incision. Channels profiles are converted to chi-elevation space, where chi is a transformed longitudinal coordinate that takes drainage area into account. The tool uses a variety of statistical tests to extract the most likely series of segments with distinct steepness in chi-elevation space. It also performs statistical tests to determine the best fit m/n ratio, where m is an area (A) exponent and n is a slope (S) exponent in the SPIM with E = K A^m S^n, where E is an erosion rate and K is an 'erodibility'. | |Extended model description=This tool is used for examining bedrock channels. The tool is based on the assumption that the stream power incision model (SPIM) adequately describes channel incision. Channels profiles are converted to chi-elevation space, where chi is a transformed longitudinal coordinate that takes drainage area into account. The tool uses a variety of statistical tests to extract the most likely series of segments with distinct steepness in chi-elevation space. It also performs statistical tests to determine the best fit m/n ratio, where m is an area (A) exponent and n is a slope (S) exponent in the SPIM with E = K A^m S^n, where E is an erosion rate and K is an 'erodibility'. | ||
}} | }} | ||
{{Start model keyword table}} | {{Start model keyword table}} | ||
{{Model keywords | {{Model keywords | ||
|Model keywords= | |Model keywords=geomorphology | ||
}} | }} | ||
{{Model keywords | {{Model keywords | ||
Line 36: | Line 36: | ||
|Email address=simon.m.mudd@ed.ac.uk | |Email address=simon.m.mudd@ed.ac.uk | ||
|Phone=+44 (0) 131 650 2535 | |Phone=+44 (0) 131 650 2535 | ||
}} | |||
{{Additional modeler information | |||
|Additional first name=David | |||
|Additional last name=Milodowski | |||
|Additional type of contact=Model developer | |||
|Additional institute / Organization=University of Edinburgh | |||
|Additional postal address 1=School of GeoSciences | |||
|Additional postal address 2=Drummond Street | |||
|Additional town / City=Edinburgh | |||
|Additional postal code=EH8 9XP | |||
|Additional state=NO STATE | |||
|Additional country=United Kingdom | |||
}} | |||
{{Additional modeler information | |||
|Additional first name=Stuart | |||
|Additional last name=Grieve | |||
|Additional type of contact=Model developer | |||
|Additional institute / Organization=University of Edinburgh | |||
|Additional postal address 1=School of GeoSciences | |||
|Additional postal address 2=Drummond Street | |||
|Additional town / City=Edinburgh | |||
|Additional postal code=EH8 9XP | |||
|Additional state=NO STATE | |||
|Additional country=United Kingdom | |||
}} | |||
{{Additional modeler information | |||
|Additional first name=Declan | |||
|Additional last name=Valters | |||
|Additional type of contact=Model developer | |||
|Additional institute / Organization=University of Edinburgh | |||
|Additional postal address 1=School of GeoSciences | |||
|Additional postal address 2=Drummond Street | |||
|Additional town / City=Edinburgh | |||
|Additional postal code=EH8 9XP | |||
|Additional state=NO STATE | |||
|Additional country=United Kingdom | |||
}} | }} | ||
{{Model technical information | {{Model technical information | ||
Line 41: | Line 77: | ||
|Programming language=C++, Python | |Programming language=C++, Python | ||
|Code optimized=Single Processor | |Code optimized=Single Processor | ||
|Start year development=2012 | |Start year development=2012 | ||
|Does model development still take place?=Yes | |Does model development still take place?=Yes | ||
|Model availability=As code | |Model availability=As code | ||
|Source code availability=Through CSDMS repository | |Source code availability=Through CSDMS repository | ||
|Source web address= | |Source csdms web address=https://github.com/csdms-contrib/chi_analysis_tools | ||
|Program license type=GPL v2 | |Program license type=GPL v2 | ||
|Memory requirements=Depends on the size of your DEM | |Memory requirements=Depends on the size of your DEM | ||
|Typical run time=A few minutes to several days depending on parameters and DEM size | |Typical run time=A few minutes to several days depending on parameters and DEM size | ||
}} | }} | ||
{{Input - Output description | {{Input - Output description | ||
|Describe input parameters=The input is a 'channel file' and a 'driver file'. The channel file contains data on channel profiles within a channel network composed of a main stem and tributaries flowing into that main stem (that is, there are no tributaries of tributaries). The driver file contains parameters for the model run. The format of these files is described in the documentation that accompanies the model source code. | |Describe input parameters=The input is a 'channel file' and a 'driver file'. The channel file contains data on channel profiles within a channel network composed of a main stem and tributaries flowing into that main stem (that is, there are no tributaries of tributaries). The driver file contains parameters for the model run. The format of these files is described in the documentation that accompanies the model source code. | ||
|Input format=ASCII | |Input format=ASCII | ||
|Describe output parameters=The outputs are | |Describe output parameters=The outputs are | ||
i) For chi_m_over_n_analysis.exe, a *.movern file that contains information about the goodness of fit of channel profiles to a series of linear segments as a function of the m/n ratio: this file is used to determine the best fit m/n ratio of a channel network. | i) For chi_m_over_n_analysis.exe, a *.movern file that contains information about the goodness of fit of channel profiles to a series of linear segments as a function of the m/n ratio: this file is used to determine the best fit m/n ratio of a channel network. | ||
ii) For chi_get_profiles.exe, a series of *.tree files which contain information about the best fit channel segments in chi-elevation space. This data can be used to infer erosion rates, tectonics, or variations in erodibility. | ii) For chi_get_profiles.exe, a series of *.tree files which contain information about the best fit channel segments in chi-elevation space. This data can be used to infer erosion rates, tectonics, or variations in erodibility. | ||
|Output format=ASCII | |Output format=ASCII | ||
|Pre-processing software needed?=Yes | |Pre-processing software needed?=Yes | ||
|Describe pre-processing software=Channel profile data must be prepared to meed the specifications of the *.chan file. The format of this file is described in the documentation provided with the source code. | |Describe pre-processing software=Channel profile data must be prepared to meed the specifications of the *.chan file. The format of this file is described in the documentation provided with the source code. | ||
|Post-processing software needed?=No | |Post-processing software needed?=No | ||
|Visualization software needed?=No | |Visualization software needed?=No | ||
}} | }} | ||
{{Process description model | {{Process description model | ||
|Describe processes represented by the model=The module performs topographic analysis but the analysis is based on the assumption that the stream power incision model is a good approximation for channel incision. | |Describe processes represented by the model=The module performs topographic analysis but the analysis is based on the assumption that the stream power incision model is a good approximation for channel incision. | ||
|Describe key physical parameters and equations=Topographic analysis so no physical parameters as such. The main parameters fro the module are: | |Describe key physical parameters and equations=Topographic analysis so no physical parameters as such. The main parameters fro the module are: | ||
1) sigma: the uncertainty (in metres, or meters for you yanks), of the topographic data. | 1) sigma: the uncertainty (in metres, or meters for you yanks), of the topographic data. | ||
2) The target nodes: These are the number of nodes you wish to use in each subset of the channel. For details see associated documentation. This should vary between 60-140. The module partitions data and the number of partitions is a highly nonlinear function of the number of nodes so target node values of >150 will lead to compute times of many months - forever. | 2) The target nodes: These are the number of nodes you wish to use in each subset of the channel. For details see associated documentation. This should vary between 60-140. The module partitions data and the number of partitions is a highly nonlinear function of the number of nodes so target node values of >150 will lead to compute times of many months - forever. | ||
3) The minimum segment length: The shortest contiguous number of nodes the user is willing to consider for statistical analysis (in testing, 8-20 performed reasonably well). If you chose 2 you will be performing linear analysis on segments with 2 data points which is clearly nonsense. | 3) The minimum segment length: The shortest contiguous number of nodes the user is willing to consider for statistical analysis (in testing, 8-20 performed reasonably well). If you chose 2 you will be performing linear analysis on segments with 2 data points which is clearly nonsense. | ||
4) Mean skip: See associated documentation, but the module uses a Monte Carlo sampling regime which skips nodes, analyses a subset of data, and then performs this skipping and analysis routine over a number of iterations. For SRTM and ASTER data this should be 1-2. For 10m data it can be 1-10, and for LiDAR data you could skip up to 100 nodes. Note that because of the iterative processes you will need to increase the number of iterations as you increase the skip value if you are to sample all of the data. | 4) Mean skip: See associated documentation, but the module uses a Monte Carlo sampling regime which skips nodes, analyses a subset of data, and then performs this skipping and analysis routine over a number of iterations. For SRTM and ASTER data this should be 1-2. For 10m data it can be 1-10, and for LiDAR data you could skip up to 100 nodes. Note that because of the iterative processes you will need to increase the number of iterations as you increase the skip value if you are to sample all of the data. | ||
|Describe length scale and resolution constraints=With reasonable parameter values channel networks with individual channel of ~1000 nodes will take a few to 10s of minutes to analyse. The longer the channel network, the longer the analysis. | |Describe length scale and resolution constraints=With reasonable parameter values channel networks with individual channel of ~1000 nodes will take a few to 10s of minutes to analyse. The longer the channel network, the longer the analysis. | ||
|Describe time scale and resolution constraints=Topographic analysis so no time scale constraints. | |Describe time scale and resolution constraints=Topographic analysis so no time scale constraints. | ||
|Describe any numerical limitations and issues=See documentation. The major limitation is computational time. This can be alleviated with sensible selection of module parameters. See documentation for guidance. | |Describe any numerical limitations and issues=See documentation. The major limitation is computational time. This can be alleviated with sensible selection of module parameters. See documentation for guidance. | ||
}} | }} | ||
{{Model testing | {{Model testing | ||
|Describe available calibration data sets=Topographic analysis: no calibration required. | |Describe available calibration data sets=Topographic analysis: no calibration required. | ||
|Describe available test data sets=There are 3 datasets within this repository. | |Describe available test data sets=There are 3 datasets within this repository. One from Southwest Pennsylvania, one from the Apennines in Italy, and one from a CHILD model run. Datasets are the same as those run with the associated manuscript (http://www.geos.ed.ac.uk/homes/smudd/Muddetal_JGRsubmit.pdf) so users can compare results with figures from the manuscript. | ||
|Describe ideal data for testing=Any bedrock channel profile. Module is intended for use on topographic data. | |Describe ideal data for testing=Any bedrock channel profile. Module is intended for use on topographic data. | ||
}} | }} | ||
{{Users groups model | {{Users groups model | ||
Line 91: | Line 119: | ||
{{Documentation model | {{Documentation model | ||
|Manual model available=Yes | |Manual model available=Yes | ||
|Model manual=Chi analysis movern and profiles.pdf,Chi analysis getting channel.pdf, | |||
}} | }} | ||
{{Additional comments model | {{Additional comments model}} | ||
}} | |||
{{CSDMS staff part | {{CSDMS staff part | ||
|OpenMI compliant=No | |OpenMI compliant=No but possible | ||
| | |IRF interface=No but possible | ||
| | |CMT component=No but possible | ||
| | |CCA component=No but possible | ||
}} | }} | ||
{{Start coupled table}} | {{Start coupled table}} | ||
{{End a table}} | {{End a table}} | ||
{{End headertab}} | {{End headertab}} | ||
{{{{PAGENAME}}_autokeywords}} | |||
<!-- PLEASE USE THE "EDIT WITH FORM" BUTTON TO EDIT ABOVE CONTENTS; CONTINUE TO EDIT BELOW THIS LINE --> | <!-- PLEASE USE THE "EDIT WITH FORM" BUTTON TO EDIT ABOVE CONTENTS; CONTINUE TO EDIT BELOW THIS LINE --> | ||
==Introduction== | ==Introduction== | ||
This code has been developed by the University of Edinburgh Land Surface Dynamics group to analyse channel longitudinal profiles using the integral method of channel analysis. The integral method is a means of comparing the relative steepness of channel profiles, normalized for drainage area. It is similar to, but has some advantages over, traditional slope area analysis. See Perron and Royden (2013) for more details on the integral method. | |||
This particular code not only transforms channels using the integral method, but it also uses statistical techniques to find the most likely series of channel segments with distinct steepness and the most likely m/n ratio. The details of these algorithms can be found in Mudd et al. (submitted manuscript) which is available here: http://www.geos.ed.ac.uk/homes/smudd/Muddetal_JGRsubmit.pdf. | |||
The code consists of objects that contain topographic data and their member functions, and driver programs perform the analyses. Full documentation is available under the ‘other’ tab in the wiki. | |||
Reference: | |||
Perron, J.T. and L. Royden (2013). An integral approach to bedrock river profile analysis. Earth Surface Processes and Landforms, 38, 570-576, doi:10.1002/esp.3302. | |||
== History == | == History == | ||
This code has been under development since 2012; Simon Mudd began the project and David Milodowski, Stuart Grieve and Declan Valters have all contributed to its development. Fiona Clubb, Declan Valters and Jen Merritt at the University of Edinburgh began beta testing in January 2013. Rahul Devrani and Vimal Singh of the University of Dehli began beta testing in April of 2013, and TC Hales of Cardiff University began beta testing in May 2013. The working version of the code used in the manuscript was finished in June 2013. | |||
== References == | == References == | ||
{{ | <br>{{AddReferenceUploadButtons}}<br><br> | ||
{{ | {{#ifexist:Template:{{PAGENAME}}-citation-indices|{{{{PAGENAME}}-citation-indices}}|}}<br> | ||
{{Include_featured_references_models_cargo}}<br> | |||
== Issues == | == Issues == | ||
If you find bugs, email simon.m.mudd _at_ ed.ac.uk | |||
== Help == | == Help == | ||
Go to the 'other' tab in this wiki: documentation and instructions are found in two .pdf documents available there. | |||
{{#ifexist:Model_help:{{PAGENAME}}|[[Model_help:{{PAGENAME}}]]|}} | {{#ifexist:Model_help:{{PAGENAME}}|[[Model_help:{{PAGENAME}}]]|}} | ||
== Input Files == | == Input Files == | ||
See documentation: go to 'other' tab in this wiki. | |||
== Output Files == | == Output Files == | ||
See documentation: go to 'other' tab in this wiki. |
Latest revision as of 20:15, 16 September 2020
Chi analysis tools
Metadata
|
|
Introduction
This code has been developed by the University of Edinburgh Land Surface Dynamics group to analyse channel longitudinal profiles using the integral method of channel analysis. The integral method is a means of comparing the relative steepness of channel profiles, normalized for drainage area. It is similar to, but has some advantages over, traditional slope area analysis. See Perron and Royden (2013) for more details on the integral method.
This particular code not only transforms channels using the integral method, but it also uses statistical techniques to find the most likely series of channel segments with distinct steepness and the most likely m/n ratio. The details of these algorithms can be found in Mudd et al. (submitted manuscript) which is available here: http://www.geos.ed.ac.uk/homes/smudd/Muddetal_JGRsubmit.pdf.
The code consists of objects that contain topographic data and their member functions, and driver programs perform the analyses. Full documentation is available under the ‘other’ tab in the wiki.
Reference: Perron, J.T. and L. Royden (2013). An integral approach to bedrock river profile analysis. Earth Surface Processes and Landforms, 38, 570-576, doi:10.1002/esp.3302.
History
This code has been under development since 2012; Simon Mudd began the project and David Milodowski, Stuart Grieve and Declan Valters have all contributed to its development. Fiona Clubb, Declan Valters and Jen Merritt at the University of Edinburgh began beta testing in January 2013. Rahul Devrani and Vimal Singh of the University of Dehli began beta testing in April of 2013, and TC Hales of Cardiff University began beta testing in May 2013. The working version of the code used in the manuscript was finished in June 2013.
References
Nr. of publications: | 1 |
Total citations: | 144 |
h-index: | 1 |
m-quotient: | 0.1 |
Featured publication(s) | Year | Model described | Type of Reference | Citations |
---|---|---|---|---|
Mudd, Simon M.; Attal, Mikaël; Milodowski, David T.; Grieve, Stuart W. D.; Valters, Declan A.; 2014. A statistical framework to quantify spatial variation in channel gradients using the integral method of channel profile analysis: CHANNEL SEGMENT FITTING. Journal of Geophysical Research: Earth Surface, 119, 138–152. 10.1002/2013JF002981 (View/edit entry) | 2014 | Chi analysis tools |
Model overview | 144 |
See more publications of Chi analysis tools |
Issues
If you find bugs, email simon.m.mudd _at_ ed.ac.uk
Help
Go to the 'other' tab in this wiki: documentation and instructions are found in two .pdf documents available there.
Input Files
See documentation: go to 'other' tab in this wiki.
Output Files
See documentation: go to 'other' tab in this wiki.