Next Article in Journal
Chromium Removal in the Presence of NOM during Fe(II) Reductive Precipitation for Drinking Water Treatment
Next Article in Special Issue
Development of Aquatic Index of Biotic Integrity and Its Driving Factors in the Diannong River, China
Previous Article in Journal
Economic Analysis of Flood Risk Applied to the Rehabilitation of Drainage Networks
Previous Article in Special Issue
Effect of Rainfall and pH on Musty Odor Produced in the Sanbe Reservoir
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biomathematical Model for Water Quality Assessment: Macroinvertebrate Population Dynamics and Dissolved Oxygen

by
Jair J. Pineda-Pineda
1,*,
Jesús Muñoz-Rojas
1,
Y. Elizabeth Morales-García
1,2,
Juan C. Hernández-Gómez
3 and
José M. Sigarreta
3
1
Ecology and Survival of Microorganisms Research Group (ESMRG), Laboratorio de Ecología Molecular Microbiana (LEMM), Centro de Investigaciones en Ciencias Microbiológicas (CICM), Instituto de Ciencias (IC), Benemérita Universidad Autónoma de Puebla (BUAP), Puebla 72570, Mexico
2
Facultad de Ciencias Biológicas, Benemérita Universidad Autónoma de Puebla, Ciudad Universitaria, Col. Jardines de San Manuel, Puebla 72570, Mexico
3
Unidad Académica de Matemáticas, Universidad Autónoma de Guerrero, Carlos E. Adame No. 54 Col. Garita, Acapulco 39650, Mexico
*
Author to whom correspondence should be addressed.
Water 2022, 14(18), 2902; https://doi.org/10.3390/w14182902
Submission received: 21 June 2022 / Revised: 5 September 2022 / Accepted: 13 September 2022 / Published: 16 September 2022
(This article belongs to the Special Issue Water Quality Assessment and Ecological Monitoring in Aquatic System)

Abstract

:
Sustainable water management is important to ensure its availability for future generations. The study of water quality is fundamental for this purpose. Assessing the health of aquatic ecosystems through bioindicators has been shown to be reliable and inexpensive. The objective of this work was to evaluate water quality through a biomathematical model that involves environmental stress indicator organisms and their close relationship with dissolved oxygen. In this direction, a system of differential equations describing the population dynamics of aquatic macroinvertebrates under the influence of dissolved oxygen is proposed. The model is validated by its application in the Coyuca Lagoon, Mexico. Likewise, population changes over time were represented, which allowed us to deduce that the increase or decrease in the aeration/oxygenation rate significantly affects the population dynamics of the bioindicator organisms. In addition, to classify water quality, a one-to-one correspondence was established between water quality and the equilibrium points of the system of differential equations. The results obtained allow inferring that the proposed techniques are useful for the study of water quality since they can predict significant changes in the ecosystem and provide researchers and water managers with tools for decision making.

1. Introduction

One of the sustainable development goals for 2030 is to guarantee the availability of water and its sustainable management [1]. To achieve this goal, the health of aquatic ecosystems must be continuously studied and monitored to know their current and future conditions, ensure better quality, and prevent degradation. In this direction, the assessment of the quality of freshwater ecosystems is fundamental, and the development of prospective methods contributes to this purpose.
Some methods for assessing water quality make use of indicators that include biological, physicochemical, and geomorphological parameters (see [2,3,4,5,6], and references therein). In addition, it has been shown that such indicators directly or inversely relate to some parameters [7]. Theoretical representations of observable phenomena in different contexts have also been developed [8,9,10]. Beyond the advantages of these methods, their main limitation is that they only provide a temporal state of the ecosystem, but do not integrate changes in water quality over time. As a solution, differential mathematical models are used, which are more representative of the state of water quality over time.
Dissolved oxygen (DO) is a very important element that has been used in the evaluation of water quality in aquatic ecosystems [11,12]. It is also fundamental for life and as an indicator of ecosystem degradation [13,14,15,16]. Likewise, DO is key to changes in aquatic communities because some organisms respond negatively or positively to different levels in their concentrations [17,18]. For example, eutrophication causes a decrease in water quality [19] and an increase in hypoxia [20], which implies the absence or presence of environmental quality indicator organisms, particularly some aquatic macroinvertebrates (AMs) taxa [21,22]. Furthermore, Galic et al. [23] in their meta-analysis of the adverse impacts of hypoxia on aquatic invertebrates infer that reproductive and respiration rates decrease when DO availability decreases and those organisms rapidly reduce their reproductive capacity in the face of environmental stress.
AMs as biological indicators have proven to be reliable in the study of water quality, habitat disturbances, and restoration of aquatic ecosystems. Likewise, AMs are sensitive or tolerant to high organic pollution loads and increased or decreased DO [24,25,26,27]. Moreover, they have different ecological requirements and degrees of tolerance to natural and anthropogenic disturbances [28,29]. In addition, AMs are used to study the ecological health of natural areas [30] and the effectiveness of the use and application of artificial barriers, which fragment the landscape, as a method of reclamation [31]. Consequently, the application of biotic indices, qualitative and quantitative, using AMs has become popular worldwide (see [32] and references therein), even though such methods are retrospective, and do not help to determine or have a close idea of what may occur in the future. On the other hand, to contribute to the reliability of AMs applications, they are included in the construction of multimetric and topological indices [33,34], and in some prospective methods based on a set of environmental parameters.
Models for assessing water quality involve one or several environmental parameters, with different methods and theories [34,35,36,37,38,39]. These models evaluate water quality in a timely manner or pose possible scenarios over time. For example, the Streeter–Phelps model relates DO deficit and biochemical oxygen demand [40], through a system of two coupled differential equations. In addition, this model is taken as a starting point in a historical review of models for water quality assessment [41]. According to the context, the models are either very simple or complex depending on the number of equations, input data, or time in the sampling effort [16]. Recently, as reported by Costa [42], AQUATOX is one of the most applied models in the last five years, because it predicts the behavior of various pollutants and their effects on the ecosystem. AQUATOX within its set of mass balance equations contains differential equations to observe the change, with respect to time, in benthic biomass and dissolved oxygen concentration, independently [43].
The objective of this work is to study water quality in freshwater ecosystems through time and population fluctuations of environmental stress indicator organisms. In this direction, a differential mathematical model describing the population dynamics of AM indicators of water quality and their relationship with DO is proposed. To this end, Section 2 describes the state variables, parameters, and the system of four coupled nonlinear differential equations that represent the dynamics between AMs and DO; Section 3 shows the main qualitative theoretical results in the ecological context. Finally, Section 4 shows the numerical simulations. In addition, we validate the model by applying it to the water quality assessment of a real case: the Coyuca Lagoon, Mexico.

2. Model Conceptualization

Disturbances in freshwater ecosystems can be studied through the population dynamics of environmental stress indicator organisms, such as AMs, and since the population fluctuations of AMs have a direct or inverse relationship with DO concentration, in water quality assessment, and restoration of aquatic ecosystems, certain taxa can be classified according to their oxygen requirements and tolerance to organic pollution. Thus, we define the sets A = { a 1 , , a n } , B = { b 1 , , b m } , and C = { c 1 , , c s } as Class I, II, and III organisms, respectively, where each class corresponds to the AM families with higher/lower, regular/regular, and lower/higher oxygen requirement and tolerance to contamination. Some images, common names, and stress tolerance classifications (Classes) of AMs can be consulted in [44,45].
The cardinal of A, denoted by # A , is the number of Class I AM families. Similarly, the cardinals for sets B and C are defined. The cardinal of a i A , for i = 1 , , n , denoted by # a i , is the abundance of the aquatic macroinvertebrate family a i . The cardinal of b j B and c k C , for j = 1 , , m and k = 1 , , s , are defined analogously to # a i . Taxonomic richness is the sum of the cardinals of sets A, B, and C. Total abundance is the sum of the cardinals of a i , b j , and c k .
The structure of aquatic macroinvertebrate communities is defined by their abundance, richness, and diversity [46,47,48]. Richness can often be assessed at higher taxonomic levels, e.g., genus, family, and order. In addition, richness measures reflect the diversity of the aquatic assemblage. On the other hand, relative abundance, defined as the ratio of a taxon to all others contained in an ecosystem, determines how rare, common, or dominant a taxon is. A population is defined as a group of organisms of the same species, each population has several characteristics, such as its population density or relative size in an available space. In this study, each population is a conglomerate of AM populations that have similar characteristics in the face of environmental impact, particularly in terms of their DO requirements.
The research on the evaluation of water quality, previously analyzed, has allowed us to infer the basic principles for the development of this research. Thus, the differential mathematical model proposed in this study considers the intraspecific competition between AM families and as a limited resource: the concentration of DO. In addition, the oxygenation/aeration of the water body, either natural or mechanical, reaches a maximum saturation level. Moreover, regarding the DO consumption and population growth of the different classes of AMs, the concentration of DO in the water is diminished by pollution associated with biotic components, particularly AMs, and the amount of energy absorbed by living organisms is greater than that which they use.
Based on the above, the state variables represent: o ( t ) the concentration of DO in the water and x 1 ( t ) , x 2 ( t ) , and x 3 ( t ) , the sum of the cardinals a i , b j , and c k , i.e., the sum of the abundances of the AM families of class I, II, and III, respectively. Note that, the state variables depend on t, but in denoting them, we will dispense with this. The flow diagram of the model is shown in Figure 1.

Mathematical Formulation

The mathematical representation of the above is described by the following system of non-linear differential equations:
o = r 1 o 2 m 2 o i = 1 3 β i x i , x i = r i x i 1 m x i k i o + γ i ( m o ) x i ,
for i = 1 , 2 , 3 , in the positive invariant region
Ω = { ( o , x 1 , x 2 , x 3 ) R 4 : 0 < o m , 0 x i k i , with i = 1 , 2 , 3 } ,
and the parameter set ω = { m , r , β i , r i , k i , γ i R : m , r , β i , r i , k i , γ i > 0 , with i = 1 , 2 , 3 } .
The differential equations associated with system (1) represent the rate of change over time in DO concentration and AM populations. The equations for the AM families make use of the logistic equation since this allows modeling phenomena with limited resources and populations of organisms that reach their maximum density, which is limited by the availability of resources. Moreover, its application in the study of populations has been widely extended [49,50]. All parameters are positive constants and their interpretation is as follows: r is the constant oxygenation/aeration rate of the water body, m is the maximum value to which the DO tends asymptotically in the absence of AMs, and the DO consumption for the different classes of AMs is given by o i = 1 3 β i x i . For the different classes of organisms: β i is the average respiration rate, r i is the average reproduction rate, k i is the average carrying capacity i.e., the average upper limit of population growth [51], γ i is the average rate of DO utilization, for x i with i = 1 , 2 , 3 . Furthermore, it is established that γ i < β i because organisms consume some amount of energy, but do not fully utilize it. Note that if o = m , then x i asymptotically approaches k i , for i = 1 , 2 , 3 . The parameter values or ranges, which will be used in the numerical simulations in Section 4, are shown in Table 1.

3. Qualitative Analysis

Although the system of differential Equation (1) has 23 equilibrium points (two without the presence of AMs, nine with the presence of only one class of AMs, nine with the presence of two classes of AMs, and three with the presence of the three classes of AMs), our study will focus only on those that are in the region of interest and that describe the absence or presence of different classes of AMs; more specifically, we focus on the equilibrium point of coexistence of the three AM classes for the information provided to us.

3.1. Absence, Presence, and Coexistence of AMS

Proposition 1. 
Given the system of the nonlinear differential Equation (1), there is o * ( 0 , m ) , such that the following points of R 4 are the equilibrium points of the system and lie in the region of interest Ω.
(i) 
AM-free equilibrium: E 0 = ( m , 0 , 0 , 0 ) .
(ii) 
Presence equilibrium: E 1 = ( o * , x 1 * , 0 , 0 ) , E 2 = ( o * , 0 , x 2 * , 0 ) , E 3 = ( o * , 0 , 0 , x 3 * ) ,
E 4 = ( o * , x 1 * , x 2 * , 0 ) , E 5 = ( o * , x 1 * , 0 , x 3 * ) and E 6 = ( o * , 0 , x 2 * , x 3 * ) .
(iii) 
Coexistence equilibrium: E 7 = ( o * , x 1 * , x 2 * , x 3 * ) with
x i * = k i m r i o * r i + γ i ( m o * ) f o r i = 1 , 2 , 3 .
Proof. 
To find the equilibrium points, we equal to zero the right-hand sides of system (1) and, thus, obtain the following homogeneous system of nonlinear algebraic equations
0 = r 1 o 2 m 2 o i = 1 3 β i x i ,
0 = r i x i 1 m x i k i o + γ i ( m o ) x i , for i = 1 , 2 , 3 .
Equation (4) can be written as follows
0 = x i r i 1 m x i k i o + γ i ( m o ) , for i = 1 , 2 , 3 .
Which are satisfied if x i = 0 or
x i = k i m r i o r i + γ i ( m o ) , for i = 1 , 2 , 3 .
When substituting x i = 0 , for all i = 1 , 2 , 3 , in Equation (3), and solve it; it is obtained that o = m . Thus, an equilibrium point of the system is E 0 = ( m , 0 , 0 , 0 ) , and as m > 0 , then E 0 Ω , so the proposition has been proven for (i). On the other hand, if x i 0 , substituting expression (5) in Equation (3) we have
0 = 1 m i = 1 3 β i γ i k i r i o 3 r m 2 + 1 m i = 1 3 β i k i + i = 1 3 β i γ i k i r i o 2 + r .
Note that Equation (6) is a polynomial of degree three in the variable o and we intend to calculate its roots. Then, we define the function
f ( o ) = 1 m i = 1 3 β i γ i k i r i o 3 r m 2 + 1 m i = 1 3 β i k i + i = 1 3 β i γ i k i r i o 2 + r .
Since f ( o ) is a polynomial, then it is a continuous function, particularly in 0 , m ; furthermore, f ( 0 ) = r > 0 and
f ( m ) = m i = 1 3 β i k i < 0 , so   that f ( 0 ) · f ( m ) < 0 ,
consequently, by Bolzano’s Theorem, there exists o * ( 0 , m ) , such that f ( o * ) = 0 ; thus, we guarantee the existence of a positive real solution in ( 0 , m ) for f ( o ) . In addition, when o + , it follows that f ( o ) + ; on the other hand, if o , then f ( o ) ; therefore, the function f ( o ) has two positive real solutions, namely o * ( 0 , m ) and o * * > m ; and one negative real solution, o * * * < 0 . Note that for roots o * * > m and o * * * < 0 , the equilibrium points with first coordinate o * * > m or o * * * < 0 , are not in Ω .
Finally, for the root o * ( 0 , m ) , substituting it in the expression (5), we obtain
x i = k i m r i o * r i + γ i ( m o * ) > 0 , for i = 1 , 2 , 3 .
Thus, any equilibrium point having o * as the first coordinate, and at least one of the coordinates x i , as in (5), will be an equilibrium point of system (1), since it satisfies (3) and (4), and will also be in Ω . If we denote as x i * the equilibrium values of the expression (8), the equilibrium points E 1 , E 2 , , E 7 Ω , and, thus, ( i i ) and ( i i i ) are proven. □
The solution E 0 and absence of all classes of AMs represent oxygen depletion or supersaturation in aquatic ecosystems, when the DO tends to zero or when it tends to m, respectively. From expression (3) we can see that, in the absence of AMs, since o = m is an equilibrium value of this scenario, the DO tends asymptotically to the value m, and in the presence of at least one AM family, the level of DO will be less than m.
According to expression (5), we have that
lim o m x i = lim o m k i m r i o r i + γ i ( m o ) = k i para i = 1 , 2 , 3 ,
that is, if o tends to m, the equilibrium values x i will tend to k i . Furthermore, if o ( 0 , m ) and x i 0 for some i = 1 , 2 , 3 , then
k i m r i o r i + γ i ( m o ) < k i .
Expression (9) indicates that the equilibrium value of each x i tends to the value k i for i = 1 , 2 , 3 , while relation (10) suggests that such equilibrium values will not be able to reach the value of k i ; summarizing these two results, we can say that if o tends to the value m, then the equilibrium values x i tend asymptotically to the values k i below, for each i = 1 , 2 , 3 .
The equilibrium point E 7 = ( o * , x 1 , x 2 , x 3 ) , having all AM classes is associated with water bodies with significant richness and abundance, because of this; water quality, as a function of E 7 , accepts the following classifications shown in Table 2.
The strict inequalities, described in Table 2, between the different classes of organisms implies that there is a dominant class; therefore, the one that determines the quality of the water. On the other hand, if the abundances are approximately equal ( x 1 x 2 x 3 ), then the water body is diverse and uniform, which implies that it tends to recover much faster than one that is not. On the other hand, the system requires the presence or absence of different classes of AMs and their relationship with the DO, i.e., the loss of diversity and abundance. Note that each equilibrium solution, explicit in the presence equilibria of the previous proposition, represents the loss of diversity of AM families, in the function of the DO concentration, which allows establishing, in a more specific way, the stress level of a freshwater ecosystem. For example, the absence of Classes II and III families in the equilibrium solution E 1 implies an ecosystem with sufficient DO to sustain the Class I organism communities, the absence of pollution, and excellent water quality. In E 2 , the absence of Class III organisms implies two cases: (a) if x 1 < x 2 , then there is regular pollution with a tendency to decrease and the water quality is good and (b) if x 2 < x 1 , then there is low pollution with a tendency to increase and the water quality is good. Note that the absence of Classes I and II AMs in the equilibrium solution E 6 allows inferring that there is high pollution, and the water quality is very poor. The degree of contamination and water quality for the equilibrium solutions E 3 , E 4 , and E 5 , can be inferred analogously.
The importance of analyzing the stability of equilibrium solutions lies in knowing under what conditions an ecosystem is altered, it returns to its original state. In addition, it is possible to infer which parameters have the greatest influence on the dynamics of the system. This analysis is developed in the following section.

3.2. Stability Analysis

Holling defines ecological resilience as the capacity of a system to absorb disturbances and reorganize while undergoing changes in order to essentially preserve its function, structure, identity, and feedback. In the literature, the main mathematical definitions of resilience are based on the dynamical systems theory, and specifically on attractors and basins of attraction [52]. In this study, the necessary conditions for system (1), after some perturbations, to return to its original state, are established in the following result:
Proposition 2.
Given the system of nonlinear differential Equation (1), the AM-free equilibrium point, E 0 = ( m , 0 , 0 , 0 ) , is unstable. In addition, the coexistence equilibrium point, E 7 = ( o * , x 1 * , x 2 * , x 3 * ) , will be locally stable if m r i γ i , for i = 1 , 2 , 3 .
Proof. 
To analyze the stability of the equilibrium points, we will use the indirect Lyapunov method. For this purpose, we return to the autonomous system of differential equations given by (1):
o = f 0 ( o , x 1 , x 2 , x 3 ) = r 1 o 2 m 2 o i = 1 3 β i x i ,
x i = f i ( o , x 1 , x 2 , x 3 ) = r i x i 1 m x i k i o + γ i ( m o ) x i , for i = 1 , 2 , 3 .
The partial derivatives of (11), with respect to o and x i , are given by:
f 0 o = 2 r m 2 o i = 1 3 β i x i
and
f 0 x i = β i o for i = 1 , 2 , 3 .
On the other hand, the partial derivatives of (12), with respect to o and x i , are given by:
f i o = m r i x i 2 k i o 2 γ i x i
and
f i x j = r i 1 2 m k i o x i + γ i ( m o ) if i = j 0 if i j
for i = 1 , 2 , 3 .
Thus, the Jacobian matrix of system (1) is shown in Appendix A.1.
Evaluating the Jacobian matrix at the equilibrium point E 0 we obtain
J E 0 = 2 r m β 1 m β 2 m β 3 m 0 r 1 0 0 0 0 r 2 0 0 0 0 r 3 .
Note that the matrix in (13) is upper triangular, so its eigenvalues are the elements of the diagonal. Since this matrix has eigenvalues λ 1 = 2 r m , λ 2 = r 1 , λ 3 = r 2 and λ 4 = r 3 , three of which are positive, then, the equilibrium point E 0 is unstable and this proves the first statement.
The Jacobian matrix evaluated at the equilibrium point E 7 is given by
J E 7 = S β 1 o * β 2 o * n β 3 o * S 1 T 1 0 0 S 2 0 T 2 0 S 3 0 0 T 3 ,
where
S = 2 r m 2 o * + i = 1 3 β i k i o * m r i T i ,
T i = r i + γ i ( m o * ) for i = 1 , 2 , 3 ,
S i = k i T i m r i T i o * γ i for i = 1 , 2 , 3 .
Computing the determinant J λ I for the matrix (14), we have
D = S λ β 1 o * β 2 o * β 3 o * S 1 T 1 λ 0 0 S 2 0 T 2 λ 0 S 3 0 0 T 3 λ
Consequently, the characteristic equation of the matrix (18) is given by
λ 4 + a λ 3 + b λ 2 + c λ + d = 0 ,
and their coefficients are denoted in Appendix A.2.
The first observations are o * < m for being the root of Equation (6) with that property, from this we have that T i > 0 for each i = 1 , 2 , 3 , and S > 0 for how it is defined. Now, it will be proven that, if for each i = 1 , 2 , 3 , m r i γ i , then S i > 0 . Suppose that for each i = 1 , 2 , 3 it is satisfied that m r i γ i , then as o * < m , you have to o * < r i γ i , but this implies that γ i o * < r i , i.e., 2 ( r i o * γ i ) > 0 , so that r i + γ i ( r i γ i o * ) o * γ i > 0 , i.e., T i o * γ i > 0 , which implies that S i > 0 for i = 1 , 2 , 3 .
With this, we have proved that the coefficients of Equation (19) satisfy a , b , c , d > 0 . Moreover, analogous to the previous statement we can easily prove that a b c c 2 a 2 d > 0 ; therefore, Expression (19) is a Hurwitz polynomial. Now, we will apply the Routh–Hurwitz criterion to the characteristic Equation (19) to demonstrate the stability of the equilibrium point E 7 . The Hurwitz matrices with the coefficients of Equation (19) are determined by:
H 1 = a , H 2 = a 1 c b , H 3 = a 1 0 c b a 0 d c and H 4 = a 1 0 0 c b a 1 0 d c b 0 0 0 d .
Hence, with a , b , c , and d given in the form of Appendix A.2, the determinants of the matrices: H 1 , H 2 , H 3 , and H 4 are greater than zero. Therefore, if m r i γ i , for i = 1 , 2 , 3 , all eigenvalues of the characteristic Equation (19) have the negative real part and, consequently, the equilibrium point E 7 is locally stable. □
Proposition 2 establishes that the stability of system (1) can be studied through the equilibrium solution E 7 . This solution allows inferring that the average utilization rate ( γ i ), oxygen saturation (m), and average reproduction rate ( r i ) are the most important parameters in the variation of the system dynamics. From an ecological point of view, if m r i γ i , for i = 1 , 2 , 3 , then the coexistence equilibrium is stable, implying that regardless of the initial conditions (abundances) in which the ecosystem is found, it will be stable, as long as this relationship is satisfied, in addition to being able to obtain information on long-term water quality through AM population fluctuations.

4. Results

In this section, the influence of oxygenation/aeration rate (r) on dissolved oxygen concentration (o) and different population densities of aquatic macroinvertebrates (state variables x 1 , x 2 , and x 3 ) are studied, as well as the effects of the coexistence equilibrium point, E 7 , on water quality. Figure 2 shows the behaviors of the state variables o , x 1 , x 2 , and x 3 for different oxygenation/aeration rates. In general, when r is less than 7 % , the AM abundances decrease (red band in b–d of Figure 2). Moreover, it takes more time for these abundances to reach their maximum and be constant. If r is between 20 and 80%, the maximum is reached in less time and takes less time to be constant. It should be noted that when the oxygenation rate is greater than 3 mgL−1, there is no significant difference in the behavior of the population densities (blue band in b–d of Figure 2).
Figure 3 shows the population fluctuations of AMs under the equilibrium conditions established in Proposition 2. Moreover, Figure 3a indicates that, for all time t, the most pollution-tolerant AM families (Class III) are the most abundant, so it can be inferred that the water quality is poor. Through the behavior of the AM populations shown in Figure 3b, it can be inferred that for 0 < t < 30 , the water quality is bad, but for t > 30 , the water quality is regular because the most abundant AM families are Class II, which means that there is a transition in two different time intervals, at the same time that under those initial conditions, the ecosystem changes when the oxygenation/aeration rate increases. Likewise, Figure 3c shows that for t < 19 the water quality is poor; however, for 9 < t < 68 the water quality is fair, and for t > 68 , the water quality is good, this is inferred because the most abundant AM families are Class I (non-tolerant to pollution).

Validation

Rosas-Acevedo et al. [53] evaluated the water quality of the Coyuca Lagoon, Mexico, through environmental stress indicator insects. They made monthly visits to 10 sampling stations to collect insects between April 2011 and April 2012, collecting 2907 organisms, of which, 1343 were AM insects that were indicators of environmental stress. AMs are distributed as follows: 82 Class I organisms (Ephemeroptera, Trichoptera, Plecoptera, and Megaloptera), 1091 Class II organisms (Notonectidae, Diptera, Odonata, Hemiptera, Coleoptera, Chironomida), and 188 Class III organisms (Diptera), in other words, x 1 < x 3 < x 2 . In addition, they reported that the minimum and maximum values observed for DO were 1.5 and 4 mg/L, respectively. From the above, they concluded that the dominant organisms were Class II, followed by Class III. While for the observed DO values, Class I organisms tended to disappear. They concluded that there was moderate contamination with a tendency to increase.
The logistic population growth model was fitted to the observed data using the least squares method, assuming an asymptotic trend response to the carrying capacity of the indicator organisms. The dashed curves, shown in Figure 4, are the logistic growth curves fitted to the observed data, using Equation (20) with k ^ 1 0 , r ^ 1 = 0.1758 , k ^ 2 98 , r ^ 2 = 0.3387 , k ^ 3 17 , and r ^ 3 = 0.1585 for Classes I, II, and III, respectively.
y ( t ) = k ^ * P 0 P 0 + ( k ^ P 0 ) * exp r ^ t
Note that x 1 < x 3 < x 2 allows inferring that there is moderate contamination, with a tendency to increase. The logistic growth curves do not consider the influence of DO. Likewise, the simulated population growth curves that consider the influence of the DO (proposed in this work) are shown in Figure 4. Observe that the population fluctuations in the observed vs. simulated data have similar behaviors but at different time intervals. For example, in the observed data, the abundance of Class I AMs is higher than that of Class III for t < 5 . Whereas for t > 5 , this relationship is reversed. On the other hand, in the simulated data, these same relationships hold for t < 15 and t 15 , respectively (see Figure 4). It should be noted that the influence of oxygenation/aeration prevents the loss of diversity and abundance in short intervals of time.

5. Discussion

Oxygenation and aeration of a water body have proven to be useful for ecosystem restoration [54], for example, in Durowskie Lake, water oxygenation by wind aerators contributed to the increase in oxygen concentration, richness, and abundance of macroinvertebrates [55], which coincides with that reported by Kail et al. [56]. From Figure 2, the changes over time of the different state variables were simulated and it is evident that, the increase or decrease, the oxygenation/aeration rate, r, significantly affects the population dynamics of the organisms and dissolved oxygen concentration.
The qualitative analysis of the model determined analytically its equilibrium solutions. Proposition 1 established a set of equilibrium points belonging to the positive invariant region Ω . The existence of the coexistence equilibrium solution, E 7 , stands out for its importance in the evaluation of water quality. As a consequence and as a function of E 7 , it was possible to classify the water quality through the dominant AM families. In the same direction, the equilibrium points described in part two of the same proposition classify the stress level of an ecosystem (water quality) as more strictly associated with the loss of diversity.
The stability of the equilibrium points demonstrated in Proposition 2 made it possible to identify those parameters that bring the system to stable or unstable behaviors, since it is essential to arrive at conclusions of value for researchers or managers since stability represents the ability to return to the state of equilibrium after temporary disturbances, which are often unpredictable. From Figure 3 it is observed that, regardless of the initial conditions of the state variables, the populations of organisms tend to stabilize when the conditions necessary for the equilibrium point, E 7 , to be stable are satisfied. Moreover, this positively or negatively affects water quality when there are dominant classes of AMs.
Model validation, shown in Figure 4, evidenced the behavior of AM populations, with and without the influence of dissolved oxygen, through the rate of oxygenation/aeration. Highlighting the benefits of long-term oxygenation/aeration for the abundance and richness of bioindicator organisms. On the other hand, to strengthen the study of water quality, in subsequent works, taxa could be considered individually and take advantage of their functional responses or qualities of AMs as bioindicators, although this implies a greater number of equations and more complex analyses. Furthermore, since the conditions established here are general, they could be extrapolated to the study of different ecosystems with bioindicators with similar dynamics, adjusting the parameters to particular contexts.

6. Conclusions

Aquatic ecosystems involve many challenges for their restoration, understanding as systems, and analyses after being disturbed. The conceptualization of the model proposed here was formulated with the natural dynamics between two important environmental variables for the evaluation of water quality—bioindicator organisms and dissolved oxygen.
In this study, the changes over time of the different state variables were simulated and it is evident that the oxygenation/aeration rate, r, significantly affected the population dynamics of the organisms. In addition, knowing the current conditions of an ecosystem makes it possible to understand its current quality and, thus, establish a one-to-one correspondence between water quality and equilibrium points. On the other hand, the simulations showed that under equilibrium conditions (and regardless of initial conditions), the carrying capacity and oxygenation/aeration rate have the greatest influence over time in determining water quality.
Therefore, the study of water quality through aquatic macroinvertebrates and their close relationship with dissolved oxygen using the predictive model proposed here made it possible to show future scenarios that help guarantee water availability and its sustainable management.

Author Contributions

The authors contributed equally to this work. J.J.P.-P., J.C.H.-G., Y.E.M.-G., J.M.-R., and J.M.S. conceived, designed, and performed the numerical experiments; J.C.H.-G., J.M.-R., J.J.P.-P., Y.E.M.-G., and J.M.S. analyzed the data and wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

J.M.S. is thankful for the support of two grants from the Ministerio de Economía y Competitividad, Agencia Estatal de Investigación (AEI) and Fondo Europeo de Desarrollo Regional (FEDER) (MTM2016-78227-C2-1-P and MTM2015-69323-REDT), Spain. J.J.P.-P. is thankful for the support from PRODEP-SEP (grant no. 511-6/2019.-16017) and Internationalization of Research BUAP, México.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

J.J.P.-P. thanks the Internationalization of Research BUAP, for the support received for the development of this work. In addition, to the Vicerrectoría de Investigación y Estudios de Posgrado (VIEP) BUAP and CONACYT for the complementary support through the Sistema Nacional de Investigadores (SNI).

Conflicts of Interest

The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Mathematical Calculation

Appendix A.1. Jacobian Matrix

The Jacobian of system (1) is calculated as follows:
J = 2 r m 2 o i = 1 3 β i x i β 1 o β 2 o β 3 o m r 1 x 1 2 k 1 o 2 γ 1 x 1 r 1 1 2 m k 1 o x 1 + γ 1 ( m o ) 0 0 m r 2 x 2 2 k 1 o 2 γ 2 x 2 0 r 2 1 2 m k 2 o x 2 + γ 2 ( m o ) 0 m r 1 x 1 2 k 1 o 2 γ 1 x 1 0 0 r 3 1 2 m k 3 o x 3 + γ 3 ( m o ) .

Appendix A.2. Coefficients

The coefficients of polynomial (19) are determined by:
a = S + T 1 + T 2 + T 3 , b = T 1 T 2 + T 3 T 2 + T 1 T 3 + S T 1 + T 2 + T 3 + o * β 1 S 1 + β 2 S 2 + β 3 S 3 , c = T 1 T 2 T 3 + S T 1 T 2 + T 3 T 2 + T 1 T 3 + o * β 1 S 1 T 2 + β 1 S 1 T 3 + β 2 S 2 T 1 + β 2 S 2 T 3 + β 3 S 3 T 1 + β 3 S 3 T 2 , d = S T 1 T 2 T 3 + o * β 1 S 1 T 2 T 3 + β 2 S 2 T 1 T 3 + β 3 S 3 T 1 T 2 .

Appendix B. Parameters

Figure A1 shows AM respiration and reproduction rates when DO availability increases or decreases.
Figure A1. Impacts of hypoxia on reproduction and respiration of aquatic invertebrates. Data digitized at different oxygen concentrations, scaled to control conditions (normoxia) and the Michaelis–Menten model fit (±95% CI). Adapted from Galic et al. [23].
Figure A1. Impacts of hypoxia on reproduction and respiration of aquatic invertebrates. Data digitized at different oxygen concentrations, scaled to control conditions (normoxia) and the Michaelis–Menten model fit (±95% CI). Adapted from Galic et al. [23].
Water 14 02902 g0a1

References

  1. United Nations. Transforming our world: The 2030 agenda for sustainable development. In General Assembley 70 Session; United Nations: New York, NY, USA, 2015. [Google Scholar]
  2. Uddin, M.G.; Nash, S.; Olbert, A.I. A review of water quality index models and their use for assessing surface water quality. Ecol. Indic. 2021, 122, 3400–3411. [Google Scholar]
  3. Almeida, C.; González, S.O.; Mallea, M.; González, P. A recreational water quality index using chemical, physical and microbiological parameters. Environ. Sci. Pollut. Res. 2012, 19, 3400–3411. [Google Scholar] [CrossRef] [PubMed]
  4. Farhat, N.; Kim, L.H.; Vrouwenvelder, J.S. Online characterization of bacterial processes in drinking water systems. Npj Clean Water 2020, 3, 1–7. [Google Scholar] [CrossRef]
  5. Jørgensen, S.E.; Xu, F.L.; Salas, F.; Marques, J.C. Handbook of Ecological Indicators for Assessment of Ecosystem Health; CRC Press: Boca Raton, FL, USA, 2016; pp. 9–64. [Google Scholar]
  6. Hernández-Mira, F.A.; Rosas-Acevedo, J.L.; Reyes-Umaña, M.; Violante-González, J.; Sigarreta-Almira, J.M.; Vakhania, N. Multimetric Index to Evaluate Water Quality in Lagoons: A Biological and Geomorphological Approach. Sustainability 2021, 13, 4631. [Google Scholar] [CrossRef]
  7. Bain, M.B.; Stevenson, N.J. (Eds.) Aquatic Habitat Assessment: Commom Methods; Asian Fisheries Society: Bethesda, MD, USA, 1999; p. 186. ISBN 1-888569-18-2. [Google Scholar]
  8. Ejigu, M.T. Overview of water quality modeling. Cogent Eng. 2021, 8, 1891711. [Google Scholar] [CrossRef]
  9. Zhen-Gang, J. Hydrodynamics and Water Quality: Modeling Rivers, Lakes, and Estuaries; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  10. Loucks, D.P.; Van Beek, E. Water Resource Systems Planning and Management: An Introduction to Methods, Models, and Applications; Springer: Cham, Switzerland, 2017. [Google Scholar]
  11. Gonenc, I.E.; Wolflin, J.P. (Eds.) Coastal Lagoons: Ecosystem Processes and Modeling for Sustainable Use and Development; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  12. Jacoby, J.; Welch, E. Pollutant Effects in Freshwater: Applied Limnology, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  13. Sánchez, E.; Colmenarejo, M.F.; Vicente, J.; Rubio, A.; García, M.G.; Travieso, L.; Borja, R. Use of the water quality index and dissolved oxygen deficit as simple indicators of watersheds pollution. Ecol. Indic. 2007, 7, 315–328. [Google Scholar] [CrossRef]
  14. Kannel, P.R.; Lee, S.; Lee, Y.S.; Kanel, S.R.; Khan, S.P. Application of water quality indices and dissolved oxygen as indicators for river water classification and urban impact assessment. Environ. Monit. Assess. 2007, 132, 93–110. [Google Scholar] [CrossRef]
  15. Best, M.A.; Wither, A.W.; Coates, S. Dissolved oxygen as a physico-chemical supporting element in the Water Framework Directive. Mar. Pollut. Bull. 2007, 55, 53–64. [Google Scholar] [CrossRef]
  16. Kannel, P.R.; Kanel, S.R.; Lee, S.; Lee, Y.S.; Gan, T.Y. A review of public domain water quality models for simulating dissolved oxygen in rivers and streams. Environ. Model. Assess. 2011, 16, 183–204. [Google Scholar] [CrossRef]
  17. Kaller, M.D.; Kelso, W.E. Association of macroinvertebrate assemblages with dissolved oxygen concentration and wood surface area in selected subtropical streams of the southeastern USA. Aquat. Ecol. 2007, 41, 95–110. [Google Scholar] [CrossRef]
  18. Wilhm, J.; McClintock, N. Dissolved oxygen concentration and diversity of benthic macroinvertebrates in an artificially destratified lake. Hydrobiologia 1978, 57, 163–166. [Google Scholar] [CrossRef]
  19. Wang, J.; Fu, Z.; Qiao, H.; Liu, F. Assessment of eutrophication and water quality in the estuarine area of Lake Wuli, Lake Taihu, China. Sci. Total Environ. 2019, 650, 1392–1402. [Google Scholar] [CrossRef] [PubMed]
  20. Meier, H.E.M.; Eilola, K.; Almroth-Rosell, E.; Schimanke, S.; Kniebusch, M.; Höglund, A.; Pemberton, P.; Liu, Y.; Väli, G.; Saraiva, S. Disentangling the impact of nutrient load and climate changes on Baltic Sea hypoxia and eutrophication since 1850. Clim. Dyn. 2019, 53, 1145–1166. [Google Scholar] [CrossRef]
  21. Zhang, J.; Wang, C.; Jiang, X.; Song, Z.; Xie, Z. Effects of human-induced eutrophication on macroinvertebrate spatiotemporal dynamics in Lake Dianchi, a large shallow plateau lake in China. Environ. Sci. Pollut. Res. 2020, 27, 1–15. [Google Scholar] [CrossRef]
  22. Bazzanti, M.; Mastrantuono, L.; Solimini, A.G. Selecting macroinvertebrate taxa and metrics to assess eutrophication in different depth zones of Mediterranean lakes. Fundam. Appl. Limnol. Hydrobiol. 2012, 180, 133–143. [Google Scholar] [CrossRef]
  23. Galic, N.; Hawkins, T.; Forbes, V.E. Adverse impacts of hypoxia on aquatic invertebrates: A meta-analysis. Sci. Total Environ. 2019, 652, 736–743. [Google Scholar] [CrossRef]
  24. Etemi, F.Z.; Bytyçi, P.; Ismaili, M.; Fetoshi, O.; Ymeri, P.; Shala–Abazi, A.; Muja-Bajraktari, N.; Czikkely, M. The use of macroinvertebrate based biotic indices and diversity indices to evaluate the water quality of Lepenci river basin in Kosovo. J. Environ. Sci. Heal. Part A 2020, 55, 748–758. [Google Scholar] [CrossRef]
  25. Slimani, N.; Sánchez-Fernández, D.; Guilbert, E.; Boumaiza, M.; Guareschi, S.; Thioulouse, J. Assessing potential surrogates of macroinvertebrate diversity in North-African Mediterranean aquatic ecosystems. Ecol. Indic. 2019, 101, 324–329. [Google Scholar] [CrossRef]
  26. Silva, D.R.O.; Herlihy, A.T.; Hughes, R.M.; Macedo, D.R.; Callisto, M. Assessing the extent and relative risk of aquatic stressors on stream macroinvertebrate assemblages in the neotropical savanna. Sci. Total Environ. 2018, 633, 179–188. [Google Scholar] [CrossRef]
  27. Croijmans, L.; de Jong, J.F.; Prins, H.H. Oxygen is a better predictor of macroinvertebrate richness than temperature—A systematic review. Environ. Res. Lett. 2020, 38, 1820–1832. [Google Scholar] [CrossRef]
  28. Su, P.; Wang, X.; Lin, Q.; Peng, J.; Song, J.; Fu, J.; Wang, S.; Cheng, D.; Bai, H.; Li, Q. Variability in macroinvertebrate community structure and its response to ecological factors of the Weihe River Basin, China. Ecol. Eng. 2019, 140, 105595. [Google Scholar] [CrossRef]
  29. Mezgebu, A.; Lakew, A.; Lemma, B. Water quality assessment using benthic macroinvertebrates as bioindicators in streams and rivers around Sebeta, Ethiopia. Afr. J. Aquat. Sci. 2019, 44, 361–367. [Google Scholar] [CrossRef]
  30. Liu, Z.; Fan, B.; Huang, Y.; Yu, P.; Li, Y.; Chen, M.; Cai, M.; Lv, W.; Jiang, Q.; Zhao, Y. Assessing the ecological health of the Chongming Dongtan Nature Reserve, China, using different benthic biotic indices. Mar. Pollut. Bull. 2018, 146, 76–84. [Google Scholar]
  31. Liu, Z.; Chen, M.; Li, Y.; Huang, Y.; Fan, B.; Lv, W.; Yu, P.; Wu, D.; Zhao, Y. Different effects of reclamation methods on macrobenthos community structure in the Yangtze Estuary, China. Mar. Pollut. Bull. 2018, 127, 429–436. [Google Scholar]
  32. Pineda-Pineda, J.J.; Rosas-Acevedo, J.L.; Sigarreta, J.M.; Hernández-Gómez, J.C.; Reyes-Umaña, M. Biotic Indices to Evaluate Water Quality: BMWP. Int. J. Environ. Ecol. Fam. Urban Stud. IJEEFUS 2018, 8, 23–36. [Google Scholar]
  33. Zhou, X.D.; Xu, M.Z.; Lei, F.K.; Zhang, J.H.; Wang, Z.Y.; Luo, Y.Y. Responses of Macroinvertebrate Assemblages to Flow in the Qinghai-Tibet Plateau: Establishment and Application of a Multi-metric Habitat Suitability Model. Water Resour. Res. 2022, 58, e2021WR030909. [Google Scholar] [CrossRef]
  34. Pineda-Pineda, J.J.; Martínez-Martínez, C.T.; Méndez-Bermúdez, J.A.; Muñoz-Rojas, J.; Sigarreta, J.M. Application of Bipartite Networks to the Study of Water Quality. Sustainability 2020, 12, 5143. [Google Scholar] [CrossRef]
  35. Schleiter, I.M.; Borchardt, D.; Wagner, R.; Dapper, T.; Schmidt, K.D.; Schmidt, H.H.; Werner, H. Modelling water quality, bioindication and population dynamics in lotic ecosystems using neural networks. Ecol. Model. 1999, 120, 271–286. [Google Scholar] [CrossRef]
  36. Maier, H.R.; Dandy, G.C. The use of artificial neural networks for the prediction of water quality parameters. Water Resour. Res. 1996, 32, 1013–1022. [Google Scholar] [CrossRef]
  37. Villamarín, C.; Rieradevall, M.; Paul, M.J.; Barbour, M.T.; Prat, N. A tool to assess the ecological condition of tropical high Andean streams in Ecuador and Peru: The IMEERA index. Ecol. Indic. 2013, 29, 79–92. [Google Scholar] [CrossRef]
  38. Karaouzas, I.; Smeti, E.; Kalogianni, E.; Skoulikidis, N.T. Ecological status monitoring and assessment in Greek rivers: Do macroinvertebrate and diatom indices indicate same responses to anthropogenic pressures? Ecol. Indic. 2019, 101, 126–132. [Google Scholar] [CrossRef]
  39. El Sayed, S.M.; Hegab, M.H.; Mola, H.R.; Ahmed, N.M.; Goher, M.E. An integrated water quality assessment of Damietta and Rosetta branches (Nile River, Egypt) using chemical and biological indices. Environ. Monit. Assess. 2020, 192, 1–16. [Google Scholar] [CrossRef] [PubMed]
  40. Streeter, H.W.; Phelps, E.B. A Study of the Pollution and Natural Purification of the Ohio River; Public Health Bulletin No 146; Public Health Service: Washington, DC, USA, 1925.
  41. Wang, Q.; Li, S.; Jia, P.; Qi, C.; Ding, F. A review of surface water quality models. Sci. World J. 2013. [Google Scholar] [CrossRef]
  42. da Silva Burigato Costa, C.M.; Leite, I.R.; Almeida, A.K.; de Almeida, I.K. Choosing an appropriate water quality model—A review. Environ. Monit. Assess. 2021, 193, 1–15. [Google Scholar] [CrossRef]
  43. Park, R.A.; Clough, J.S.; Wellman, M.C. AQUATOX: Modeling environmental fate and ecological effects in aquatic ecosystems. Ecol. Model. 2008, 213, 1–15. [Google Scholar] [CrossRef]
  44. West Virginia Department of Environmental Protection. Available online: https://dep.wv.gov/WWE/getinvolved/sos/Documents/Benthic/VisualMacroGuide.pdf (accessed on 1 September 2022).
  45. Researchgate. Available online: https://www.researchgate.net/profile/Pablo-Gutierrez-Fonseca/publication/295854904_Guia_fotografica_de_familias_de_macroinvertebrados_acuaticos_de_Puerto_Rico/links/56ce23e508aeb52500c36b4f/Guia-fotografica-de-familias-de-macroinvertebrados-acuaticos-de-Puerto-Rico.pdf (accessed on 1 September 2022).
  46. Everaert, G.; De Neve, J.; Boets, P.; Dominguez-Granda, L.; Mereta, S.T.; Ambelu, A.; Hoang, T.H.; Goethals, P.L.M.; Thas, O. Comparison of the abiotic preferences of macroinvertebrates in tropical river basins. PLoS ONE 2014, 9, e108898. [Google Scholar] [CrossRef] [PubMed]
  47. Weiwei, L.; Youhui, H.; Zhiquan, L.; Yang, Y.; Bin, F.; Yunlong, Z. Application of macrobenthic diversity to estimate ecological health of artificial oyster reef in Yangtze Estuary, China. Mar. Pollut. Bull. 2016, 103, 137–143. [Google Scholar]
  48. Weiwei, L.; Zhiquan, L.; Yang, Y.; Youhui, H.; Bin, F.; Qichen, J.; Yunlong, Z. Loss and self-restoration of macrobenthic diversity in reclamation habitats of estuarine islands in Yangtze Estuary, China. Mar. Pollut. Bull. 2016, 103, 128–136. [Google Scholar]
  49. Pearl, R.; Reed, L.J. On the rate of growth of the population of the United States since 1790 and its mathematical representation. Proc. Natl. Acad. Sci. USA 1920, 6, 275. [Google Scholar] [CrossRef]
  50. Hutchinson, G.E. An Introduction to Population Ecology; Yale University Press: New Haven, CT, USA, 1978. [Google Scholar]
  51. Chapman, E.J.; Byron, C.J. The flexible application of carrying capacity in ecology. Glob. Ecol. Conserv. 2018, 13, e00365. [Google Scholar]
  52. Gunderson, L.H. Ecological resilience-in theory and application. Annu. Rev. Ecol. Syst. 2000, 31, 425–439. [Google Scholar] [CrossRef]
  53. Rosas-Acevedo, J.L.; Ávila-Pérez, H.; Sánchez-Infante, A.; Rosas-Acevedo, Y.; García-Ibañez, S.; Sampedro-Rosas, L.; Granados-Ramírez, J.G.; Juárez-López, A.L. índice BMWP, FBI y EPT para determinar la calidad del agua en la laguna de Coyuca de Benítez, Guerrero, México. Rev. Iberoam. Cienc. 2014, 1, 81–88. [Google Scholar]
  54. Rosińska, J.; Kozak, A.; Dondajewska, R.; Kowalczewska-Madura, K.; Gołdyn, R. Water quality response to sustainable restoration measures–Case study of urban Swarzędzkie Lake. Ecol. Indic. 2018, 84, 437–449. [Google Scholar] [CrossRef]
  55. Gołdyn, R.; Podsiadłowski, S.; Dondajewska, R.; Kozak, A. The sustainable restoration of lakes—Towards the challenges of the Water Framework Directive. Ecohydrol. Hydrobiol. 2014, 14, 68–74. [Google Scholar] [CrossRef]
  56. Kail, J.; Brabec, K.; Poppe, M.; Januschke, K. The effect of river restoration on fish, macroinvertebrates and aquatic macrophytes: A meta-analysis. Ecol. Indic. 2015, 58, 311–321. [Google Scholar] [CrossRef]
Figure 1. Conceptual diagram of the mathematical model. The state variables o, x 1 , x 2 , and x 3 represent the DO concentrations in the water and the AM populations—intolerant, regularly tolerant, and tolerant to organic contamination, respectively. The processes are: (a) oxygenation/aeration, (b) DO consumption, (c) AM natality, (d) DO utilization, and (e) interspecific competition.
Figure 1. Conceptual diagram of the mathematical model. The state variables o, x 1 , x 2 , and x 3 represent the DO concentrations in the water and the AM populations—intolerant, regularly tolerant, and tolerant to organic contamination, respectively. The processes are: (a) oxygenation/aeration, (b) DO consumption, (c) AM natality, (d) DO utilization, and (e) interspecific competition.
Water 14 02902 g001
Figure 2. Influence of oxygenation/aeration rate r on AM population densities. (a) Variation of DO through time for different values of r. The population fluctuations of the different classes of AMs are shown in (b), (c), and (d) for x 1 , x 2 , and x 3 , respectively. Note that, in all cases, if r is smaller, then AM reproduction is smaller and slower than when r is larger.
Figure 2. Influence of oxygenation/aeration rate r on AM population densities. (a) Variation of DO through time for different values of r. The population fluctuations of the different classes of AMs are shown in (b), (c), and (d) for x 1 , x 2 , and x 3 , respectively. Note that, in all cases, if r is smaller, then AM reproduction is smaller and slower than when r is larger.
Water 14 02902 g002
Figure 3. Behavior of AM population densities. In all cases; m = 14.5 , r 1 = 0.1 , r 2 = 0.2 , r 3 = 0.3 , β 1 = 1 * 10 3 , β 2 = 1 * 10 4 , β 3 = 1 * 10 5 , γ 1 = 1 * 10 4 , γ 2 = 1 * 10 5 , γ 3 = 1 * 10 6 , and the initial conditions are the same, except (a) k 1 < k 2 < k 3 and r = 0.2 (poor water quality), (b) k 1 < k 3 < k 2 and r = 2.5 (regular water quality), and (c) k 3 < k 2 < k 1 and r = 10.5 (good water quality).
Figure 3. Behavior of AM population densities. In all cases; m = 14.5 , r 1 = 0.1 , r 2 = 0.2 , r 3 = 0.3 , β 1 = 1 * 10 3 , β 2 = 1 * 10 4 , β 3 = 1 * 10 5 , γ 1 = 1 * 10 4 , γ 2 = 1 * 10 5 , γ 3 = 1 * 10 6 , and the initial conditions are the same, except (a) k 1 < k 2 < k 3 and r = 0.2 (poor water quality), (b) k 1 < k 3 < k 2 and r = 2.5 (regular water quality), and (c) k 3 < k 2 < k 1 and r = 10.5 (good water quality).
Water 14 02902 g003
Figure 4. Population growth curves of observed vs. simulated data. The dotted lines do not consider the influence of the DO, while the continuous line curves do. Data from [53].
Figure 4. Population growth curves of observed vs. simulated data. The dotted lines do not consider the influence of the DO, while the continuous line curves do. Data from [53].
Water 14 02902 g004
Table 1. Description and parameter values of system (1).
Table 1. Description and parameter values of system (1).
ParameterDescriptionValueReference
roxygenation/aeration rate0–14.5 mgL−1Galic et al. [23]
moxygen saturation constant a0 < m ≤ 14.5 mgL−1
βiaverage respiration rate of AMs a0–0.918Galic et al. [23]
riaverage reproduction rate of AMs a0–1.02Galic et al. [23]
kiaverage carrying capacity of AMsvaries
γiaverage DO utilization ratevaries
a—Assuming an oxygen saturation response for each of the rates (see Figure A1 in Appendix B).
Table 2. Relationships between the cardinalities of different AM classes, meaning, and classification of water quality.
Table 2. Relationships between the cardinalities of different AM classes, meaning, and classification of water quality.
RelationMeaningWater Quality
x3 < x2 < x1Low pollution with a tendency to increasegood
x3 < x1 < x2Regular pollution with a tendency to decreasemoderate
x2 < x3 < x1Low pollution with a tendency to increasegood
x2 < x1 < x3High pollution with a tendency to decreasepoor
x1 < x3 < x2Regular pollution with a tendency to increasemoderate
x1 < x2 < x3High pollution with a tendency to decreasepoor
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pineda-Pineda, J.J.; Muñoz-Rojas, J.; Morales-García, Y.E.; Hernández-Gómez, J.C.; Sigarreta, J.M. Biomathematical Model for Water Quality Assessment: Macroinvertebrate Population Dynamics and Dissolved Oxygen. Water 2022, 14, 2902. https://doi.org/10.3390/w14182902

AMA Style

Pineda-Pineda JJ, Muñoz-Rojas J, Morales-García YE, Hernández-Gómez JC, Sigarreta JM. Biomathematical Model for Water Quality Assessment: Macroinvertebrate Population Dynamics and Dissolved Oxygen. Water. 2022; 14(18):2902. https://doi.org/10.3390/w14182902

Chicago/Turabian Style

Pineda-Pineda, Jair J., Jesús Muñoz-Rojas, Y. Elizabeth Morales-García, Juan C. Hernández-Gómez, and José M. Sigarreta. 2022. "Biomathematical Model for Water Quality Assessment: Macroinvertebrate Population Dynamics and Dissolved Oxygen" Water 14, no. 18: 2902. https://doi.org/10.3390/w14182902

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop