Next Article in Journal
Understanding the Spatial and Temporal Distribution and Environmental Characteristics of Polychaete Assemblages in the Coastal Waters of Ulleungdo, East Sea of Korea
Next Article in Special Issue
Making of Massoia Lactone-Loaded and Food-Grade Nanoemulsions and Their Bioactivities against a Pathogenic Yeast
Previous Article in Journal
Sustainable Development of Oil and Gas Resources: A System of Environmental, Socio-Economic, and Innovation Indicators
Previous Article in Special Issue
Vertebrate Palaeoecology of the Pisco Formation (Miocene, Peru): Glimpses into the Ancient Humboldt Current Ecosystem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pumping Rate and Size of Demosponges—Towards an Understanding Using Modeling

by
Poul S. Larsen
1,* and
Hans Ulrik Riisgård
2
1
Department of Mechanical Engineering, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
2
Marine Biological Research Centre, University of Southern Denmark, 5300 Kerteminde, Denmark
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(11), 1308; https://doi.org/10.3390/jmse9111308
Submission received: 15 October 2021 / Revised: 11 November 2021 / Accepted: 16 November 2021 / Published: 22 November 2021
(This article belongs to the Special Issue Feature Papers in Marine Biology)

Abstract

:
Filter-feeding sponges pump large amounts of water and contribute significantly to grazing impact, matter transport and nutrient cycling in many marine benthic communities. For ecological studies it is therefore of interest to be able to estimate the pumping rate of different species from their volume size or osculum cross-sectional area by means of experimentally determined allometric correlations. To help understand allometric data correlations and observed large variations of volume-specific pumping rate among species we developed a model that determines the pumping rate as a function of the size (volume) of a tubular-type demosponge described by 4 geometric length scales. The model relies on a choanocyte-pump model and standard pressure loss relations for flow through the aquiferous system, and density and pumping rate per choanocyte is assumed to be constant. By selecting different possibilities for increase of the length scales, which may also simulate different growth forms, we demonstrate that the model can imitate the experimental allometric correlations. It is concluded that the observed dependence of pumping rate on size is primarily governed by the hydraulics of pump performance and pressure losses of the aquiferous system rather than, e.g., decreasing density of choanocytes with increasing sponge size.

1. Introduction

Filter-feeding sponges are one of the major components of marine benthic communities with a worldwide distribution [1]. Data on sponge size and pumping rates (filtration rates) are important for ecological studies in determining the grazing impact, matter transport and nutrient cycling of observed populations of sponges [2,3,4,5,6]. While sponge volume (V) may be more difficult to measure or estimate in situ than the osculum cross-sectional area (OSA), the latter has also been used as a correlation parameter [7], who used the following allometric relationships to characterize sponges:
Q/V = aVb; Q = a1OSAb1; Q = a2Vb2; U = a3Vb3; OSA = a4Vb4
where Q/V denotes the volume-specific pumping rate, Q the pumping rate, V the sponge volume, OSA the osculum cross-sectional area, and U the exhalant jet speed at osculum.
Demosponges are modular filter feeders that consist of one or several units, known as “aquiferous modules” [8,9,10]. An aquiferous module is a functional unit that draws ambient water through numerous inhalant openings (ostia) into an incurrent canal system by means of pumping units (choanocyte chambers, CC) that filter the water for nutrition and then eject it via an excurrent canal system to an atrial cavity and through a single exhalant opening (osculum) as a jet to the surrounding water.
Temporal variation in OSA and pumping activity is a common feature among sponges [11,12,13]. Laboratory observation of the pumping activity of a multi-oscula sponge colony indicated that the sponge is sensitive to changes in water motion resulting in contractions of oscula, and further that a constriction of an individual osculum may cause closure of adjacent oscula and thus become a behavioral response of the entire modular sponge colony [14]. More recently, [15] compared relations between sponge size, OSA, exhalant jet speed, and pumping rate (Q = OSA × U) of single-osculum explants and of individual aquiferous modules of multi-oscula explants of the demosponge Halichondria panicea. The latter modules were identified by observing from which osculum surface administered dye emerged, and there was a fair agreement in results between the two types of modules. For both types, the pumping rate was a linear function of the modules’ size (volume), and it was suggested [15] that this finding along with data from the literature on five demosponges [16] indicated that the density of choanocytes may be of the same order, independent of sponge size, although it is still unknown whether the choanocyte chamber density and its pumping rate vary through the ontology of long-lived species, especially for those that reach considerable size. Suggesting a constant density of pumping units (choanocytes) throughout any given sponge size, and assuming simple geometric (isometric) scaling V~OSA3/2, [17] suggested a theoretical scaling of Q~OSA3/2, which has been experimentally supported and confirmed in other studies, see Table 4 in [15]. However, [6,7] found in recent in situ studies on 20 demosponges that Q on average increased allometrically Q = a1OSAb1 with scaling exponent b1 ≈ 0.75, 1.07 and 0.67 (instead of the theoretical value of b1 = 3/2), thus suggesting that the choanocyte density may not be constant but on the contrary decreases with sponge size.
A consequence of possible constant density and pumping rate of choanocytes is that Q as a function of size V is linear, i.e., Q = aV, and further that the volume-specific pumping rate is constant Q/V = aVb, b = 0, regardless of sponge size. This is supported by [13] and see Figure 5 in [15], but not by Table 5 in [6] that give negative b values ranging from −0.2 to −0.7. Further, [7] recently reported b = −0.52 for seven Red Sea sponge species and b = −0.43 for eight Caribbean species and therefore suggested that the theoretical considerations that explained the relationships observed by [15,17] for relatively small sponge explants may not be “expected to hold for much larger and more static sponges”. Thus, the scaling b < 1 of power-law exponents may, according to [7], be explained by different aquiferous system structures and by choanocyte chamber densities, which may vary among species and with sponge size, and that the distribution of choanocyte chambers may vary in different parts of the sponge body. Further, according to [7], “as the sponge grows, the aquiferous system may become more complex” and thus a “constraint on the growth”. Obviously, the relationships between sponge dimensions and pumping rate are not simple and vary among species. Specifically, the increase in pressure losses with increasing size in the aquiferous system structure may account for the observed scaling b < 1 as discussed in the present study.
Here, as a supplement to the recent detailed summary of allometric data correlations by Table 4 in [7], we present available experimental literature data on volume-specific pumping rate (Q/V = aVb) of demosponges characterized by values of the exponent b. Furthermore, we present possible modelling of this relationship. Previous studies [16,18] have determined the pressure losses for a specified pumping rate while the present approach, also using a model for the choanocyte pumps, gives the pumping rate. The way in which morphological length scales change with size increase affects pumping rate through changing pressure losses and may be an important factor and thus some help to explain observed trends. Thus, a key issue of the present study has been to establish trends in Q/V versus V based on use of pump-system characteristics.

2. Available Data

Over the years, there has been a considerable interest in measuring the pumping rate (Q) of sponges as a function of their size (volume V). Results are often expressed as volume-specific pumping rate (Q/V, in units of mL water pumped per min per mL sponge = mL min−1 mL−1= min−1) versus volume (V, in units of cm3), see Table 1 and Figure 1 for representative data. Reiswig Figure 1 in [2] found Q/V to be essentially constant for increasing size of three specimens, likely because of constant density of choanocytes delivering the same volumetric flow rate, but specimens also showed a reduction of activity in largest sized individuals. For Halichondria panicea and Haliclona urceolus, Tables 1 and 2 in [19] measured a declining Q/V with increasing size characterized by the negative power-law exponents of b = −0.6 and −0.39 of the correlation Q/V = aVb, similar to b = −0.35 for large Xestospongia muta by [20]. The average of data from 14 species [21] and five species [16], on the other hand, showed Q/V to be nearly constant (b = −0.045 and 0.051, respectively). Similar nearly constant Q/V trends are seen in the data [17] for singe-osculum explants of H. panicea (b = −0.038, b1 = 1.45 and b2 = 0.974), and in the data [22] for Cinachyrella cf. cavernose (b = −0.071, b2 = 0.929) where b2-values (Equation (1)) near unity suggest constant pumping rate and density of choanocytes. Smaller single- and multi-oscular explants [15] showed different behaviors (b = −0.19 and 0.726, respectively, and b1 = 1.22), while the data of five Mediterranean species [6] showed consistent negative power-law exponents (b = −0.19 to −0.49) over large ranges of increasing size and b1 = 0.75. For 15 tropical species Figures S7 and 3 in [7], average values are suggested to be b = −0.52 and −0.43, and b1 = 1.07 and 0.67, respectively. While it appears that other investigators based their data on total sponge volume, [13] explicitly used sponge tissue volume in arriving at an average volume-specific pumping rate of 0.06 ± 0.04 l (s l tissue)−1 = 3.6 ± 2.4 min−1 for 274 large X. muta specimens, and the linear increase of pumping rate with size (b2 = 1.1) suggests b~0.
Data on other allometric scaling correlations in Equation (1), as well as growth forms of species, are given by Table 4 in [7].
From the available data, it appears that the most fundamental feature of scaling laws of pumping rate is associated with the performance of the choanocyte pumps driving the flow against the pressure losses determined by the morphological structure of the aquiferous system of canals and restrictions, which are determined by several characteristic length scales [16]. Data on choanocyte density (Figure 2) show no clear trend for increasing sponge size, but all appear to be of the order of 106 mm−3 within a factor of about 2. Therefore, the modelling in next section provisionally assumes density of choanocytes to be constant.

3. Modelling

The purpose of the present modelling is to determine how pumping rate changes when pressure losses in the aquiferous system changes for increasing size of a demosponge due to various changes of morphological length scales (diameter, height, etc.) on the assumption of constant choanocyte density and to compare the results to available data. For a given configuration, the total pressure loss versus flow rate defines the so-called system characteristic. Furthermore, the pumping rate depends on how the basic pump unit, the choanocyte, or the group of these in a choanocyte chamber, responds to a given pressure to be delivered, which is given by the so-called pump characteristic [26]. Intersection of the system and pump characteristics determines the actual pumping rate and pump pressure.
Demosponges exhibit, according to [24], two extremes of external morphology: A low, solid, or encrusting form and an upright, hollow, or tubular form, and “all other types of habitus can be viewed as either slight modifications of these two dominant patterns or intermediates between these extremes.” In the present study, we focus on the tubular sponge form to illustrate basic features and hydraulic consequences of increasing various morphological length scales. However, by suitable choice of model length scales, a tubular form may degenerate to a low and wide form with long incurrent and excurrent canals and thus imitate a massive or encrusting form.
For the model of a tubular sponge (Figure 3), whose morphology and size are given by four geometric length scales, model calculations give pumping rate as a function of sponge volume. Size increase in a sponge is determined by a genetically controlled ground plan (Bauplan) of morphological features common to each sponge species for how the individual length scales increase. One feature of this plan, for example, is revealed in the recorded data on the correlation Table 4 in [7] OSA = a4Vb4, which relates geometric dimensions. Later we shall use this as a constraint in the study of model performance. For each ground plan, the model calculations will show a trend of the volume-specific pumping rate versus size. Among possible scaling trends, we seek to identify the most likely ones to match a given observed trend.
Recognizing the structural differences in the geometrical form of different sponges, we consider three basic variants defined by the same four length scales, H, DO, Da, Dosc: (i) The slender tubular form (Figure 3A) where increasing size is due to increase primarily in height and less in diameter; (ii) the massive or encrusted form Figure S1H–L in [7], represented by a cylindrical form where increasing size is due to increase primarily in diameter and less in height, e.g., due to new aquiferous modules growing next to existing ones in multi-oscula sponges; (iii) the small cone-shaped single-osculum explants Figure 1 in [17] where increasing size is probably due to increase of both height and diameter. The individual modules of explants Figure 1 in [15] may be treated as the slender tubular form. For the tubular cases, volume and canal length are defined in Figure 3 as V = (π/4) DO2 H and L = (DODa)/2, respectively. For the cone-shaped case, V = (π/12) DO2 H and the average canal length is approximated by L = (DODa)/4 and the outer surface area by As = π (H2 + DO2/4)1/2 DO/2.

3.1. Materials and Methods

The model calculations are based on use of the computer modeled pump characteristic for a leucon-type demosponge Equation (2.13) for C1 = 0 in [27] which leads to the nonlinear characteristic equation,
Δppmax = (1 − Q/Qmax)1/3
where Q denotes the pumping rate for a given delivered pressure rise Δp against and equal to the sum of pressure losses of the system represented by the flow through the sponge. The constants Δpmax and Qmax were determined from the data for morphological dimensions and pumping rate of the reference sponge Haliclona urceolus [19] as explained in Appendix A where Table A1 and Table A2 summarize sponge parameters.
The pressure losses are calculated according to a previous model [19] now accounting also for the pressure loss of the atrium flow. In the model, choanocyte chambers (CCs) are imbedded in wall structures between incurrent (IC) and excurrent (EC) canals such that these pump units act in parallel (Figure 3C). In the annular tissue structure of the tubular sponge, ambient water is drawn through ostia restriction into an IC canal by a row of CCs that discharge the flow into an EC canal that empties into the atrial cavity where the upward volume flow Qa cumulates to leave ultimately through the osculum as an exhalant jet. There are several inlets (prosopyles) to a CC but only one outlet (apopyle). As part of the model, we assume that the observed decrease of diameter of ICs with length and with branching of ICs implies an approximately constant pressure gradient along the canal (see Appendix A.2 for a discussion). If CCs of the same pumping rate draw flow from points uniformly distributed over an IC inner surface (dQ/dx~−D) and the pressure gradient due to friction is constant (Q~D4), the diameter of the IC would decrease as D/D0 = (1 − x/L)1/3, where L is the total length. Note also that a constant pressure gradient along IC and EC would imply that each one of the CCs in a row (Figure 3C) would face the same length of flow path, hence same total pressure drop and thus provide the same volume flow. The equations of the various pressure losses are detailed in Appendix A.

3.2. Results

First, let us ignore the effect that pressure losses have on reducing the pumping rate. In the search of scaling relations for volume-specific pumping rate versus volume (size) of a tubular sponge (Figure 3) we note that V~DO2H and Q~(DO2Da2)H provided the pumping rate and density of choanocytes were constant in the structural annular volume. This scaling leads to the expression Q/V~1 − (Da/DO)2, which for increasing V implies a decrease, an increase, or no change depending on how Da/DO changes with increasing V, irrespective of how H changes. However, including the effect of flow reduction due to pressure losses will change the outcome, in most cases leading to decreasing pumping rate with increasing volume, as illustrated in the following examples. First, however, a simple matching of a linear pump characteristic with the pressure loss in a canal [28] leads to the relation Q/V~1/(c + V2/3) which gives Q/V~V−2/3 for c << V and approaches Q/V~constant for large values of the constant c (see Appendix A.4 for the analysis). These results correspond to the range of trends shown by much of the data reported in Table 1.
Using the flow model described in Appendix A we first study isometric scaling of the geometry given by the four length scales H, DO, Da, Dosc defined in Figure 3. Introducing a scale factor sf that takes values from 0.25 to 4, we calculate how the volume-specific pumping rate Q/V changes with the relative sponge volume (V/Vref) as the four length scales are changed by the factor sf (Figure 4 upper curve). The volume given in Figure 3 will increase as sf3, but the pumping rate will increase less (even though the choanocyte density is constant) because increased length of various canals increases pressure losses and thus reduce flow rate as dictated by the pump characteristic. Here, this implies a negative power-law exponent b = −0.29 which is similar to values displayed by some data in Table 1. Next, reducing the scaling of the osculum diameter to Dosc~sf1/2 changes the power-law exponent to b = −0.57 and b1 = 1.3 (Figure 4 lower curve). These values should be compared to b = −0.6 and −0.39 for Halichondria panicea and Haliclona urceolis, respectively, of Table 1 (#4 and #5) and b1 = 1.45 of the correlations of Q = a1 OSAb1 (not shown).
To see the effect of keeping some length scales constant, Figure 5A,B shows in detail how constant values of height or diameter, respectively, affects pumping rate for increasing sponge size when the remaining parameters are scaled. The figures show the existence of regions of increase, no change, and decrease in volume-specific pumping rate, corresponding to positive, zero, and negative values of the power-law exponent b. Increase occurs, for example, when diameters of atrium and osculum of a smaller sponge increase despite increasing size. Decrease occurs when length of canals and size increase.
To model single-osculum explants, we choose as reference ID#21 of Table 1 in [17] for which Q/V = 1.77 min−1, b = −0.38 (see Table 1), and b1 = 1.45. For this reference, values of the four length scales are H =13 mm, DO = 15.6 mm, Da = 1.1 mm, and Dosc = 15.6 mm. Model results in Figure 6 suggest that isometric scaling, giving b = −0.04 and b1 = 1.4, is in fair agreement with measured values. Single-osculum explants #9a of Table 1 in [15] are of about the same size but with lower Q/V and show somewhat stronger decrease with increasing size, b = −0.19 and b1 = 1.22.
The individual modules of multi-oscula explants of Halichondria panicea are tubular Figure 1 in [15] rather than cone-shaped, and the data show a strong increase in pumping rate versus size (#9b of Table 1, b = 0.726) that is difficult to interpret. Size increase is probably associated with increase in both height and outer diameter, and if canal length (L = (DODa)/2) does not increase appreciably, pumping rate would tend to increase with size.
For demosponges of massive or encrusting form with several oscula, hence composed of several aquiferous modules, such as many of the species #10–14 of Table 1, the structure and canal system is not well known. During growth, we imagine a first module can be treated as tubular until it reaches its maximal size and a new module starts to grow adjacent to the first one. Ref. [10] states that “Among encrusting demosponges new modules arise at the periphery and often result in the formation of very large sheets bearing more than 50 oscula”. As new modules are added, the free surface area available to each module for inflow through ostia and incurrent canals is steadily being reduced which may imply that canals become longer which in turn causes a relative reduction in pumping rate of all modules such that b < 0. Such volume increase during growth may be approximated by the tubular model in a phenomenological way by a scaling with little increase in height, large increase of outer diameter and canal lengths, and with an increase in osculum diameter, that reflects increasing number of oscula and pumping rate.
Finally, it is of interest to see how the model performs when using some experimentally determined scaling data relating to morphology. Here we consider the correlation OSA = a4Vb4 with a4 = 9.42 and b4 = 0.6 from Figure 5A in [7] for low-microbial-abundance, LMA, as constraint on the model parameters. Using the OSA-relation and the relations V = (π/4)DO2H and Dosc = (4OSA/π)1/2, determines Dosc as a function of DO and H, leaving the three remaining length scales H, DO, Da free to be scaled. As shown in Appendix A Figure A1A, sample results for three scaling cases (scaling H, or diameters DO, Da) or all three parameters give resulting values of exponents b, b1, and b2. Values of b (=−0.16, −0.2 and −0.3) are comparable to those of Table 5 in [6] for LMA species, also listed as #10–11 in Table 1 of the current paper. Similar calculations (Appendix A, Figure A1B) for high microbial-abundance HMA species (a4 = 3.81 and b4 = 0.56) give b = −0.2, −0.26 and −0.37. These results suggest that some or all observed trends of reduced Q/V with increased V could be explained almost entirely by the hydrodynamic effects of increasing pressure losses in aquiferous structures of increasing size. For more detail, Figure A2A–C shows the different contributions to the total pressure loss for the three cases of Figure A1A. As expected, according the model, as canal length increases for scaled diameters at constant height, pressure losses in incurrent and excurrent canals is the dominant contribution (by 50–85% of the total) which could represent encrusted and massive growth forms, but also for isometric scaling of height and diameter (by 20–75%). For scaled height at constant diameter, on the other hand, pressure losses in the atrium flow dominate (by 26–75%), which could represent the slender tubular growth form.

4. Discussion

The available literature data on pumping rate versus size of demosponges covers six decades of volumes from small single- and multi-oscula explants to large temperate and tropical sponges of various forms (Table 1). The morphological scaling is presented in terms of power-law regressions of volume-specific pumping rate versus volume (Q/V = aVb), which shows a variety of trends of positive, near-zero and—most often—negative exponents b (Figure 1). Alternative scaling involving the parameters pumping rate Q, sponge volume V, osculum cross-sectional area OSA, and exhalant jet speed at osculum U (Equation (1)) are summarized in [7].
Available data show that scaling laws of pumping rate appear to depend on species, temperate or tropical habitat, high or low microbial abundance (HMA or LMA) and growth form (tubular, vase, encrusting, massive). Possible other factors include density of choanocytes, seasonal variation in food availability and thus variation of the condition index (ratio of organic to inorganic matter [29], and height of sponge modules above the sea floor, as well as the morphological structure. The latter aspect is the focus of the present study because it is obvious that, for example, increased canal lengths with increased sponge size leads to increased pressure losses, hence to reduced pumping rate as dictated by the characteristic of the flagellum pump of choanocytes. To illustrate such effects, we developed a computation model that calculates the pumping rate as function of sponge size as specified by four length scales for a tubular model (Figure 3).
The model assumes a constant density of choanocytes because there is no clear trend in available data for its possible size dependence (Figure 2), but it could readily be modified to include any specified variation of density with size. The model (see Appendix A for details) also accounts for the pressure losses of the upward flow in the atrial cavity, which implies a decreasing pressure with increasing height facing the choanocyte pumps. However, to facilitate an algebraic solution we approximate the pressure in the atrium by its mean value such that all pumps deliver the same pressure rise and same volume flow. It is estimated that the error of this approximation amounts to a few percent for short atrial canals.
For realistic results, the model is based on the data for the “standard sponge” Haliclona urceolus [19] which is used as the reference sponge and for which the equation of pump characteristic was developed [27]. Reference values appearing in results are Vref = 1.86 cm3, (Q/V)ref = 3.41 min−1 and Qref = 6.34 mL min−1, and an isometric scaling with sf = 0.25 to 4 will cover a range of sponge sizes from 0.029 to 119 cm3. The model results in Figure 4 for isometric scaling (H, DO, Da, Dosc~sf) and a modified scaling (H, DO, Da, ~sf, Dosc~sf1/2, i.e., OSA~sf) are compared to the scarce experimental results for Halichondria panicea and H. urceolus (b = −0.6 and −0.39, respectively, of #4 and #5 of Table 1). The case of modified scaling implies a reduced OSA, hence increased exit pressure loss and thus reduced pumping rates that seem to fit the data better, hence supporting the credibility of the model.
To help understand the response of model prediction to changes of the four length scales, Figure 5A shows results for four values of constant height (H) as diameters are scaled (DO, Da, Dosc~sf), and Figure 5B shows results for four values of constant outer diameter (DO) as other parameters are scaled (H, Da, Dosc~sf). The results show regions in which scaled increase of sponge volume may involve increase (b > 0), essentially no change (b~0), or decrease (b < 0) of volume-specific pumping rate. First, for some constant value of height H (Figure 5A), scaled increase of all diameters implies V~DO2H and Q~(DO2Da2)H, hence Q/V~1 − (Da/DO)2~constant if no pressure losses. This is consistent with an assumption of constant pumping rate and density of choanocyte pumps. However, for small diameters and large heights, such pressure losses are large primarily due to the flow along the atrial cavity and through osculum and these contributions decrease with increasing diameters such that Q/V increases at first. Later, for larger diameters, the pressure loss in canals of increasing length, L = (DODa)/2, becomes important and Q/V levels off and begins to decrease, particularly for the smaller values of height. Such situations may represent encrusting or massive growth forms. Secondly, for some constant outer diameter DO and constant width of the annular structure of a tubular sponge, which equals the canal length L = (DODa)/2 (Figure 5B), scaled increase of height H implies again Q/V~constant, if there are no pressure losses. However, for small height, the pressure loss is primarily due to flow through the small osculum and this contribution decreases with increasing height of sponge and diameter of the osculum such that Q/V increases at first. Later, for increasing height the increasing pressure loss in the atrium flow causes Q/V to level off and begin to decrease, particularly at the smaller values of outer diameter. Such situations may represent slender tubular growth forms.
To study single-osculum explants, the tubular model was changed such that the four length scales could represent a cone-shaped geometry. The results in Figure 6 for three strategies show that scaled height at constant diameters leads to decreasing pumping rate during growth (b < 0), while scaled diameters at constant height gives increasing pumping rate (b > 0). Isometric scaling of height and diameters gives a weakly decreasing, almost constant pumping rate (b = −0.041) that is close to the observation [17] (b = −0.038). The correlation of the data on OSA versus volume [17] shows Dosc~V0.34 which is in good agreement with the isometric scaling, which implies Dosc~fs and V~fs3, hence Dosc~V1/3. For the correlation Q = a2Vb2, the experiment and the model give almost identical trends, b2 = 0.97 and 0.96, respectively, supporting the hypothesis of constant choanocyte density. However, for magnitudes we should compare the complete correlations, which for the pumping rate of the experiment and the model are Qexp = 2.3V0.97 and Qmodel = 4.5V0.96, respectively, and for the volume-specific pumping rate, (Q/V)exp = 2.3V−0.038 and (Q/V)model = 1.79V−0.041. The discrepancies in magnitude could be due to underestimation of the model pressure losses, overestimation of the choanocyte pump performance, and the power-law fit to model results that show significant variation over the scaling range. Nevertheless, results strongly suggest that observed dependence of pumping rate on size of demosponges might primarily be governed by the hydraulics of pump and pressure losses of the aquiferous system.
We also explored the effect of constraining the scaling strategy by the geometrical correlation OSA = a4Vb4 obtained from experiments, leaving the three length scales H, DO, Da free to be scaled (see Figure A1 in Appendix A). For both an LMA-type (e.g., Dysidea avara) and an HMA-type (e.g., Agelas oriodes), the best agreement of trend in terms of b-values appeared to be for scaled height at constant diameters, which would be expected for a tubular growth form. The two specimens are classified as encrusting and massive, but images in the literature (e.g. Figure S1 in [7]) show specimens and colonies with tubular growth form.

5. Conclusions

The present model determines the pumping rate as a function of size of a tubular demosponge described by four geometric length scales. The model relies on the availability of a choanocyte pump model and standard pressure loss relations for flow through restrictions and a canal system. Each choice of how the four length scales increase leads to a correlation between volume-specific pumping rate and size in terms of volume of a sponge which is compared to the published data on such allometric correlations that are summarized in Table 1. The trends of data, e.g., decrease or increase of pumping rate with increasing size (Figure 1), vary a great deal among species, their temperate or tropical habitat, high or low microbial abundance (HMA or LMA), and possibly other factors such as density of choanocytes and growth form, which affects the morphological structure. By selecting different possibilities for increase of the four length scales it has been demonstrated that the model imitates the experimental allometric correlations qualitatively (through the sign and magnitude of the exponent b) and to some extent also quantitatively (through the magnitude of the coefficient a). It is concluded that observed dependence of pumping rate on size of sponges might primarily be governed by the hydraulics of pump and pressure losses of the aquiferous system and not by, e.g., a reducing density of choanocytes with increasing size.

Author Contributions

P.S.L. developed the model and did the data analysis, H.U.R. contributed with biological input and text writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Pump Characteristic of a Demosponge

To quantify the pumping performance of a sponge we need a model for the basic pump unit, the choanocyte or the choanocyte chamber, or for the whole sponge, in terms of the pump characteristic Δpp{Q}, and the relation between volume flow Q and pressure losses of the flow through the sponge, in terms of the system characteristic Δps{Q}. Equating Δpp to Δps, determines the operating point Δpop, Qop.
For the pump, we use the results of the computational modeling of the leaky leucon pump Equation (2.13) for C1 = 0 [27] that lead to the nonlinear characteristic equation,
Δppmax = (1 − Q/Qmax)1/3
where Δpmax and Qmax were determined from the data of the demosponge Haliclona urceolus [19], Table 1 and Figure 1c in [18]. Here, the figure gives the measured back-pressure characteristic while the table gives the following parameters: size of specimen, length L = 32 mm; outer and atrium diameters DO = 8.6 mm, Da = 2.3 mm, respectively, hence volume Vsp = 1726 mm3; pumping rate and pressure drop at the operating point Qop = 6 mL min−1, Δpop = 0.673 mm H2O (=6.6 Pa), respectively; estimated CC density n = 12,000 mm−3; estimated number of choanocytes per CC, nch = 80. The measured value for the sponge, Δpmax = 2.69 mm H2O (=26.4 Pa) ([18], Table 1 and Figure 1c therein) applies equally to the choanocyte chamber and the individual choanocyte when they all act in parallel. Because Equation (A1) implies a vertical tangent at Qmax we approximately set Qmax,sp ≈ 6 mL min−1 for the sponge. For Equation (A1) to apply to a choanocyte chamber, we calculate the corresponding pumping rate per CC to Qmax = 6 × 1012/(60 × 1726 × 12000) = 4828 µm3 s−1, and per choanocyte Qmax,ch = 60.3 µm3 s−1.

Appendix A.2. Simplified Model of Flow System

Despite several studies of the dimensions of incurrent (IC) and excurrent (EC) canals and maps of morphometric data [16,24,25], the details of actual branching aquiferous canal systems are complex and not fully understood. The branching, say from diameter D1 to D2 and D3 at a bifurcation, should follow Murray’s law Equation (14.6) in [30] to ensure constant wall shear stress, D1e = D2e + D3e with e = 3. On the other hand, assuming the pressure gradient of Poiseuille flow to be constant (Q~D4 and Q1 = Q2 + Q3) leads to a diameter exponent of e = 4. The detailed study by Figure 4 in [31] found that the exponent e varied considerably, with a spectral distribution that peaked between 1 and 5. For the present flow model, we assume e = 4 such that the pressure gradient remains constant in the IC canal system that branches and thins because of outflow to choanocyte chambers as well as in EC canals that thicken due to inflow from chambers (Figure 3C). In this way it is likely that the total pressure loss through IC and EC canals experienced by flow through a given choanocyte chamber would be the same, hence the parallel coupled CC-pumps would all operate at the same pumping rate.

Appendix A.3. Model Equations

The present flow model for the sponge is based on a model [19] as derived for a specimen of the tubular demosponge Haliclona urceolus of dimensions: outer diameter DO, inner atrium diameter Da, height H and osculum diameter Dosc (Figure 3B,C). The volume of sponge is Vsp = (π/4)DO2H, density of choanocyte chambers (CCs) is denoted n and total pumping rate Q. Central to the flow model is the size and number of incurrent canals (IC) and the number of choanocyte chambers (CC) associated with each IC. An IC of mean diameter D is surrounded by a wall structure of thickness Dcc in which choanocyte chambers (CCs) of diameter Dcc are imbedded at a spacing l~Dcc in a rectangular array over the inner IC surface area. For an IC canal of length L this implies a number of CC pumps ncc = πDL/l2. The mean diameter D is then determined by equating ncc to the product of expected sponge average CC number density n and the volume of the IC canal plus its wall, VIC = (π/4) (D + 2Dcc)2L, ncc = n VIC. We obtain the number of IC canals from the count of CCs, nIC = n Vsp/ncc, and the pumping rate per CC, which is the variable Q in Equation (A1), becomes Qcc = Q/(nIC ncc). To complete the model of the aquiferous canal system, we consider an IC canal of varying diameter D(x) due to outflow to the CCs along its wall, i.e., dQ/dx = −(πD/l2) Qcc, and in which the pressure gradient of Poiseuille flow −dp/dx = (128µD4) Q(x) is furthermore assumed to be constant. This leads to a diameter variation D(x)/D0 = (1 − x/L)1/3, where D0 denotes the inlet diameter at x = 0 in the canal of length L, and a mean diameter of D = ¾ D0. It is convenient to express the constant pressure gradient by its value at the inlet where the total flow rate is Q0 = nccQcc = πDL/l2 Qcc and using also D = ¾ D0 leads to −dp/dx = (96µ/D03) (L/l2)Qcc, and a total pressure loss in the IC canal of
ΔpIC = (96µ/D03) (L/l)2 Qcc
where µ (=0.00127 Pa s at 15 °C) is the dynamic viscosity.
We now turn to the other pressure losses of flow through the sponge. According to Figure 2 in [24] there appears to be several, say nost, ostia apertures of diameter Dost leading the flow Q0 = ncc Qcc into one IC canal. Treating each as an orifice Equation (4-29.2) in [32] gives,
Δpost = (24µ/Dost3) ncc Qcc/nost
Given the fraction fost of total outer surface area taken up by ostia, their number per IC inlet may be estimated from nost = fost (D0/Dost)2. (For example, Table 2 in [24] gives fost = 30% and fIC = 10% for Haliclona permollis which for D0 = 100 µm and Dost = 20 µm gives nost ≈ 8). Similarly, for flow through npro prosopyles of diameter Dpro at the entrance to CCs,
Δppro = (24µ/Dpro3) ncc Qcc/npro
Given nch choanocytes per CC, the mean velocity of approach to the microvilli collar filter is u = (Qcc/nch)/(πDcolLcol), and the pressure drop across the array of cylindrical fibrils (microvilli) is estimated as in [19], using the theory of [33],
Δpcol = Kµu/dmv
where K = 8τ/(1 − 2lnτ + τ2/6), τ = π dmv/bmv, dmv (=0.14 µm) is the fibril diameter and bmv (=0.25 µm) the spacing.
From the apopyle of the CCs, flow enters EC canals, in which—depending on the architecture of the aquiferous canal system—there may be a pressure loss equal to that of Equation (A2) for the IC canal. Flow from many EC canals cumulate in the atrial cavity to give the upward flow rate Qa that varies with height z and reaches the maximum Qa(H) = Q = nIC ncc Qcc at the osculum at height H. The number of EC canals per unit of sponge height is nIC/H, each of which provides the volume flow ncc Qcc, hence the continuity equation of flow in atrium that describes the variation of Qa becomes
dQa/dz = (nIC/H) ncc Qcc
Because of the upward pressure decrease due to friction of the atrium flow the choanocyte chamber pumps deliver increasing flow with increasing height. This spatial variation gives rise to a complicated analysis, so to simplify the problem we approximate the total pressure loss in the atrial cavity flow by that of half the total flow acting over the full height,
Δpa = ½(128µHDa4) nIC ncc Qcc
It follows that the flow from EC canals, hence from the CC pumps, is approximated by a constant mean value Qcc. The average pump units at height z = ½H face a contribution that is one-half that of Equation (A7) while those below and above face more and less, respectively. Finally, the exit pressure loss at osculum is given by the kinetic energy ½ρUex2 of the exhalant jet, where the Uex = nIC ncc Qcc/(πDosc2/4), hence
Δpex = ½ρ[nIC ncc Qcc/(πDosc2/4)]2
The formulation of the problem is completed by inserting the sum of pressure losses, Equations (A2)–(A5), (A7) and (A8), into Equation (A1),
ΣΔp = Δpmax (1 − Qcc/Qmax)1/3
which is solved for Qcc by a Newton iteration, and then all other quantities can be calculated, such as total and volume-specific pumping rates, Q = nIC ncc Qcc and Q/V.

Appendix A.4. Simple Scaling Model

Consider the simple linear pump characteristic in terms of flow per choanocyte, Qch,
Qch = Qmax,ch (1 − Δpchpmax,ch)
where Qmax,ch and Δpmax,ch are constants. We assume constant density of choanocyte pumps acting in parallel and facing the same total pressure loss, Δpch, of the flow through the aquiferous system. Expressing this pressure loss as that of flow rate Qcan in a canal of length Lcan~L of scale independent diameter gives the scaling Δpch~QcanL. The number of canals is proportional to the sponge surface area ~L2 hence Qcan~Q/L2 and Δpch~ Q/L, where the total pumping rate is proportional to the volume Q~ QchV, V~L3 because of the constant choanocyte density. Using these relations to eliminate Qch and Δpch from Equation (A10) gives Q/V~c1 (1−c2 Q/L) or
Q/V~c3/(c + V2/3)
where c’s are constants. Equation (A11) leads to Q/V~V−2/3 for c << V and approaches Q/V~constant for large values of the constant c, i.e., b < 0 approaching b~0, a range of trends shown by much of the data reported in Table 1.
Table A1. Nomenclature and sample values of the ‘standard’ sponge, Haliclona urceolis that serves as reference case (CC = choanocyte chamber, IC = incurrent canal, EC = excurrent canal).
Table A1. Nomenclature and sample values of the ‘standard’ sponge, Haliclona urceolis that serves as reference case (CC = choanocyte chamber, IC = incurrent canal, EC = excurrent canal).
SymbolDescriptionStd. SpongeUnits
D0 = (4/3)DIC inlet diameter0.136–0.313mm
DIC mean diameter0.102–0.235mm
DaAtrium diameter2.3mm
DccCC diameter0.030mm
DoscOsculum diameter1.5mm
DostOstia diameter0.0206mm
DproProsopyle diameter0.005mm
DOSponge outer diameter8.6mm
HSponge height32mm
lCC spacing in IC canal0.036–0.030mm
L = ½(DODa)IC canal length3.15mm
nDensity of CCs in sponge12000mm−3
nchNo. of choanocytes per CC chamber80-
ncc = πDL/l2No. of CC along one IC canal778–2585-
nIC = n Vsp/nccNo. of IC (or EC) canals in sponge26,628–8013-
nostNo. of ostia per IC canal inlet8–13-
nproNo. of prosopyles in CC3-
Q0Flow rate at IC inlet
QaFlow rate in atrial cavity
QccPumping rate of CC4760µm3 s−1
Q = nIC ncc QccTotal pumping rate of sponge5.91mL min−1
Q/VVolume-specific pumping rate3.43min−1
VICVolume of IC canal with wall
VspVolume of sponge1726mm3
xIC canal coordinate
zVertical coordinate of atrium
QmaxCC pump characteristic parameter4828µm3 s−1
ΔpmaxCC pump characteristic parameter26.4Pa
Table A2. Pressure losses and pumping rate of present Haliclona urceolus model for the reference case with choanocyte chamber spacing l = 0.036 mm and counting ΔpIC once, compared to those of ([19] Table 3 therein]).
Table A2. Pressure losses and pumping rate of present Haliclona urceolus model for the reference case with choanocyte chamber spacing l = 0.036 mm and counting ΔpIC once, compared to those of ([19] Table 3 therein]).
Present ModelReference [19]
Pa%Pa%
OstiaΔpost0.78140.376
IC/EC canal flowΔpIC1.40262.3636
ProsopylesΔppro0.3061.1317
Collar filterΔpcol0.2551.2018
Atrium flowΔpa1.152100
Exahalent jetΔpex1.56291.5523
Total pressure lossΔptot5.431006.60100
Pumping rate (mL min−1)Qsp5.9 6.0
Figure A1. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref) for 3 cases constrained by the experimental correlation OSA = a4Vb4. (A): a4 = 9.42 and b4 = 0.6 and (B): a4 = 3.81 and b4 = 0.56 from Figure 5A in [7]. Isometric increase of H, DO, Da by scale factor sf = 0.5 to 4 for 4 (dashed), H constant, and diameters DO, Da constant, respectively, (solid). Values of exponents b, b1 and b2 (see Equation (1)) of power-law fit for sf > 1.5 shown.
Figure A1. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref) for 3 cases constrained by the experimental correlation OSA = a4Vb4. (A): a4 = 9.42 and b4 = 0.6 and (B): a4 = 3.81 and b4 = 0.56 from Figure 5A in [7]. Isometric increase of H, DO, Da by scale factor sf = 0.5 to 4 for 4 (dashed), H constant, and diameters DO, Da constant, respectively, (solid). Values of exponents b, b1 and b2 (see Equation (1)) of power-law fit for sf > 1.5 shown.
Jmse 09 01308 g0a1
Figure A2. Tubular sponge flow model. Variation of total pressure loss and individual contributions as a function of scaled size increase for the 3 cases of Figure A1A. (A) H constant, DO, Da scaled by factor sf = 0.5 to 4, (B) H, DO, Da scaled (isometric case), (C) DO, Da constant, H scaled. Pressure losses in incurrent and excurrent canals is the dominant contribution in (A,B) (by 50–85% and 20–75%, respectively, of the total pressure loss), while it is that of the atrium flow in (C) (by 26–75% of the total pressure loss). Pressure loss in ostia (ost), prosopyles (pro), incurrent and excurrent canals (IC, EC), collar (col), atrium (atr), exhalant jet (ex), and the sum (tot).
Figure A2. Tubular sponge flow model. Variation of total pressure loss and individual contributions as a function of scaled size increase for the 3 cases of Figure A1A. (A) H constant, DO, Da scaled by factor sf = 0.5 to 4, (B) H, DO, Da scaled (isometric case), (C) DO, Da constant, H scaled. Pressure losses in incurrent and excurrent canals is the dominant contribution in (A,B) (by 50–85% and 20–75%, respectively, of the total pressure loss), while it is that of the atrium flow in (C) (by 26–75% of the total pressure loss). Pressure loss in ostia (ost), prosopyles (pro), incurrent and excurrent canals (IC, EC), collar (col), atrium (atr), exhalant jet (ex), and the sum (tot).
Jmse 09 01308 g0a2aJmse 09 01308 g0a2b

References

  1. de Voogd, N.J.; Alvarez, B.; Boury-Esnault, N.; Carballo, J.L.; Cárdenas, P.; Díaz, M.-C.; Dohrmann, M.; Downey, R.; Hajdu, E.; Hooper, J.N.A.; et al. World Porifera Database. 2021. Available online: http://www.marinespecies.org/porifera (accessed on 10 November 2021).
  2. Reiswig, H.M. Water transport, respiration and energetics of three tropical marine sponges. J. Exp. Mar. Biol. Ecol. 1974, 14, 231–249. [Google Scholar] [CrossRef]
  3. Pile, A.J.; Patterson, M.R.; Witman, J.D. In situ grazing on plankton <10 μm by the boreal sponge. Mar. Ecol. Prog. Ser. 1996, 141, 95–102. [Google Scholar] [CrossRef] [Green Version]
  4. Ribes, M.; Coma, R.; Gili, J.-M. Natural diet and grazing rate of the temperate sponge Dysidea avara (Demospongiae, Dendroceratida) throughout an annual cycle. Mar. Ecol. Prog. Ser. 1999, 176, 179–190. [Google Scholar] [CrossRef]
  5. Coppari, M.; Gori, A.; Viladrich, N.; Saponari, L.; Canepa, A.; Grinyó, J.; Olariaga, A.; Rossi, S. The role of Mediterranean sponges in benthic–pelagic coupling processes: Aplysina aerophoba and Axinella polypoides case studies. J. Exp. Mar. Biol. Ecol. 2016, 477, 57–68. [Google Scholar] [CrossRef]
  6. Morganti, T.; Ribes, M.; Yahel, G.; Coma, R. Size Is the Major Determinant of Pumping Rates in Marine Sponges. Front. Physiol. 2019, 10, 1474. [Google Scholar] [CrossRef] [Green Version]
  7. Morganti, T.M.; Ribes, M.; Moskovich, R.; Weisz, J.B.; Yahel, G.; Coma, R. In situ Pumping Rate of 20 Marine Demosponges Is a Function of Osculum Area. Front. Mar. Sci. 2021, 8, 1–17. [Google Scholar] [CrossRef]
  8. Ereskovskii, A.V. Problems of coloniality, modularity, and individuality in sponges and special features of their mor-phogeneses during growth and asexual reproduction. Russ. J. Mar. Biol. 2003, 29, 46–56. [Google Scholar] [CrossRef]
  9. Fry, W.G. The sponge as a population: A biometric approach. Symp. Zool. Soc. Lond. 1970, 25, 135–162. [Google Scholar]
  10. Fry, W.G. Taxonomy, the individual and the sponge. Biology and Systematics of Colonial organisms. Syst. Ass. Spec. 1979, 11, 49–80. [Google Scholar]
  11. Reiswig, H.M. In situ pumping activities of tropical Demospongiae. Mar. Biol. 1971, 9, 38–50. [Google Scholar] [CrossRef]
  12. Savarese, M.; Patterson, M.R.; Chernykh, V.I.; Fialkov, V.A. Trophic effects of sponge feeding within Lake Baikal’s littoral zone. 1. In situ pumping rates. Limnol. Oceanogr. 1997, 42, 171–178. [Google Scholar]
  13. McMurray, S.E.; Pawlik, J.R.; Finelli, C.M. Trait-mediated ecosystem impacts: How morphology and size affect pumping rates of the Caribbean giant barrel sponge. Aquat. Biol. 2014, 23, 1–13. [Google Scholar] [CrossRef] [Green Version]
  14. Riisgård, H.U.; Kumala, L.; Charitonidou, K. Using the F/R-ratio for an evaluation of the ability of the demosponge Halichondria panicea to nourish solely on phytoplankton versus free-living bacteria in the sea. Mar. Biol. Res. 2016, 12, 907–916. [Google Scholar] [CrossRef]
  15. Kealy, R.A.; Busk, T.; Goldstein, J.; Larsen, P.S.; Riisgård, H.U. Hydrodynamic characteristics of aquiferous modules in the demosponge Halichondria panicea. Mar. Biol. Res. 2019, 15, 531–540. [Google Scholar] [CrossRef]
  16. Ludeman, D.A.; Reidenbach, M.A.; Leys, S.P. The energetic cost of filtration by demosponges and their behavioural re-sponse to ambient currents. J. Exp. Biol. 2017, 220, 995–1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Goldstein, J.; Riisgård, H.U.; Larsen, P.S. Exhalant jet speed of single-osculum explants of the demosponge Halichondria panicea and basic properties of the sponge-pump. J. Exp. Mar. Biol. Ecol. 2019, 511, 82–90. [Google Scholar] [CrossRef]
  18. Larsen, P.S.; Riisgård, H.U. The Sponge Pump. J. Theor. Biol. 1994, 168, 53–63. [Google Scholar] [CrossRef]
  19. Riisgård, H.U.; Thomassen, S.; Jakobsen, H.; Weeks, J.; Larsen, P.S. Suspension feeding in marine sponges Halichondria panicea and Haliclona urceolus:effects of temperature on filtration rate and energy cost of pumping. Mar. Ecol. Prog. Ser. 1993, 96, 177–188. [Google Scholar] [CrossRef]
  20. Fiore, C.L.; Baker, D.M.; Lesser, M.P. Nitrogen Biogeochemistry in the Caribbean Sponge, Xestospongia muta: A Source or Sink of Dissolved Inorganic Nitrogen? PLoS ONE 2013, 8, e72961. [Google Scholar] [CrossRef] [Green Version]
  21. Southwell, M.W.; Weisz, J.B.; Martens, C.S.; Lindquist, N. In situ fluxes of dissolved inorganic nitrogen from the sponge community on Conch Reef, Key Largo, Florida. Limnol. Oceanogr. 2008, 53, 986–996. [Google Scholar] [CrossRef] [Green Version]
  22. Dahihande, A.S.; Thakur, N.L. Temperature- and size-associated differences in the skeletal structures and osculum cross-sectional area influence the pumping rate of contractile sponge Cinachyrella cf. cavernosa. Mar. Ecol. 2019, 40. [Google Scholar] [CrossRef]
  23. Thomassen, S.T.; Riisgård, H.U. Growth and energetics of the sponge Halichondria panicea. Mar. Ecol. Prog. Ser. 1995, 128, 239–246. [Google Scholar] [CrossRef]
  24. Reiswig, H.M. The aquiferous systems of three marine demospongiae. J. Morphol. 1975, 145, 493–502. [Google Scholar] [CrossRef] [PubMed]
  25. Leys, S.P.; Yahel, G.; Reidenbach, A.; Tunnicliffe, V.; Shavit, U.; Reiswig, H. The sponge pump: The role of current induced flow in the design of the sponge body plan. PLoS ONE 2011, 6, e27787. [Google Scholar] [CrossRef] [Green Version]
  26. Riisgård, H.U.; Larsen, P.S. Filter-Feeding in Marine Macro-Invertebrates: Pump Characteristics, Modelling and Energy Cost. Biol. Rev. 1995, 70, 67–106. [Google Scholar] [CrossRef]
  27. Asadzadeh, S.S.; Larsen, P.S.; Riisgård, H.U.; Walther, J.H. Hydrodynamics of the leucon sponge pump. J. R. Soc. Interface 2019, 16, 20180630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Larsen, P.S.; Riisgård, H.U.; Asadzadeh, S.S.; Walther, J.H. Scaling in Filter-Feeding Sponges. In Proceedings of the ALSO Aquatic Science Virtual Meeting, 22–27 June 2021. Abstract 2899. [Google Scholar]
  29. Lüskow, F.; Riisgård, H.; Solovyeva, V.; Brewer, J. Seasonal changes in bacteria and phytoplankton biomass control the condition index of the demosponge Halichondria panicea in temperate Danish waters. Mar. Ecol. Prog. Ser. 2019, 608, 119–132. [Google Scholar] [CrossRef]
  30. Vogel, S. Life in Moving Fluids, 2nd ed.; Princeton University Press: Princeton, NJ, USA, 1994. [Google Scholar]
  31. Hammel, J.U.; Filatov, M.V.; Herzen, J.; Beckmann, F.; Jaap, A.; Kaandorp, J.A.; Nickel, M. The non-hierarchical, non-uniformly branching topology of a leuconoid sponge aquiferous system revealed by 3D reconstruction and morpho-metrics using corrosion casting and X-ray microtomography. Acta Zool. 2012, 93, 160–170. [Google Scholar] [CrossRef]
  32. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics; Springer: Singapore, 1981. [Google Scholar]
  33. Tamada, K.; Fujikawa, H. The steady two-dimensional flow of viscous fluid at low reynolds numbers passing through an infinite row of equal parallel circular cylinders. Q. J. Mech. Appl. Math. 1957, 10, 425–432. [Google Scholar] [CrossRef]
Figure 1. Power-law regression lines of experimental data showing trends. Numbers refer to the data of Table 1.
Figure 1. Power-law regression lines of experimental data showing trends. Numbers refer to the data of Table 1.
Jmse 09 01308 g001
Figure 2. Density of choanocytes in various demosponge species studied by some authors [24,25].
Figure 2. Density of choanocytes in various demosponge species studied by some authors [24,25].
Jmse 09 01308 g002
Figure 3. Model of tubular-type demosponge. (A) Haliclona permollis (redrawn from Figure 2 in [24]. (B) Schematic vertical section of model defined by four length scales, H, DO, Da, Dosc, which also define canal length L = (DODa)/2 and volume V = (π/4)DO2H. The model shows flow rate Q0 entering an incurrent canal (IC), driven by choanocyte pump units (CC), and leaving an excurrent canal (EC) to the atrial cavity in which flow rate Qa cumulates to the sponge pumping rate Q that leaves the osculum as an exhalant jet with velocity Uex. (C) Schematic detail showing the pressure variation through body of the sponge with inlet pressure loss Δpost at ostia, pressure rise Δpp at choanocyte chamber pump, and prevailing pressure Δpa at atrium at given height z above the bottom of the sponge of height H.
Figure 3. Model of tubular-type demosponge. (A) Haliclona permollis (redrawn from Figure 2 in [24]. (B) Schematic vertical section of model defined by four length scales, H, DO, Da, Dosc, which also define canal length L = (DODa)/2 and volume V = (π/4)DO2H. The model shows flow rate Q0 entering an incurrent canal (IC), driven by choanocyte pump units (CC), and leaving an excurrent canal (EC) to the atrial cavity in which flow rate Qa cumulates to the sponge pumping rate Q that leaves the osculum as an exhalant jet with velocity Uex. (C) Schematic detail showing the pressure variation through body of the sponge with inlet pressure loss Δpost at ostia, pressure rise Δpp at choanocyte chamber pump, and prevailing pressure Δpa at atrium at given height z above the bottom of the sponge of height H.
Jmse 09 01308 g003
Figure 4. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref) for the case of isometric scaling of H, DO, Da, Dosc~sf (upper curve) and for reduced scaling of osculum Dosc~sf0.5. (lower curve). Values of scale factor sf = 0.5, 1, 2, 3, and 4 give points marked by a star, and the reference sponge is marked by a circle. Shown values of power-law exponents may be compared to values b = −0.6 and −0.39 in Table 1. Data for Halichondria panicea (#4 of Table 1) are shown by squares and Haliclona urceolus (#5) by circles. Shown values of power-law exponents may be compared to values b = −0.6 and −0.39 in Table 1.
Figure 4. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref) for the case of isometric scaling of H, DO, Da, Dosc~sf (upper curve) and for reduced scaling of osculum Dosc~sf0.5. (lower curve). Values of scale factor sf = 0.5, 1, 2, 3, and 4 give points marked by a star, and the reference sponge is marked by a circle. Shown values of power-law exponents may be compared to values b = −0.6 and −0.39 in Table 1. Data for Halichondria panicea (#4 of Table 1) are shown by squares and Haliclona urceolus (#5) by circles. Shown values of power-law exponents may be compared to values b = −0.6 and −0.39 in Table 1.
Jmse 09 01308 g004
Figure 5. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref). (A) Effect of increasing diameters by scale factor sf = 0.5 to 4 for 4 values of fixed height. (B) Effect of increasing height by scale factor sf = 0.5 to 4 for 4 values of diameters for fixed annular structure thickness L = (DODa)/2 at its reference values. Values at sf = 0.5, 1, 2, 3, and 4 are marked on each curve by a star. The reference case is marked by a circle on the full-line curve.
Figure 5. Tubular sponge flow model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref). (A) Effect of increasing diameters by scale factor sf = 0.5 to 4 for 4 values of fixed height. (B) Effect of increasing height by scale factor sf = 0.5 to 4 for 4 values of diameters for fixed annular structure thickness L = (DODa)/2 at its reference values. Values at sf = 0.5, 1, 2, 3, and 4 are marked on each curve by a star. The reference case is marked by a circle on the full-line curve.
Jmse 09 01308 g005
Figure 6. Sponge explant cone model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref). The reference value (circle) corresponds to the geometry of ID#21 of Table 1 in [17]. The scaling factor range sf = 0.1 to 1.5 applied to height for constant diameters leads to a decreasing trend with increasing size, while applied to diameters at constant height it leads to an increasing trend (solid curves). Isometric scaling (dashed curve) leads to a weakly decreasing trend (b = −0.04).
Figure 6. Sponge explant cone model. Volume-specific pumping rate (Q/V) versus normalized size (V/Vref). The reference value (circle) corresponds to the geometry of ID#21 of Table 1 in [17]. The scaling factor range sf = 0.1 to 1.5 applied to height for constant diameters leads to a decreasing trend with increasing size, while applied to diameters at constant height it leads to an increasing trend (solid curves). Isometric scaling (dashed curve) leads to a weakly decreasing trend (b = −0.04).
Jmse 09 01308 g006
Table 1. Data on volume-specific pumping rate of demosponges (in mL = cm3 per min per cm3 of sponge = min−1) versus sponge volume (in cm3) approximated by the power-law, Q/V = aVb (min−1). V-range gives Vmin and Vmax and <Q/V> the algebraic average of data. Number (#) refers to the plot of regression lines in Figure 1; #8 and #9a are single-osculum explants; #9b individual modules of multi-osculum explants; #15 and #16 tropical sponges.
Table 1. Data on volume-specific pumping rate of demosponges (in mL = cm3 per min per cm3 of sponge = min−1) versus sponge volume (in cm3) approximated by the power-law, Q/V = aVb (min−1). V-range gives Vmin and Vmax and <Q/V> the algebraic average of data. Number (#) refers to the plot of regression lines in Figure 1; #8 and #9a are single-osculum explants; #9b individual modules of multi-osculum explants; #15 and #16 tropical sponges.
#ReferenceSpeciesabV-Range (cm3)<Q/V>
1[2]Mycale sp.15.1~ 0100 to 230014
2-“ -Verongia gigantea6.9~ 0600 to 10,0003.3
3-“ -Tethya crypta9~ 0100 to 170011
4[19]Halichondria panicea8.18−0.64 to 172.7
5-“ -Haliclona urceolus3.96−0.390.6 to 452.5
[23]H. panicea3.12−0.1191.3 to 962.7
[21]14 species6.42−0.045300 to 109,0006.3
[20]Xestospongia muta148−0.34812,000 to 110,0005.5
[13]X. muta ~017.8 to 451,6493.6
6[14]H. panicea2−0.821 to 12.56.1
7[16]5 species3.770.0511.9 to 4496.7
8[17]H. panicea2.3−0.0380.42 to 15.63.6
9a[15]H. panicea1.31−0.190.09 to 0.611.8
9b-“ -H. panicea0.9850.7260.1 to 3.81.24
[22]Cinachyrella cf. cavernose0.944−0.0710.1 to 1001.1
10[6]Dycidea avara17.88−0.221.5 to 148
11-“ -Crambe crambe26.24−0.290.8 to 25
12-“ -Petrosia ficiformis5.01−0.193.12 to 421
13-“ -Chondrosia reniformis26.31−0.492.72 to 181
14-“ -Agelas oroides30.39−0.702.98 to 440
15[7]7 species, Red Sea66.3−0.521 to 120
16-“ -8 species, Caribbean176−0.4350 to 35,500
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Larsen, P.S.; Riisgård, H.U. Pumping Rate and Size of Demosponges—Towards an Understanding Using Modeling. J. Mar. Sci. Eng. 2021, 9, 1308. https://doi.org/10.3390/jmse9111308

AMA Style

Larsen PS, Riisgård HU. Pumping Rate and Size of Demosponges—Towards an Understanding Using Modeling. Journal of Marine Science and Engineering. 2021; 9(11):1308. https://doi.org/10.3390/jmse9111308

Chicago/Turabian Style

Larsen, Poul S., and Hans Ulrik Riisgård. 2021. "Pumping Rate and Size of Demosponges—Towards an Understanding Using Modeling" Journal of Marine Science and Engineering 9, no. 11: 1308. https://doi.org/10.3390/jmse9111308

APA Style

Larsen, P. S., & Riisgård, H. U. (2021). Pumping Rate and Size of Demosponges—Towards an Understanding Using Modeling. Journal of Marine Science and Engineering, 9(11), 1308. https://doi.org/10.3390/jmse9111308

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