(Superscript numbers refer to the appended annotations)
Note:
Regions are non-overlapping major ocean areas. For species found in or migrating to higher latitudes, these will normally be the Arctic and adjacent waters, the North Atlantic and adjacent waters, the North Pacific and adjacent waters, and the Southern Hemisphere. For species confined to lower latitudes, the Regions will normally be the Atlantic, Pacific and Indian Oceans. Regions can be combined for species where the interchange is not negligible.
Small Areas are disjoint areas small enough to contain whales from only one biological stock, or be such that if whales from different biological stocks are present in the Small Area, catching operations would not be able to harvest them in proportions substantially2 different to their proportions in the Small Area.
Medium Areas3 correspond to known or suspected ranges of distinct biological stocks.
Large Areas4 coincide with Regions, unless evidence exists to support the selection of one or more areas smaller than a Region which fully covers the range of some biological stocks of a species and definitely excludes whales from all other biological stocks of that species in the Region.
Residual Areas5 are all geographical areas in a Region which are outside any Small Areas. Medium Areas comprise unions of Small and, where identified, Residual Areas. Large Areas comprise unions of Medium and, where identified, Residual Areas.
Combination Areas are disjoint unions of Small Areas to which the Catch Limit Algorithm is applied when Catch-cascading is used.
Management Area is a generic term denoting a Small, Medium, Large, Residual or Combination Area.
Catch Limit Algorithm is the process (described in Section 4) that is used to calculate a catch limit for a Management Area.
Years6 are consecutive periods of 12 months used for the compilation of time series of catches and abundance data for application of the Catch Limit Algorithm.
Catch-cascading7 is the process by which a catch limit calculated for a Combination Area is distributed among the Small Areas that make up the Combination Area in proportion to the calculated relative abundances in those Small Areas. When Catch-cascading occurs, the relative abundances for Small Areas within the Combination Area shall normally be calculated from the same estimates of absolute abundance as were used for the application of the Catch Limit Algorithm to the Combination Area. The calculated relative abundance in a Small Area shall be an appropriate form of weighted average of the available abundance indices for that Small Area, with the statistically appropriate weighting, except that each estimate shall also be further weighted by the factor 0.9,n where n is the number of years that have elapsed between the Year to which the estimate refers and the Year of the Catch Limit Calculation.
Catch-capping8 is the process by which Catch Limits calculated for Small Areas are adjusted by reference to those calculated for either Medium or Large Areas containing those Small Areas. It consists of the following rules. If the sum of the catch limits calculated for those Small Areas that make up a Medium (or Large) Area exceeds the catch limit calculated for the Medium (or Large) Area, then both the Small and Medium (or Large) Area catch limits shall apply in such a way that the maximum catch allowed in each Small Area is the appropriate Small Area catch limit and the maximum catch allowed in the Medium (or Large) Area is the Medium (or Large) Area catch limit. This definition does not preclude the possibility of applying Catch-capping to overlapping Medium Areas.
An Implementation involves the designation of the Management Areas and their boundaries and the selection of Catch-cascading and/or Catch-capping options for a particular species and Region. These designations and/or selections may be changed in a subsequent Implementation Review.
A Catch Limit Calculation is the process by which catch limits for a species in a Region are calculated for all Small (and where appropriate Medium or Large) Areas within that Region, as specified in Sections 3.3, 3.4 and 3.5, by application of the Catch Limit Algorithm as described in Section 4. This algorithm uses historic catch data and estimates of absolute abundance for each Management Area that meet the requirements of Section 3.2.
Implementations and Implementation Reviews are conducted by the Scientific Committee on a Regional basis. They involve the delineation of Small Areas and, where appropriate, Medium and Large Areas. A selection between possible options for the Catch-cascading and/or Catch-capping is made during an Implementation (Review), which includes the designation of Combination Areas as may be appropriate. This process is described as an Implementation on the first occasion it takes place for a species in a Region; subsequent revisions are termed Implementation Reviews.9 An Implementation (Review) shall take account of the available biological and operational data, including in particular those data pertaining to stock identity. An Implementation (Review) is conducted by species or other suitable taxonomic unit below the species level.10 Such taxonomic units should be treated separately for the purpose of Catch Limit Calculations (see Section 3) where the extent of geographical separation is sufficient to make this feasible. In the following text, 'species' should be taken to refer to taxonomic units below the species level where appropriate.
3.1 Scope and period of validity
Catch limits pertain to the first Year commencing after their
calculation by the Scientific Committee, and for each of the following four
Years.11
A catch limit is calculated for each Small Area in a Region for
each of these five Years.
The five catch limits calculated for each Management Area shall be
equal, except where adjustments are made under the phaseout rule specified in
section 3.4. A Catch Limit Calculation involves the (re)calculation of
catch limits for all Small Areas and, where appropriate, Medium
or Large Areas in the Region.
At the request of the Commission, the first of these catch limits calculated
may alternatively refer to the Year in which the calculation takes
place, and for each of the following four Years.
3.2 Data requirements12
3.2.1 Catch history
Time series of catches by sex shall be compiled for each of the
Management Areas specified within the region, using the best available
information.
These catch histories shall cover a period beginning not later than the
Year of the first recorded or estimated13 catch and ending
with the Year preceding
the first Year for which catch limits are to be calculated.14
If there are catches known to have occurred in the Region, but the Small Area in which they were taken is not known, they shall be assigned to the Small Area in which they are considered most likely to have been taken. Pro rata allocations are allowed. Where the sex ratio of catches is not accurately known, the best available estimate of the sex ratio shall be used to divide the catches; in the absence of any information, a 50:50 sex ratio shall be assumed. Unspecified catches of whales shall be allocated to species using the best available information on the species composition of the catch.15 Known or estimated numbers of whales struck and lost shall be added to the catches. If the timing of catches is uncertain, they shall be assigned to Years according to the best available information. No catches known to have occurred in the Region shall be omitted from the Catch Limit Calculation on grounds of uncertainty over their location, timing, sex ratio or other details. All known removals16 from a Region shall be included in the catch series.
3.2.2 Absolute abundance estimates
Absolute abundance data to be used in the calculation of catch limits shall
have been obtained by direct methods,17 such as sightings surveys,
and collected and analysed using methods approved by the Scientific Committee.
Management Areas to which the Catch Limit Algorithm is applied
should normally be surveyed at intervals not exceeding six years.
The methods shall be such as to provide estimates of whale abundance that have
acceptable levels of bias and precision.
They shall also permit estimation of the variance of each estimate and of their
variance-covariance matrix, or alternative variance-related statistics where
appropriate.
Data for any sightings survey18 to be used to calculate abundance estimates for the purposes of conducting a Catch Limit Calculation shall be documented and provided to the Secretariat in computer readable data files before a specified time in advance of the Scientific Committee meeting during which the data are to be used. All such data should be archived by the Secretariat in an appropriate database such that abundance estimates can be calculated for any specified Small Area. Data should be in a fully disaggregated form so that estimates can be recalculated appropriately if the boundaries of Management Areas are altered. Once lodged with the Secretariat, these data shall be available to accredited scientists as defined in the Scientific Committee's Rules of Procedure.
Estimates of absolute abundance are required for each Management Area to which the Catch Limit Algorithm is to be applied under the procedures described in Section 3.3.19 For each such Management Area, a time series of absolute abundance estimates shall be calculated, along with an estimate of their variance-covariance matrix, or alternative variance-related statistics where appropriate. The approximate distributional properties of the abundance estimates shall also be determined. Care should be taken to avoid substantially underestimating the variance (or alternative variance-related statistic) of each abundance estimate used for input into the Catch Limit Algorithm.20
The absolute abundance estimate for a given Year should ideally be calculated from data collected in that Year. Data collected in different Years may be used, for example to account for parts of the area that were not covered in that Year, to pool results from surveys conducted over consecutive or nearly consecutive Years in order to reduce variance, or to provide estimates of calibration factors, provided that appropriate statistical methods are used.21
Data from surveys conducted in different Years or at different times of year may only contribute to a single abundance estimate if adequate precautions are taken to avoid substantial double counting of whales due to migration or other factors. In the calculation of an absolute abundance estimate for a Management Area in a given Year, parts of the Area for which there are no absolute abundance estimates available at any time meeting the above specifications shall be treated as having an absolute abundance of zero.
The absolute abundance estimates should pertain to the total number of whales aged one year and above in the Management Area, regardless of any size limits that may be in force or the selectivity or otherwise of any past or present exploitation.22 Animals aged less than one year shall be excluded where possible.
3.3 Options for determination of catch limits
Catch limits shall always be set at the Small Area level and they shall
be set for each Small Area in a Region.
In addition, where Catch-capping is invoked at the Medium or
Large Area level, corresponding catch limits will be set for those
Medium or Large Areas.
Catch limits for all Residual Areas within a Region shall be set
at zero.
Catch limits for the total number of whales that may be taken in a season in each Small Area will be calculated by:
Catch limits for the total number of whales that may be taken in a Year in Medium or Large Areas, as required when Catch-capping is invoked, will be calculated by application of the Catch Limit Algorithm to those Medium or Large Areas.
The decision for any particular species or Region on whether or not Catch-capping is to be applied, and if so whether it should be applied at the Medium or Large Area level, and whether or not Small Areas are to be combined for the purposes of Catch-cascading, will be made on the basis of biological evidence available to the Scientific Committee, and, where necessary, the results of computer simulation trials23 conducted by the Scientific Committee. Where computer simulation trials are carried out, they shall, as far as possible, encompass the full range of plausible hypotheses (regarding, for example, stock identity) consistent with existing biological data.
3.4 Phaseout rule
The catch limits for a Small Area calculated under Section 3.3 shall be
adjusted downwards when the time series of absolute abundance estimates used
for the application of the Catch Limit Algorithm to the
Small Area (or, if Catch-cascading has been applied, to the
Combination Area containing it) does not include an absolute abundance
estimate pertaining to a Year not more than eight years24
prior to the Year to which the catch limit pertains.
Under these circumstances, the catch limit for the Small Area shall be
reduced by 20% of the unadjusted catch limit for that Small Area and
Year for each year in excess of eight years that has or will have
elapsed since the Year of the most recent such abundance
estimate.25
This rule shall also be invoked in a Small Area included in a
Combination Area for Catch-cascading if the data used for the
derivation of absolute abundance estimates for input to the
Catch Limit Algorithm do not contain any survey effort in that
Small Area within this eight year period.
3.5 Adjustments for recent sex ratios in the catch
If the proportion, Pf, of female whales in the total catch
taken
from a Small Area in the most recent five Years prior to the
Catch Limit Calculation for which the catch data are available exceeds
50% , the catch limits for the Small Area calculated according to the
procedure described in sections 3.3 and 3.4 shall be adjusted downwards by the
ratio 0.5/Pf.26
However, should the Scientific Committee decide it to be more appropriate,
this adjustment ratio shall be determined from the proportion of females in
the total catch taken from a union of Small Areas, and applied to the
catch limit for each Small Area in the union.
Further, the sex ratio adjustment shall be waived if the Scientific Committee
agrees that the catches taken in the most recent five Years for which
the catch data are available are too few to provide a useful indication of the
expected future sex ratio of the catch.
The nominal catch limit for a Management Area shall be calculated using the algorithm defined below if at least one estimate of absolute abundance as defined in Section 3.2 is available for the Area in question. Otherwise, the nominal catch limit for the Management Area shall be zero.
4.1 Input data
The input data for application of the Catch Limit Algorithm for any
Management Area shall include the time series of annual catches as
detailed in Section 3.2.1 and the time series of absolute abundance estimates,
along with their variance-covariance matrix or other appropriate
variance-related statistics and a specification of the distributional form of
the absolute abundance estimates, as specified in Section 2.2.2.
4.2 Population model
The following population dynamics model27 shall be used:
Provided there have been at least some catches, the population dynamics model
is fully determined when the catch series and the values of
DT and 'mu' are specified.
If there have been no catches, a nominal catch of one whale in Year
zero is assumed.
4.3 Fitting of the model
The annual absolute abundance estimate (if there is one) for each
Year t, is assumed to have expectation bPt where
b is the bias parameter.
The joint likelihood function of the parameters b, DT
and 'mu' is determined using the absolute abundance estimates, the
variance-covariance matrix of the absolute abundance estimates (or alternative
variance-related statistics where appropriate) and information on their
distributional form.
Unless there are specific indications to the contrary,29 the
absolute abundance estimates shall be assumed to be lognormally distributed
with a variance-covariance matrix of the log estimates to be estimated from
the data using methods judged appropriate by the Scientific Committee.
In this case, the formula for the likelihood is:
The stock depletion parameter DT is assigned a prior
probability distribution30 that is uniform from zero to one, and
zero outside this range.
The productivity parameter 'mu' is assigned a prior probability distribution30 that is uniform from zero to 0.05, and zero outside this range.
The bias parameter b is assigned a prior probability distribution30 that is uniform from zero to 5/3, and zero outside this range.
The above three prior distributions are treated as independent and combined accordingly to determine the joint prior distribution of the parameters DT, 'mu' and b.
The joint 'posterior' distribution of the parameters DT,
'mu' and b is defined as follows:
where s, the scale parameter, is set equal to 1/16.
The presence of the scale parameter represents an intended deviation from a
strictly Bayesian approach.
4.4 The catch control law
The internal catch limit, LT, is the following function of
DT, 'mu' and PT:
The marginal posterior distribution of LT is obtained by
integration of the joint posterior distribution of (DT,
'mu', b).
This requires that, for each value of LT, the joint
posterior distribution of (DT, 'mu', b) is to be
integrated over the subset of parameter space that corresponds to that value
of LT.
The nominal catch limit is equal to the lower 41.02 percentile of the marginal
posterior distribution of LT.31
4.5 Computation
All steps in the above algorithm for the calculation of the nominal catch
limit shall be performed using a computer program validated by the IWC
Secretariat and with sufficient numerical accuracy that the calculated nominal
catch limit is numerically accurate to within one whale.
The development of the RMP has been a long and difficult task, involving a wide range of scientific and technical issues and a thorough and extensive testing process. The Scientific Committee has recommended a protocol for evaluating amendments to the RMP which is given in Rep. int. Whal. Commn 44:47-8 (Item 7.1.4).
An example of the calculation involved is as follows. If the absolute abundance estimates are treated as being lognormally distributed, then the relative abundance for a Small Area would normally be calculated using the following formula.
Let:
a = vector of log abundance estimates in the Small Area;
ti = difference between the current Year and the Year
of the ith estimate;
F = information matrix of a.
If F is non-singular, F = V-1 where V is the variance-covariance matrix of a. G is the matrix such that Gij = Fij(0.9) -(ti+tj)/2
The relative abundance in the Small Area is given by:
Catch-capping Can be applied together with Catch-cascading. In this case, after the Small Area catch limits have been calculated under Catch-cascading, the capping option is invoked.
In some cases, it may be appropriate to carry out an Implementation Review less than 5 years after the previous Implementation (Review). This would be warranted, for example, if important new evidence on stock identity becomes available, if major advances are made in methodology of calculating absolute abundance estimates, if major changes occur in the areas covered by the abundance surveys, or if other evidence becomes available to the Scientific Committee suggesting that the premises on which the previous Implementation (Review) was conducted are no longer appropriate.
The RMP has been demonstrated to be robust to considerable uncertainties in catch histories in single stock robustness trials (Rep. int. Whal. Commn 42:272).
Note that the above does not preclude the use of estimates of relative abundance during Catch-cascading (see annotation 7) or in analysing abundance data collected in different Years (see Section 3.2 and annotation 2l).
In some circumstances, the best available estimates of absolute abundance may come from mark-recapture analyses, e.g. those resulting from photo-identification studies. The properties of such estimates, and the implications of these with respect to possible uncertain stock identity and migration patterns need to be evaluated before estimates of abundance based on them may be used when implementing the RMP for a particular species and Region. Until this is done, sightings surveys or other direct methods of estimation with similar statistical properties remain the primary tools for obtaining suitable estimates of absolute abundance for Catch Limit Calculations.
Note:
Process error reflects the extent to which abundance estimates from repeat
surveys of the same area in successive years will vary more than would be
expected on the basis of the observation error alone, for example due to
variations in the numbers of whales moving into or out of the survey area.
Observation error is the sampling error arising from the survey methods and
design.
The level of observation error is inversely related to the amount of survey
effort, provided that the survey is well designed.
Where simulation trials are used during implementation to evaluate the appropriateness or otherwise of Catch-cascading and/or Catch-capping, and in the latter case whether at Medium or Large Area level, judgements will be based on comparisons of performance of the different options against a base case where catch limits are calculated and set by Small Area only. The addition of Catch-capping to other options leads to the setting of catch limits lower than or equal to those which would be set in the absence of Catch-capping. Where the performance of suitable simulation trials of the base case option for setting catch limits is satisfactory in terms of statistics related to lowest and final depletion levels, it would not normally be judged necessary to invoke Catch-capping ('depletion' is defined in Section 4.2). However, where the performance of the base case option is judged unsatisfactory in terms of the depletion statistics, and this is rectified when one of the Catch-capping options is used, Catch-capping at the relevant level shall be invoked.
Catch-cascading normally leads to higher catch limits than the base case option. Accordingly, Catch-cascading may only be invoked when simulation trials show that it does not lead to unsatisfactory performance on depletion statistics related to lowest and final depletions.
Examples of examination of these issues in the context of potential implementation of the RMP to Southern Hemisphere and North Atlantic minke whales are given in Annexes E and F of the 1992 Report of the Scientific Committee (Rep. int. Whal. Commn 43:104-14 and 115-29).
132 x 0.5/Pf = 132 x 0.5/0.6 = 110 whales per annum
Note that the aim of the Catch Limit Algorithm in setting the pre-adjustment catch limit is that this comprise equal numbers (66 in this case) of males and females. The intent of the adjustment is that no more than 66 females will be caught: if the female proportion remains at 0.6, this will be achieved exactly by the adjustment process because 0.6 x 110 = 66.
Although several factors contribute to the variance of an estimate of absolute abundance, the variance is dominated by the variance in the number seen when the number of sightings is very low. The variance of the number of sightings will be at least as high as the variance of a random variable with a Poisson distribution with expectation equal to the expectation of the number of sightings. The number of sightings refers to the number of schools or groups, rather than to individual animals.
The expected number of sightings, E(n), is proportional to the true
absolute abundance, P:
The parameter 'alpha' represents the estimate of absolute abundance that would
have been obtained had there been exactly one sighting.
This will be a function of the survey effort, the size of the area, and survey
parameters that may need to be estimated by adopting values from similar
surveys.
Ignoring the variance of 'alpha', the likelihood of the zero estimate of
absolute abundance is the following function of the true absolute abundance:
Since the only covariance between the absolute abundance estimate and other
absolute abundance estimates is that due to the 'alpha' parameter, whose
variance is being ignored, the joint likelihood function of the zero estimate
of absolute abundance and the remaining estimates is taken to be the product
of the respective likelihood functions.
The information about the zero estimate of absolute abundance that needs to be supplied to the Catch Limit Algorithm is: (i) the Year of the zero estimate; (ii) the fact that it is a zero estimate; and (iii) the value of the 'alpha' parameter. The computer program implementing the Catch Limit Algorithm that has been validated by the Secretariat has the facility to handle zero estimates of absolute abundance in this manner. P is identified with the simulated population size generated by the Catch Limit Algorithm's internal calculations.
Since the treatment above ignores some contributions to the variance of a zero estimate of absolute abundance, it assigns more weight to a zero estimate than is strictly warranted.
_