A mathematical model of the commercial harvest of <italic>Palmaria palmata</italic> (Palmariales, Rhodophyta) on Digby Neck, Nova Scotia, Canada
A mathematical model of the commercial harvest of Palmaria palmata (Palmariales, Rhodophyta) on Digby Neck, Nova Scotia, Canada
ALGAE. 2012. Mar, 27(1): 43-54
Copyright ©2012, The Korean Society of Phycology
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://cre-ativecommons. org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : January 26, 2012
  • Accepted : February 29, 2012
  • Published : March 15, 2012
Export by style
Cited by
About the Authors
Ryan J. Lukeman
Department of Mathematics, Statistics and Computer Science, St. Francis Xavier University, Antigonish, NS B2G 2W5, Canada
Leah F. Beveridge
Department of Biology, St. Francis Xavier University, Antigonish, NS B2G 2W5, Canada
Andrea D. Flynn
Department of Biology, St. Francis Xavier University, Antigonish, NS B2G 2W5, Canada
David J. Garbary
Department of Biology, St. Francis Xavier University, Antigonish, NS B2G 2W5, Canada

A mathematical model of the commercial harvest of Palmaria palmata (Dulse) is presented based on a logistic model and field data collected on Digby Neck, Nova Scotia from 14 harvested shores during May to August, 2010. Field observations used to estimate model parameters included cover of Dulse before and after harvest from Dulse dominated boulders for which surface area was estimated, and from which fresh biomass of harvested Dulse was weighed. Over all the surveys the average harvest fraction was about 50%, and the total resource was about 1,600 g m -2 . With growth rates in excess of 4% per day and a 50% harvest of the standing crop each month, the model suggests that the Dulse resource is sustainable at current harvest levels.
Palmaria palmata (Linnaeus) Weber et Mohr (common name Dulse) is an abundant low intertidal and shallow subtidal zone red seaweed in the Bay of Fundy (Taylor 1957, Edelstein et al. 1970, Ffrench 1974, Wilson et al. 1979, South et al. 1988, Vadas et al. 2004). The species is commercially harvested in eastern Canada (MacFarlane 1964, 1966, Ffrench 1974, Bird and Van der Meer 1993, Chopin and Ugarte 2006). The harvesting of Dulse is a significant cottage industry in Nova Scotia (MacFarlane 1964, 1966, Ffrench 1974). The edible alga can be eaten dried or wet, and can be used in cooking (e.g., Rhatigan 2009), and its antioxidant properties may have medicinal benefits (e.g., Yuan et al. 2009, Cornish and Garbary 2010). There have been no significant studies of the Dulse harvest in this area since MacFarlane (1966) and Ffrench (1974), and few ecological studies of Dulse have been conducted in this part of the Bay of Fundy. The ecological studies by Vadas et al. (2004) of algal biomass and production in northern Maine that included extensive data on P. palmata is an important exception.
It is important to examine the abundance of Dulse and its growth patterns in order to evaluate whether the current harvesting activities are sustainable. This is in addition to larger scale management issues that relate to macrophyte management in general (e.g., Chambers et al. 1999). Policies for sustainable development, integrated management, and precautionary approaches are in place under the Canadian Oceans Act (Sharp and Bodiguel 2001); however, there are no legal or management restrictions on Dulse harvesting. Thus, Dulse may be considered an open resource. As a consequence of the lack of regulation of Dulse, the sustainability of commercially harvestable populations of P. palmata on Digby Neck requires study.
The harvest of Dulse on Digby Neck would appear to be sustainable as it has existed on a commercial scale for at least a century (French 1974, Chopin and Ugarte 2006). In the Fundy region, Dulse harvesting is governed by the semidiurnal tidal cycle; harvesters periodically pick the resource during the approximately 7 days during the lowest (spring) tides of the 28-day tidal cycle. The local people and businesses of Digby Neck harvest Dulse during low tides throughout the summer months (Chopin and Ugarte 2006). The method of harvest has low disturbance because it is only done by hand and the rocks are not picked completely clean; only the longer, easily picked fronds are removed, leaving most of the holdfasts and smaller fronds to regenerate. The same beaches are picked repeatedly over the harvest season, and the Dulse has time to regrow between successive spring tides. Unlike Ascophyllum nodosum (Linnaeus) Le Jolis, Dulse is not harvested in vast quantities; the average harvest on the Nova Scotia coast of the Bay of Fundy from 1982-1991 was 2 dry tonnes (Chopin and Ugarte 2006). This is much less compared to the New Brunswick shores of the Bay of Fundy, which, during the same time frame, had an average yield of 38 dry tonnes (Chopin and Ugarte 2006). This difference between the two sides of the Bay Fundy suggests that either the resource in Nova Scotia is much less, or it is not being harvested to its maximal potential. To address the latter issue we carried out field studies on which to base a mathematical model of the harvest. Our study aims to determine the optimal harvest quantities that will maximize yields while maintaining the sustainability of the resource. While this is largely a theoretical study, it is based on actual harvest data conducted during commercial harvesting.
- Introduction to resource modelling in a seaweed context
Biological growth processes are modeled in a variety of ways, depending on the characteristics of the system, the availability of data for parameter estimation, and on choices of the modeller (see Edelstein-Keshet 1988 for an excellent overview). In the context of seaweed growth, a wide variety of models arise, depending on the above factors. Where the temporal domain spans years or generations, discrete models are typically used. Discrete models are composed of difference equations, which are recursive relations usually in the form At+1 = f ( At ), such that the state of the system (e.g., population) at the next time period depends on some function of the current state of the system. Conversely, for temporal domains lacking a generational structure (e.g., growth within one season), continuous models are appropriate, and lead to differential equations in the form dp / dt = f ( p , t ).
The dimension of the model is governed by the structure of the seaweed population modeled. Total biomass (i.e., production) models (e.g., Seip 1980, Lee and Ang 1991) represent the entire population by one variable, whereas populations that contain a class structure lead to higher-dimensional models (Ang and De Wreede 1990, Nyman et al. 1990, Åberg 1992 a , 1992 b , Duarte and Ferreira 1997, Santos and Nyman 1998). Common classes include age (‘age-structured models’), size, or developmental stage. The models that result from this classification generally include a system of difference or differential equations, with one equation per class.
For growth of a population subject to environmental limitation, the most basic and widely used model is the logistic model (Tsoularis and Wallace 2002), originally developed by Verhulst (1838). The central assumption of this model is that the relative growth rate (percentage growth per unit time) decreases as the population approaches its environmental carrying capacity; the simplest assumption which leads to the classical logistic model is a linear decrease. Evidence of logistic-type growth in seaweeds include reports of increased growth rate attributed to reduced standing crop post-harvest for kelp (Smith 1986), and growth curves of Ascophyllum nodosum (Linnaeus) Le Jolis that are well-approximated by a logistic curve (i.e., of sigmoid type) (Baardseth 1955). More recently, in cultivation bags of Palmaria palmata , frond growth decreased with increased frond density (Martínez et al. 2006). Both continuous (e.g., Seip 1980, Lee and Ang 1991) and discrete (Scrosati 2004) logistic growth models have been applied to seaweeds.
Logistic models are parameterized by the intrinsic growth rate, r , and carrying capacity, K , the values of which vary according to species and environment. In one tank-cultivation study, P. palmata , cultivated at a density of 4.2 kg m -2 , had maximal yield of 609 g m -2 per day, and a maximal growth rate of 12% per day (Pang and Lüning 2004). Lower growth rates of 2-8% per day were obtained by Morgan and Simpson (1981 a , 1981 b ). In an open-sea cultivation study, the maximal growth rate of Dulse was 14% of the initial fresh weight per day (Martínez et al. 2006). On the other hand, in the natural environment, maximal specific growth rates reported range from 0.89% (Northern Ireland unpublished) to 1.5% (Spain) (Faes and Viejo 2003).
We construct a one-year, continuous logistic growth model for the growth of P. palmata on Digby Neck. We
PPT Slide
Lager Image
Map of Nova Scotia (A) with boxed area enlarged (B) indicating Digby Neck and adjacent Long and Brier Islands with double headed arrow indicating general area of most collection beaches and arrows indicating other beaches harvested in this study.
then investigate the effects of harvesting on the resource, both in the context of the current harvest, and the optimal harvest possible.
Dulse collection beaches from northeast to southwest along the coast of Digby Neck, Long Island and Brier Island used for developing the logistic modelNote: sites 1-9, Zone 20T; sites 10-14, Zone 19T.aLong Island.bBrier Island.
PPT Slide
Lager Image
Dulse collection beaches from northeast to southwest along the coast of Digby Neck, Long Island and Brier Island used for developing the logistic model Note: sites 1-9, Zone 20T; sites 10-14, Zone 19T. aLong Island. bBrier Island.
Data was collected during the one-week period of the spring tides of each month from May to August along the western coast of Digby Neck ( Table 1 , Fig. 1 ). These months span the primary period of the Dulse harvest with reduced harvesting extending into September and early October depending upon weather conditions. During each low tide series, beaches were accessed by foot or with local Dulse harvesters by boat. Beach and transect locations were recorded using a hand-held GPS (Garmin etrex, Olathe, Kansas). Arrival time at the beach was approximately one hour before the documented time for low tide ( ) using Saint John (NB), Digby (NS) or Westport (NS) as reference ports. Twenty-five rocks on each beach were picked by hand to the extent demonstrated by a local harvester. Before- and after-harvest percent cover of Dulse on each rock was recorded ( cs , cf , respectively) as well as the size of the rock being picked. As the Dulse was picked, it was put into a mesh onion bag, swung to remove excess water, and weighed using a 2,000 g spring scale (Stansi Scientific Com, Chicago, IL, USA) with a resolution of 25 g. Resource density, R was obtained by taking the quantity of resource harvested divided by area ( w / A ) and then scaling by the inverse of the relative fraction of Dulse cover picked, via
PPT Slide
Lager Image
The median value of all 1 m-2 density measurements per shore; n = 25. Note: the two low values from August are from nonharvested shores.
PPT Slide
Lager Image
Thus R represents the approximate total available resource (how much would be harvested should all the Dulse be removed from a 1 m 2 patch). R was calculated for each rock measurement, and the median value of these measurements was taken for each shore ( Fig. 2 ). By averaging all the data, a carrying capacity for P. palmata on Digby Neck was estimated.
- The model
To understand the growth dynamics and potential harvest of P. palmata , we construct a model for one-year growth in a domain of size 1 m 2 , which is meant to represent an average patch of the resource in the harvest area. As a starting point, we assume that the growth dynamics are well approximated by the classical logistic model
PPT Slide
Lager Image
where x (grams) is the amount of the resource, r (percent/day) is the fixed per-capita growth rate throughout time, and K (grams) is the carrying capacity; in this context, the amount of the resource found on an unharvested 1 m 2 patch is assumed to have grown to maturity.
To decouple the model results with the estimated carrying capacity, the logistic model is non-dimensionalized via the substitution
PPT Slide
Lager Image
, in which case Eq. 2 reduces to
PPT Slide
Lager Image
PPT Slide
Lager Image
Graphical representation of terms in Eq. 5. The growth rate r(t) (Eq. 4) is modeled using a piecewise quadratic function with growth occurring for 214 days, March through September inclusively, with a maximal growth rate of rmax.
where ^ notation has been removed for simplification, and x is now interpreted as the unitless fraction of the carrying capacity K for P. palmata , and therefore, x ranges between 0 and 1, assuming the initial resource density is below K .
- Growth rate
Environmental factors such as light level, water temperature, average wave height, etc., result in an intrinsic growth rate that is time-dependent (Faes and Viejo 2003). For instance, unpublished data from Northern Ireland showed that Palmaria palmata exhibited dominating growth in the summer months, whereas dominating decay through frond breakage occurred from August to March. Therefore, the assumption of fixed per-capita growth rate r must be modified, and a suitable r ( t ) constructed. Our model does not require growth between October and February inclusively as this is largely a period of large frond reproduction and senescence, and initiation of small fronds. For the remaining months, we approximate r ( t ) with the quadratic function ( Fig. 3 ), leading to
PPT Slide
Lager Image
where rmax is the maximal per-capita growth rate, and ts and tf are the day of year for which positive new growth starts and finishes, respectively. We set ts = 59 days and tf = 273 days (a total of 214 days) to correspond to positive growth during March through September inclusively. Vadas et al. (2004) estimated a frond generation time
PPT Slide
Lager Image
Graphical representation of terms in Eq. 5. The harvest function h(t), comprised of four (instantaneous) annual harvests that occur once a month on Digby Neck and Islands on a 28-day cycle between May and August inclusively.
of 30 days, and observed growth for 7 months (roughly 210 days). We cannot estimate per-capita growth rate for the resource in the Fundy region because data were not collected at precisely the same location across monthly harvests. Thus, we let rmax vary over the range [0.04, 0.15], which corresponds to a range of 4 to 15% daily growth. The minimum rmax value (4%) was chosen because it corresponds to approximately 50% logistic regrowth of the resource over the 28 day harvest cycle, assuming half the resource has been previously harvested, which is lower than observed regrowth in the data collected. Fifteen percent was chosen as an upper limit, somewhat higher than the maximal recorded growth rate of 14% (Martínez et al. 2006).
- Harvesting
Next, Eq. 3 is modified to include the effect of a harvest, giving
PPT Slide
Lager Image
where h ( t ) (per day) is a time-dependent harvesting policy. To model the harvesting cycle that takes place in the Fundy region, a periodic proportional impact-harvest function h ( t ) is used such that a constant fraction f of the resource is harvested at regular 28-day intervals. The fraction h is collected during one time-step of the model, as opposed to harvesting over a 7-day period because we are modeling a hypothetical 1 m 2 patch, which can be effectively harvested instantly, relative to model timescales. Fig. 4 demonstrates the four annual harvests that occur once a month on a 28-day cycle between May and August inclusively. Although a positive per-capita growth of the resource is observed through September, there is typically little harvest in this month due predominantly to cultural and customary reasons and vagaries of weather.
The total annual yield per square meter of resource x ( t ) is defined as y = ∫ h ( t ) x ( t ) dt . To obtain actual yield in grams, values can simply be scaled by the carrying capacity; i.e., Ky .
- Optimal harvesting
For a given value of rmax , total yield is the result of a trade-off between the frequency of harvesting and the harvest fraction h . To investigate the relative productivity of the current Fundy Dulse harvest, optimal harvesting strategies h ( t ) are investigated under two assumptions: a fixed growth rate [ r ( t ) = r , as in Eq. 2], for which analytical results are available, and the more realistic time-dependent growth rate r ( t ) as in Eq. 4. It is assumed that n harvests occur at equally spaced time intervals throughout [ ts , tf ] and for each, a constant fraction f of the resource x ( t ) is harvested. Harvest frequency n = 1, 2, . . . , 20, and the harvest fraction h is varied over the interval [0, 1]; no harvest to complete harvest.
To study optimality in the context of the Dulse harvest, the harvest period is modified to the typical summer harvest season May-August, in which it is assumed that n harvests of constant fraction h occur at equally spaced time intervals throughout the period. As done previously, rmax is varied as a model parameter.
Mean densities for each site time sample are given in Table 2 . In May, the three sampled shores had a mean density of 740 ± 75 g m -2 (mean ± standard deviation [SD] with values round to nearest 5 g). In June, six beaches were visited, two of which were along the same shoreline (Point Prim A and Point Prim B). Density dropped slightly to 685 ± 115 g m -2 , but there was no significant difference between months (p > 0.05). Six beaches were also surveyed in July, with a mean density of 640 ± 87 g m -2 . This, again, was lower than the previous month, but was most likely an underestimate because partial harvesting occurred on some beaches before data collection. There was no significant difference between the densities of beaches in June and July (p > 0.05). Eight beaches were sampled in August, two of which were from Brier Island and are non-harvested shores with low densities. Thus mean density of 755 ± 340 g m -2 is an underestimate of
PPT Slide
Lager Image
Example growth curve for resource x(t) over one year with four harvests. Parameters are rmax = 0.1, h = 0.478. The growth follows an initially logistic (sigmoidal) profile, after which harvesting periodically reduces the biomass. After each harvest, the rate of recovery is governed by the time-dependent growth rate r(t).
biomass available for harvest on harvest shores. Again, this value is not significantly different than in July (p > 0.05) even when the two Brier Island samples are omitted.
The pre- and post-harvest cover of P. palmata were calculated at cs = 0.795, cf = 0.419 respectively. Over all the surveys, the average harvest fraction was h = 0.478, and the average total resource was R = 1,601 g.
- Annual growth and harvesting
Density of Palmaria palmata picked from harvest boulders at all beaches and surveys used in construction of the logistic modelValues indicate mean ± standard deviation, n = 25.aLong Island.bBrier Island.
PPT Slide
Lager Image
Density of Palmaria palmata picked from harvest boulders at all beaches and surveys used in construction of the logistic model Values indicate mean ± standard deviation, n = 25. aLong Island. bBrier Island.
Fig. 5 demonstrates annual growth with rmax = 0.10. The total yield y , as a function of rmax is asymptotic to 4 h = 4(0.478) = 1.912 for increasing values or rmax . Using the carrying capacity estimate K = 1,601 g, we show actual yields (g m -2 ) in Fig. 6 .
- Optimal harvesting
For r ( t ) = r , analytical conditions of optimality have been derived in Zhang et al. (2003). Specifically, for a given value of r , the optimal harvest fraction depends on the harvest frequency n via
  • h* = 1 – e-rT/2(6)
where T = ( tf ts )/ n is the period of the harvest. Thus, as n increases, optimal harvest fraction is reduced, as intuition suggests. Although there are inherent difficulties in applying this optimality condition to the Fundy harvest, as the assumption of fixed growth rate is not realistic, nevertheless choosing T = 28 gives an optimal harvest fraction of h * that increases from 0.429 to 0.877 as r increases from 0.04 to 0.15. The long-term (sustainable) yield per unit time for this model is (Zhang et al. 2003):
PPT Slide
Lager Image
To obtain an approximate total seasonal yield (neglecting effects of initial quantity of the resource), the quantity
PPT Slide
Lager Image
The total annual yield (grams) as a function of the maximal growth rate rmax assuming a carrying capacity which results in an average density of 1,600 g m-2.
in Eq. 7 is multiplied by tf ts . Again, limitations arise in applying this analytical result to the Fundy harvest, yet for 4 harvests of period T = 28 days, tf ts = 112, and total yield ranges from 1.043 ( r = 0.04) to 2.98 ( r = 0.15).
For a clearer picture of how optimal yield depends on harvest fraction and number of harvests, simulations were performed over a range of values of the parameters h , n , and r . Yield is plotted as a function of both harvest frequency n and fraction h for r = 0.025, 0.05, 0.075, 0.10, and 0.125 in Fig. 7 A-E. The yield value at a given ( n , h ) is represented by the color at that point. Values of parameters n and h corresponding to the optimal yield are denoted by ( n *, h *). These plots demonstrate the trade-off between n and h : highest yields occur at frequent small-fraction harvests, infrequent large-fraction harvests, or some intermediate combination. Contour lines indicate that this trade-off is sharpest at low n values, in keeping with the negative exponential trade-off function in Eq. 6. For smaller values of r , yield is more sensitive to changes in ( n *, h *) ( Fig. 8 A & B). The non-smooth nature of these plots stems from the discrete variation of the parameter n , which must take on integer values. Clearly, higher r leads to optimal harvests of increasing frequency ( n ) and decreasing fraction per harvest ( h ).
When the time-dependent growth rate r ( t ) (Eq. 4) is incorporated, analytical results are no longer available, although optimal harvest fraction / frequency trade-off is qualitatively similar ( Fig. 9 A-E). Overall yields are reduced due to the lower average growth rate, and optimal harvest for intermediate and high values of rmax are obtained by more frequent, low-fraction harvests (limited in these plots by chosen range n ). For both fixed r and time-dependent r ( t ), Figs 7 and 9 indicate that yields very near the optimal are possible for a broad combination of n and h values, so long as the exponential-type trade-off rela-tionship is obeyed.
Optimal yields are greater than yields on Digby Neck and Islands (as expected), but only slightly for low values of rmax ( Fig. 10 ); e.g., for rmax < 0.05, the Fundy harvest is sub-optimal by no more than 0.175 K per m -2 , or approximately 280 g m -2 . As rmax increases, the harvest yield diverges from the optimal yield due to the increased benefit of more frequent harvesting due to the higher growth rate; e.g., for rmax = 0.10, the Fundy harvest yield is 0.888 K per m -2 less than the optimal yield, or approximately 1420 g m -2 less.
Palmaria palmata is a perennial, clonal species in which thalli can survive for more than one year, and form fronds that successively grow, reproduce and senesce from a holdfast that survives multiple years (i.e., a pseudo-perennial) (Garbary 1976, Vadas et al. 2004). The basic phenology of the species on Digby Neck is one where thalli grow vegetatively through the winter, spring and summer producing numerous fronds from a discoid holdfast. During late summer and fall fronds become increasing reproductive with spermatagia or tetrasporangia. Over the winter many large reproductive fronds may remain in the population, but as these fronds senesce the population becomes increasingly dominated by the new fronds developing from the holdfasts (Garbary et al. 2012, Katelyn White unpublished data). We do not have sufficient information on the relative importance of vegetative regeneration from holdfasts (vs. reproduction through tetraspores) to effectively study the population from a multi-year perspective. Also, at established locations, it appears that growth is predominantly from exist-
PPT Slide
Lager Image
Annual harvest yield as a function of harvest fraction h and frequency n for fixed r with (A) rmax = 0.025, (B) rmax = 0.05, (C) rmax = 0.075, (D) rmax = 0.10, and (E) rmax = 0.125. For each plot, the asterisk represents the maximal yield.
PPT Slide
Lager Image
Optimal harvest (A) frequency n* and (B) fraction h* as a function of the per-capita growth rate r.
PPT Slide
Lager Image
Annual harvest yield as a function of harvest fraction h and frequency n for r(t) with (A) rmax = 0.025, (B) rmax = 0.05, (C) rmax = 0.075, (D) rmax = 0.10, and (E) rmax = 0.125. For each plot, the asterisk represents the maximal yield.
PPT Slide
Lager Image
Total annual yield as a function of rmax for Fundy harvest (solid) and optimal harvest (n*, h*) (dotted).
ing holdfasts; accordingly, the ‘spread’ of the resource in space across generations is not modelled, and our time period is limited to one season. Furthermore, the specialized life history of P. palmata with its microscopic female gametophytes and absence of carposporophytes (Van der Meer and Todd 1980, Gabrielson and Garbary 1986) have not been studied in terms of perennation or colonization. While these are important constraints on the basic ecology of the species, they do not alter the key elements of our model.
Within one growth season, there is not a clear discretization of growth phases; we, therefore, opt for a continuous model over a discrete model. Because we have neither growth rate data nor frond size distributions throughout the growing season, we lack sufficient data to incorporate a class structure into our model. If we were to do so, the obvious choice would be frond size-classes; then harvest strategies that are based on picking particular-sized fronds could be investigated.
- Growth rates
The chosen range of maximal growth rate is most likely larger than the actual range; 4% is lower because, according to the data, almost complete regrowth takes place over the harvest cycle, while the maximal growth rate used in this model (15%) is greater than that found in the literature for both cultivated and natural growth. The range of growth rates is higher than that recorded from some non-aquaculture field samples, e.g., Faes and Viejo (2003), suggesting the growth rates for P. palmata in the Fundy region significantly exceed those in other regions (e.g., Spain, N. Ireland).
The most detailed evaluation of in situ productivity of P. palmata comes from Cobscook Bay in Maine (Vadas et al. 2004). Their values for biomass at three sites of about 25-120 g dry wt m -2 from June to August converted into 34-240 g dry wt m -2 y -1 at their low site (mean minimum and maximum), to 119-830 g dry wt m -2 y -1 (mean minimum to maximum) for the most productive site. Using a 1 : 10 ratio of dry to wet mass, the productivity of Cobscook Bay P. palmata (up to 8,300 g wet wt m -2 y -1 ) seems much higher than for Digby Neck (i.e., between 1,250 and 3,060 g m -2 y -1 , depending on growth rate). The reconciliation of the two studies is derived from the complete replacement of frond biomass each month over a seven month period (Maine) rather than our four months of our harvesting in which only 50% of the material is removed at each sampling. We also consider our value of 1,600 g wet wt m -2 y -1 an underestimate based on the fact that some of our beaches had previously been harvested when our quantification had occurred. Furthermore, the values from Vadas et al. (2004) are for average m 2 on the shore whereas our data is from harvested boulders that do not comprise 100% of the shore. Incidental observations of P. palmata populations on Digby Neck in 2011 suggest that large biomasses continue into September and October that remain unharvested (Garbary unpublished observations).
- Harvesting
Total harvest yield for a range of growth rates were obtained using the model based on the estimated carrying capacity K = 1,601 g. However, this value is an underestimate of the actual carrying capacity of the resource because, especially in July and August, some beaches were surveyed after partial harvesting had taken place, because of the time restrictions placed upon the harvesting by the tides, and because several of the beaches included in the model are not currently harvested.
- Optimal harvesting
In order for the Dulse industry of Digby Neck to expand, the growth rate of the Palmaria palmata in the region must first be determined as the combination of harvest frequency and fraction, and total optimal yield are dependent on rmax . Our model is based on quadratic growth of P. palmata from March to September inclusively, therefore, in the spring and fall, the maximal growth rate is less than in mid-summer. We have demonstrated that at low maximal growth rates, the current harvest yields are slightly less than the optimal yields. During mid-summer, however, when growth rates are higher, harvest yields are much less than the optimal yields. Our model predicts that in order to maximize yields when maximal growth rates are high, harvesting should occur at very high frequencies, and only small fractions of the standing population should be removed. However, this is unrealistic based on the time constraints of harvesting at low tides, and the 7-day duration of the spring tides. This could be compensated for by adding harvesters, but, if there are more people, each removing a small quantity, it may not be profitable for the individuals. However, our model also showed that larger fractions, taken at only moderately higher frequencies than exist currently, result in yields not significantly lower than the optimal yield.
- Management of intertidal seaweed resources in Nova Scotia
The two intertidal seaweeds with commercial harvests in Nova Scotia and New Brunswick are Ascophyllum nodosum and P. palmata . These seaweeds have radically different patterns of government regulation and resource management. The A. nodosum industry is highly managed with 90% of the harvest licences granted to Acadian Seaplants Ltd., and the resource base is carefully monitored to ensure that only 17 and 25% of biomass is harvested annually in New Brunswick and Nova Scotia, respectively (Ugarte et al. 2009). The Dulse harvest is unregulated, without harvest licences or official monitoring. Given the extreme difference in management strategies, it is of interest that both resources appear to be equally sustainable. This is especially surprising for P. palmata given the paucity of study of the regional ecology in the context of resource management. Indeed, the current study, and the related paper by Garbary et al. (2012) are the first detailed studies of Dulse ecology in Nova Scotia, and the first to document the state of the resource rather than simply harvest amounts (e.g., MacFarlane 1966, Ffrench 1974) or growth rates based on tank cultivation (e.g., Morgan et al. 1980, Morgan and Simpson 1981 a , 1981 b , Corey et al. 2011).
We thank Wanda and Calvin Vantassel of Fundy Dulse Ltd. for their generosity and logistic support during this research, and Eun Ju Kang and Katelyn White for assistance with preparation of figures. Grants from the University Council for Research of St. Francis Xavier University and the Nova Scotia Department of Fisheries and Aquaculture to DG, and the Natural Sciences and Engineering Research Council of Canada to DG and RL are gratefully acknowledged.
?berg P. 1992a A demographic study of two populations of the seaweed Ascophyllum nodosum. Ecology 73 1473 - 1487    DOI : 10.2307/1940691
?berg P. 1992b Size-based demography of the seaweed Ascophyllum nodosum in stochastic environments. Ecology 73 1488 - 1501    DOI : 10.2307/1940692
Ang P. O. , De Wreede R. E. Jr. 1990 Matrix models for algal life history stages. Mar. Ecol. Progr. Ser. 59 171 - 181    DOI : 10.3354/meps059171
Baardseth E. 1955 Regrowth of Ascophyllum nodosum after harversting. Institute for Industrial Research and Standards Dublin 67 -
Bird C. J. , Van der Meer J. P. 1993 Systematics of economically important marine algae: a Canadian perspective. Can. J. Bot. 71 361 - 369    DOI : 10.1139/b93-040
Chambers P. A. , DeWreede R. E. , Irlandi E. A. , Vandermeulen H. 1999 Management issues in aquatic macrophyte ecology: a Canadian perspective. Can. J. Bot. 77 471 - 487
Chopin T. , Ugarte R. 2006 The seaweed resources of eastern Canada. In Critchley, M. & Largo, D. B. (Eds.) World Seaweed References. An Authoritative Reference System. DVD ROM. ETI Bioinformatics Amsterdam.
Corey P. , Kim J. K. , Garbary D. J. , Prithiviraj B. , Duston J. 2011 Bioremdiation potential of Chondrus crispus (Basin Head) and Palmaria palmata: effect of temperature and high nitrate on nutrient removal. J. Appl. Phycol.    DOI : 10.1007/s10811-011-9734-8
Cornish M. L. , Garbary D. J. 2010 Antioxidants from macroalgae: potential applications in human health and nutrition. Algae 25 155 - 171    DOI : 10.4490/algae.2010.25.4.155
Duarte P. , Ferreira J. G. 1997 A model for the simulation of macroalgal population dynamics and productivity. Ecol. Model. 98 199 - 214    DOI : 10.1016/S0304-3800(96)01915-1
Edelstein T. , Chen L. , McLachlan J. 1970 Investigations of the marine algae of Nova Scotia. VIII. The flora of Digby Neck Peninsula, Bay of Fundy. Can. J. Bot. 48 621 - 629    DOI : 10.1139/b70-086
Edelstein-Keshet L. 1988 Mathematical models in biology. Random House New York, NY 586 -
Faes V. A. , Viejo R. M. 2003 Structure and dynamics of a population of Palmaria palmata (Rhodophyta) in Northern Spain. J. Phycol. 39 1038 - 1049    DOI : 10.1111/j.0022-3646.2003.02-142.x
Ffrench R. A. 1974 Rhodymenia palmata: an appraisal of the Dulse industry. Laboratory technical report. Atlantic Regional Laboratory, National Research Council of Canada Ottawa 49 -
Gabrielson P. W. , Garbary D. 1986 Systematics of red algae (Rhodophyta). CRC Crit. Rev. Plant Sci. 3 325 - 366    DOI : 10.1080/07352688609382215
Garbary D. 1976 Life-forms of algae and their distribution. Bot. Mar. 19 97 - 106    DOI : 10.1515/botm.1976.19.2.97
Garbary D. J. , Beveridge L. F. , Flynn A. D. , White K. L. 2012 Population ecology of Palmaria palmata (Palmariales) from harvested and non-harvested beds on Digby Neck, Nova Scotia. Algae 27 33 - 42    DOI : 10.4490/algae.2012.27.1.033
Lee C. S. , Ang P. Jr. 1991 A simple model for seaweed growth and optimal harvesting strategy. Ecol. Model. 55 67 - 74    DOI : 10.1016/0304-3800(91)90065-9
MacFarlane C. I. 1964 The seaweed industry of the Maritime Provinces. In De Virville, A. D. & Feldman, J. (Eds.) Proc. 4th Int. Seaweed Symp. Pergamon Press Oxford 414 - 419
MacFarlane C. I. 1966 A report on some aspects of the seaweed industry in the Maritime Provinces of Canada. Industrial Development Service, Department of Fisheries of Canada Ottawa 24 -
Martínez B. , Viejo R. M. , Rico J. M. , Rødde R. H. , Faes V. A. , Oliveros J. , Álvarez D. 2006 Open sea cultivation of Palmaria palmata (Rhodophyta) on the northern Spanish coast. Aquaculture 254 376 - 387    DOI : 10.1016/j.aquaculture.2005.10.025
Morgan K. C. , Shacklock P. F. , Simpson F. J. 1980 Some aspects of the culture of Palmaria palmata in greenhouse tanks. Bot. Mar. 23 765 - 770
Morgan K. C. , Simpson F. J. 1981a The cultivation of Palmaria palmata: effect of light intensity and nitrate supply on growth and chemical composition. Bot. Mar. 24 273 - 277    DOI : 10.1515/botm.1981.24.5.273
Morgan K.C. , Simpson F. J. 1981b The cultivation of Palmaria palmata: effect of light intensity and temperature on growth and chemical composition. Bot. Mar. 24 547 - 552    DOI : 10.1515/botm.1981.24.10.547
Nyman M. A. , Brown M. T. , Neushul M. , Keogh J. A. 1990 Macrocystis pyrifera in New Zealand: testing two mathematical models for whole plant growth. J. Appl. Phycol. 2 249 - 257    DOI : 10.1007/BF02179782
Pang S. , L?ning K. 2004 Tank cultivation of the red alga Palmaria palmata: effects of intermittent light on growth rate, yield and growth kinetics. J. Appl. Phycol. 16 93 - 99    DOI : 10.1023/B:JAPH.0000044779.30182.d8
Rhatigan P. 2009 Irish seaweed kitchen: the comprehensive guide to healthy everyday cooking with seaweeds. Booklink Hollywood, LA 288 -
Santos R. , Nyman M. 1998 Population modelling of Gelidium sesquipedale (Rhodophyta, Gelidiales). J. Appl. Phycol. 10 261 - 272    DOI : 10.1023/A:1008083926349
Scrosati R. 2004 A discrete-time logistic model of frond dynamics for Mazzaella parksii (Rhodophyta, Gigartinales). J. Appl. Phycol. 16 69 - 72    DOI : 10.1023/B:JAPH.0000019072.58963.4c
Seip K. L. 1980 A computational model for growth and harvesting of the marine alga Ascophyllum nodosum. Ecol. Model. 8 189 - 199    DOI : 10.1016/0304-3800(80)90037-X
Sharp G. , Bodiguel C. 1957 Introducing integrated management, ecosystem and precautionary approaches in seaweed management: the Ascophyllum nodusum (rockweed) harvest in New Brunswick, Canada and implications for industry. In Chapman, A. R. O., Anderson, R. J., Vreeland, V. J. & Davison, I. (Eds.) Oxford University Press New York Proc. 17th Int. Seaweed Symp. 107 - 114
Smith B. D. 1986 Implications of population dynamics and interspecific competition for harvest management of the seaweed Laminaria. Mar. Ecol. Prog. Ser. 33 7 - 18    DOI : 10.3354/meps033007
South G. R. , Tittley I. , Farnham W. F. , Keats D. W. 1988 A survey of the benthic marine algae of southwestern New Brunswick, Canada. Rhodora 90 419 - 451
Taylor W. R. Marine algae of the northeastern coast of North America University of Michigan Press Ann Arbor, MI 509 -
Tsoularis A. , Wallace J. 2002 Analysis of logistic growth models. Math. Biosci. 179 21 - 55    DOI : 10.1016/S0025-5564(02)00096-2
Ugarte R. A. , Critchley A. , Serdynska A. R. , Deveau J. P. 2009 Changes in composition of rockweed (Ascophyllum nodosum) beds due to possible recent increase in sea temperature in eastern Canada. J. Appl. Phycol. 21 591 - 598    DOI : 10.1007/s10811-008-9397-2
Vadas R. L. , Beal B. F. , Wright W. A. , Emerson S. , Nickl S. 2004 Biomass and productivity of red and green algae in Cobscook Bay, Maine. Northeast. Nat. 11 (Special issue 2) 163 - 196
Van der Meer J. P. , Todd E. R. 1980 The life history of Palmaria palmata in culture: a new type for the Rhodophyta. Can. J. Bot. 58 1250 - 1256    DOI : 10.1139/b80-155
Verhulst P. -F. 1838 Notice sur la loi que la population poursuit dans son accroissement. Corresp. Math. Phys. 10 113 - 121
Wilson J. S. , Bird C. J. , McLachlan J. , Taylor A. R. A. 1979 An annotated checklist and distribution of benthic marine algae of the Bay of Fundy. Meml. Univ. Newfoundland Occas. Pap. Biol. 2 1 - 65
Yuan Y. V. , Westcott N. D. , Hu C. , Kitts D. D. 2009 Mycosporine-like amino acid composition of the edible red alga, Palmaria palmata (Dulse) harvested from the west and east coasts of Grand Manan Island, New Brunswick. Food Chem. 112 321 - 328    DOI : 10.1016/j.foodchem.2008.05.066
Zhang X. , Shuai Z. , Wang K. 2003 Optimal impulsive harvesting policy for single population. Nonlinear Anal. Real World Appl. 4 639 - 651    DOI : 10.1016/S1468-1218(02)00084-6