How to fill sinks in DEM with specific amount of water

I identified the sinks in my terrain based on the DEM (grid) by using the 'fill'-tool and substracting the original DEM from the filled DEM. I also identified watersheds as well as flow paths.

Now i want to model which water level each sink will have in case of a specific rain event. Can any of you recommend a understandable tutorial how to do that in ArcGIS 10.2 or explain the procedure for me?

• Do you plan to model in ArcGIS which water level will be the result of a specific rain event or do you get the level from another software and only want to visualize it in ArcGIS? – Martin Nov 20 '14 at 10:04
• The idea was to model the water level resulting from a specific rain event directly in ArcGIS. – chris3o12 Nov 20 '14 at 13:46
• I don't know of any functions in ArcGIS that can model water volumes, only waterways, catchments etc (which you have done already). For that you need something like Mike Flood. You can, however, from a given volume find out which water level it corresponds to in a specific sink. From that you can visualize the water depth in the area and a couple of more things, but none of that is of much use unless you have done the modelling first. – Martin Nov 20 '14 at 13:58

This is actually an exceedingly complex problem and not one that you're likely going to be able to solve using conventional ArcGIS tools. To do this, you'll need to develop for each depression in your landscape the relation between depth and volume, i.e. the depth vs volume curve, which will be uniquely defined for each depression based on it's form. To calculate this curve you would essentially need to enter the height of each grid cell contained within a depression into a sorted list. Then visit each entry, from lowest to highest. For each new entry i of height zi calculate the sum of (zi - zj) x g2, where j is a grid cell within the depression of lower elevation than zi. This is an estimate of the volume of the depressional feature at a depth d = zi - z0 where z0 is the lowest elevation in the feature.

Now, since you've mentioned that you've defined the flowpaths and watersheds of each depression I assume that you plan on relating storm events which result in a specified depth of rainfall (p) over your defined watershed of area a, such that you can get a volume (p x a) and then use the volume-depth relation described above to fill the depression to a certain level, or depth. That's great and I've certainly seen this approach applied before, but you're missing an important complexity here (several actually). First, because a depression may well have several other upslope depressions within their watersheds, the watershed itself must be dynamically defined, e.g. if you have a storm of 15 mm then this upslope depression will become filled and its watershed will then contribute to and expand the watershed of this downslope depression. This quickly becomes an exceedingly complex thing to calculate although I have seen it done before. The next main thing that you may be overlooking is that the watershed areas that we define using DEMs in GIS are not really hydrologically realistic, i.e. the entire watershed for a downslope point will rarely if ever contribute flow to that point at any given time. Instead the catchment area that drains to a point of interest in a catchment itself will expand and contract based on moisture conditions much like the dynamic surface-saturated area concept. It's possible of course to try to model this dynamic expansion/contraction of the contributing areas to each depression using a hydrological model like TOPMODEL, but this complexity combined with the fact that you have fill-and-spill of upslope depressions resulting in thresholds or step-changes will mean that you are likely well outside the scope of what can be accomplished within a traditional GIS like ArcGIS and likely will need to program a custom solution.

It's an exciting problem and one that I have spent a great deal of time thinking about in the past. I recall a researcher at the Centre for Hydrology, the University of Saskatchewan named Dr Kevin Shook was looking at something very similar to this. You may want to contact him directly. I wish you the best of luck in pursuing this interesting problem.

In response to Jakob's suggestion to use Cut Fill, you could do a comparison of the original DEM and a depression-less DEM (a DEM resulting from using the Fill tool) and that should give you a volume for any given area showing as a sink on the original DEM.

Design a flat raster surface for each height then run the Cut Fill tool. This will also give you volumes if needed.

The Cut Fill tool summarizes the areas and volumes of change from a cut-and-fill operation. By taking surfaces of a given location at two different time periods, it identifies regions of surface material removal, surface material addition, and areas where the surface has not changed. Source: ArcGIS 10.2.1 Help