Difference between revisions of "IOTC Area Predictive analysis"

From D4Science Wiki
Jump to: navigation, search
(Description)
(Analysis)
Line 32: Line 32:
 
The bechmark for our analysis can be found [http://goo.gl/ykwIyI here]. The dataset was extracted from the reports by IOTC published on their [http://www.iotc.org/sites/default/files/documents/2014/01/CEPSBBFinal_0.zip website]. The dataset reports fishing activity in the Indian Ocean from January 1997 to December 2012. From this dataset, we selected only the long trail purse seine activities and aggregated the corresponding Yellowfin tuna catch reports. The following analyses have been obtained by using the D4Science [https://i-marine.d4science.org/group/biodiversitylab Statistical Manager].
 
The bechmark for our analysis can be found [http://goo.gl/ykwIyI here]. The dataset was extracted from the reports by IOTC published on their [http://www.iotc.org/sites/default/files/documents/2014/01/CEPSBBFinal_0.zip website]. The dataset reports fishing activity in the Indian Ocean from January 1997 to December 2012. From this dataset, we selected only the long trail purse seine activities and aggregated the corresponding Yellowfin tuna catch reports. The following analyses have been obtained by using the D4Science [https://i-marine.d4science.org/group/biodiversitylab Statistical Manager].
  
[[Image:barycentereffort.png|frame|center|Fishing hours aggregated per month along the entire period January 1997-December 2012]]
+
[[Image:barycentereffort.png|frame|center|Fishing hours aggregated per month along the entire period, from January 1997 to December 2012]]
  
 
Our first analysis aimed at discovering hidden periodicities in the exploitation of the ocean area. Thus, we extracted the trend of the fishing hours effort in the period between January 1997 and December 2012. We summed the monthly hours spent by superposing reports and eventually obtained the above trend. We associated the total fishing hours for one month to the barycenter of the areas interested by the fishing activity in that month. The monthly barycenters were calculated by weighting each 1 degree cell in the original dataset, with the fishing hours spent in that month inside the cell. The complete resulting dataset can be found [http://goo.gl/VzExKB here]. This contains unfolded information about the barycenters coordinates and their distances from the origin of the reference system.
 
Our first analysis aimed at discovering hidden periodicities in the exploitation of the ocean area. Thus, we extracted the trend of the fishing hours effort in the period between January 1997 and December 2012. We summed the monthly hours spent by superposing reports and eventually obtained the above trend. We associated the total fishing hours for one month to the barycenter of the areas interested by the fishing activity in that month. The monthly barycenters were calculated by weighting each 1 degree cell in the original dataset, with the fishing hours spent in that month inside the cell. The complete resulting dataset can be found [http://goo.gl/VzExKB here]. This contains unfolded information about the barycenters coordinates and their distances from the origin of the reference system.
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:barycenterlongitude.png|thumb|center|upright=2.5|Distribution of the vessels longitudes for Yellofin Tuna by Purse Seine technique, in the entire period January 1997-December 2012]]
+
| [[Image:barycenterlongitude.png|thumb|center|upright=2.5|Distribution of the variations of the vessels longitudes for Yellowfin tuna catch by purse seining]]
| [[Image:barycenterlatitude.png|thumb|center|upright=2.5|Distribution of the vessels latitudes for Yellofin Tuna by Purse Seine technique, in the entire period January 1997-December 2012]]
+
| [[Image:barycenterlatitude.png|thumb|center|upright=2.5|Distribution of the variations of the vessels latitudes for Yellowfin tuna catch by purse seining]]
 
|}
 
|}
  
Line 44: Line 44:
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:barycentercatch.png|thumb|center|upright=2.5|Catch quantity in tonnes for Yellofin Tuna, in the entire period January 1997-December 2012]]
+
| [[Image:barycentercatch.png|thumb|center|upright=2.5|Catch quantity trend in tonnes for Yellowfin tuna catch.]]
| [[Image:barycentercpue.png|thumb|center|upright=2.5|CPUE for Yellofin Tuna, in the entire period January 1997-December 2012]]
+
| [[Image:barycentercpue.png|thumb|center|upright=2.5|CPUE trend for Yellowfin tuna catch.]]
 
|}
 
|}
  
Line 51: Line 51:
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:barycenterspectrogramlongitude.png|thumb|center|upright=2.5|Spectrogram for the longitude trend.]]
+
| [[Image:barycenterspectrogramlongitude.png|thumb|center|upright=2.5|Spectrogram of the longitude trend.]]
 
| [[Image:Barycenterlongitudesamples.png|thumb|center|upright=2.5|Longitude trend in samples.]]
 
| [[Image:Barycenterlongitudesamples.png|thumb|center|upright=2.5|Longitude trend in samples.]]
 
|-
 
|-
| [[Image:barycenterspectrogramlatitude.png|thumb|center|upright=2.5|Spectrogram for the latitude trend.]]
+
| [[Image:barycenterspectrogramlatitude.png|thumb|center|upright=2.5|Spectrogram of the latitude trend.]]
 
| [[Image:Barycenterlatitudesamples.png|thumb|center|upright=2.5|Latitude trend in samples.]]
 
| [[Image:Barycenterlatitudesamples.png|thumb|center|upright=2.5|Latitude trend in samples.]]
 
|}
 
|}
Line 62: Line 62:
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
 
| [[Image:barycentereffortsamples.png|thumb|center|upright=2.5|Fishing hours trend in samples.]]
 
| [[Image:barycentereffortsamples.png|thumb|center|upright=2.5|Fishing hours trend in samples.]]
| [[Image:barycenterspectrogrameffort.png|thumb|center|upright=2.5|Spectrogram for the fishing hours trend.]]
+
| [[Image:barycenterspectrogrameffort.png|thumb|center|upright=2.5|Spectrogram of the fishing hours trend with 3 half-superposing windows.]]
 
|-
 
|-
 
|}
 
|}
Line 69: Line 69:
 
The analysis detected a weak overall periodicity of 12.91 samples with indecision interval between 13.25 and 12.59 samples. This indicates an annual periodicity in the variation of the fishing effort.
 
The analysis detected a weak overall periodicity of 12.91 samples with indecision interval between 13.25 and 12.59 samples. This indicates an annual periodicity in the variation of the fishing effort.
  
[[Image:barycentercsquares.png|frame|center|Distribution of the presence of vessels on 1 degree csquares, in the entire period January 1997-December 2012]]
+
[[Image:barycentercsquares.png|frame|center|Distribution of the presence of vessels on 1 degree square areas.]]
  
 
We also report the trend of the [http://en.wikipedia.org/wiki/C-squares Csquare] codes associated to the barycenters. Each code was assigned to an index and the above chart displays the distribution of such indexes in time. We adopted this representation in order to report a bidimensional phenomenon on one dimension. The trend of the indexes is periodic and this is due to the periodic trend in longitudes and latitudes.
 
We also report the trend of the [http://en.wikipedia.org/wiki/C-squares Csquare] codes associated to the barycenters. Each code was assigned to an index and the above chart displays the distribution of such indexes in time. We adopted this representation in order to report a bidimensional phenomenon on one dimension. The trend of the indexes is periodic and this is due to the periodic trend in longitudes and latitudes.
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:barycenterforecastlatitude.png|thumb|center|upright=2.5|Forecasted trend for latitude.]]
+
| [[Image:barycenterforecastlatitude.png|thumb|center|upright=2.5|Forecasted trend for vessels latitudes on 12 further months.]]
| [[Image:barycenterforecastlongitude.png|thumb|center|upright=2.5|Forecasted trend for longitude.]]
+
| [[Image:barycenterforecastlongitude.png|thumb|center|upright=2.5|Forecasted trend for vessels longitudes on 12 further months.]]
 
|-
 
|-
| [[Image:barycentereffortforecast.png|thumb|center|upright=2.5|Forecasted trend for fishing hours.]]
+
| [[Image:barycentereffortforecast.png|thumb|center|upright=2.5|Forecasted trend for fishing hours on 24 further months.]]
 
|}
 
|}
  
Line 86: Line 86:
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:month12effort.png|thumb|center|upright=2.5|Aggregation on 12 months of the number of Purse Seine fishing hours.]]
+
| [[Image:month12effort.png|thumb|center|upright=2.5|Aggregation of fishing hours over 12 months.]]
| [[Image:month12catch.png|thumb|center|upright=2.5|Aggregation on 12 months of the catch quatities.]]
+
| [[Image:month12catch.png|thumb|center|upright=2.5|Aggregation of Yellowfin tuna catch quantities in tonnes over 12 months.]]
 
|-
 
|-
| [[Image:month12cpue.png|thumb|center|upright=2.5|Aggregation on 12 months of the CPUE for Yellfin Tuna.]]
+
| [[Image:month12cpue.png|thumb|center|upright=2.5|Aggregation of Yellowfin tuna CPUE over 12 months.]]
| [[Image:month12latitudes.png|thumb|center|upright=2.5|Aggregation on 12 months of the barycenters latitudes.]]
+
| [[Image:month12longitudes.png|thumb|center|upright=2.5|Report of barycenter longitudes over 12 months.]]
 
|-
 
|-
| [[Image:month12longitudes.png|thumb|center|upright=2.5|Aggregation on 12 months of the barycenters longitudes.]]
+
| [[Image:month12latitudes.png|thumb|center|upright=2.5|Report of barycenter latitudes over 12 months.]]
 
|}
 
|}
  
Line 131: Line 131:
 
| 33038.7
 
| 33038.7
 
|-
 
|-
|+ Most exploited areas at 10 degrees resolution
+
|+ Most exploited areas at 10 degrees resolution, with respect to the total amount of fishing hours.
 
|}
 
|}
  
Line 137: Line 137:
  
 
{|style="margin: 0 auto;"
 
{|style="margin: 0 auto;"
| [[Image:month12_3005cpue.png|thumb|center|upright=2.5|Aggregation on 12 months of the CPUE for the most exploited area, with csquare code 3005.]]
+
| [[Image:month12_3005catch.png|thumb|center|upright=2.5|Aggregation of Yellowfin tuna catch in tonnes over 12 months, in the csquare with code 3005.]]
| [[Image:month12_3005fishinghours.png|thumb|center|upright=2.5|Aggregation on 12 months of the fishing hours for the most exploited area, with csquare code 3005.]]
+
| [[Image:month12_3005fishinghours.png|thumb|center|upright=2.5|Aggregation of fishing hours over 12 months, in the csquare with code 3005.]]
 
|-
 
|-
| [[Image:month12_3005catch.png|thumb|center|upright=2.5|Aggregation on 12 months of the catch in tonnes for the most exploited area, with csquare code 3005.]]
+
| [[Image:month12_3005cpue.png|thumb|center|upright=2.5|Aggregation of Yellowfin tuna CPUE over 12 months, in the csquare with code 3005.]]
 
|-
 
|-
 
|}
 
|}
  
 
In the case of the firstmost exploited area, a periodic trend of 3\4 months is visible in the CPUE chart. The chart also shows that the catch per effort is lower and lower. This scenario indicates that higher catch has been obtained with less hours of fishing, during the last years.
 
In the case of the firstmost exploited area, a periodic trend of 3\4 months is visible in the CPUE chart. The chart also shows that the catch per effort is lower and lower. This scenario indicates that higher catch has been obtained with less hours of fishing, during the last years.

Revision as of 22:08, 21 May 2014

Description

The Indian Ocean Tuna Commission (IOTC) is an intergovernmental organization that coordinates the regulation and management of tuna in the Indian Ocean. The commission publishes data about fishing activity, to be used by fishery managers for monitoring catch effort and ocean exploitation.

The experiment described in this page, focuses on the purse seine fishing activity, a very efficient way of catching fish that can produce large amounts of revenue in a short period of time. On the other side, this kind of fishing practice can also cause damage in a very short time. Unregulated purse seine activity would result in non-sustainable fisheries resources. Purse seining, in fact, is a non-selective fishing method that captures everything that it surrounds, including protected species. Examples are sea turtles and marine mammals that can be captured by catch.

Thus, monitoring and regulating purse seining is important and the example in this page reports an analysis to detect, monitor and predict such activity.

Hypothesis and Thesis

Premise

  • The analysis benchmark is made up of a report about purse seine fishing activity in the Indian Ocean involving Yellowfin tuna.

Hypothesis

  • The exploitation trends and the locations involved in fishing activity could be periodic and could be forseen.

Theses

  • Basic statistical analysis, trends displaying and signal processing can detect periodic properties;
  • Forecasting techniques can predict future exploitation.

Outcome

Starting from IOTC reports about purse seine fishing activity, in the Indian Ocean from 1997 to 2012, we detected the following phenomena:

1 an overall periodic, yet decreasing, number of Purse Seine fishing hours in the Indian Ocean;
2 a periodic trend in the visitation of fishing locations by vessels, with periodicity around 51 months;
3 an overall periodic variation of 12 months in the latitudes of the purse seine activity;
4 forecasted locations for the fishing areas in the following years;
5 a list of the most exploited areas at 10 degrees resolution;
6 an overall periodic reduction of the Catch per Unit Effort (CPUE) in the most exploited 10-degrees area.

Analysis

The bechmark for our analysis can be found here. The dataset was extracted from the reports by IOTC published on their website. The dataset reports fishing activity in the Indian Ocean from January 1997 to December 2012. From this dataset, we selected only the long trail purse seine activities and aggregated the corresponding Yellowfin tuna catch reports. The following analyses have been obtained by using the D4ScienceAn e-Infrastructure operated by the D4Science.org initiative. Statistical Manager.

Fishing hours aggregated per month along the entire period, from January 1997 to December 2012

Our first analysis aimed at discovering hidden periodicities in the exploitation of the ocean area. Thus, we extracted the trend of the fishing hours effort in the period between January 1997 and December 2012. We summed the monthly hours spent by superposing reports and eventually obtained the above trend. We associated the total fishing hours for one month to the barycenter of the areas interested by the fishing activity in that month. The monthly barycenters were calculated by weighting each 1 degree cell in the original dataset, with the fishing hours spent in that month inside the cell. The complete resulting dataset can be found here. This contains unfolded information about the barycenters coordinates and their distances from the origin of the reference system.

Distribution of the variations of the vessels longitudes for Yellowfin tuna catch by purse seining
Distribution of the variations of the vessels latitudes for Yellowfin tuna catch by purse seining

The figures above display the trends for the longitudes and latitudes of the monthly barycenters. The signals show a periodic behaviour, which we investigate more accurately later in this page.

Catch quantity trend in tonnes for Yellowfin tuna catch.
CPUE trend for Yellowfin tuna catch.

The above figures, instead, report the trends for the amount of catch in tonnes for the Yellowfin Tuna, and the related catch per unit of effort (CPUE) calculated as the ratio between the catch and the fishing hours.

Spectrogram of the longitude trend.
Longitude trend in samples.
Spectrogram of the latitude trend.
Latitude trend in samples.

In order to inspect the periodic behaviour of the trends of the coordinates, we applied Fourier analysis to the signals and traced spectrograms. For both longitude and latitude, the analysis detected a period of 51.0 samples with indecision interval [56.64 ; 46.38]. This corresponds to a periodic exploitation of areas after approximately 51 months.

Fishing hours trend in samples.
Spectrogram of the fishing hours trend with 3 half-superposing windows.

Fishing effort in hours, instead, is not a stationary signal, thus we tried to inspect its hidden periodicities by means of a Short-Time Fourier analysis. The analysis was made on 3 chunks of the signal presenting quite stationary means and variances. The analysis detected a weak overall periodicity of 12.91 samples with indecision interval between 13.25 and 12.59 samples. This indicates an annual periodicity in the variation of the fishing effort.

Distribution of the presence of vessels on 1 degree square areas.

We also report the trend of the Csquare codes associated to the barycenters. Each code was assigned to an index and the above chart displays the distribution of such indexes in time. We adopted this representation in order to report a bidimensional phenomenon on one dimension. The trend of the indexes is periodic and this is due to the periodic trend in longitudes and latitudes.

Forecasted trend for vessels latitudes on 12 further months.
Forecasted trend for vessels longitudes on 12 further months.
Forecasted trend for fishing hours on 24 further months.

Finally, we tried to forecast the coordinates trends in order to understand which places could be interested by exploitation in the 12 following months. To such aim, we used Singular Spectrum Analysis. The forecasts confirm the periodic nature of the coordinates signals. Furthermore, the forecast of the fishing hours effort, suggests that the amount of hours per month is going to oscillate around lower and lower values in average. Thus, the effort in purse seine is likely to be reduced.

Aggregation on 12 Months

Aggregation of fishing hours over 12 months.
Aggregation of Yellowfin tuna catch quantities in tonnes over 12 months.
Aggregation of Yellowfin tuna CPUE over 12 months.
Report of barycenter longitudes over 12 months.
Report of barycenter latitudes over 12 months.

We repeated the previous analysis, this time aggregating the fishing hours on 12 months. The figures above display the trends we extracted after the aggregation. The periodic nature of the fishing activity is no more visible, except for the latitudes of the barycenters. For this signal, a periodic behaviour of 12 months is visible, which demonstrates that the latitude periodic behaviour is stronger for latitude than for longitude.

Analysis on the Most Exploited Csquare

Csquare code Exploitation (fishing hours)
3005 614135.3
1005 709254.3
3004 513798.8
3104 211621.1
3006 258669.6
3105 24718.3
1004 12235.1
3004 513798.8
3007 33038.7
Most exploited areas at 10 degrees resolution, with respect to the total amount of fishing hours.

As last analysis, we investigated the exploitation of the areas involved in the catch. We selected a coarse granularity of 10 degrees and ranked the most exploited square areas at this resolution. The trends and the data for each single square can be found here. The above table lists the most exploited squares, in Csquare encoding.

Aggregation of Yellowfin tuna catch in tonnes over 12 months, in the csquare with code 3005.
Aggregation of fishing hours over 12 months, in the csquare with code 3005.
Aggregation of Yellowfin tuna CPUE over 12 months, in the csquare with code 3005.

In the case of the firstmost exploited area, a periodic trend of 3\4 months is visible in the CPUE chart. The chart also shows that the catch per effort is lower and lower. This scenario indicates that higher catch has been obtained with less hours of fishing, during the last years.