Estimating Nitrogen and Phosphorus Cycles in a Timber Reef Deployment Area

: In an oligotrophic bay, Mitsu Bay, Japan, artiﬁcial timber reefs (ATRs) are deployed to increase ﬁsh production. In such man-made ecosystems, the biological activities of other organisms as well as the physical structures of ATRs could inﬂuence nutrient cycling. A pelagic–benthic coupling model expressing both phosphorus and nitrogen cycling was developed to investigate seasonal variation in the associated nutrients and their annual budget in the ATR areas and the entire bay system. The model consists of equations representing all the relevant physical and biological processes. The model reproduced the observed seasonal variations in dissolved inorganic P, ammonium, and nitrate concentrations that were low in spring and summer and high in autumn and winter. The internal regeneration rates of the nutrients were two times higher in the ATRs than in the bay area, so that ﬁsh production was predicted to be higher in the ATRs than in the bay area. the inﬂows from the land and precipitation are quite low, nutrient regeneration is an important source of nutrients for the water in Mitsu Bay. ATR deployment could be an important local nutrient source in an oligotrophic bay, and could increase ﬁsh production. and T.U.; validation, J.F.A.; formal analysis, J.F.A.; investigation, S.N. and K.H.; data curation, T.U., S.N. and K.H.; writing—original preparation, J.F.A.; writing—review and T.Y.; supervision, T.Y. and T.U.; T.Y.;


Introduction
The Seto Inland Sea, Japan, has been a eutrophicated coastal area owing to industrialization and urbanization activities in the 1960s. However, currently, the coastal marine environment is oligotrophic because of the strict implementation of nutrient discharge reduction measures in line with the Law Concerning Special Measure for Conservation of the Environment of the Seto Inland Sea, enacted in 1979 [1]. The condition of Mitsu Bay (MB), in the western part of the Seto Inland Sea, is also oligotrophic, and oyster farming in the bay is hampered by poor phytoplankton concentrations. In addition, fish catch in the bay is low [2]. Consequently, the deployment of artificial timber reefs (ATRs) ( Figure S1) was initiated by the local Higashi-Hiroshima City government with the support of the central government, in a project named Regional Revitalization and Local Residents Life, such as Emergency Assistance Grant (Local Creation Proactive) 2015, Cabinet Office, Government of Japan. From the second year, Higashi-Hiroshima City has been overseeing successive research projects, from November 2015 to 2018 [3][4][5][6].
ATRs facilitate the development of habitats for various marine organisms in the deployment area of MB [7]. The encrusting of algae on the ATR surfaces followed by the invasion of other organisms sustains newly established food webs, which eventually attract fish and enhance fish catch.
Several factors influence the habitat development processes, including nutrient supply, which regulates primary production.
However, the material cycles through food webs established among fish and feed organisms have not been investigated quantitatively. Organisms attached to ATRs and fish that gather can serve as nutrient sources through their physiological activities such as excretion and egestion. In the former process, ammonium and urea are excreted, and in the latter process, organic matter and feces are produced and are finally decomposed and released into water as nutrients, such as ammonium and phosphate. In addition, natural mortality may be high during spawning by fish and other organisms at ATRs. The deposition of organic matter through such physiological activities and dead matter that settle on the sea floor also facilitate the organic enrichment of benthic environments. In addition to facilitating biological processes, ATR structures act as traps that collect organic particles from neighboring areas by decreasing the speed of currents [8,9]. Furthermore, other physical processes, reflected in water movement, which facilitate the constant exchange between water in the ATR deployment areas and the surrounding area in MB, are key factors that influence the nutrient status in the deployment area.
Several studies on the modeling of nutrient cycling in natural coastal ecosystems have been conducted [10,11]. Similar attempts to evaluate the success of environmental management activities have been made in man-made ecosystems such sea aquaculture systems [12,13]. In the present study, we estimated nutrient cycling in a new man-made ecosystem established within the ATR systems and bay area based on a pelagic-benthic coupling model. The model consists of two major elements, nitrogen (N) and phosphorus (P), whose dissolved inorganic forms are called "nutrients" and influence the productivity of aquatic ecosystems. Either element can be a limiting factor for the growth of phytoplankton, and they sometimes alternate as limiting factors seasonally and spatially [14]. Therefore, biogeochemical processes are promoted by either element seasonally or spatially based on a benthic-pelagic ecosystem model. However, a numerical model could offer more comprehensive insights on ATR functions.

Materials and Methods
The present study was conducted in MB, Hiroshima, Japan ( Figure 1). The bay is not geographically enclosed, with an open mouth to the Seto Inland Sea. Tidal flow is the major factor driving material transport in the bay, with 0.5-1.2 cm/s in the northern half and up to 2.0 cm/s at the bay mouth [15]. The major components are the principal lunar semi-diurnal (M2) and the solar semi-diurnal (S2). The residual flow in the bay is clockwise, and the water enters from the straight between Oshibajima Island and the mainland in the west, going up to the north, and then down to the southeast passing through the mouth. The average water retention time in the bay has been calculated to be 2.5 days [15]. In the model, the water exchange rate was adjusted to reflect all the observed nutrient concentrations, with an average water exchange rate of 0.4/day. The present model illustrates the nutrient cycles that occur in the pelagic and benthic systems in the ATRs and surrounding waters in the bay, in addition to the food web models of organisms in the area described in the model of Alam et al. [16]. The volume of the entire MB is 144,000,000 m 3 , with a surface area of 16,800,000 m 2 and an average depth of 8.6 m [15]. The total area of the two sites where the ATRs have been established was 320 m 2 , at a depth of 10 m. In the model, we set the total ATR volume to be 3200 m 3 . The structure of the model was illustrated in Figure 2 and the food web structure was illustrated in detail in Figure S2. Three compartments for different forms of P (detritus P; DET-P, dissolved organic P; DOP, and dissolved inorganic P; DIP), and four compartments for different forms of N (detritus N; DET-N, dissolved organic N; DON, ammonium; NH4 + , and nitrate; NO3 − ) were applied for both the pelagic and benthic systems. Subsequently, in the MB area, which does not include the ATR area, three P compartments (DET-P_MB, DOP_MB, DIP_MB) and four N compartments (DET-N_MB, DON_MB, NH4 + _MB, and NO3 − _MB) for each pelagic and benthic system were input in the model. DET-N was calculated from DET-P using the Redfield ratio (molar N/P = 16) as a conversion factor [17]. Subsequently, in the MB water, three P compartments (DET-P_MB, DOP_MB, DIP_MB) and four N compartments (DET-N_MB, DON_MB, NH4 + _MB, and NO3 − _MB) for each pelagic and benthic system were also set to the model to reflect nutrient exchange over time with the materials in the ATR deployment area. Furthermore, there were 21 compartments of living organisms in the ATR area, and eight compartments of living organisms in MB. Their mass balance, parameters, and equations are based on our previously published paper [16]. The present model illustrates the nutrient cycles that occur in the pelagic and benthic systems in the ATRs and surrounding waters in the bay, in addition to the food web models of organisms in the area described in the model of Alam et al. [16]. The volume of the entire MB is 144,000,000 m 3 , with a surface area of 16,800,000 m 2 and an average depth of 8.6 m [15]. The total area of the two sites where the ATRs have been established was 320 m 2 , at a depth of 10 m. In the model, we set the total ATR volume to be 3200 m 3 . The structure of the model was illustrated in Figure 2 and the food web structure was illustrated in detail in Figure S2. Three compartments for different forms of P (detritus P; DET-P, dissolved organic P; DOP, and dissolved inorganic P; DIP), and four compartments for different forms of N (detritus N; DET-N, dissolved organic N; DON, ammonium; NH + 4 , and nitrate; NO − 3 ) were applied for both the pelagic and benthic systems. Subsequently, in the MB area, which does not include the ATR area, three P compartments (DET-P_MB, DOP_MB, DIP_MB) and four N compartments (DET-N_MB, DON_MB, NH + 4 _MB, and NO − 3 _MB) for each pelagic and benthic system were input in the model. DET-N was calculated from DET-P using the Redfield ratio (molar N/P = 16) as a conversion factor [17]. Subsequently, in the MB water, three P compartments (DET-P_MB, DOP_MB, DIP_MB) and four N compartments (DET-N_MB, DON_MB, NH + 4 _MB, and NO − 3 _MB) for each pelagic and benthic system were also set to the model to reflect nutrient exchange over time with the materials in the ATR deployment area. Furthermore, there were 21 compartments of living organisms in the ATR area, and eight compartments of living organisms in MB. Their mass balance, parameters, and equations are based on our previously published paper [16].  [18]. The loading amounts of nutrients from rainfall (precNH + 4 , NO − 3 , DIP) obtained by multiplying the rainwater volume [19] with the reported DIP, NH4 + , and NO3 -concentrations in the rainwater [20]. Furthermore, water exchange (inwaterexch, outwaterexch DET-N, DON, NH + 4 , NO − 3 , DET-P, DOP, DIP) between the ATR area and the MB area was obtained by multiplying the concentrations in either the ATRs or MB with the water exchange rate estimated from the diffusion and advection processes calculated based on the 3D-flow dynamic model (14). In the N and P cycles in the system, excretion (inexcNH + 4 , DIP) from the animals and fish and decomposed products of detritus matter from the processes of mortality (Ntotalmortality, Ptotalmortality) and egestion (Nf + Pf) were considered as recycled nutrient sources that enhance primary production.
The general mass balance equations for the N and P cycles are as follows: In the ATR pelagic system:   [18]. The loading amounts of nutrients from rainfall ( prec NH + 4 , NO − 3 , DIP) obtained by multiplying the rainwater volume [19] with the reported DIP, NH + 4 , and NO − 3 concentrations in the rainwater [20]. Furthermore, water exchange ( in waterexch, out waterexch DET-N, DON, NH + 4 , NO − 3 , DET-P, DOP, DIP) between the ATR area and the MB area was obtained by multiplying the concentrations in either the ATRs or MB with the water exchange rate estimated from the diffusion and advection processes calculated based on the 3D-flow dynamic model (14). In the N and P cycles in the system, excretion (inexc NH + 4 , DIP) from the animals and fish and decomposed products of detritus matter from the processes of mortality (N totalmortality , P totalmortality ) and egestion (N f + P f ) were considered as recycled nutrient sources that enhance primary production.
The general mass balance equations for the N and P cycles are as follows: In the ATR pelagic system: In the Mitsu Bay pelagic system: In the Mitsu Bay benthic system: Most of the parameters for the processes in the model are based on literature, and some parameters are tuned (tuning) to obtain an optimal fit to the observed values (Table 1). Environmental parameters such as temperature were obtained using chlorophyll-turbid meters (INFINITY-CLW, JFE Advantech, Co. Ltd., Nishinomiya, Japan), while light intensity and dissolved oxygen (DO) data were obtained from government website [21]. excOYS Input from excretion of oyster cultivation mgP m −3 day −1 0.012 [30] faeOYS Input from egestion of oyster cultivation mgP m −3 day −1 0.003 [30] mortalOYS Input from mortality of oyster cultivation mgP m −3 day −1 0.5 [30] * Tuning: the parameter is tuned to fit the outputs to the observed data. The parameters are tuned based on the range value reported from several related papers. In addition, the parameters were tuned by slightly increasing or decreasing (in a reasonable range) closer to the reported value from the related paper.
The numerical model was developed using STELLA Architect v1.6.0 (https://www.iseesystems. com). We ran the model for 1020 days with a time step of 0.36 h. The calculations were performed using the fourth-order Runge-Kutta method. The outputs were compared with the values obtained from field observations [20] and the accuracy of validation was evaluated with the percent bias (PBIAS). PBIAS is a measure of over-and under-estimation of bias for predicted and measured values, and is expressed as a percentage [31]. The performance rating in fitting of the model output to the observed value was evaluated with the following criteria [32]: very good (< ±25%), good (±25 PBIAS < ±40), satisfactory (±40 PBIAS < ±70), and non-satisfactory (≥ ±71).

Reproducibility of Calculation Outputs to the Observed Seasonal Variations
The calculated concentrations of DIP, NH + 4 , and NO − 3 nutrients in the ATR deployment area and MB are illustrated in Figure 3. In the ATR deployment area, the calculated DIP and NO − 3 concentrations were consistent with the observed values, which were generally high in autumn and winter, and low in spring and summer. According to the PBIAS values, the calculation outputs for DIP were "very good," with a value of 5%. Conversely, the calculated NH + 4 and NO − 3 concentrations were judged "good," with PBIAS values of −67% and −26%, respectively. The calculated DIP, NH + 4 , and NO − 3 concentrations in the ATR deployment area were approximately two-fold higher than the concentrations in the water in MB. The seasonal trends of the nutrients in the MB area, which were high in autumn and winter, and low in spring and summer, were similar to the trends observed in the ATR area. The observed values of nutrient concentrations outside of the bay (boundary) seem to be higher than the calculated concentrations in the bay.   The calculated P and N concentrations in the sediment are illustrated in Figure 5. The calculated DET-P and DET-N were high in summer and low in winter, respectively. Such trends were also observed in DOP and DON. The seasonal DIP, NH + 4 , and NO − 3 trends in the sediment pore water were similar to those in the pelagic system, being higher in winter and lower in summer.   The calculated P and N concentrations in the sediment are illustrated in Figure 5. The calculated DET-P and DET-N were high in summer and low in winter, respectively. Such trends were also observed in DOP and DON. The seasonal DIP, NH + 4 , and NO − 3 trends in the sediment pore water were similar to those in the pelagic system, being higher in winter and lower in summer. The calculated P and N concentrations in the sediment are illustrated in Figure 5. The calculated DET-P and DET-N were high in summer and low in winter, respectively. Such trends were also observed in DOP and DON. The seasonal DIP, NH + 4 , and NO − 3 trends in the sediment pore water were similar to those in the pelagic system, being higher in winter and lower in summer.

Balances in Inputs and Outputs
The inputs to and outputs from DIP, NH + 4 , and NO − 3 stocks in the ATR area and the MB area are presented in Table 2. In the ATR area, the major inputs to the DIP stock were organism excretions (27%), followed by decomposition processes (11%). Conversely, the major outputs from the DIP stock were the phytoplankton uptake (35%) and physical transport to the surrounding water via water exchange (16%). With regard to NH + 4 , organism excretions (24%) was also the highest contributor to the stock. In addition to organism excretions, the contribution of benthic flux was high (20%).

Balances in Inputs and Outputs
The inputs to and outputs from DIP, NH + 4 , and NO − 3 stocks in the ATR area and the MB area are presented in Table 2. In the ATR area, the major inputs to the DIP stock were organism excretions (27%), followed by decomposition processes (11%). Conversely, the major outputs from the DIP stock were the phytoplankton uptake (35%) and physical transport to the surrounding water via water exchange (16%). With regard to NH + 4 , organism excretions (24%) was also the highest contributor to the stock.
In addition to organism excretions, the contribution of benthic flux was high (20%). Nitrification, the process of oxidation of ammonium into nitrate, was the major pathway for the decrease in NH + 4 concentration, accounting for 42%. Other outputs, including phytoplankton uptake (7%) and transport to the surrounding water (3.3%), were also pathways of reduction of NH + 4 concentrations. Nitrification was a major pathway of the increase in NO − 3 stock, accounting for 37%. For the loss of NO − 3 stock, uptake by producers accounted for 31% and 12% was attributed to denitrification.
The major inputs to DIP in the MB area were transport from the ATR areas (25%), transport from outside the MB (13%), and excretion by organisms (11%). Other sources such as input from rivers, sewage discharge, precipitation, and benthic flux were very low (<1% in total). With regard to DIP outputs, phytoplankton uptake represented the greatest pathway (35%). In the case of NH + 4 , the major input to the MB water was excretion (32%), although the rate was two times lower in the MB water, while the highest output was the nitrification process (33%). In the case of NO − 3 in MB, the major input was transport from the ATR areas (23%). Other inputs were nitrification from ammonium and transport from outside of the MB area, which accounted for 16% and 11%, respectively. In addition, the major NO − 3 outputs were phytoplankton uptake (23%) and transport to the ATR areas (14%).

Internal Regeneration
Internal regeneration was the major source of P and N inputs in both the ATR systems and the MB area (Figures 6 and 7). In the ATR area, internal regeneration accounted for 0.04 kg P day −1 (14.9 mg P m −3 day −1 ), while in MB, it accounted for 449 kg P day −1 (3.1 mg P m −3 day −1 ). In the case of N, internal regeneration accounted for 0.1 kg N day −1 (31.7 mg N m −3 day −1 ) in the ATR area and 1695 kg N day −1 (11.8 mg N m −3 day −1 ) in MB. The ATR area was a source of P, NH + 4 , and NO − 3 for MB, at 0.01 kg P day −1 , 0.01 kg N day −1 , and 0.05 kg N day −1 , respectively.

Internal Regeneration
Internal regeneration was the major source of P and N inputs in both the ATR systems and the MB area (Figures 6 and 7). In the ATR area, internal regeneration accounted for 0.04 kg P day −1 (14.9 mg P m −3 day −1 ), while in MB, it accounted for 449 kg P day −1 (3.1 mg P m −3 day −1 ). In the case of N, internal regeneration accounted for 0.1 kg N day −1 (31.7 mg N m −3 day −1 ) in the ATR area and 1,695 kg N day −1 (11.8 mg N m −3 day −1 ) in MB. The ATR area was a source of P, NH + 4 , and NO − 3 for MB, at 0.01 kg P day −1 , 0.01 kg N day −1 , and 0.05 kg N day −1 , respectively.

Internal Regeneration
Internal regeneration was the major source of P and N inputs in both the ATR systems and the MB area (Figures 6 and 7). In the ATR area, internal regeneration accounted for 0.04 kg P day −1 (14.9 mg P m −3 day −1 ), while in MB, it accounted for 449 kg P day −1 (3.1 mg P m −3 day −1 ). In the case of N, internal regeneration accounted for 0.1 kg N day −1 (31.7 mg N m −3 day −1 ) in the ATR area and 1,695 kg N day −1 (11.8 mg N m −3 day −1 ) in MB. The ATR area was a source of P, NH + 4 , and NO − 3 for MB, at 0.01 kg P day −1 , 0.01 kg N day −1 , and 0.05 kg N day −1 , respectively.

Seasonal Variation in Nutrient Concentrations and the Potential Causes
The dissolved inorganic forms of P and N, such as DIP, NH + 4 , and NO − 3 , which are referred to as "nutrients", are required by marine algae for their growth and they regulate primary production in the seas [33]. The calculated nutrient concentrations in the ATR pelagic system were consistent with the observed values, with similar seasonal trends, being high in autumn and winter, and low in spring and summer (Figure 3). These seasonal trends were also similar to the results observed in Hiroshima Bay, which is next to MB [11], and to results reported in artificial reefs in southern Portugal [33]. As represented in Table 2, the major causes of increase in DIP and NH + 4 in the pelagic system were excretion by organisms, whose biomass peaked in winter.
The calculated concentrations of DOP, DON, DET-P, and DET-N stocks in the ATRs were in the range of those in Hiroshima Bay, which is a more eutrophicated bay than MB, and is located approximately 50 km west of MB [11]. Among these parameters, the calculated average concentrations of DET-P and DET-N in MB were higher than those reported for Hiroshima Bay. They were 14 mg P m −3 and 44 mg N m −3 , compared to approximately 7 mg P m −3 and 40 mg N m −3 for DET-P and DET-N, respectively, in Hiroshima Bay.
The nutrient concentrations in MB were lower than those in the ATR areas and those outside the bay. The nutrient concentrations outside the bay used as the boundary conditions for MB were not as low as the concentrations in the MB areas ( Figure 3). The average data were taken from three layers in the water column, at 0, 5, and 1 m above the bottom. Regarding the bottom topography, it goes deeper probably sharply at the mouth of the bay. The outside of the bay, which had a deeper layer, had potentially higher concentrations when compared to the relatively shallow MB. Deeper layers often have higher nutrient concentrations as a result of the decomposition of organic matter. Another factor is the cooling of the water column in winter. The overturning of the water column in the cooler season is another cause of increase in nutrient concentration, as observed in the larger gap during winter in the observed data outside of MB.

Internal Regeneration of Nutrients in the ATR Deployment Area
Higher nutrient concentrations were estimated in the ATR area compared to the waters in MB. In the ATRs, internal nutrient regeneration was the major process compared to the other sources. In the entire MB, the internal regeneration also played a key role as a major source of nutrient input when compared to the relatively low inputs of other sources such as river inflows, sewerage, etc., ( Table 2). As shown in Table 2, the internal regeneration contributed considerably to the bay primary production. In the ATR system, the biomass of organisms and fish was high compared to those in the surrounding area. As they produce organic matter as feces, urine, and sometimes dead matter, such organic matter becomes nutrient supply sources, which are utilized by algae again [34]. The estimates summarized in Table 2 show that the DIP and NH + 4 inputs from excretion in the ATR areas were two to three times higher than the inputs in MB.
ATRs also increase biodeposits, which are not only feces but also dead attached animal bodies peeled off from the ATRs. The deposition rate was calculated to be 47.7 mg P m −3 day −1 and 62.6 mg N m −3 day −1 for P and N, respectively. The deposition, in turn, increases the organic matter load to the bottom, and can lead the anoxic sediment status and reduced quality status overall [8]. Enhanced deposition to the benthic systems was also observed in an AR study in southern Portugal [35]. The materials used in the present study, including timber, branches, and leaves, are also sources of organic matter. Such materials decay following degradation by ship worms or the effects of physical activities such as wake waves, which were observed by divers during the observations in the third year after ATR deployment. Overall, the deployment of ATRs had positive effects by serving as sources of nutrients both through increases in the associated organisms on the surfaces of ATRs and the gathering of fish at the sites to feed and hide among the frames and structures.

Denitrification
Denitrification is a major pathway through which nitrogen is removed from a system, and it often occurs in surface anoxic sediment under nitrate supply. The N2 removed in this process was released to the air through the water surface. The sediment quality in MB is partially anoxic due to feces and pseudofeces loads from intensive oyster aquacultures in the area [15]. As calculated in the present study, 14.8 mg N m −2 day −1 or approximately 0.04 kg N day −1 was removed in the ATR area through denitrification, while 5.6 mg N m −2 day −1 was removed in the MB area (Table 2). Compared to the increase in feces sedimentation, the removal of N by denitrification was up to 12% in ATRs and 8% in MB. The denitrification rate estimated in the present study was in the range reported in previous studies, including at 4.2-5.88 mg N m −2 day −1 [11] and 0.0-23.2 mg N m −2 day −1 [36], which was reported in studies conducted in Hiroshima Bay, which is adjacent to MB. The range reported in the present study was also within the range reported in Kansas, United States (0-35 mg N m −2 day −1 ) [37], and that reported in an estuary in Portugal (2-84 mg N m −2 day −1 ) [38]. Therefore, denitrification can alleviate the deterioration of sediment quality below the ATRs. However, P would definitely be accumulated in the sediment.

Potential Oligotrophication Countermeasures
Our model reveals that ATRs function as a source of nutrients for the surrounding water (Table 2). Because other nutrient sources such as inputs from land, precipitation, and benthic flux are quite low, the nutrient regeneration in the ATR area is considered an important nutrient source for MB. This means that ATRs could be one of the promising ways for increasing fish production in an oligotrophic bay such as MB.
In the Seto Inland Sea, cultural oligotrophication is widespread due to the strict implementation of reduction measures for nutrient loads from land as mentioned in the introduction. In several prefectures facing the Seto Inland Sea, sewer discharge regulations have been relaxed. For example, Hyogo Prefecture began relaxing sewer discharge regulations more than five years ago, and they set the target levels of N and P nutrient concentrations in seawater required for the maintenance of the local fisheries represented by Nori farming and sand eel [39]. Hiroshima Prefecture, in which MB is found, has initiated discussions on the measures of counteracting the oligotrophication of the coastal seas. However, they are very careful about relaxing sewer operation regulations, because they could inadvertently generate phytoplankton blooms.
In addition, the deployment of ATRs, as demonstrated in the present study, can be a superior strategy for addressing oligotrophication to the relaxation of sewerage disposal regulations. No net increase in nutrient loads have been observed from out of the system in the case of the deployment of timber reefs, which would otherwise be observed following sewage discharge. To implement either of the measures or both, their effects on the alleviation of oligotrophication need to be evaluated based on scientific analyses. Currently, Higashi-Hiroshima City is exploring strategies for increasing fish production in the future, and the optimal strategy could involve a mixture of the two approaches, based on the results of quantitative scientific evaluations. We believe that the results of the present study provide a scientific basis for decision-making by the stakeholders in Higashi-Hiroshima City.

Conclusions
Our model demonstrates that the deployment of ATRs could enhance the internal regeneration of nutrients (DIP, NH + 4 , and NO − 3 ) in ATR areas by increasing physiological activities, as exemplified by the excretion and feces production by the associated biological communities, including fish. ATRs, by transporting nutrients to the surrounding water, could also be an important local nutrient source in the oligotrophic area. In several neighboring provinces, local authorities have begun relaxing sewer discharge regulations. In Hiroshima Prefecture, the local government is also planning to allow the discharge of sewer into the marine ecosystem. The model we developed in the present study can guide administrators to select appropriate measures of increasing productivity in fisheries. In conclusion, the deployment of ATRs is a potential measure of increasing fish production, by accelerating nutrient cycling without any extra external nutrient inputs. In addition, it facilitates the establishment of a recycling-oriented society, in the case of Higashi-Hiroshima City, which has large forest areas that need to be maintained.