Model:Plume: Difference between revisions

From CSDMS
No edit summary
m (removed old questionnaire information. Replaced with form)
Line 75: Line 75:
<!-- Edit the part above to update info on other papers -->
<!-- Edit the part above to update info on other papers -->


== Plume ==
==Introduction==
__TOC__
== PLUME <ref>Hutton and Syvitski, 2008. Sedflux-2.0: An advanced process-response model that generates three-dimensional stratigraphy.  Computers and Geosciences, v. 34. [http://dx.doi.org/10.1016/j.cageo.2008.02.013 doi:10.1016/j.cageo.2008.02.013]</ref> <ref>Syvitski et al., 1998. PLUME1.1: Deposition of sediment from a fluvial plume  ([http://dx.doi.org/10.1016/S0098-3004(97)00084-8 doi:10.1016/S0098-3004(97)00084-8] </ref><ref> Peckham, S.D., 2008. A new method for estimating suspended sediment concentrations and deposition rates from satellite imagery based on the physics of plumes. Computer & Geosciences, 34, 1198-1222. [http://dx.doi.org/10.1016/doi:10.1016/j.cageo.2008.02.009 doi:10.1016/j.cageo.2008.02.009]</ref> ==
 
===Introduction===
=== PLUME <ref>Hutton and Syvitski, 2008. Sedflux-2.0: An advanced process-response model that generates three-dimensional stratigraphy.  Computers and Geosciences, v. 34. [http://dx.doi.org/10.1016/j.cageo.2008.02.013 doi:10.1016/j.cageo.2008.02.013]</ref> <ref>Syvitski et al., 1998. PLUME1.1: Deposition of sediment from a fluvial plume  ([http://dx.doi.org/10.1016/S0098-3004(97)00084-8 doi:10.1016/S0098-3004(97)00084-8] </ref><ref> Peckham, S.D., 2008. A new method for estimating suspended sediment concentrations and deposition rates from satellite imagery based on the physics of plumes. Computer & Geosciences, 34, 1198-1222. [http://dx.doi.org/10.1016/doi:10.1016/j.cageo.2008.02.009 doi:10.1016/j.cageo.2008.02.009]</ref> ===


Plume simulates hypopycnal plumes generated by a river draining its suspended sediment load into a receiving basin.
Plume simulates hypopycnal plumes generated by a river draining its suspended sediment load into a receiving basin.
Line 102: Line 99:




===References to plume theory===
==References to plume theory==


Albertson, M.L., Dai, Y.B., Jensen, R.A., Hunter, R., 1950, Diffusion of submerged jets. American Society Civil Engineers Trans, v. 115, p. 639-697.
Albertson, M.L., Dai, Y.B., Jensen, R.A., Hunter, R., 1950, Diffusion of submerged jets. American Society Civil Engineers Trans, v. 115, p. 639-697.
Line 109: Line 106:
<br>
<br>


===Plume model papers===
==Plume model papers==
  <references/>
  <references/>


=== Plume Questionnaire ===
== Issues ==


==== Contact Information ====
== Help ==
 
=== Help on Model Output ===
{| class="wikitable"
| class="model_col1"| Model:
| class="model_col2"| Plume
|-
| class="model_col1"| Contact person:
| class="model_col2"| Eric Hutton
|-
| class="model_col1"| Institute:
| class="model_col2"| University of Colorado
|-
| class="model_col1"| City:
| class="model_col2"| Boulder, CO
|-
| class="model_col1"| Country:
| class="model_col2"| USA
|-
| class="model_col1"| Email:
| class="model_col2"| --
|}
 
==== Model description ====
 
{| class="wikitable"
| class="model_col1"| Model type:
| class="model_col2"| Single model
|-
| class="model_col1"| Description:
| class="model_col2"| Run a hypopycnal sediment plume
|}
 
==== Technical information ====
 
{| class="wikitable"
| class="model_col1"| Supported platforms:
| class="model_col2"| UNIX, Linux, Mac OSX, Windows
|-
| class="model_col1"| Programming language:
| class="model_col2"|  C
|-
| class="model_col1"| Model was developed started from:
| class="model_col2"| 1997 and development still takes place
|-
| class="model_col1"| To what degree will the model become available:
| class="model_col2"| As code
|-
| class="model_col1"| Current license type:
| class="model_col2"| Apache public license
|-
| class="model_col1"| Memory requirements:
| class="model_col2"| Minimal
|-
| class="model_col1"| Typical run time:
| class="model_col2"| seconds
|}
 
==== Input / Output description ====
 
{| class="wikitable"
| class="model_col1"| Input parameters:
| class="model_col2"| River velocity, width, depth; Sediment concentrations
|-
| class="model_col1"| Input format:
| class="model_col2"| ASCII
|-
| class="model_col1"| Output parameters:
| class="model_col2"| Grid of Sediment rate in m/day for specified grain size classes
|-
| class="model_col1"| Output format:
| class="model_col2"| ASCII
|-
| class="model_col1"| Post-processing software (if needed):
| class="model_col2"| no
|-
| class="model_col1"| Visualization software (if needed):
| class="model_col2"| no
|}
 
==== Process description ====
 
{| class="wikitable"
| class="model_col1"| Processes represented by model:
| class="model_col2"| Steady-state river generated hypopycnal sediment plume
|-
| class="model_col1"| Key physical parameters & equations:
| class="model_col2"| 2D advection-diffusion equation
|-
| class="model_col1"| Length scale & resolution constraints:
| class="model_col2"| kilometers to tens of kilometers; resolution typically 10 to 100s of meters
|-
| class="model_col1"| Time scale & resolution constraints:
| class="model_col2"| Daily; Steady-state
|-
| class="model_col1"| Numerical limitations and issues :
| class="model_col2"|
|}
 
==== Testing ====
 
{| class="wikitable"
| class="model_col1"| Available calibration data sets:
| class="model_col2"| Eel River (California), Knight and Bute Inlet (British Columbia)
|-
| class="model_col1"| Available test data sets:
| class="model_col2"| --
|-
| class="model_col1"| Ideal data for testing:
| class="model_col2"| --
|}
 
==== User groups ====
 
{| class="wikitable"
| class="model_col1"| Currently or plans for collaborating with:
| class="model_col2"| None
|}
 
==== Documentation ====
 
{| class="wikitable"
| class="model_col1"| Key papers of the model:
| class="model_col2"| Syvitski et al., 1998 ([http://dx.doi.org/10.1016/S0098-3004(97)00084-8 doi:10.1016/S0098-3004(97)00084-8])
|-
| class="model_col1"| Is there a manual available:
| class="model_col2"| no
|-
| class="model_col1"| Model website if any:
| class="model_col2"| --
|}
 
==== Additional comments ====
 
{| class="wikitable"
| class="model_col1"| Comments:
| class="model_col2"| --
|}
 
=== Issues ===
 
=== Help ===
== Help on Model Output ==


PLUME generates a comma separated (*.csv) file which shows sedimentation rates per specified grain size class in m/day for the entire model grid.
PLUME generates a comma separated (*.csv) file which shows sedimentation rates per specified grain size class in m/day for the entire model grid.
Line 279: Line 136:


To create a planview map the plumes:
To create a planview map the plumes:
<geshi lang=matlab>
<syntaxhighlight lang=matlab>
> imagesc(c);
> imagesc(c);
</geshi>
</syntaxhighlight>
Usually the first grain size is the coarser fraction traveling in suspension, but it is dependent on your input file. If you just want to create a planview map of a single grainsize plume:
Usually the first grain size is the coarser fraction traveling in suspension, but it is dependent on your input file. If you just want to create a planview map of a single grainsize plume:


<geshi lang=matlab>
<syntaxhighlight lang=matlab>
> cg1=c(101:200,:);
> cg1=c(101:200,:);
> cg2=c(301:400,:);
> cg2=c(301:400,:);
Line 290: Line 147:
> figure
> figure
> imagesc(cg2)
> imagesc(cg2)
</geshi>
</syntaxhighlight>


If you want the plume area of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:
If you want the plume area of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:


<geshi lang=matlab>
<syntaxhighlight lang=matlab>
> n=sum(sum(cg1>0.0005);
> n=sum(sum(cg1>0.0005);
> plume_area=n*10*10;
> plume_area=n*10*10;
</geshi>
</syntaxhighlight>


Similarly, if you want the plume volume of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:
Similarly, if you want the plume volume of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:


<geshi lang=matlab>
<syntaxhighlight lang=matlab>
> v=sum(cg1>0.0005);
> v=sum(cg1>0.0005);
> plume_volume=v*10*10;
> plume_volume=v*10*10;
</geshi>
</syntaxhighlight>


It may be usefull to compare X-sections for the different grain size classes:
It may be usefull to compare X-sections for the different grain size classes:
<geshi lang=matlab>
<syntaxhighlight lang=matlab>
> plot(cg1(25, :), 'r');
> plot(cg1(25, :), 'r');
> hold
> hold
> plot(cg2(25,:));
> plot(cg2(25,:));
</geshi>
</syntaxhighlight>
 
=== Input Files ===
 


== Input Files ==


=== Download ===
== Download ==
{{Download Model|plume}}
{{Download Model|plume}}
<br>
<br>


==== Simple Input Files for Testing ====
=== Simple Input Files for Testing ===


A very simple scenario for the plume model has been posted here. It models a bankfull flood event for a small meandering river, draining into a lake. These files may be useful for testing whether the model is running.
A very simple scenario for the plume model has been posted here. It models a bankfull flood event for a small meandering river, draining into a lake. These files may be useful for testing whether the model is running.
Line 331: Line 186:


<br>
<br>
==== Simple Output File for Testing ====
=== Simple Output File for Testing ===


The associated output file, with grids of sedimentation rates per day, for both grain sizes can be found here:
The associated output file, with grids of sedimentation rates per day, for both grain sizes can be found here:
Line 337: Line 192:
* [[media:MC8.csv.zip| plume output file]] (January 2009, csv format)
* [[media:MC8.csv.zip| plume output file]] (January 2009, csv format)


=== Source ===
== Source ==
''plume'' is part of the ''sedflux'' model.  Although it is also available as a seperate distribution, its source code is contained within the sedflux repository.
''plume'' is part of the ''sedflux'' model.  Although it is also available as a seperate distribution, its source code is contained within the sedflux repository.


To browse the repository, point your browser to: [http://csdms.colorado.edu/viewvc/sedflux/?root=sedflux http://csdms.colorado.edu/viewvc/sedflux/?root=sedflux]
To browse the repository, point your browser to: [http://csdms.colorado.edu/viewvc/sedflux/?root=sedflux http://csdms.colorado.edu/viewvc/sedflux/?root=sedflux]


==== Command-Line Access ====
=== Command-Line Access ===


If you plan to make changes, use this command to check out the code as yourself using HTTPS:
If you plan to make changes, use this command to check out the code as yourself using HTTPS:


<geshi lang=bash>
<syntaxhighlight lang=bash>
# Project members authenticate over HTTPS to allow committing changes.
# Project members authenticate over HTTPS to allow committing changes.
svn checkout https://csdms.colorado.edu/svn/sedflux/
svn checkout https://csdms.colorado.edu/svn/sedflux/
</geshi>
</syntaxhighlight>


When prompted, enter your CSDMS Subversion password.
When prompted, enter your CSDMS Subversion password.
Line 357: Line 212:
To obtain a CSDMS Subversion account or to become a member of this project, please email [mailto:csdms@colorado.edu csdms@colorado.edu].
To obtain a CSDMS Subversion account or to become a member of this project, please email [mailto:csdms@colorado.edu csdms@colorado.edu].


==== GUI and IDE Access ====
=== GUI and IDE Access ===


This project's Subversion repository may be accessed using many different client programs and plug-ins. See your client's documentation for more information.
This project's Subversion repository may be accessed using many different client programs and plug-ins. See your client's documentation for more information.


==== Subversion Help ====
=== Subversion Help ===


For help on how to use Subversion, an excellent manual is available online at [http://svnbook.red-bean.com/ http://svnbook.red-bean.com/]
For help on how to use Subversion, an excellent manual is available online at [http://svnbook.red-bean.com/ http://svnbook.red-bean.com/]

Revision as of 11:34, 17 September 2009

Contact

Name Eric Hutton
Type of contact Model developer
Institute / Organization CSDMS, INSTAAR, University of Colorado
Postal address 1 1560 30th street
Postal address 2
Town / City Boulder
Postal code 80305
State Colorado
Country USA"USA" is not in the list (Afghanistan, Albania, Algeria, Andorra, Angola, Antigua and Barbuda, Argentina, Armenia, Australia, Austria, ...) of allowed values for the "Country" property.
Email address huttone@colorado.edu
Phone
Fax



Plume


Metadata

Summary

Also known as
Model type Single
Model part of larger framework
Note on status model
Date note status model

Technical specs

Supported platforms
Unix, Linux, Mac OS, Windows
Other platform
Programming language

C

Other program language
Code optimized Single Processor
Multiple processors implemented
Nr of distributed processors
Nr of shared processors
Start year development 1997
Does model development still take place? Yes
If above answer is no, provide end year model development
Code development status
When did you indicate the 'code development status'?
Model availability As code
Source code availability
(Or provide future intension)
Through CSDMS repository
Source web address
Source csdms web address
Program license type Apache public license
Program license type other
Memory requirements Minimal
Typical run time Seconds


In/Output

Describe input parameters River velocity, width, depth; Sediment concentrations
Input format ASCII
Other input format
Describe output parameters Grid of Sediment rate in m/day for specified grain size classes
Output format ASCII
Other output format
Pre-processing software needed? No
Describe pre-processing software
Post-processing software needed? No
Describe post-processing software
Visualization software needed? No
If above answer is yes
Other visualization software


Process

Describe processes represented by the model Steady-state river generated hypopycnal sediment plume
Describe key physical parameters and equations 2D advection-diffusion equation
Describe length scale and resolution constraints kilometers to tens of kilometers; resolution typically 10 to 100s of meters
Describe time scale and resolution constraints Daily; Steady-state
Describe any numerical limitations and issues


Testing

Describe available calibration data sets Eel River (California), Knight and Bute Inlet (British Columbia)
Upload calibration data sets if available:
Describe available test data sets --
Upload test data sets if available:
Describe ideal data for testing --


Other

Do you have current or future plans for collaborating with other researchers? None
Is there a manual available? No
Upload manual if available:
Model website if any --
Model forum / discussion board
Comments --


Introduction

PLUME [1] [2][3]

Plume simulates hypopycnal plumes generated by a river draining its suspended sediment load into a receiving basin. Satellite images of any river-delta emphasize the importance of river plumes. A plume’s behavior is dependent on the density contrast between the river water and the standing water (Albertson, 1950; Bates, 1953). Ocean water has a high density, and the plumes often flow buoyantly on the surface (hypopycnal). Another complementary model that deals with more rare hyperpycnal flows is [[Model:Sakura|Sakura]. The river’s sediment concentration adds density to the freshwater, but usually the effluent remains below the density of seawater. The shape that a hypopycnal plume will have, depends on a variety of factors:

  1. Angle between the river course and the coastline
  2. Strength and direction of the coastal current
  3. Wind direction influencing local upwelling or downwelling
  4. Mixing tidal or storm energy near the river mouth
  5. Latitude and thus the strength of the Coriolis effect.

The plume equations follow those of Albertson (1950) developed for a jet flowing into a steady receiving basin. Plumes of similar shape but differing concentrations result for each grain size in the model. Fine sand will generally settle rapidly, whereas clay can travel much larger distances. Naturally, this affects the geometry of the deposited sediments on the basin floor.

Plumetheory coarseandfine.jpg

River dimensions, i.e. the channel width, depth and velocity at the river mouth are input conditions. In addition, river sediment concentration and settling velocities for specific grain size classes are input parameters as well. Plume is a steady-state model, meaning that it simulates constant input conditions, representative of a 'unit' event.

River dimensions for plume range over orders of magnitude, small streams of only a few meters wide have been run, as well as large continental scale rivers (for example the Ganges-Brahmaputra). Consequently, the spatial resolution of the grid is highly variable depending on the modeling objective. If plume is used in stand-alone mode, it runs events of a single day. If you are interested in exploring deposits of changing plumes over time you will need to use the PLUME model within the framework of the stratigraphic model Sedflux.

MC8bothGSD.jpg


References to plume theory

Albertson, M.L., Dai, Y.B., Jensen, R.A., Hunter, R., 1950, Diffusion of submerged jets. American Society Civil Engineers Trans, v. 115, p. 639-697.

Bates, C.C., 1953, Rational theory of delta formation. AAPG Bulletin, v. 37, p. 2119-2162.

Plume model papers

  1. Hutton and Syvitski, 2008. Sedflux-2.0: An advanced process-response model that generates three-dimensional stratigraphy. Computers and Geosciences, v. 34. doi:10.1016/j.cageo.2008.02.013
  2. Syvitski et al., 1998. PLUME1.1: Deposition of sediment from a fluvial plume (doi:10.1016/S0098-3004(97)00084-8
  3. Peckham, S.D., 2008. A new method for estimating suspended sediment concentrations and deposition rates from satellite imagery based on the physics of plumes. Computer & Geosciences, 34, 1198-1222. doi:10.1016/j.cageo.2008.02.009

Issues

Help

Help on Model Output

PLUME generates a comma separated (*.csv) file which shows sedimentation rates per specified grain size class in m/day for the entire model grid.

Half of the model grid is land, the other half is the receiving marine or lake basin as shown in the accompanying figure. The sedimentation rate for the first specified grainsize is listed for every gridcell. The small plume in the example is visible in the middle of the grid.

Then the grid repeats itself for the next grainsize, of which the plume has a different shape and sedimentation rate. So a grid of 2000m basin width and length, with gridcells of 10 by 10 m will have 400 rows by 200 columns, if the simulation was only run for two grain size classes.

This would repeat on to match the total number of simulated grain size classes.

GriddimensionsFIG.jpg

Help on Simple Matlab Visualization of Output

If you want to use PLUME output with Matlab you will have to cut off the header lines of the file in a text editor. These are just a couple of commands to get started with analyzing the output.

In Matlab use the following commands:

To import your generated output file: <geshi lang=matlab> > c=dlmread('*.csv'); </geshi>

To create a planview map the plumes:

> imagesc(c);

Usually the first grain size is the coarser fraction traveling in suspension, but it is dependent on your input file. If you just want to create a planview map of a single grainsize plume:

> cg1=c(101:200,:);
> cg2=c(301:400,:);
> imagesc(cg1)
> figure
> imagesc(cg2)

If you want the plume area of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:

> n=sum(sum(cg1>0.0005);
> plume_area=n*10*10;

Similarly, if you want the plume volume of all cells with more than a certain threshold of sediment deposited, in this example all area with more than 0.005m deposits, you can get the no-of-gridcells that contain layers thicker than a certain threshold:

> v=sum(cg1>0.0005);
> plume_volume=v*10*10;

It may be usefull to compare X-sections for the different grain size classes:

> plot(cg1(25, :), 'r');
> hold
> plot(cg2(25,:));

Input Files

Download

Template:Download Model

Simple Input Files for Testing

A very simple scenario for the plume model has been posted here. It models a bankfull flood event for a small meandering river, draining into a lake. These files may be useful for testing whether the model is running.



Simple Output File for Testing

The associated output file, with grids of sedimentation rates per day, for both grain sizes can be found here:

Source

plume is part of the sedflux model. Although it is also available as a seperate distribution, its source code is contained within the sedflux repository.

To browse the repository, point your browser to: http://csdms.colorado.edu/viewvc/sedflux/?root=sedflux

Command-Line Access

If you plan to make changes, use this command to check out the code as yourself using HTTPS:

# Project members authenticate over HTTPS to allow committing changes.
svn checkout https://csdms.colorado.edu/svn/sedflux/

When prompted, enter your CSDMS Subversion password.

Non-members may only check out a read-only working copy of the project source.

To obtain a CSDMS Subversion account or to become a member of this project, please email csdms@colorado.edu.

GUI and IDE Access

This project's Subversion repository may be accessed using many different client programs and plug-ins. See your client's documentation for more information.

Subversion Help

For help on how to use Subversion, an excellent manual is available online at http://svnbook.red-bean.com/