Numerical simulation of water renewal timescales in the Mahakam delta, Indonesia

Water renewal timescales, namely age, residence time, and exposure time, which are defined in accordance with the Constituent-oriented Age and Residence time Theory (CART), are computed by means of the unstructured-mesh, finite element model Second-generation Louvain-la-Neuve Ice-ocean Model (SLIM) in the Mahakam Delta (Borneo Island, Indonesia). Two renewing water types, i.e., water from the upstream boundary of the delta and water from both the upstream and the downstream boundaries, are considered, and their age is calculated as the time elapsed since entering the delta. The residence time of the water originally in the domain (i.e., the time needed to hit an open boundary for the first time) and the exposure time (i.e., the total time spent in the domain of interest) are then computed. Simulations are performed for both low and high flow conditions, revealing that (i) age, residence time, and exposure time are clearly related to the river volumetric flow rate, and (ii) those timescales are of the order of one spring-neap tidal cycle. In the main deltaic channels, the variation of the diagnostic timescales caused by the tide is about 35% of their averaged value. The age of renewing water from the upstream boundary of the delta monotonically increases from the river mouth to the delta front, while the age of renewing water from both the upstream and the downstream boundaries monotonically increases from the river mouth and the delta front to the middle delta. Variations of the residence and the exposure times coincide with the changes of the flow velocity, and these timescales are more sensitive to the change of flow dynamics than the age. The return coefficient, which measures the propensity of water to re-enter the domain of interest after leaving it for the first time, is of about 0.3 in the middle region of the delta.