the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Water Ecosystems Tool (WET) 1.0 – a new generation of flexible aquatic ecosystem model
Nicolas Azaña Schnedler-Meyer
Tobias Kuhlmann Andersen
Fenjuan Rose Schmidt Hu
Karsten Bolding
Anders Nielsen
Dennis Trolle
Download
- Final revised paper (published on 13 May 2022)
- Preprint (discussion started on 21 Jan 2022)
Interactive discussion
Status: closed
-
RC1: 'Comment on gmd-2021-366', Anonymous Referee #1, 07 Feb 2022
Schnedler-Meyer et al. introduced an open source, highly customizable aquatic ecosystem model, the Water Ecosystems Tool (WET), in this work. The WET model is developed based on the FABM-PCLake aquatic ecosystem model framework but has several advantages, including that this new tool is extensively modularized, empowering users and with flexibility of food web configurations and has new features of nitrogen fixation and vertical migration. The manuscript is well writen and easy to follow. I want to thank the authors for it. In my view, this work deserves to be published because it has the potential to help fill the gap between scientific findings and pratical applications in aquatic sicence. However, there are several outstanding questions to be addressed before acceptance.
First, to make the tool useful for broad users, it should be easy to deploy and configure. However, for many open-source tools, the developers focus on more on development but documentation. It is understandable for those new developments. But for a tool that is built on a mature model framework, I did not expect it. I downloaded the source code and test case from the provided link. However, there are none of documentations in the folder that teach the users how to compile and configure. I typed "make" in the folder but cannot compile WET successfully. For most of the users, they will give up after this first attempt. It is unfortunate for a tool that can benefit the society.
Second, the authors stated that WET can be used to test for the optimal food web configuration in a specific case. But they have very limited discussions about it. For example, how should we distinguish the effect of model calibration and module settings in the situation that both can improve the model performance? As a complex model presented here, there may be tens of different modules. Each module may have more than ten parameters. Under such complex situations, the optimal food web configuration is easily said by done. So I want to hear some insights from the authors.
Third, does WET have any unit testing features? As the tool is extensively modulized and supposed to be under community development, unit testing would be a key procesure to ensure software quality.
Specific Comments:
L53: supports
L122-124: Please split this sentence which is too complex.
L166: should be "never limit"
L197-199: please split this sentence. It is too complex.
L245: remove "most"
L268: what does "Each water layer included a sediment layer of 10 cm" mean? Please illustrate.
L279: Are boundary conditions only set at surface layer or all water layers? If the latter, how the inflow are distributed across different layers?
L284: In Chen et al. (2020), ACPy was said to be used for calibration. What is the relationship between ACPy and parsac?
L398: "try to" not "try"
L306: runs at a lower resolution?
L309: that contains
L336: Section 5 needs to be divided into several sub-sections, for example model performance, model limitation and future work.
L381: "improve" not "improving"
Citation: https://doi.org/10.5194/gmd-2021-366-RC1 -
AC1: 'Reply on RC1', Nicolas Schnedler-Meyer, 16 Feb 2022
Dear Referee #1
Thank you for your support of the manuscript and your insights into how it may be improved!
We have considered your questions and will reply to these individually below. With regards to smaller comments not specifically mentioned below, we have generally followed your suggestions.
First comment:
While we realize that this was not clearly stated anywhere in the manuscript, comprehensive documentation on how to download and compile WET is available from the WET website at wet.au.dk under ‘for developers’. We have highlighted this fact in the code availability section of the manuscript, and will refer to this guide in the README.txt file as well. Regarding model usability, we have revised the description of QWET (i.e. the graphical user interface for the GOTM-WET model complex) to reflect that a tool is available for model users to ease model configuration and execution (see line 66).
Second comment:
We might have presented the wrong impression here, by implying that there is a tested workflow for doing this. Unfortunately, this is indeed an area where the researcher must make use of their own judgement. We tried to allude to these thoughts in section 5, lines 354-358. We have changed line 94-95 to be more precise, and improved on the discussion of this aspect in section 5, lines 369-376.
Third comment:
There is as of yet no comprehensive unit testing suite available for WET. WET is currently maintained and developed by a small team of researchers, and our resources do unfortunately not cover this. Thus, users must be prepared to check their results for unexpected or erroneous behaviors, and are encouraged to post any concerns as e.g. support questions on the gitlab website. However, as WET is built for the FABM framework, it does benefit from the excellent error handling within FABM. In addition, some tools for e.g. stress-testing models are available for the FABM framework, but their description is unfortunately beyond the scope of this manuscript.
L268: what does "Each water layer included a sediment layer of 10 cm" mean? Please illustrate.
In this case, this is a feature of the lake GOTM’s hypsograph setup. In order to capture lake sediment-water column interactions at all depths, the bottom is effectively split up, such that each model layer in the 1D setup has an attached bottom layer. Interactions between the water column of a layer, its attached bottom, and the water column layer below is governed by the hypsograph, which specifies the fraction of the bottom area to total layer area. We have added this short description to the section, and included a reference which contains a further description of the hypsograph.
We would however like to avoid tying the mind of the reader too tightly to the chosen physical model used, as WET can in principle be utilized with a variety of physical models, and so we would prefer to avoid a figure of the hypsograph setup in this manuscript.
L279: Are boundary conditions only set at surface layer or all water layers? If the latter, how the inflow are distributed across different layers?
We have specified that both in- and outflow was applied in the top layer.
L284: In Chen et al. (2020), ACPy was said to be used for calibration. What is the relationship between ACPy and parsac?
ACPy is identical to the Parsac package, except for the name, which has been changed in the meantime.
L306: runs at a lower resolution?
We are unsure what you mean here, but to clarify, we ran the model at higher resolution, in order to increase detail on e.g. the vertical movement dynamics.
L336: Section 5 needs to be divided into several sub-sections, for example model performance, model limitation and future work.
We have separated this section into four sub-sections.
Citation: https://doi.org/10.5194/gmd-2021-366-AC1
-
AC1: 'Reply on RC1', Nicolas Schnedler-Meyer, 16 Feb 2022
-
RC2: 'Comment on gmd-2021-366', Anonymous Referee #2, 28 Feb 2022
Review report for Manuscript ID gmd-2021-366 by Schnedler-Meyer et al. entitled:
Water Ecosystems Tool (WET) 0.1.0 - a new generation of flexible aquatic ecosystem model
General comment:
The present study elaborated to introduce a new aquatic ecosystem model - WET, based on the previous model FABM-PCLake. A number of major new features are added to the model such as the flexible food web configuration, vertical migration of organisms, and N fixation, which all make the new model more powerful in addressing key ecosystem principles in lakes. I find this working quite interesting and believe it will become a major step forward of the WET model family, and serve as a milestone for the field of aquatic ecosystem modeling. Nevertheless, I do have several major concerns as I have explained below in detail, and based on these concerns, I suggest that the manuscript need a major revision before acceptance by the journal of GMD.
- My major concern is a lack of sufficient validation for the new functions to the model in the current manuscript. Clearly these new features, such as vertical migration and N fixation of cyanobacteria, are crucial for modeling the food web interaction and ecosystem dynamics, which jointly improve the applicability of the WET model. However, the test case of Lake Bryrup and the available dataset are insufficient to explicitly test the adequacy of the new functions. The lake is also a bit too shallow to reflect the features of the deep lakes with pronounced vertical migration of organisms. More importantly, the main advance of the ‘flexible food web’ is also not addressed in the test case. Nevertheless, I can understand that field data to validate these new features are rare, particularly the spatial distributions of the organism along the vertical dimension. Therefore, I recommend to use the case of Lake Bryrup, but rather than validating the model, the main focus here should be on providing a more detailed evaluation of the new features. In additional to the existing analyses, the authors may consider the following (but not limited to): 1) compare the modeled cyanobacteria with and without N fixing function (ideally with two types of cyanobacteria), and illustrate how the N limitation regulates the growth of the cyanobacteria and potentially other algae groups; 2) design a gradient of the vertical migration parameters for zooplankton, as shown currently in Fig. 5 as an nice example (but now only with one parameter set). Such additional evaluations will help readers to better understand the importance of the new features for modeling and unravel the sensitivity of the key parameters, thereby opening up the opportunities to test the new model in other lake cases with more comprehensive datasets in the future.
- I found that some details of the new features are not clearly described. This could be fine for current users but may puzzle the new users in the future, thus obscure the application and further development of the model. The authors may consider using more Tables to list the new features of some modules with the explanation of the related parameters and meaning of the values in the configuration files. Below I list some specific points in the minor comments, but overall I would appreciate taking advantages of Tables in such model description papers as the key document.
- It remains unclear how to practically change the structure of the model by adding or deleting organism components. Likely it will be done in the configuration file, but it is not described (if I did not miss anything). Tutorial cases and/or a detailed manual for such configuration procedure should be provided because such flexible feature is one of the major advantages of the new WET over the previous version and also many other aquatic ecosystem models. Besides, I recommend to briefly discuss the possibility to implement the flexible food web in a more user-friendly way (such as the GUI of ECOPATH or AQUATOX), if such function would be envisioned for further model development.
- It is uncommon to number the first version of the WET by ‘0.1.0’. The first release of software or model is usually numbered as ‘1.0’ or ‘1.00’ etc. Current version of ‘0.1.0’ gives the impression that the model is still at the very beginning stage and far from completed for first release. Please consider changing the version numbering, if appropriate.
Below please find additional minor points for the author to consider:
- Line 57: -coding -> recoding
- Line 66-67, this is a duplication of the information in line 59-61. Recommend removing one of them.
- Line 84: I am not sure if all the biogeochemical processes equation are ‘unchanged’ in WET compared to PCLake. If so, it is contradictory to the previous statement that processes such as resuspension are changed compared to previous model. Better add ‘mostly unchanged’, if proper.
- Figure 2: suggest to change the subtitle in vertical ‘FOOD WEB MODULES’ to ‘FLEXIBLE FOOD WEB MODULES’ to highlight the features of the model.
- Line 105, this paragraph lacks the description whether new phytoplankton groups can be added to the community, and whether they can be at least pre-parameterized with readily available parameter values (like a database). In addition, the new algal properties, e.g. Si dependency, N fixation, are better to list in the Table with details of how they can be switched on/off in the configuration file. This would be highly useful for the users.
- Line 110, the same issues here: rather than the descriptive text, a Table summarizing the new features and method to implementation in the configuration file would be helpful.
- Equation 1: reverse the super- and subscript of the ∑ (i.e. i=1 and nPISC), and delete the first product sign (×) in the inner equation, both in the numerator and denominator.
- Line 138: need to clarify how to calculate the terms such as aFunVegAss and uFunTmFish, or at least refer to related references.
- Line 142: population structure, or age structure? Please clarify this statement.
- One comment about the N fixation modeling in section 2.2: removing the ‘aNLim’ term in the growth modeling of algae with N fixing may be not fully realistic. Recent field observations imply that N fixation cannot compensate for the N depletion in lakes with reduced external N laoding over a long-term scale, and phytoplankton still exhibits N limitation in summer (see Shatwell et al. 2019 L&O for the long-term study in Mueggelsee, Germany). No need to change anything here, but maybe keeping this in mind and slightly discuss in the end of the manuscript.
- Line 185-195, This paragraph fits better to the introduction rather than the model description section here.
- Line 213 and 217, better describe each qTrans option in one individual paragraph. It is confusing to read such as (qTrans = 2 or 4) or (qTrans = 3 or 4) at the first glance. The same applies to the zooplankton and fish migration section below.
- Section 4.2, I feel that the results demonstration in this section can be improved by (1) adding the location of the thermocline during summer and (2) making parallel panels showing the model results of original FABM-PCLake without the newly added features such as the vertical migrations. This can be used as part of the additional evaluation of the model, as I mentioned in the major comments above.
Citation: https://doi.org/10.5194/gmd-2021-366-RC2 -
AC2: 'Reply on RC2', Nicolas Schnedler-Meyer, 17 Mar 2022
Dear Referee #2
Thank you for your very valid and constructive criticisms and suggestions. We have considered your questions and will reply to these individually below. With regards to smaller comments not specifically mentioned below, we have generally followed your suggestions.
First major comment:
- My major concern is a lack of sufficient validation for the new functions to the model in the current manuscript. Clearly these new features, such as vertical migration and N fixation of cyanobacteria, are crucial for modeling the food web interaction and ecosystem dynamics, which jointly improve the applicability of the WET model. However, the test case of Lake Bryrup and the available dataset are insufficient to explicitly test the adequacy of the new functions. The lake is also a bit too shallow to reflect the features of the deep lakes with pronounced vertical migration of organisms. More importantly, the main advance of the ‘flexible food web’ is also not addressed in the test case. Nevertheless, I can understand that field data to validate these new features are rare, particularly the spatial distributions of the organism along the vertical dimension. Therefore, I recommend to use the case of Lake Bryrup, but rather than validating the model, the main focus here should be on providing a more detailed evaluation of the new features. In additional to the existing analyses, the authors may consider the following (but not limited to): 1) compare the modeled cyanobacteria with and without N fixing function (ideally with two types of cyanobacteria), and illustrate how the N limitation regulates the growth of the cyanobacteria and potentially other algae groups; 2) design a gradient of the vertical migration parameters for zooplankton, as shown currently in Fig. 5 as an nice example (but now only with one parameter set). Such additional evaluations will help readers to better understand the importance of the new features for modeling and unravel the sensitivity of the key parameters, thereby opening up the opportunities to test the new model in other lake cases with more comprehensive datasets in the future.
Our reply:
We appreciate that the reviewer recognizes the challenges involved, but agree that the model tests can be improved. We are thankful for the excellent suggestions and recommendations on how to improve this aspect of the manuscript. We have gone with the suggestion of the reviewer and provided a more thorough presentation of the new features, in accordance with suggestion 1) and 2) of the reviewer’s comment.
- As suggested by the reviewer, we have added a second, N-fixing cyanobacteria to the model, and compare the dynamics of these in Lake Bryrup in a new figure (attached to this reply). As Lake Bryrup is mainly P-limited during the growth season (as are most temperate lakes), we have also run scenarios with N-inflow to the lake disabled, and report on the dynamics observed.
- We interpret this suggestion to mean that the reviewer would like us to compare the dynamics of the different vertical migration options in the model. To this end, we have introduced two new panels in figure 5 (now 6, attached to this reply), illustrating the dynamics of passive transport, and hypoxia avoidance by zooplankton, and comment on the results in the results and discussion.
Second major comment:
- I found that some details of the new features are not clearly described. This could be fine for current users but may puzzle the new users in the future, thus obscure the application and further development of the model. The authors may consider using more Tables to list the new features of some modules with the explanation of the related parameters and meaning of the values in the configuration files. Below I list some specific points in the minor comments, but overall I would appreciate taking advantages of Tables in such model description papers as the key document.
Our reply:
We have done as suggested, and added a table detailing the parameters and options related to each new feature.
Third major comment:
- It remains unclear how to practically change the structure of the model by adding or deleting organism components. Likely it will be done in the configuration file, but it is not described (if I did not miss anything). Tutorial cases and/or a detailed manual for such configuration procedure should be provided because such flexible feature is one of the major advantages of the new WET over the previous version and also many other aquatic ecosystem models. Besides, I recommend to briefly discuss the possibility to implement the flexible food web in a more user-friendly way (such as the GUI of ECOPATH or AQUATOX), if such function would be envisioned for further model development.
Our reply:
This comment is in line with reviewer #1, who also requested more documentation on model setup and configuration. Although we failed to adequately cover this in the manuscript, there is in fact extensive documentation on how to setup and configure the model on the WET homepage. We have emphazised this more clearly in the manuscript. We have also included a short description on how to add or remove food web modules to the manuscript.
With regards to the possibility of running the model with a more user-friendly interface, the QWET plugin for the (freeware) GIS software QGIS does exactly this. There is also extensive documentation in the QWET section of the WET homepage on how to set up QWET. We have further emphasized these facts in the manuscript.
Fourth major comment:
- It is uncommon to number the first version of the WET by ‘0.1.0’. The first release of software or model is usually numbered as ‘1.0’ or ‘1.00’ etc. Current version of ‘0.1.0’ gives the impression that the model is still at the very beginning stage and far from completed for first release. Please consider changing the version numbering, if appropriate.
Our reply:
In this case, we simply followed the examples of some other model software, but upon reflection, we agree with your comment. We will make another release of the model under the version number 1.0, when and if our manuscript reaches the accepted status.
Peer review completion
re-inventions of the wheelthat are common to our scientific modeling community.