Bayesian inference of environmental effects on seaweed production in Japan via a production-environmental suitability model (2024)

As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsem*nt of, or agreement with, the contents by NLM or the National Institutes of Health.
Learn more: PMC Disclaimer | PMC Copyright Notice

Bot Stud. 2019 Dec; 60: 2.

Published online 2019 Feb 1. doi:10.1186/s40529-018-0250-x

PMCID: PMC6358630

PMID: 30707346

Hungyen Chen

Author information Article notes Copyright and License information PMC Disclaimer

Associated Data

Data Availability Statement

Abstract

Background

Both natural and human-induced disturbances are commonly responsible for an overall decrease of the world’s seaweed. Along Japan’s coastal areas, edible seaweed production has been decreasing for decades. In this study, a production-environmental suitability model to estimate the impacts of environmental factors on seaweed production was developed. The developed model not only estimates human-induced disturbances but also quantifies the impacts of environmental factors responsible for the decline of annual seaweed production. The model estimated the temporal variation in human-induced disturbances and the effects of environmental factors (i.e., rainfall, CO₂ concentrations, temperature, typhoons, solar radiation, water nutrient levels, and water quality) on edible seaweeds in Japan from 1985 to 2012.

Results

The environmental suitability for seaweed production in Japan was about 4.6 times greater in 1992 than in 2011, meanwhile as a result of human activities, human-induced disturbances of seaweed increased at a rate of 4.9 times faster during the period of 1998–2012 than the period of 1985–1997. The ratio of decreased production to decreased environmental suitability for seaweed production in Japan increased by 15.2% during the study years, which means that seaweed production has become more sensitive to environmental disturbances, including climatic factors and human activities in 1998–2012.

Conclusions

The results are novel in demonstrating temporal variations in the level of environmental suitability to seaweed production by using a simple mathematical model. The production-environmental suitability model successfully predicted seaweed production by reflecting the 28-year temporal variation of the observed seaweed production in Japan.

Keywords: Bayesian estimation, Climatic change, Seaweed growth, Yield prediction

Background

Seaweed, bathed in nutrient rich seawater, is a family of algae and is one of the most important aquatic living resources of the world’s oceans (Dhargalkar and Pereira 2005). Since 300 BC, seaweeds have been a source of food, feed, fertilizer, chemicals, medicine, and storm protection (Chapman 1970; Verkleij 1992; Smit 2004; Rönnbäck et al. 2007). The majority of the world’s seaweed cultivation, production, and consumption occurs particularly in East and Southeast Asian, e.g., China and Japan (Dhargalkar and Pereira 2005).

Seaweeds play a central ecological role in coastal habitats; as ecosystem engineers, they assist in supplying oxygen and absorbing CO₂ to the sea (Hemminga and Duarte 2000; Graham 2004). Seaweeds act as one of the primary producers and competitors in the marine food chain (Graham 2004; Norderhaug et al. 2005), with some having the capacity to remove heavy metals from the seawater and potentially being used as bio-indicators for heavy metal pollution (Johansen et al. 1991; Phillips and Segar 1986; Campanella et al. 2001; Serfor-Armah et al. 2006) and can remove nutrients from the water column and inhibit eutrophication (Harlin and Thorne-Miller 1981). They also provide an important niche, e.g., create habitat for living, eating, parenting, and hiding, for plenty of marine species, and thus maintain community diversity (Stachowicz et al. 2008). However, although seaweed in general has good survival strategies, enabling the algae to withstand environmental stresses, the impacts of climatic change and human-induced disturbances on seaweeds remain poorly understood (Chan et al. 2006).

Seaweeds depend on an adequate degree of water quality and clarity to sustain productivity in their growth environment (Lüning 1990; Lobban and Harrison 1994). Seaweed loss, however, is increasing due to (1) natural disturbances, such as hurricanes, earthquakes, disease, and grazing by herbivores; and (2) human activities, including nutrient and sediment loading from runoff and sewage disposal, dredging and filling, pollution, upland development, and certain fishing practices (Sousa 1979; Chapman and Johnson 1990; Rainbow 1995; Short and Burdick 1996; Airoldi 1998; Coelho et al. 2000). The loss of seaweed influences the biodiversity of the ecosystem’s communities (Graham 2004; Norderhaug et al. 2005) and should be of concern to society.

Japan is over 3000km long from the Okhotsk Sea to the Philippine Sea in the Pacific Ocean and has 29,751km long coastline, the sixth longest in the world (MHLW 2018). As an island nation, Japan harvests numerous types of seaweed from all around the country (MAFF 2018). Seaweed farming began in Japan as early as 1670 in Tokyo Bay and seaweeds and seaweed-derived products have been central ingredients of Japanese cuisine for centuries (Borgese 1980). In 2012, Japan dedicated about 100,000 tonnes to seaweed production around the country (MAFF 2018).

Previous models for seaweed have used to simulate the growth and optimal annual harvesting strategy (Lee and Ang 1991), study the relationships between environmental factors and annual growth pattern (Friedlander et al. 1990), and explore the critical environmental parameters influencing the growth and establishment (Murphy et al. 2016). Cabral et al. (2016) proposed an approach to present an assessment of the potential impact of the installation of seaweed farms on ecosystem services and the induced compensation costs. Préat et al. (2018) used a method to model the seaweed production and develop indicators to assess the potential reductions in fisheries yield due to seaweed farming.

The crop yield-fertility model, modified from Michaelis–Menten kinetics (Michaelis and Menten 1913), was developed to estimate the soil fertility level and to quantify the contributions of fertilizer to the improvements in crop yields (Chen et al. 2014). It was used to interpret the effects of soil fertility and fertilizers on the crop yield and the extent of saturation of additional inputs. The model can be used to predict crop yields given the total fertility level of the soil and to describe the variation in crop yield and the effect of fertilizer treatment over years for crops using the estimated soil fertility. Bayesian inference of environmental effects on seaweed production.

In this study, a production-environmental suitability model, which was modified from the crop yield-fertility model, to estimate the effect of environmental suitability on seaweed production using Bayesian algorithm was developed. One of the major advantages of the Bayesian approach is to incorporate prior information (Box and Tiao 1973). The Bayesian approach forces the analyst to canvass expert knowledge by setting priors to determine what is known about the biological or environmental parameters and processes (Punt and Hilborn 1997). The developed model not only estimates human-induced disturbance but also quantifies the disturbance impacts of environmental factors (i.e., rainfall, CO₂ concentrations, temperature, typhoons, solar radiation, water nutrient levels, and water quality) contributing to the decline of annual seaweed production in Japan. This paper defines human-induced disturbance as the combination of pollution, upland development, fishing practices, and other human activities that disrupt seaweed production. I used the model to estimate the annual variation in human-induced disturbances and the effects of environmental factors on edible seaweeds along the coastal areas around Japan from 1985 to 2012.

Methods

Production-environmental suitability model

Michaelis–Menten kinetics (Michaelis and Menten 1913) and the crop yield-fertility model (Chen et al. 2014) were taken and modified to produce this study’s seaweed model.

$\begin{matrix} Y_{t} & ≅ \frac{V}{1 + \frac{K}{E S_{t}}} \\ E S_{t} & = & \frac{1}{H D_{t}} + \frac{a}{R_{t}} + \frac{b}{C_{t}} + \frac{c}{T_{t}} + \frac{d}{W_{t}} + e \times S_{t} + f \times N_{t} + g \times Q_{t} \end{matrix}$

The model was used to interpret the effects of human-induced disturbances, rainfall, CO₂ concentrations, temperature, typhoons, solar radiation, water nutrient levels, and water quality on seaweed production and the extent of saturation of additional inputs of environmental suitability for seaweed production. The model related productions (Y_t) to the environmental suitability (ES_t) of the seaweed at time t. V is the maximum production in response to the maximum environmental suitability; K is the suitability level before disturbances of human and environmental factors that is required to produce half of the maximum production. Human-induced disturbance (HD_t) affecting the seaweed at time t is assumed to gradually vary over time. As previously mentioned, this paper defines human-induced disturbance as the combination of pollution, upland development, fishing practices, and other human activities. In other words, human-induced disturbance represents all the factors which relate to human activities and may influence the seaweed production. R_t is the value of rainfall at time t, C_t is the value of CO₂ concentrations at time t, T_t is the value of sea surface temperature at time t, W_t is the value of number of typhoons at time t, S_t is the value of solar radiation at time t, N_t is the value of water nutrient levels at time t and Q_t is the value of water quality at time t. a, b, c, d, e, f, and g represent the effects of R_t, C_t, T_t, W_t, S_t, N_t, and Q_t relative to the seaweed production, respectively. The values of R_t, C_t, T_t, W_t, S_t, N_t, and Q_t were transferred to the relative value to the mean value (from 1985 to 2012). In the model, rainfall, CO₂ concentrations, sea surface temperatures, and typhoons were assumed to be the factors of disturbance for seaweed production, decreasing the environmental suitability for seaweed production. On the other hand, solar radiation, water nutrient levels, and water quality increased the environmental suitability for seaweed.

Likelihood and priors

A Bayesian framework was adopted for parameter estimation. Here, a gradual change in the human disturbance was assumed for analysis. However, we have to keep in mind that human disturbance is sometimes more often of abrupt or exponential change. The likelihood of the production at time t followed a normal distribution with the mean $\frac{V}{1 + \frac{K}{\frac{1}{H D_{t}} + \frac{a}{R_{t}} + \frac{b}{C_{t}} + \frac{c}{T_{t}} + \frac{d}{W_{t}} + e \times S_{t} + f \times N_{t} + g \times Q_{t}}}$ and the variance δ:

$Y_{t} \sim N (\frac{V}{1 + \frac{K}{\frac{1}{H D_{t}} + \frac{a}{R_{t}} + \frac{b}{C_{t}} + \frac{c}{T_{t}} + \frac{d}{W_{t}} + e \times S_{t} + f \times N_{t} + g \times Q_{t}}}, δ) .$

Setting the HD in the first year to 1 normalized this value. The smoothness priors of the HD from the second year followed a normal distribution with the mean equal to the value of the HD in the preceding year and the variance τ:

$H D_{t} \sim N (H D_{t - 1}, τ), t > 1 .$

The inverse of δ followed a gamma distribution with a shape parameter of 0.1 and a scale parameter of 10. The inverse of τ followed a gamma distribution with a shape parameter of 10 and a scale parameter of 0.1. The prior of V followed a normal distribution with a mean of 0 and a standard deviation of 10⁴. The priors of K, a, b, c, d, e, f, and g followed a gamma distribution with a shape parameter of 1 and a scale parameter of 1. The priors of the estimates were designed to be as non-informative as possible within a realistic range of the parameter values. Runs of estimates in the model were performed using a Bayesian Markov chain Monte Carlo (MCMC) method as implemented in the software program WinBUGS (Lunn et al. 2000) and ran from R (R Development Core Team 2013) using the package R2WinBUGS (Sturtz et al. 2005). The number of iteration of the MCMC operation was set to 1,000,000 and thin was set to 100. All calculations and data analyses were performed using R v 3.0.2 (R Development Core Team 2013).

Seaweed production in Japan

The data of yearly edible seaweed (i.e., kelp from the family Laminariaceae, Undaria pinnatifida, Sargassum fusiforme, and other edible seaweeds not harvested from aquiculture) production in the coastal areas around Japan from 1985 to 2012 were downloaded from the official website of the Ministry of Agriculture, Forestry and Fisheries (MAFF 2018). In this study, the fishery records were assumed to reflect the seaweed production. The size of growth area for seaweed can also respond to environmental change well, however it may be difficult to measure and evaluate annually for a long period around Japan, especially the records for the early years could be hard to obtain.

Data of environmental factors

Annual data for precipitation, CO₂ concentrations in the seawater, sea surface temperatures, number of typhoons, and solar radiation in Japan were downloaded from the Japan Meteorological Agency official website (available online from http://www.jma.go.jp/jma/menu/menureport.html). Data for precipitation was collected using annual values from around Japan. To measure CO₂ concentrations, data were collected using values along coastal areas of central Japan (137ºE, 34ºN) in winter. Sea surface temperatures were collected using the annual mean in the seas around Japan. The number of recorded typhoons was based on the number of occurrences around Japan. Solar radiation data was based on the annual value of sunlight duration around Japan. Seawater nutrient levels and seawater quality were downloaded from the Japanese Ministry of Land, Infrastructure, Transport and Tourism official website (available online from http://www.pa.cgr.mlit.go.jp/chiki/suish*tu/). The data for seawater nutrient levels and seawater quality were collected using the annual value in the Seto Inland Sea in Japan, with nutrient levels calculated as the sum of concentrations of discharged loads of total phosphorus (TP), total nitrogen (TN), dissolved inorganic phosphate (PO₄), and nitrogen (NO₃, NO₂ and NH₄). Seawater quality was assessed according to water transparency values. To integrate the collected data, the value in each year was transferred to the ratio of the mean value among the 28years for each environmental factor (Fig.1). The mean values during the growth season of environmental parameters are reasonable for the analyses of growth in plants and algae. The annual means were used here because this analysis included different kinds of seaweeds at large-ranged areas (latitude) with variable suitable growth seasons and it may be difficult to select an appropriate period for growth season. Although annual mean may not reflect the exact situation for the seaweed growth, it could reflect the annual variation among different years in long-term temporal analyses.

Bayesian inference of environmental effects on seaweed production in Japan via a production-environmental suitability model (3)

Open in a separate window

Fig.1

Temporal variation of the environmental factors in Japan from 1985 to 2012. a Rainfall. b Sea water CO₂ concentration. c Sea surface temperature. d Number of typhoon. e Solar radiation. f Sea water nutrient level. g Sea water quality. The value in each year was transferred to the ratio of the mean value among the 28years for each environmental factor. See “Methods” for details

Results

Decreasing trend of seaweed production

According to Fig.2, the annual temporal variation in seaweed production in Japan ranged between 87,779 (in 2011) and 209,141 (in 1992) tonnes with a mean of 139,119 and standard deviation 36,808 during the 28years from 1985 to 2012. A decreasing trend of temporal variation was found by using simple linear regression of the production in response to the year during 1985–2012 (slope = − 039; p value < 0.001).

Bayesian inference of environmental effects on seaweed production in Japan via a production-environmental suitability model (4)

Open in a separate window

Fig.2

Temporal variations in observed seaweed production (1000 tonnes) in Japan from 1985 to 2012. The gray line represents the simple linear regression line

Maximum production and environmental suitability

Table1 shows the results of the Bayesian estimates (V, K, a, b, c, d, e, f, and g) of the model for seaweed production in Japan. The traces of the MCMC samples show good mixing and convergence to the posterior distributions (Fig.3). The largest observed seaweed production in Japan for 1year was 209,141 (1992), which is 68.7% of the estimated maximum production (V, 304,508 tonnes). The environmental suitability (ES_t) for each year was estimated by using environmental variables (R_t, C_t, T_t, W_t, S_t, N_t, and Q_t) given the estimates HD_t, a, b, c, d, e, f, and g. The estimated environmental suitability for seaweed in Japan had a mean of 0.925 and standard deviation of 0.518 during the 28years. The estimated environmental suitability ranged between 0.459 (in 2011) and 2.148 (in 1992), which are 0.46 × half-saturated suitability (K, 1.008) and 2.13 × half-saturated suitability for seaweed in Japan. This means that the environmental suitability for seaweed in Japan decreased from more than 2 times (2.13 × K) to about half of half-saturated suitability (0.46 × K) over the years. Therefore, the environmental suitability for seaweed production in Japan was about 4.6 times (2.13/0.46) greater in 1992 than in 2011. The environmental suitability was first estimated below half-saturated suitability in 1994 (0.833) and has never reached half-saturated suitability again since then.

Table1

The posterior mean and standard deviation (SD) of the Bayesian estimates

Factor	Estimate	Mean	SD	Percentage of contribution to all the environmental factors (%)
Maximum production (1000 tonnes)	V	304.508	43.085	–
Half-saturated suitability	K	1.008	0.421	–
Rainfall	a	0.054	0.053	13.2
CO₂ concentration	b	0.065	0.072	15.9
Temperature	c	0.055	0.065	13.5
Typhoon	d	0.018	0.022	4.4
Solar radiation	e	0.053	0.060	13.0
Nutrient level	f	0.137	0.119	33.6
Water quality	g	0.026	0.039	6.4

Open in a separate window

Effects of environmental factors

The Bayesian estimates for the contribution of environmental factors to seaweed and their percentages of contribution in the studied period are showed in Table1. The impact from nutrient levels (f, 0.137; 33.6%) was the largest among the seven environmental factors; the second largest was CO₂ concentrations (b, 0.065; 15.9%); third was temperature (c, 0.055; 13.5%); fourth was rainfall (a, 0.054; 13.2%); solar radiation (e, 0.053; 13.0%) was fifth; second last was water quality (g, 0.026; 6.4%); and last was number of typhoons (d, 0.018; 4.4%). This result suggests that nutrient levels may play the most influential role, of the seven environmental factors, regarding seaweed production in Japan. In fact, the impact was twice as large as seawater CO₂ concentrations, the second largest impact. Temperature, rainfall, and solar radiation impacts were all similar in value. Water quality had the smallest impact among the selected environmental factors, which had an effect of about one-fifth the significance of seawater nutrient levels.

Increasing trend of human-induced disturbance

The estimated human-induced disturbance (HD_t) to seaweed in Japan was highly correlated with year (correlation coefficient = 0.957; p-value < 0.001) and has been increasing over the 28-year period (slope of simple linear regression = 0.368; p-value < 0.001). Setting the initial year (1985) to 1 normalized the value of human-induced disturbance. Figure4 shows the long-term temporal variation of the human-induced disturbance to seaweed and its standard deviation band in Japan.

Bayesian inference of environmental effects on seaweed production in Japan via a production-environmental suitability model (6)

Open in a separate window

Fig.4

Temporal variations in the posterior mean and the gray band that corresponds to the standard deviation (± SD) of the Bayesian estimates of human-induced disturbance to seaweed in Japan from 1985 to 2012. The straight lines represent the simple linear regression lines. The value in the initial year was set to 1

The trend of human-induced disturbance has been increasing since 1998 (Fig.4). The correlation coefficient for human-induced disturbance to year was 0.636 (p-value = 0.02), and slope of simple linear regression was 0.099 (p-value = 0.02) during 1985–1997 (Fig.4). However, the slope of simple linear regression became 0.483 (p-value < 0.001) with a correlation coefficient, 0.980 (p-value < 0.001) during 1998–2012 (Fig.4). The results suggest a much faster rate of increase of the human-induced disturbance to seaweed in Japan during the period of 1998–2012 than the period of 1985–1997 (about 4.9 times faster).

Predicted crop yields reflected the trend of the observed yields

Using the estimated maximum production, half-saturated suitability and environmental suitability, the range of predicted production was estimated by the production-environmental suitability model for seaweed in Japan from 1985 to 2012. The estimated environmental suitability was transformed to the seaweed production by the production-environmental suitability curve (Fig.5). The predicted seaweed production had a range of 95,321–207,241 tonnes with a mean of 136,666 and standard deviation of 36,031 during the 28years. The predicted production was highly correlated with the observed production for seaweed in Japan (correlation coefficient = 0.995; p-value < 0.001; Fig.6).

Open in a separate window

Fig.5

Production-environmental suitability curve for seaweed in Japan. The gray band represents the band that corresponds to the standard deviation (± 0.5 SD) of the curve. V, maximum production; K, half-saturated suitability

Open in a separate window

Fig.6

The relationship between the observed and predicted seaweed productions in Japan. The gray line represents diagonal 1:1 line

The environmental suitability was estimated at 2.103 in 1990, 0.875 in 1997, 0.555 in 2004, and 0.459 in 2011, yielding 205,857, 141,475, 108,156, and 95,321 tonnes of seaweed, respectively (Fig.5). The environmental suitability decreased by 1.229, then 0.319, and then 0.096 in these three 7-year periods (1990–1997, 1997–2004, 2004–2011), respectively, which resulted in an overall decrease of 64,382, 33,319, and 12,835 tonnes of seaweed production, respectively (Table2). Compared with the value in the initial year of each period, the environmental suitability decreased by 58.4%, 36.5%, and 17.3%, respectively, resulting in a decrease in seaweed production by 31.3%, 23.6%, and 11.9%, respectively (Table2). The ratios of the decreased production to the decreased environmental suitability were 53.5%, 64.5%, and 68.7% for these three periods, respectively (Table2). The ratio of the decreased production to the decreased environmental suitability for the seaweed production in Japan increased by 15.2% (68.7%–53.5%) during the years. Figure5 shows the increasing ratio of the decreased production to the decreased environmental suitability. This result suggests that in 1998–2012, seaweed production in Japan has become more and more sensitive to environmental disturbances, including climatic factors and human activities.

Table2

Environmental suitability and production for seaweed in Japan estimated using the production-environmental suitability model

Year	Environmental suitability			Production			Ratio of decreased percentage of production to decreased percentage of environmental suitability (%)
Year	Estimated value	Decreased value	Decreased percentage (%)	Estimated value (ton)	Decreased value (ton)	Decreased percentage (%)
1990	2.103	–	–	205,857	–	–	–
1997	0.875	1.229	58.4	141,475	64,382	31.3	53.5
2004	0.555	0.319	36.5	108,156	33,319	23.6	64.5
2011	0.459	0.096	17.3	95,321	12,835	11.9	68.7

Open in a separate window

See the production-environmental suitability curve in Fig.5

Discussion

Both large-scale and local losses of seaweed have been frequently reported around the world for decades (Short and Burdick 1996; Coelho et al. 2000; Heck et al. 2003; Orth et al. 2006; Airoldi and Beck 2007). Natural events, such as extreme climatic events, have been identified as causes of seaweed loss in temperate and tropical regions (Orth et al. 2006). Reported anthropogenic threats to seaweed health include global warming, shifts in water quality, and increased loading of sediment, contaminants, and nutrients (Kemp et al. 2005; Orth et al. 2006; Walker et al. 2006). Extreme climatic events (e.g., frequency and intensity of storms, hurricanes, typhoons, and associated surges and swells) also can have large-scale impacts on seaweed (Orth et al. 2006). However, Short and Wyllie-Echeverria’s (1996) findings suggest that, although natural events have been responsible for loss of seaweed habitat, the most serious cause of seaweed loss is attributable to human population expansion and increasing anthropogenic inputs to the coastal oceans.

The environmental factors analyzed in the study were selected because: the amount of solar radiation are related to seaweed communities (Bischof et al. 2006); the degrees of water temperature and CO₂ concentrations are positively correlated and they reflect the degrees of global warming and ocean acidification which affect seaweed communities (Pörtner 2008; Harley et al. 2012); the amount of precipitation and the number of typhoons can be considered as the degree of environmental disturbances (Short and Wyllie-Echeverria 1996); water nutrient levels and water quality are related to the growth of seaweed (Msuya and Neori 2008). Among these selected environmental factors, the results suggest nutrient levels in seawater may have the largest impact on seaweed production in Japan. Various nutrients in different concentrations come from industrial pollutants, agricultural runoff, municipal sewage systems, and other human activities. Seawater CO₂ concentrations are also significant environmental factors, suggesting that pH value could be another important factor in the production of seaweed in Japan and could also be correlated with the human activities.

Generally, the reasons behind the global reduction of seaweed could be due to changing sea currents, sea temperatures, nutrient levels, and the abundance of seaweed-eating species, such as fish, shellfish, and sea urchins (Burkepile and Hay 2006; Filbee-dexter and Scheibling 2014; Ling et al. 2014). At the local scale, the loss of seaweed could be due to the varied coastal topography, seafloors, seawater transparency, and eutrophication (Lüning 1990; Schramm 1999; Walker et al. 2006; Yang et al. 2015). The loss of seaweed has been a serious problem in Japan since early 2000s (Fujita et al. 2006; Fujita 2009). Besides the reasons mentioned above, these losses are most likely due to a combination of many natural factors such as increased sea surface temperatures, decreased nutrient levels, consumption by herbivores, inflow of fresh water from rivers, extreme climatic events, and human-induced factors such as coastal destruction, industrial and agricultural pollution, drainage from power plants, fisheries, aquaculture, and overfishing (Nakatsu 2005; Fujita 2009; Yamaguchi 2010; Fisheries Agency 2015). The increasing rate of human-induced disturbance to seaweed production suggested by the results of this study may be due to a high increase in human activities such as, pollution, upland development, and fishing practices (Short and Wyllie-Echeverria 1996; Dulvy et al. 2003; Crain et al. 2009). The other serious problem in some of the shallow seas along the Japanese coast have been suffering from a phenomenon known as barren ground, sea desert or denudation of rocks (Fujita 2010). This phenomenon was known and observed in Japan, however it may be difficult to measure and evaluate annually for a long period and a large-scaled area.

The objectives of this study are to estimate the human-induced disturbance and to quantify the disturbance impacts of environmental factors contributing to the decline of annual seaweed loss in Japan by using a mathematical model. However, one point concerns the temporal data of seaweed production may not reflect the exact temporal variation of seaweed abundance. The quality and reliability of this analysis could be improved by including the information of the efforts involved to get the seaweed production, market demand of seaweed, and the number of fisherman every year. Because it is very difficult to get this information for the data of large-scaled investigation in several years, we may assume these variables were consistent during the studying period. We have to keep this point in mind and interpret the results carefully. The model used in this study can be applied at other areas in the world, as long as the information of environmental variables is provided. In addition, the model can be revised easily depending on what kind of information or type of data we have.

Conclusions

By assuming that the fishery records can reflect the seaweed production, the results succeeded in demonstrating a mathematical method to estimate the increasing effect of human-induced disturbance to the seaweed community. The environmental suitability for seaweed in Japan was about 4 times greater in 1992 than in 2011. The proposed model was used to estimate the size of the effects of different environmental factors. The results show that human activity related factors, i.e., seawater nutrient levels and CO₂ concentrations, had larger impacts on the variation of seaweed production in Japan. Human-induced disturbance of seaweed increased at a rate 4.9 times faster during the period of 1998–2012 than the period of 1985–1997. The ratio of the decreased production to the decreased environmental suitability for seaweed production in Japan increased by 15.2% during the years, which means that the seaweed production has become more sensitive to environmental disturbances in 1998–2012, including climatic factors and human activities. The results are novel in demonstrating temporal variations in the level of environmental suitability to seaweed production by using a simple mathematical model. The production-environmental suitability model successfully predicted seaweed production by reflecting the 28-year temporal variation of the observed seaweed production in Japan. Thereby this study has provided a method to estimate the effects of human-induced disturbances and climatic changes on aquatic plants.

Authors’ contributions

HC designed and conducted the analyses and wrote the manuscript. The author read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The author declares no competing interests.

Availability of data and materials

Requests for data and materials should be addressed to HC.

Consent for publication

The author agrees to the terms of the Springer Open Copyright and License Agreement.

Ethics approval and consent to participate

Not applicable, the study involves no human participants.

Funding

This work was supported by funding from the Ministry of Science and Technology, Taiwan (Grant Number: MOST 107-2313-B-002-010-MY3).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Abbreviations

ES	environmental suitability
HD	human-induced disturbance

References

Agency Fisheries. Isoyake Taisaku Guideline Re-ed. Japan: Tokyo; 2015. [Google Scholar]
Airoldi L. Roles of disturbance, sediment stress, and substratum retention on spatial dominance in algal turf. Ecology. 1998;79:2759–2770. doi:10.1890/0012-9658(1998)079[2759:RODSSA]2.0.CO;2. [CrossRef] [Google Scholar]
Airoldi L, Beck MW. Loss, status and trends for coastal marine habitats of Europe. Oceanogr Mar Biol Annu Rev. 2007;45:345–405. [Google Scholar]
Bischof K, Gómez I, Molis M, Hanelt D, Karsten U, Lüder U, Roleda MY, Zacher K, Wiencke C. Ultraviolet radiation shapes seaweed communities. Rev Environ Sci Bio. 2006;5:141–166. doi:10.1007/s11157-006-0002-3. [CrossRef] [Google Scholar]
Borgese EM. Seafarm: the story of aquaculture. New York: Abrams HN Incorporated; 1980. [Google Scholar]
Box GEP, Tiao GC. Bayesian inference in statistical analysis. New York: Wiley; 1973. [Google Scholar]
Burkepile DE, Hay ME. Hervibore vs. nutrient control of marine primary producers: context-dependent effects. Ecology. 2006;87:4128. [PubMed] [Google Scholar]
Cabral P, Levrel H, Viard F, Frangoudes K, Girard S, Scemama P. Ecosystem services assessment and compensation costs for installing seaweed farms. Mar Policy. 2016;71:157–165. doi:10.1016/j.marpol.2016.05.031. [CrossRef] [Google Scholar]
Campanella L, Conti ME, Cubadda F, Sucapane C. Trace metals in seagrass, algae and mollusks from an uncontaminated area in the Mediterranean. Environ Pollut. 2001;111:117–126. doi:10.1016/S0269-7491(99)00327-9. [PubMed] [CrossRef] [Google Scholar]
Chan C-X, Ho C-L, Phang S-M. Trends in seaweed research. Trends Plant Sci. 2006;11:165–166. doi:10.1016/j.tplants.2006.02.003. [PubMed] [CrossRef] [Google Scholar]
Chapman VJ. Seaweeds and their uses. 2. London: The Camelot Press Ltd.; 1970. [Google Scholar]
Chapman ARO, Johnson CR. Disturbance and organization of macroalgal assemblages in the Northwestern Atlantic. Hydrobiologia. 1990;192:77–121. doi:10.1007/BF00006228. [CrossRef] [Google Scholar]
Chen H, Yamagishi J, Kishino H. Bayesian inference of baseline fertility and treatment effects via a crop yield-fertility model. PLoS ONE. 2014;9:e112785. doi:10.1371/journal.pone.0112785. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
Coelho SM, Rijstenbil JW, Brown MT. Impacts of anthropogenic stress on the early development stages of seaweeds. J Aquat Ecosys Stress Recov. 2000;7:317–333. doi:10.1023/A:1009916129009. [CrossRef] [Google Scholar]
Crain CM, Halpern BS, Beck MW. Understanding and managing human threats to the coastal marine environment. The year in ecology and conservation biology, 2009. Ann NY Acad Sci. 2009;1162:39–62. doi:10.1111/j.1749-6632.2009.04496.x. [PubMed] [CrossRef] [Google Scholar]
Dhargalkar VK, Pereira N. Seaweed: promising plant of the millennium. Sci Cul. 2005;71:60–66. [Google Scholar]
Dulvy NK, Sadovy Y, Reynolds JD. Extinction vulnerability in marine populations. Fish Fish. 2003;4:25–64. doi:10.1046/j.1467-2979.2003.00105.x. [CrossRef] [Google Scholar]
Filbee-dexter K, Scheibling RE. Sea urchin barrens as alternative stable states of collapsed kelp ecosystems. Mar Ecol Prog Ser. 2014;495:1–25. doi:10.3354/meps10573. [CrossRef] [Google Scholar]
Friedlander M, Galai N, Farbstein H. A model of seaweed growth in an outdoor culture in Israel. Hydrobiologia. 1990;204(205):367–373. doi:10.1007/BF00040258. [CrossRef] [Google Scholar]
Fujita D (2009) Decline and restoration of seaweed beds. Tokyo Fisheries Promotion Foundation
Fujita D. Current status and problems of Isoyake in Japan. Bull Fish Res Agen. 2010;32:33–42. [Google Scholar]
Fujita D, Ishikawa T, Kodama S, Kato Y, Notoya M. Distribution and recent reduction of Gelidium beds in Toyama Bay, Japan. J Appl Phycol. 2006;18:591–598. doi:10.1007/s10811-006-9060-8. [CrossRef] [Google Scholar]
Graham MH. Effects of local deforestation on the diversity and structure of southern California giant kelp forest food webs. Ecosystems. 2004;7:341–357. doi:10.1007/s10021-003-0245-6. [CrossRef] [Google Scholar]
Harley CDG, Anderson KM, Demes KW, Jorve JP, Kordas RL, Coyle TA, Graham MH. Effects of climate change on global seaweed communities. J Phycol. 2012;48:1064–1078. doi:10.1111/j.1529-8817.2012.01224.x. [PubMed] [CrossRef] [Google Scholar]
Harlin MM, Thorne-Miller B. Nutrient enrichment of seagrass beds in a Rhode Island coastal lagoon. Mar Biol. 1981;65:221–229. doi:10.1007/BF00397115. [CrossRef] [Google Scholar]
Heck KL, Jr, Hays G, Orth RJ. Critical evaluation of the nursery role hypothesis for seagrass meadows. Mar Ecol Prog Ser. 2003;253:123–136. doi:10.3354/meps253123. [CrossRef] [Google Scholar]
Hemminga M, Duarte CM. Seagrass Ecology. Cambridge: Cambridge University Press; 2000. [Google Scholar]
Johansen P, Hansen MM, Asmund G, Nielsen PB. Marine organisms as indicators of heavy metal pollution—experience from 16years of monitoring at a lead zinc mine in Greenland. Chem Ecol. 1991;5:35–55. doi:10.1080/02757549108035243. [CrossRef] [Google Scholar]
Kemp WM, Boynton WR, Adolf JE, Boesch DF, Boicourt WC, Brush G, Cornwell JC, Fisher TR, Glibert PM, Hagy JD, Harding LW, Houde ED, Kimmel DG, Miller WD, Newell RIE, Roman MR, Smith EM, Stevenson JC. Eutrophication of Chesapeake Bay: historical trends and ecological interactions. Mar Ecol Progr Ser. 2005;303:1–19. doi:10.3354/meps303001. [CrossRef] [Google Scholar]
Lee CS, Ang P., Jr A simple model for seaweed growth and optimal harvesting strategy. Ecol Model. 1991;55:67–74. doi:10.1016/0304-3800(91)90065-9. [CrossRef] [Google Scholar]
Ling SD, Scheibling RE, Rassweiler A, Johnson CR, Shears N, Connell SD, Salomon AK, Norderhaug KM, Pérez-Matus A, Hernández JC, Clemente S, Blamey LK, Hereu B, Ballesteros E, Sala E, Garrabou J, Cebrian E, Zabala M, Fujita D, Johnson JE. Global regime shift dynamics of catastrophic sea urchin overgrazing. Phil Trans R Soc B. 2014;370:20130269. doi:10.1098/rstb.2013.0269. [CrossRef] [Google Scholar]
Lobban CS, Harrison PJ. Seaweed ecology and physiology. Cambridge: Cambridge University Press; 1994. [Google Scholar]
Lüning K. Seaweeds, their environment, biogeography and ecophysiology. New York: Wiley; 1990. [Google Scholar]
Lunn DJ, Thomas A, Best N, Spiegelhalter D. WinBUGS—a Bayesian modelling framework: concepts, structure, and extensibility. Stat Comput. 2000;10:325–337. doi:10.1023/A:1008929526011. [CrossRef] [Google Scholar]
MAFF (2018) Ministry of Agriculture, Forestry and Fisheries. http://www.maff.go.jp/j/tokei/index.html
MHLW (2018) Water Supply in Japan. Ministry of Health, Labour and Welfare. https://www.mhlw.go.jp/english/policy/health/water_supply/1.html
Michaelis L, Menten ML. Die kinetik der invertinwirkung. Biochem Z. 1913;49:333–369. [Google Scholar]
Msuya FE, Neori A. Effect of water aeration and nutrient load level on biomass yield, N uptake and protein content of the seaweed Ulva lactuca cultured in seawater tanks. J Appl Phycol. 2008;20:1021–1031. doi:10.1007/s10811-007-9300-6. [CrossRef] [Google Scholar]
Murphy JT, Johnson MP, Viard F. A modelling approach to explore the critical environmental parameters influencing the growth and establishment of the invasive seaweed Undaria pinnatifida in Europe. J Theor Biol. 2016;396:105–115. doi:10.1016/j.jtbi.2016.01.038. [PubMed] [CrossRef] [Google Scholar]
Nakatsu T. Promotion of seaweed bed creation in fisheries infrastructure improvement. Fish Engine. 2005;42:185–187. [Google Scholar]
Norderhaug KN, Christie H, Fossa JH, Fredriksen S. Fish-macrofauna interactions in a kelp (Laminaria hyperborea) forest. J Mar Biol Assoc UK. 2005;85:1279–1286. doi:10.1017/S0025315405012439. [CrossRef] [Google Scholar]
Orth RJ, Carruthers TJB, Dennison WC, Duarte CM, Fourqurean JW, Heck KL, Hughes AR, Kendrick GA, Kenworthy WJ, Olyarnik S, Short FT, Waycott M, Williams SL. A global crisis for seagrass ecosystems. Bioscience. 2006;56:987–996. doi:10.1641/0006-3568(2006)56[987:AGCFSE]2.0.CO;2. [CrossRef] [Google Scholar]
Phillips DJH, Segar DA. Use of bio-indicators in monitoring conservative contaminants: programme design imperatives. Mar Pollut Bull. 1986;17:10–17. doi:10.1016/0025-326X(86)90797-6. [CrossRef] [Google Scholar]
Pörtner H-O. Ecosystem effects of ocean acidification in times of ocean warming: a physiologist’s view. Mar Ecol Progr Ser. 2008;373:203–218. doi:10.3354/meps07768. [CrossRef] [Google Scholar]
Préat N, De Troch M, van Leeuwen S, Taelman SE, De Meester S, Allais F, Dewulf J. Development of potential yield loss indicators to assess the effect of seaweed farming on fish landings. Algal Res. 2018;35:194–205. doi:10.1016/j.algal.2018.08.030. [CrossRef] [Google Scholar]
Punt AE, Hilborn R. Fisheries stock assessment and decision analysis: the Bayesian approach. Rev Fish Biol Fish. 1997;7:35–63. doi:10.1023/A:1018419207494. [CrossRef] [Google Scholar]
Rainbow PS. Biomonitoring of heavy metal availability in the marine environment. Mar Pollut Bull. 1995;31:183–192. doi:10.1016/0025-326X(95)00116-5. [CrossRef] [Google Scholar]
R Core Team (2013) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria http://www.R-project.org/
Rönnbäck P, Kautsky N, Pihl L, Troell M, Söderqvist T, Wennhage H. Ecosystem goods and services from Swedish coastal habitats: identification, valuation, and implications of ecosystem shifts. Ambio. 2007;36:534–544. doi:10.1579/0044-7447(2007)36[534:EGASFS]2.0.CO;2. [PubMed] [CrossRef] [Google Scholar]
Schramm W. Factors influencing seaweed responses to eutrophication: some results from EU-project EUMAC. J Appl Phycol. 1999;11:69–78. doi:10.1023/A:1008076026792. [CrossRef] [Google Scholar]
Serfor-Armah Y, Carboo D, Akuamoah RK, Chatt A. Determination of selected elements in red, brown and green seaweed species for monitoring pollution in the coastal environment of Ghana. J Radioanal Nuclear Chem. 2006;269:711–718. doi:10.1007/s10967-006-0291-6. [CrossRef] [Google Scholar]
Short FT, Burdick DM. Quantifying seagrass habitat loss in relation to housing development and nitrogen loading in Waquoit Bay, Massachusetts. Estuaries. 1996;19:730–739. doi:10.2307/1352532. [CrossRef] [Google Scholar]
Short FT, Wyllie-Echeverria S. Natural and human-induced disturbance of seagrasses. Environ Conser. 1996;23:17–28. doi:10.1017/S0376892900038212. [CrossRef] [Google Scholar]
Smit AJ. Medicinal and pharmaceutical uses of seaweed natural products: a review. J Appl Phycol. 2004;16:245–262. doi:10.1023/B:JAPH.0000047783.36600.ef. [CrossRef] [Google Scholar]
Sousa WP. Experimental investigation of disturbance and ecological succession in a rocky intertidal algal community. Ecol Monogr. 1979;49:227–254. doi:10.2307/1942484. [CrossRef] [Google Scholar]
Stachowicz JJ, Graham M, Bracken MES, Szoboszlai AI. Diversity enhances cover and stability of seaweed assemblages: the role of heterogeneity and time. Ecology. 2008;89:3008–3019. doi:10.1890/07-1873.1. [PubMed] [CrossRef] [Google Scholar]
Sturtz S, Ligges U, Gelman A. R2WinBUGS: a package for running WinBUGS from R. J Stat Softw. 2005;12:1–16. doi:10.18637/jss.v012.i03. [CrossRef] [Google Scholar]
Verkleij FN. Seaweed extracts in agriculture and horticulture: a review. Biol Agric Hortic. 1992;8:309–324. doi:10.1080/01448765.1992.9754608. [CrossRef] [Google Scholar]
Walker DI, Kendrick GA, McComb AJ (2006) Decline and recovery of seagrass ecosystems—the dynamics of change. Orth RJ, Duarte CM, eds. Seagrasses: Biology, Ecology and Conservation. Dordrecht (The Netherlands). Springer
Yamaguchi A. Biological aspects of herbivorous fishes in the coastal areas of western Japan. Bull Fish Res Agen. 2010;32:89–94. [Google Scholar]
Yang Y, Chai Z, Wang Q, Chen W, He Z, Jiang S. Cultivation of seaweed Gracilaria in Chinese coastal waters and its contribution to environmental improvements. Algal Res. 2015;9:236–244. doi:10.1016/j.algal.2015.03.017. [CrossRef] [Google Scholar]

Articles from Botanical Studies are provided here courtesy of Springer