Modeling the E ﬃ cacy of Three Types of Artiﬁcial Timber Reefs in Mitsu Bay, Japan

: With the aim of simultaneously improving ﬁshery production and utilizing forestry and oyster ﬁshery wastes, three types of artiﬁcial timber reefs (ATRs)—constructed from simple timbers, timbers with oyster shells from local oyster farms, and timbers with leaves / branches from forest thinning—were deployed in Mitsu Bay, Japan. We developed a food web model to investigate the relative e ﬃ cacies of these ATR types compared with the bare, sandy seaﬂoor. The model described the material ﬂow through the food webs formed in each ATR type and their potential to increase ﬁsheries production. The model outputs were validated with observational data over three years. The model ﬁt the observed biomass of both prey animals and ﬁsh predators. The simulation results highlighted that ATRs, particularly those with additional materials, had two to three times higher feeding ﬂow than the sandy seaﬂoor and resulted in increased ﬁsh biomass. Fish catch doubled in the ATR areas compared to the bare seaﬂoor. Aside from providing a feeding ground, the complexity of the ATRs with additional materials likely acts to provide shelter for juvenile ﬁsh. ATR deployment using by-products such as those mentioned above may not only enhance ﬁsh stock but also help foster the establishment of a recycling-oriented society. Subscripts abbreviations: micro algae biomass (MIA), macro algae biomass (MAA), benthic micro algae biomass (MIAb), zooplankton biomass (ZOO), shrimp biomass (SHR), bivalves biomass (BIV), benthic bivalves biomass (BIVb), sedentary worm biomass (SWO), benthic sedentary worm biomass (SWOb), errant worm biomass (EWO), bryozoan biomass (BRY), sea squirt (SSQ), gammarid biomass


Introduction
Artificial reefs (ARs) are human-made structures deployed in the aquatic ecosystem for various services and benefits. The deployment of ARs is believed to benefit fishery yield by providing a habitat for fish and their feeds [1,2]. The type of structures and materials used for the construction of ARs is key for meeting the expected goals. The materials used for AR construction vary around the world, ranging from either natural materials such as rocks, shells, and wood or artificial materials such as concrete, fiberglass, tires, vessels, and vehicles [3].
In Japan, fishery production has been decreasing over the past few decades [4]. The current estimated domestic fishery production is approximately 60% and the remaining 40% is imported. Ensuring the long-term sustainability of fisheries is essential to reduce dependency on imported seafoods. The national government attempted to increase fish stocks by deploying ARs in the coastal waters [5]. During the 1950s, with considerable government funding, ARs were primarily constructed of concrete [6]. Recently, with less funding from the central government, the trend of AR material has shifted from concrete to natural materials such as timber. The government of Higashi-Hiroshima City, in coordination with the local fishery association, deployed artificial timber reefs (ATRs) in Mitsu Bay, Japan in late 2015. Timber is plentiful, as a by-product of forest thinning in the forested areas surrounding this city. It is also reported that timbers attract fish more efficiently than other materials, and hence, they are ideal as an AR construction material [7]. In Mitsu Bay, the main fishery activity is oyster cultivation. This industry generates large quantities of oyster shells. Thus, we tested oyster shells as supplementary materials for ATRs, as well as the leaves and branches of timbers from thinning.
The use of timbers for ARs is not uncommon in Japan [8] and elsewhere [9][10][11]. Nevertheless, there have been few numerical studies in terms of the efficacy of ARs in structuring food webs in new ecosystems and their implications in increasing fishery production [12,13]. The dynamics of predator-prey interactions within the food web structure of an aquatic ecosystem are multifactorial and thus highly complex. ATR structures provide settlement locations for sessile organisms on the reefs and for benthic organisms beneath the reefs. Further, ATRs may enhance the production of higher trophic level animals such as fish through changes in prey composition, resulting in the formation of a new ecosystem [14,15].
In the present study, we focused on the difference in efficacy of the three types of ATRs-simple timbers, timbers with oyster shells, and timbers with leaves/branches. The efficacies, in terms of food web structure, were tested by applying a numerical model. Thus, the objectives of the study were as follows: (1) to evaluate the increase in fish production as a result of the presence of ATRs, particularly with the supplementary materials such as oyster shells and leaves/branches, and (2) to understand the difference in material flows induced by the different animal and fish compositions created in the three different types of ATRs by numerical model analyses.

Artificial Timber Reefs and Their Deployment
ATRs (1.5 m × 1.5 m × 1.5 m) were constructed with Japanese cypress timbers in a cube-shaped parallel cross, following the method of Alam et al. [1] ( Figure S1). A total of 15 ATRs (three simple timber reefs; six timber reefs with oyster shells, ATRsOS; and six timber reefs with leaves and branches, ATRsLB) were deployed in November 2015 at 10 m depth and arranged as shown in Figure 1

Sampling
Field observations were conducted eight times for three years after ATR deployment: January and July 2016; January, July, and November 2017; and January, March, and August 2018. Prey (zooplankton, attached and benthic) organisms were sampled from the ATRs and bare seafloor. In the ATR deployment area, zooplankton samples were collected using a 100-µm mesh Kitahara plankton net, 22.5 cm in diameter (0.04 m −2 ), vertically towed from the bottom to the surface just beside the ATR constructions. The zooplankton samples were fixed with 5% formaldehyde in the final concentration; then, they were sorted, counted, and weighed in the laboratory. Test pieces of timber and net baskets containing supplementary materials were collected by two SCUBA divers. All samples of test pieces of timber and the net basket were packed in plastic bags to prevent the organisms on and in them from escaping. Benthic organisms were also collected from the top 5 cm of sediment at the edge of each ATR using a 25 cm × 25 cm (625 cm 2 ) grab sampler. All samples were sieved with a 1.0-mm mesh, fixed in 10% formaldehyde, and transported to the laboratory for sorting, counting, and weighing.

Sampling
Field observations were conducted eight times for three years after ATR deployment: January and July 2016; January, July, and November 2017; and January, March, and August 2018. Prey (zooplankton, attached and benthic) organisms were sampled from the ATRs and bare seafloor. In the ATR deployment area, zooplankton samples were collected using a 100-m mesh Kitahara plankton net, 22.5 cm in diameter (0.04 m −2 ), vertically towed from the bottom to the surface just beside the ATR constructions. The zooplankton samples were fixed with 5% formaldehyde in the final concentration; then, they were sorted, counted, and weighed in the laboratory. Test pieces of timber and net baskets containing supplementary materials were collected by two SCUBA divers. All samples of test pieces of timber and the net basket were packed in plastic bags to prevent the organisms on and in them from escaping. Benthic organisms were also collected from the top 5 cm of sediment at the edge of each ATR using a 25 cm × 25 cm (625 cm 2 ) grab sampler. All samples were Fish abundance was evaluated using video recordings. This method was adopted as it provides permanent recordings of data, permits observations for fixed time periods, and enables data collection under conditions unsafe for diver-based assessments [16,17]. Furthermore, video recording avoids miscalculations that can occur when diver movements cause fish escape reactions, resulting in the loss of fish to observation. Neither method works at night, so nocturnal fish biomass was underestimated. In total, eight underwater cameras (Hero Session 4, GoPro Co. Ltd., San Mateo, CA, USA) were placed 1.5 m in front of each ATR and bare seafloor, i.e., one camera each for the different types of ATRs and the bare seafloor in the two sites. Each camera had a field of view covering the associated single ATR structure plus 1 m to each side. The total volume observed by each camera was estimated to be 11.67 m 3 . Camera recordings were operated for a net of one hour during 9:00-11:00 in the morning. With this way of sampling, we missed nocturnal fish. However, the operation of the camera at night would have required additional light to catch fish in video recording, which may have affected the results, particularly by selecting fish with positive phototaxis; as such, it was decided that nocturnal fish would not be observed, and the recording would only be done during daytime. All of the fish observed in the footage were classified to the species or genus level and enumerated. The fish abundance determined in the present study resembled the fish abundance (MaxN) metric of Scott et al. [18]. Fish abundance is the maximum number of fish species observed in a single frame of footage during one deployment [19]. As fish size measurements based on the video scenes were inaccurate, size ranges (0-15 cm, 15-30 cm, and >30 cm) were applied, and the biomass was estimated using the published length-weight relationship for each fish species [20]. The averages of the prey and predator biomasses from the two sites for each treatment were then used for further analysis in the model.
In addition to the samplings of biotic parameters, abiotic parameters that affect the growth of primary producers were measured. Chlorophyll-turbid meters (INFINITY-CLW, JFE Advantech, Co. Ltd., Kobe, Japan) were used to record temperature, relative chlorophyll fluorescence, and turbidity, and PAR loggers (DEFI2-L, JFE Advantec, Co. Ltd., Kobe, Japan) were used to record the underwater photon density; the chlorophyll-turbid meters and photosynthetic active radiation (PAR) loggers were attached on the ATR structures at each site and regularly replaced at every sampling event. The sampling interval was set to every hour, and the sample number was 10 at each sampling event. We averaged the 10 data points in each hour for use.

Structure of the Numerical Model
The numerical model developed in this study described the food webs created in the ATRs and differences in the material flow seen in the different types of ATRs. In each case, simplified frameworks were composed that illustrated the dynamics of the predator-prey relationships ( Figure 2). Each ATR had 21 compartments and each bare seafloor had 8 compartments, reflecting its scarcity of residing organisms as a result of not having ATR structures. In all 21 compartments in the ATRs, biomass changes over time were governed by biotic processes and by each compartment's relationship with the other compartments, as shown in Table 1.      Changes in the biomass of the primary producers, planktonic microalgae (MIA), benthic microalgae (MIAb), and macroalgae (MAA) were functions of nutrient uptake, temperature, and light intensity. Their loss processes were excretion, temperature-dependent mortality, and grazing by herbivores. The changes in the biomasses of zooplanktons, attached animals residing in the ATRs, and benthic animals were functions of food availability, competition with other animals at the same trophic level, egestion, mortality, and predation.

Symbols Expression
For the changes in fish biomass, swimming behavior was taken into consideration as described below, along with physiological functions such as grazing, excretion, and mortality applied as described above. Unlike attached and benthic organisms, which are mostly sessile, some fish are mobile and swim in and around the observation area. The change in fish biomass was affected by fish mobilization in addition to the changes governed by feeding, egestion, excretion, natural death, and catch described by Equation (1) [20].
where N is the fish abundance (inds m −3 ) and W is the individual wet weight (g ind −1 ) estimated from the length (TL) and weight (BW) relationship of each fish species in Equations (2)-(7) [21].
The dominant fish species incorporated into the present model were previously categorized as resident fish [22,23]. These fish tended to remain in the AR area. Other species were designated as visitor fish and were observed only once in the AR area. Transient fish such as Japanese mackerel Trachaurus japonicus were observed infrequently and would have contributed a disproportionate amount to the total biomass. Hence, they were excluded from the model calculation to avoid introducing bias into the analysis of the AR ecology. Eliminating the influences of visitor and transient fish from the analysis reveal more pronounced seasonal variations in fish abundance and biomass [24].
The feeding preference of each fish was also incorporated in the model. The values used in this model lie within the ranges reported in published papers ( Table 2). Fish egestion and excretion were set at the same values for all fish, at 20% and 11% of feeding, respectively [25,26]. When the fish reached minimum catch size, loss to fish catch became an exit flow from the relevant biomass compartment, with the rate depending on the catch. The catch size of each fish was determined from information available from the internet, fishermen, and market sources. The daily catch rates were adjusted so that the temporal changes in fish biomass fit the observed values, and also to obtain minimum percent bias (PBIAS). The equations used in all the processes are as described below. The uptake of algae (e.g., MIA) is as expressed in Equation (8), dependent on temperature and light intensity.
The function of growth of attached, benthic organisms and fish (e.g., SHR) is shown in Equation (9). Furthermore, the function of excretion, egestion, and mortality for organisms are shown in Equations (10)- (12). These equations are also applied to the other organisms in this study.
For fish, an addition function of catch (e.g., BSE) was used, as shown in Equation (13).
The parameters used in this study are shown in Table 3. Although most of the model parameter values were obtained from published literature, some of them were tuned to fit the model output to the observed data. The numerical model was developed using STELLA Architect for Macintosh (v. 1.6.0) and run with a time step of 0.015 day. Calculations were conducted for 1020 days by a 4th-order Runge-Kutta method.

Statistical Analysis
The data obtained from the observations from the two sites were compared to investigate the efficacy of the different types of ATRs in terms of increase in prey animals and fish. The difference in attached and benthic organisms and fish were tested using one-way ANOVA after verifying the homogeneity of variances. The average values of the two sites, along with their variations, were used to validate the model, and checked with the percent bias (PBIAS), which is a measure of overand under-estimation of bias for predicted values from the measured values, and it is expressed as a percentage [33]. The performance rating of the model outputs was judged using the criterion as follows: very good (PBIAS < ±25%), good (±25 PBIAS < ±40), satisfactory (±40 PBIAS < ±70), and non-satisfactory (≥±70) [34]. Table 3. List of parameters used in the model.

Results
Organisms began to encrust on the ATR surfaces soon after ATR deployment. The calculations from the numerical model of the attached and benthic animal biomass showed good and satisfactory fits to all observed values; the PBIAS values were in the 25-70% range for all, except for bivalves (BIV), errant worm (EWO), caprellid (CAP), and other amphipods (OAM) in the artificial timber reefs with leaves and branches (ATRsLB) ( Table S1). The biomass values for all attached organisms in the artificial timber reefs with oyster shells (ATRsOS) and ATRsLB were higher than those for the simple ATRs ( Figure 3). CAP, SHR, BIV, EWO, and SWO were significantly higher in ATRsOS and ATRsLB compared to ATRs (p < 0.05) ( Table 4). Among the benthic organisms, the biomass values for two major groups, benthic bivalves (BIVb) and sedentary worms (SWOb), did not differ significantly between the ATRs and the bare seafloor areas during early deployment, although biomass in the ATRs was slightly higher in the third year ( Figure 4).      Figure 3. Relative temporal changes in attached organisms on simple artificial timber reefs and those with either supplementary oyster shells or supplementary leaf and branch materials. The paired superscripts, a and b, indicate statistically significantly differences (p < 0.05), while a and a, or b and b are not significantly different.  Animals, regardless of whether attached or benthic, were the prey of fish, varying substantially, as summarized in Table 2. GAM, CAP, OAM, and SHR were generally eaten by most of the dominant fish, so their biomass profiles were shaped by heavy predation (Table 2). SWO was less eaten by dominant fish or could escape from predation, and consequently it had higher biomass values, especially in ATRsOS and ATRsLB. Their biomass was shaped mainly by the losses on egestion and natural mortality. Bryozoan (BRY) and sea squirt (SSQ) were also not fed on by fish as much, and     Figure 4. Relative temporal changes in benthic organisms in artificial timber reefs and bare seafloor areas.
The annual average stocks and flows within the food webs in the three ATRs and bare seafloor areas are shown in Figures S2-S5. The main inputs for attached organisms were feeding on producers: phytoplankton (PHY) in the water column, and microalgae (MIA) and macroalgae (MAA) on the timber surfaces. Zooplankton (ZOO) and other organisms such as gammarids (GAM), caprellids (CAP), and other amphipods (OAM) were eaten by shrimp (SHR) and errant worms (EWO). In the sediment, the main input for BIVb and SWOb growth was feeding on benthic microalgae (MIAb) and PHY. There was a larger feeding flow from producers in ATRsOS and ATRsLB than was seen in simple ATRs. The total grazing rates on algae (PHY, MIA, MAA, and MIAb) by attached organisms and benthic animals were an order of magnitude larger in ATRsOS and ATRsLB reefs (69.2 and 49.7 g m −3 day −1 , respectively) than that in simple ATRs (4.3 g m −3 day −1 ). At the bare seafloor, the rate was also low, with a value of 2.8 g m −3 day −1 (Figures S2-S5).
Animals, regardless of whether attached or benthic, were the prey of fish, varying substantially, as summarized in Table 2. GAM, CAP, OAM, and SHR were generally eaten by most of the dominant fish, so their biomass profiles were shaped by heavy predation (Table 2). SWO was less eaten by dominant fish or could escape from predation, and consequently it had higher biomass values, especially in ATRsOS and ATRsLB. Their biomass was shaped mainly by the losses on egestion and natural mortality. Bryozoan (BRY) and sea squirt (SSQ) were also not fed on by fish as much, and their biomass profile was shaped by natural mortality. In the first year of the ATRs deployment, the BRY biomass was highest on the simple ATRs and ATRsLB reefs, and it was second highest on the ATRsOS reefs ( Figure 3). However, it steadily decreased from the second year onwards, while SSQ biomass increased.
Twenty-nine fish species were observed at the study sites during the three years of field observations. Six species were found most frequently and were targeted for calculation in the model as dominant fish. Black rockfish Sebastes inermis (BRO) and Japanese surfperch Ditrema temmincki (JSU) were highly associated with all ATRs at all sampling dates. There were no clear associations of these fish species with the bare seafloor, as the fish only roamed over them.
The model reproduced the observed fish biomass ( Figure 5). The PBIAS values showed very good, good, and satisfactory fits to all observed values, i.e., the PBIAS values were in the 25-70% range (Table S1), except for black scraper Thamnaconus modestus (BSC), which was categorized as non-satisfactory. However, the trend of the calculated value followed the observed data except at the point in the third year of deployment where there was a large gap between the observed and calculated values. The average observed black seabream Acanthopagrus schlegelii (BSE) biomass was 62, 31, 59, and 48 g m −3 in ATRs, ATRsOS, ATRsLB, and bare seafloor, respectively. While the BSE observed around ATRs were mainly juveniles of smaller size, as recorded in the video, those observed in bare seafloor areas were mostly adult size. The highest average BRO biomass was recorded around the ATRsLB (223 g m −3 ), followed by that around the simple ATRs (173 g m −3 ) and that around the ATRsOS (140 g m −3 ). Meanwhile, no BRO was observed in bare seafloor sites at any time during the observation period. The BRO biomass increased over time in all ATRs, peaking during the colder season and troughing during the warmer season.
The JSU biomass was highest in the colder season and lowest in the warmer season. The average JSU biomass were 14, 8, and 17 g m −3 in ATRs, ATRsOS, and ATRsLB, respectively. In contrast, the average JSU biomass in bare seafloor areas was only 0.5 g m −3 . The average multicolorfin rainbowfish Parajulis poecilopterus (MRA) biomass was relatively equal across all ATRs and was highest in the warmer season and lowest in the colder season. The average MRA biomass for the ATRs and ATRsOS was 9 g m −3 in each case. This was somewhat lower than the average MRA biomass values for the ATRsLB, which was 15 g m −3 . The average MRA biomass in bare seafloor areas was only 0.9 g m −3 . For all ATRs, the average BSC biomass had increased by summer 2017, two years after deployment. The average BSC biomass values was fairly similar among ATRs, whereas in bare seafloor areas, no BSC were observed. The average red seabream Pagrus major (RSE) biomass were 19, 45, and 29 g m −3 in ATRS, ATRsOS, and ATRsLB, respectively. In the bare seafloor areas, it was 55 g m −3 .
The model calculation of total consumption of attached and benthic animals by fish was highest in the ATRsLB, at a rate of 37.9 g m −3 day −1 , followed by the simple ATRs and ATRsOS, at 27.6 g m −3 day −1 and 27.5 g m −3 day −1 , respectively, while in the bare seafloor areas, it was only 12.5 g m −3 day −1 . About 10% of outflow losses from fish stocks were due to excretion in both ATR and bare seafloor areas; a further 22-24% and 20% were due to egestion in the ATRs and bare seafloor areas, respectively. The natural mortality of attached and benthic animals constituted 25-28% of the loss in the ATRs and 17% in bare seafloor areas.
The calculated fish catch made further reductions to these natural losses. Its impact differed among ATRs. The loss values were approximately 3.7 kg m −3 year −1 , 4.2 kg m −3 year −1 , 5.4 kg m −3 year −1 , and 2.3 kg m −3 year −1 in the simple ATRs, ATRsOS reefs, ATRsLB reefs, and bare seafloor areas, respectively. BRO was the most frequently caught fish in all of the ATRs, at 50%, 33%, and 51% in the simple ATRs, ATRsOS, and ATRsLB, respectively. In bare seafloor areas, the fish catch accounted for 52% of total fish loss, and the most frequently caught fish there was RSE (90%). The JSU biomass was highest in the colder season and lowest in the warmer season. The average JSU biomass were 14, 8, and 17 g m −3 in ATRs, ATRsOS, and ATRsLB, respectively. In contrast, the average JSU biomass in bare seafloor areas was only 0.5 g m −3 . The average multicolorfin rainbowfish Parajulis poecilopterus (MRA) biomass was relatively equal across all ATRs and was highest in the warmer season and lowest in the colder season. The average MRA biomass for the ATRs and ATRsOS was 9 g m −3 in each case. This was somewhat lower than the average MRA biomass values for the ATRsLB, which was 15 g m −3 . The average MRA biomass in bare seafloor areas was only 0.9 g m −3 . For all ATRs, the average BSC biomass had increased by summer 2017, two years after deployment. The average BSC biomass values was fairly similar among ATRs, whereas in bare seafloor areas, no BSC were observed. The average red seabream Pagrus major (RSE) biomass were 19, 45, and 29 g m −3 in ATRS, ATRsOS, and ATRsLB, respectively. In the bare seafloor areas, it was 55 g m −3 .
The model calculation of total consumption of attached and benthic animals by fish was highest in the ATRsLB, at a rate of 37.9 g m −3 day −1 , followed by the simple ATRs and ATRsOS, at 27.6 g m −3

Discussion
Novel community assemblies developed rapidly (within three months) following the deployment of the ATRs, with the attachment of organisms to the ATR structures. ATRs also supply prey for higher trophic level animals by enabling prey animals to reside inside the construction materials. Japanese cypress trees, the timbers from which were used for the ATR frameworks in this study, emit volatile phytoncides that protect the trees from wood-boring invertebrates. For this reason, this tree has superior durability and is commonly used as a construction material in traditional Japanese buildings [74,75]. However, this repellent effect is relatively weak in seawater, as reported by Masuda et al. [7]. Nevertheless, the rough surface of the cypress wood supports organism attachment and growth. This is an advantage of using timbers, unlike the concrete materials that are commonly used for ARs around Japan.
The relatively higher prey biomass on the modified ATRs (ATRsOS and ATRsLB) was observed comparing to that on the simple ATRs. This suggests that the complex structures formed inside the ATRs with the supplemental materials may have contributed to the increase in the number of organisms attached in the space. The simulation results reproduced the higher attached animal biomass in ATRsOS and ATRsLB compared to the simple ATRs. To some extent, this can be attributed to the dominant presence of filter feeders on the additional settlement materials. This domination of filter feeders in the biomass has also been commonly reported in other studies, as in epibiotic communities in the ARs in Hong Kong [76] and in the designated artificial reef (DAR) off the coast of Sydney [77].
After the deployment of ATRs, the first settler was bryozoan as represented on ATRsLB. Bryozoans are reported to be the pioneer species in the newly formed ecosystem of this type [78]. Bryozoan biomass was succeeded by sea squirt in the second year, which was possibly because of space exclusion and food source competition. A succession from bryozoan to sea squirt on wood panels was also observed in Hong Kong waters [79]. Sea squirt is often observed as an invasive species in marine structures [80]. It may penetrate mature assemblages and subsequently become the dominant species [81,82]. Bryozoans and sea squirt are not likely the animals fed preferentially by fish. However, these filter feeders may act to collect particulate organic matter including phytoplankton and supply organic matter to the benthic animals [83,84].
Bivalve and sedentary worm biomass increased over time against the early occupation by bryozoans and sea squirt. Bivalves and sedentary worms are qualified food for fish [85]. Our study emphasizes the importance of ATR structures in that they may contribute to increasing levels of prey animals for fish. The attached and benthic prey organisms on all ATRs seemed to support far more fish biomass than in bare seafloor areas, creating and shaping new food webs. AR deployment in the Caribbean significantly altered the faunal biomass [86]. Thus, ATR deployment exerts an impact by activating secondary production and further energy transfer to higher trophic levels. As seen in Figures S2-S5, secondary production and energy transfer were higher in all ARs than that in the bare seafloor areas. In the benthic system, feeding flow to benthic bivalves and benthic sedentary worm in the ATRs was almost the same as that in bare seafloor areas at the beginning, as they were located in what could be, in principle, considered a fairly homogenous environment at their spatial scale. Nevertheless, slightly larger flow was simulated in the ATRs in the following years, when biomass in the ATRs was higher than in the bare seafloor areas.
Overall, organism biomass on the ATRsLB initially increased, but one year after post-deployment, it declined as the structure decayed. Photographs taken by divers confirmed that the leaves and branches had undergone degradation within approximately 18 months ( Figure S6). Such loss of epifaunal communities after the decay of substrate materials was also reported elsewhere [85,87]. In contrast, the ATRsOS were able to maintain their shape and the biomass of the organisms attached on them, although the fish biomass was slightly lower than that of the ATRsLB. Although the major occupants in fish species were the same consisting of black rockfish and Japanese surfperch in all ATRs, the size of them were small in the ATRsOS and ATRsLB compared to that in the simple ATRs. It is said that the number of living spaces in ARs has a positive correlation to fish production [88]. Our ATRs consisted of timber logs stacked in parallel cross formation and furnished to provide shelters for juvenile fish. Thus, ATRsOS and ATRsLB, with their additional materials, are more likely to encourage the proliferation of smaller fish due to the many small, interstitial spaces provided [89], and these refuge sites promote the aggregation of smaller fish in ARs [90].
Wintertime increases in feeding pressure due to the increase of fish biomass may have depleted the prey animal biomass, causing the observed low animal biomass in this compartment during winter. The model corroborated the observations in terms of the seasonal variations in the dominant fish species at all of the ATRs and bare seafloor areas. Both black rockfish and black scraper increased in numbers in winter and contributed to the decrease in prey animal biomass in all ATRs. The increase in feeding pressure by black seabream and multicolor rainbowfish resulted in the loss of shrimp and other amphipods in the simple ATRs, gammarid, bivalves, and other amphipods in the ATRsOS, and gammarid and other amphipods in the ATRsLB (Figure 3).
The dominant fish incorporated into the model were commercial species that are commonly captured by local fishermen in Mitsu Bay. The estimated annual fish catch in the ATR areas (4.4 kg m −3 year −1 ) was twice as high as compared to that in the bare seafloor areas (2.3 kg m −3 year −1 ). Our calculation was lower than those of large ARs (16 kg m −3 year −1 ) and regular-size ARs with 1-2-m −3 chamber modules (20 kg m −3 year −1 ) reported in previous studies [91,92]. The relatively low estimation in the present study may be explained by the short observation period after the settlement (the ecosystem has not matured) and other several reasons such as the omission of transient fish from the calculations and oligotrophic conditions in the bay.
The results of the present study indicated that the ATRs deployed in Mitsu Bay promoted the creation of food webs that established a much greater availability of prey organisms for a wide range of commercially important fish species. This finding was confirmed by the output of a numerical model using predator-prey relationship inputs. The ATR structure serves various purposes as a new habitat, especially for those with additional materials where the diversity indices of attached organisms were higher than those in the simple ATRs [1]. This indicates that leaves and branches create extra space, which could work as shelter for juvenile fish. Oyster shells function in the same manner but last much longer than leaves and branches. Oyster shells are by-products of local aquaculture. Their use, along with timber from forest thinning, for ATR construction is a good way to encourage a recycling-oriented society, especially in Japan, which has large forests and relatively low fishery production.
In our study, only 30 ATRs were deployed, which are not likely to increase fish stock dramatically in the bay. In practice, it will be necessary to increase the number of ATRs in order to achieve fish stock recovery [93]. This requires further calculation (i.e., sensitivity analysis) on how many ATRs will be needed to improve fisheries production significantly in the deployment area. A regular repeated deployment of ATRs will be needed because of the short lifespan of the materials. In our regular observations conducted 48 months after the deployments, some of the frames forming the ATR were broken down, as shown in Figure S7. As natural products, the decomposed timbers may act as a source of organic matter, which would work positively in the aspect of material cycling in oligotrophic waters as in Mitsu Bay. Increasing the area of ATR deployment may help increase fishery production in Mitsu Bay. This idea can be applicable to the whole Seto Inland Sea area where the fishery production has been decreasing because of oligotrophication [94]. As we mentioned in an earlier paragraph, the usage of coastal areas is strictly controlled by the central and local government in addition to fishermen's unions. It would be difficult to change the social system in a short amount of time. The fisheries law, which is the basic law describing the rules of fisheries in Japan, is planned to be revised in December 2020. Under the new law, considering the decrease in the number of fishermen, non-profit organizations (NPOs), non-governmental organizations (NGOs), and private companies will be allowed to join fishery activities under the condition of conserving the fisheries environments. At this stage, we will require proposals based on scientific knowledge, as in the present study, in collaboration with all involved sectors such as fishermen, forest keepers, local government, NPOs, NGOs, and volunteers, to increase fishery production while maintaining environmental conditions. Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/12/7/2013/s1, Figure S1: Artificial timber reefs (ATRs) used in the present study. (a) The appearance before deployment in Mitsu Bay, (b) at 1 month later, and (c) at 25 months later after deployment in the seawater. The net basket placed in the middle of frame was filled with additional material, oyster shells or leaves and branches. Table S1: Results of PBIAS calculated for the deviations of the calculation outputs from the observed values for each attached, benthic organism and fish. Figure S2: Annual average biomass stock (g m −3 ) (box) and flow (g m −3 d −1 ) (arrow) on simple artificial timber reefs for three years of 2016-2018. Abbreviations of stocks are same as those in Table 1. Figure S3: Annual average biomass stock (g m −3 ) (box) and flow (g m −3 d −1 ) (arrow) on artificial timber reefs with supplementary oyster shell material for three years of 2016-2018. Abbreviations of stocks are same as those in Table 1. Figure S4: Annual average biomass stock (g m −3 ) (box) and flow (g m −3 d −1 ) (arrow) on artificial timber reefs with supplementary leaf and branch material for three years of 2016-2018. Abbreviations of stocks are same as those in Table 1. Figure S5: Annual average biomass stock (g m −3 ) (box) and flow (g m −3 d −1 ) (arrow) in bare seafloor areas for three years of 2016-2018. Abbreviations of stocks are same as those in Table 1. Figure  S6: Decayed supplemented leaves and branches packed in a net basket after one and a half years of deployment. Figure S7: The condition of ATR that was broken down after 48 months of deployment.