Technical Comments

Response to Comment on “Designing river flows to improve food security futures in the Lower Mekong Basin”

See allHide authors and affiliations

Science  07 Jun 2019:
Vol. 364, Issue 6444, eaav9887
DOI: 10.1126/science.aav9887


Williams et al. claim that the data used in Sabo et al. were improperly scaled to account for fishing effort, thereby invalidating the analysis. Here, we reanalyze the data rescaled per Williams et al. and following the methods in Sabo et al. Our original conclusions are robust to rescaling, thereby invalidating the assertion that our original analysis is invalid.

Williams et al. (1) claim that the monitoring data we used to design flows that maximize harvest in the Dai fishery of Cambodia’s Tonle Sap River (2) were not properly scaled to account for both sampling and fishing effort, and hence the results are invalid. First, we note that our analysis was multifaceted, including historical analysis of the effects of dams on hydrologic variation, a historical analysis of the links between hydrology and harvest at the Dai fishery, and a forward-looking thought experiment about how dams could be managed to sustain or enhance harvest as they come online. Williams et al. address only the last two portions of the analysis. With regard to scaling of the monitoring data, Williams et al. contend that our analysis should have been performed following the methods of a technical report by Halls et al. (3). The implication is that, because monitoring inconsistencies and known biological variation were not accounted for, the entirety of the analysis is flawed.

We agree with Williams et al. that scaling of the Dai fishery monitoring data is an important consideration. However, it is impossible to say anything about the effect of scaling decisions on the results without performing the analysis presented in Sabo et al. (2) with the addition of scaling. Here, we present a reanalysis of Sabo et al. (2) using a scaling method that addresses the main considerations identified in Williams et al. This reanalysis demonstrates that the key results presented in (2) are robust.

The Tonle Sap Dai fishery consists of roughly 64 individual bagnets distributed among 14 “Dai rows” that are placed longitudinally in the Tonle Sap River. Catch at the Dais is not temporally or spatially homogeneous. Fish that are targeted by the Dai fishery migrate on a lunar cycle, with peak migration and peak catch typically about 6 days before the full moon. Downstream Dais should have lower catch per unit effort (CPUE) than upstream Dais as migrating fish are removed by the fishery upstream. Previous work has also used a classification system that identifies certain Dais as “high-yield” as a result of their position in the river. Finally, it is conceivable that fish abundance and catches would vary systematically across the course of a fishing season.

Catch at the Dai fishery has been monitored heterogeneously for a variety of reasons influenced by political conditions and local capacity. The Halls et al. method (3) attempts to scale the Dai monitoring data to correct for inconsistent spatial and temporal sampling in order to estimate absolute or total catch (kg) among two political districts in each fishing season. Thus, scaling is used to extrapolate from the subset (sample) to the total fishery yield over the full open season. Our analysis differs from this in two important ways: (i) We are interested in CPUE as an effort-independent metric of fish abundance. This involves fewer steps and assumptions than when estimating total catch as per Halls et al. (3). (ii) The subsequent multivariate autoregressive state space (MARSS) time-series analysis uses each Dai row as spatially explicit replicates, and thus we needed annual data by Dai row (n = 14). Therefore, we calculated a new average annual CPUE dataset that meets the constraints of our spatially replicated MARSS model and is adjusted for potential temporal biases due to migration timing.

This new approach considers the effect of lunar cycle but does not make any assumptions regarding productivity (high- versus low-yield classification of Dais) or the effect of period within the fishing season. Hence, we refer to the new dataset as the “lunar adjusted CPUE” dataset, and our former dataset presented in Sabo et al. (2)—catch summed across Dai rows within years, scaled by days of effort—as the “simple CPUE” dataset.

Survey data obtained from the Mekong River Commission included information on the single-day catch biomass for an individual Dai by species and identified whether a given sample occurred during the peak-migration period or low-migration period. That is, the data were previously scaled from catch per haul to catch per day by multiplying the average catch per haul by number of hauls, observed for the day of the survey. The Dai fishery operates continually 24 hours a day; thus, kg Dai–1 day–1 is a standardized metric of CPUE. In scaling the data to average CPUE of all fish for the season, the first step was to sum the total daily catch for each species into total single daily fish catch with units of kg Dai–1 day–1. Next, all samples from a Dai row (j) during each migration period (p) in a given season (i) were averaged irrespective of the month of the sample or whether the sampled Dai was categorized as high- or low-yield. The effect of migration period was accounted for by taking the peak- and low-migration CPUE within each Dai row and season and calculating a weighted average based on proportion of days in each season that are peak and low:average CPUEi,j=Pp=peak,iCPUEp=peak,i,j+Pp=low,iCPUEp=low,i,j(1)The proportion of fishing days that were peak-migration was determined on the basis of a lunar model (package lunar in R) that provides the moon phase for any location on Earth past and future. We used number of days when the moon phase was full in each season as the number of days of peak migration. This does not correspond to the same exact days of the peak migration, as this is typically 6 days prior to the full moon. Nonetheless, both the full-moon period and the peak-migration period last about a week, so we considered this a reasonable approximation of total number of days in a given season. Total fishing days was the number of days between 1 October and 31 March in each season (182 days; 183 days in leap years). The proportion of days that were peak-migration was roughly 25% of the fishing season. Because we are not scaling to total catch, it is not critical to know the exact number of peak-migration days, only the relative proportions. This method is a relatively unbiased way to estimate the effect of lunar migration that is evenly applied across years.

Plotting lunar adjusted CPUE across years by Dai row (Fig. 1) reveals the expected decline moving upstream (Dai 15) to downstream (Dai 2) mentioned by Williams et al. (1) that reflects the removal of fish by the fishery. The average annual lunar adjusted CPUE across Dai rows (Fig. 2) reflects substantial interannual variation. Annotated R code for the lunar adjusted CPUE calculation is provided via Figshare (4).

Fig. 1 CPUE among years 1997–2013 by Dai row.

Light gray lines are CPUE for individual years; the purple line is the average among years.

Fig. 2 Average CPUE across Dai rows for each year 1997–2013.

Year is reported as the year when the fishing season concluded (e.g., the 1996–1997 fishing season is coded as 1997).

The correlation between lunar adjusted CPUE and the data used in Sabo et al. (2) (simple CPUE) is high. Spearman correlation coefficients (based on ranked data within each Dai row) between the lunar adjusted and simple CPUE datasets averaged 0.68 (range 0.40 to 0.94; SD = 0.15) across 14 possible time-series pairs (one per Dai row). As in Sabo et al. (2), we fitted a density-independent MARSS model using historical hydrology as covariates [high-level drivers: flood pulse extent (FPExt) and net annual anomaly (NAA)] and the new, lunar adjusted CPUE data as variates. The bootstrapped MARSS model coefficients preserved their sign (i.e., direction) from the original analysis and increased slightly in effect size, and the confidence intervals (α = 0.05) continued not to overlap with zero (Table 1). Interestingly, the effects of the flood pulse (high flows) and the net annual anomaly (here, indicative of low flows) were both higher in magnitude, which suggests that these covariates may be identified as even more important with the lunar adjusted CPUE data. This variation and the FPExt and NAA covariate data were the key structural components of design flows in Sabo et al. (2). Hence, reconstructing the hypothetical future scenario presented in (2) with lunar adjusted CPUE data would result in a design flow with an even greater transition between high and low flow than with the simple CPUE dataset used in the original analysis. This expectation remains untested but underscores the need for a broader sensitivity analysis of methods (CPUE estimation) and shapes of design flows [as suggested originally in (2)].

Table 1 Parameter estimates for density-independent MARSS models on historical hydrology and CPUE data.

Results using the two datasets (Sabo et al. 2017 versus current approach) are compared.

View this table:

As noted by Williams et al. (1) and Halls et al. (3), the monitoring effort at the Dai fishery was inconsistent (i.e., variable frequency and location of sampling across Dais, seasons, and years). Counts of monitoring observations by year and Dai row for the low- and high-migration periods (Tables 2 and 3) show that the median number of observations in each year ranged from 62 to 488 (median = 217) and 37 to 422 (median = 229) for the low- and high-migration periods, respectively. There was relatively low overall sampling effort in the first 2 years of the time series and among Dai rows 13 to 15 (rows farthest from Phnom Penh). We therefore investigated the effect of this limited sampling effort on the results by repeating our MARSS analysis factorially, removing the first 2 years and the last three Dai rows. As before, the results remain robust whether or not these data are included (Table 4).

Table 2 Count of sampling events during the low migration period.

Year is reported as the year when the fishing season concluded (e.g., 1996–1997 fishing season is listed as 1997).

View this table:
Table 3 Count of sampling events during the high migration period.

Year is reported as the year when the fishing season concluded (e.g., 1996–1997 fishing season is listed as 1997).

View this table:
Table 4 MARSS reanalyses to test influence of years and Dai rows with low sampling effort.
View this table:

Our reanalysis with lunar adjusted CPUE data, as suggested by Williams et al. (1), preserves, and in fact strengthens, the central findings of Sabo et al. (2). Differences in how the Dai fishery monitoring data are adjusted for CPUE or total catch are ultimately a judgment call. The critical question, however, is whether these decisions influence the result of the time-series analysis in a meaningful way. The calculation of CPUE and lunar adjustment and subsequent reanalysis presented here found no substantive differences in the main results of (2), which suggests that trends in the data are strong relative to variation in scaling.

Finally, Williams and colleagues have repeatedly asserted (1, 5, 6) that Sabo et al. (2) advocate for dams to be built in the Mekong River. Holtgrieve et al. (7) rebutted this claim and we do so again here. The intent of our research is not to advocate for dam construction in the basin, which has been happening at a rapid pace in recent years regardless (8). Rather, we seek to understand the potential effects of dam reoperation for multiple purposes, including fisheries, and to explore opportunities to minimize negative impacts of existing dams. Given the ongoing development pressures, exploring options to mitigate the impacts of infrastructure development on food security is a pressing challenge in need of scientifically driven approaches like the one presented by Sabo et al. (2).


View Abstract


Navigate This Article