LULC Default Parameter Basis

July 21, 2026 ยท View on GitHub

LULC_parameters.xlsx provides a consistent starting table for a new HydroPol2D case. It does not provide universal physical constants. Terrain resolution, vegetation structure, season, water depth, soil depth, and local management all affect the effective parameters used by a distributed model.

Surface Routing, Rooting, and ET

LULC classManning n [s m^-1/3]Root depth [m]Kc [-]Basis for the generic value
Tree cover0.1001.501.05Forest roughness within standard 2D floodplain guidance; deep woody root zone.
Shrubland0.0801.000.70Shrub/scrub roughness range; intermediate woody rooting.
Grassland0.0600.600.85Conservative shallow-flow roughness; shallow herbaceous rooting.
Cropland0.0501.001.00Upper end of a generic crop roughness range; mid-season crop ET starting value.
Built-up0.0300.000.00Smooth sealed-surface value. Use a separate urban class with larger n where buildings or obstructions must be represented.
Bare / sparse vegetation0.0300.100.30Barren land roughness range; minimal vegetation water uptake.
Snow and ice0.0200.000.00Initial smooth snow/ice runoff surface; revise for rough ice, debris, or crevassed terrain.
Permanent water bodies0.0350.000.00Within the standard open-water range. Ponded cells evaporate at Ep, not through Kc.
Herbaceous wetland0.0700.401.00Emergent-wetland roughness range and shallow, water-table-limited rooting.
Mangroves0.1501.001.05High-end woody-wetland roughness; revise using local vegetation density and tidal/channel geometry.
Moss and lichen0.0800.050.40Conservative shallow-flow resistance and near-surface water uptake.

The roughness values follow the class ranges in the USACE HEC-RAS 2D land-cover guidance. That guidance also notes that its tabulated values are intended for appreciable flow depths; shallow overland flow can require larger effective values. Root-depth values are simple representative depths consistent with global land-model parameterizations, rather than estimates of the deepest roots at a site. Kc is applied only to internally computed reference ET and is an initial, time-invariant value. Seasonal crop curves, remotely sensed ET, or locally calibrated values are preferable when available.

h0 and d0 both default to zero. h0 is used by the cellular-automata routing option only; it is not a substitute for canopy storage or soil infiltration. d0 represents the initial surface-water state and should be set from a restart or an initial-depth raster when the event begins on a wet surface.

Water Quality

The generic table sets C1 = C2 = C3 = 0 and C4 = 1. This deliberately represents no initial pollutant buildup and no washoff. There is no defensible universal set of buildup/washoff coefficients by land-cover class alone. In the current mass-based HydroPol2D formulation,

B0 = C1 [1 - exp(-C2 ADD)] A / 10000
W = C3 Q^C4 B

where B0 is buildup mass, ADD is antecedent dry days, A is cell area in m2, Q is runoff in m3 s^-1, and B is the available pollutant mass. C1 is in kg ha^-1, C2 is in d^-1, C4 is dimensionless, and the units of C3 depend on C4. Define these values for the pollutant, measurement units, and case study, then calibrate them against concentration or load observations.

Snow

The snow columns use broadly plausible initial values: snow albedo of 0.55-0.80, emissivity of 0.98-0.99, a linear rain/snow transition of -1 to 2 degC, fresh-snow density of 100 kg m^-3, and maximum seasonal-snow density of 450 kg m^-3 (550 kg m^-3 for the snow-and-ice class). Degree-day melt, sublimation, and compaction coefficients are empirical rate parameters. They are retained as explicit inputs because they should be fitted or evaluated against local SWE, snow-depth, or runoff data; they should not be interpreted as transferable LULC constants.

Sources