米歇尔(2004年)的先前研究展示了如何通过长期的3H记录估算来自上游密苏里河流域的水的平均滞留时间。米歇尔(2004年)使用了一个二元混合模型来模拟河流中的3H浓度,其中一个组分代表近期降水(即1年内的即时流出),第二个组分代表来自流域长期地下水库的水。这个例子展示了如何使用TracerLPM进行类似的分析。Michel(2004年)输入了1963年至1997年之间密苏里河中的氚浓度。
通过对来自Lincoln, Nebraska和Bismarck, North Dakota的Michel(1989年更新)的月度降水3H记录求平均值,估算了降水中的氚浓度(图35)。此处构建的降水中的3H记录的高峰浓度比Michel(2004年)报告的要高,尽管使用了类似的近似方法。
如图35所示,密苏里河中的3H浓度在核爆高峰年份约为降水中的四分之一。然而,河水中的3H浓度下降速度比降水中的下降速度慢,到1966年,河水中的3H浓度超过了降水中的浓度。流入河流的水是来自流域内各种水文系统的近期降水和水的混合物;因此,河水中的3H浓度代表了近期降水的滞留时间(可能为零或接近零)和通过地下渗透并在排入河流之前移动的3H的较长滞留时间的混合物。
图35. 氚浓度在降水(蓝线)和内布拉斯加州内布拉斯加市密苏里河(黑色圆圈)的比较。降水中的氚浓度是通过对内布拉斯加州林肯和北达科他州俾斯麦的月平均氚降水记录进行估算得出的。图表修改自Michel(2004)。
Michel(2004年)使用的二元混合模型中,最近降水中的3H是通过具有平均年龄为0年的PFM进行模拟的,而地下水储层中的3H则是通过EMM进行模拟的。TracerLPM被配置为仅能优化BMM的一个组成部分,而这个组成部分是模型名称中给出的第一个组成部分。
这要求第二个组成部分的平均年龄(和模型参数)被指定。为了适应PFM组分具有零或接近零的平均年龄的假设(示踪剂的瞬时排放),在TimeSeriesFits工作表上选择进行优化的模型是BMM-EMM-PFM。通过选择平均年龄在2到8年之间变化以及混合分数在0.5到1之间变化的约束,对密苏里河模型进行了优化。这些约束产生了EMM的平均年龄为4.3年,混合分数(老化分数)为0.84(或84%的地下水;图36)。 Michel的分析得出了90%的地下水,平均年龄为4年,与TracerLPM非常接近。结果之间的差异很可能与Michel使用的3H降水记录和本分析中使用的记录之间的细微差异有关,正如前面所指出的,以及TracerLPM能够将搜索域离散化为比Michel手动校准更精细的分辨率。
Michel使用的氚浓度在4400 TU左右达到峰值,而此处使用的3H输入记录的峰值浓度约为4900 TU。因此,老地下水中的3H浓度可能比Michel预测的要高。因此,为了匹配河流中测得的3H,需要较低的老水贡献。尽管如此,这个例子说明了TracerLPM如何在时间序列模式下使用,以及一些二元混合模型如何对流域滞留时间分析有用。BMM-EMM-PFM模型可以描述降水迅速通过径流或浅层渗透层传递到河流的情况,并且下面是具有较长滞留时间的含有混合的老地下水的含水层。
图36. 上游密苏里河长期氚记录与二元混合模型结果比较。二元混合模型由EMM和PFM组成:BMM-EMM-PFM。EMM的平均年龄为4.3年,占样本的84%,PFM的平均年龄为0。
图37。模拟得出的密苏里河流域硝酸盐浓度响应的假设结果,针对两种不同情景的虚构硝酸盐输入历史。蓝线表示充沛期的硝酸盐输入历史,红线表示河流中没有反硝化作用的硝酸盐响应,而紫线则表示地下水具有每年0.1的衰减率时的硝酸盐响应。
图38. 上密苏里河的模拟年龄分布。