The blackfin flounder
기름가자미(
가자미류는 동해안에서 자망, 연승, 저인망 및 트롤 등 여러 어업에서 어획되고 있고, 가자미류의 어획량은 최근 증가하는 경향을 보이고 있다. 국립수산과학원 동해수산연구소에서 수행한 어업자원조사 결과에서도 동해안에 기름가자미의 분포밀도가 높은 것으로 나타났고, 동해안에서 어획되는 가자미류 중 기름가자미의 어획비율이 54.9%에 이른다고 보고한 바 있다 (NFRDI, 2010).
기름가자미는 주로 동해안 동해구외끌이중형저인망과 자망어업에 의해서 주로 어획되고 있으며, 우리나라 동해안의 연도별 어획량은 1990년대 초 1,600 mt의 수준에서 점차 증가하기 시작하여 1999년 4,050 mt으로 가장 높은 어획량을 보인 후 2009년까지 2,800 mt의 수준에서 변동을 보이다가 2010-2012년에는 3,500-3,700 mt의 어획량을 보여 최근 점차 증가하는 경향을 나타내었다. 동해구외끌이중형저인망 어업의 CPUE는 2007년 이후 점차 증가하여 최근에는 척당 80-90 mt의 수준을 나타내고 있다(Fig. 1).
가자미류의 어획량 증가는 기름가자미의 어획량증가를 의미하며 기름가자미 주 어업인 동해구외끌이중형저인망과 자망 어업의 특성상 항상 남획의 가능성을 가지고 있으며, 최근 자원량이 서서히 감소하는 경향을 보이고 있다. 따라서 기름가자미 자원보호와 관리를 위해 수산자원회복사업 대상종으로 선정하였고, 현재 본 종의 자원회복을 위한 생태학적 연구 등 많은 연구가 수행되고 있다.
기름가자미에 관한 연구는 국외에서 연령과 성장(Hashimoto, 1953; Ishida and Kotakata, 1953), 자치어의 변태(Okiyama, 1963), 산란생태(Ivankova, 1974), 섭식특성(Hayase and Hamai, 1974; Pushchina, 2000), 분포특성(Shvydkii and Vdovin, 2001; Tokranov, 2008) 등의 보고가 있고, 국내에서는 성숙과 산란(Cha et al., 2008), 연령과 성장(Lee et al., 2008; Yang et al., 2012), 자망과 삼중자망의 선택성 연구(Park et al., 2011) 등 연구가 최근 활발히 이루어지고 있다. 국내에서 수행된 가자미목에 속하는 종들에 관한 연구는 참가자미의 연령과 성장 (Choi et al., 1986; Lee et al., 2006), 물가자미의 생물학적 연구(Hwnag et al., 1979), 범가자미의 연령과 성장(Jeon et al., 1996), 문치가자미의 연령과 성장(Park, 1997), 용가자미의 성장과 성숙(Choi et al., 1999), 돌가자미의 연령과 성장(Jun and Im, 2004) 등이 있으며, 대부분의 연구가 상업적 가치가 높은 종들이거나 종묘생산을 위한 선행연구로 수행된 바 있다. 따라서 동해안의 주요 우점종인 기름가자미 자원의 지속적인 생산을 위한 기름가자미의 생태학적 특성치 및 자원량 추정 등 자원평가 및 관리방안 연구가 필요한 실정이다.
본 연구는 가장 최근에 추정된 성장계수(Yang et al., 2012)를 초기자료로 사용하여 추정된 생잔율, 사망계수, 어획개시연령 등의 자원생태학적 특성치와 어획물의 연령조성 및 어획량 자료를 이용하여 자원평가를 실시하고 적정어획량을 추정하여, 우리나라 동해안에 서식하는 주요 상업어종인 기름가자미의 효율적 관리방안 수립을 위한 과학적 근거를 제공하는데 목적이 있다.
기름가자미의 연령사정과 성장매개변수 추정은 2007년 강원도 해역에서 동해구외끌이중형저인망에 의해 채집된 표본을 대상으로 하였으며, 성장매개변수인 이론적 최대체장(
The von Bertalanffy growth parameters of the blackfin flounder Glyptocephalus stelleri in the East Sea, Korea (from Yang et al., 2012)
순간전사망계수는 최근 3년(2010-2012년) 동안 동해구외끌이중형저인망에서 어획된 기름가자미의 체장조성(TL, 전장)자료와 Yang et al. (2012)에 의한 von Bertalanffy 성장매개변수(Table 1)를 적용하여 어획물곡선법(Pauly, 1984)으로 추정하였다. 즉, 최근 3년간의 체장자료를 성장매개변수에 의해 연령조성자료로 변환하고, 완전가입연령 이상에 해당하는 연령조성자료를 다음 식에 대입하여 직선회귀시켜 순간전사망계수를 추정하였다.(C / Δt) = c-Z(t + Δt / 2)
여기서,
일반적인 어류의 사망은 자연사망과 어획사망으로 구분되는데, 여기서 자연사망은 어획을 제외한 다른 요인에 의한 사망을 의미하며, 이입과 이출은 고려하지 않았다. 기름가자미의 순간 자연사망계수를 추정하기 위하여 Zhang and Megrey (2006) 방법을 적용하였으며, 그 모델식은 아래와 같다.
여기서
일반적으로 어떤 자원의 어획개시연령은 망목시험을 통하여 추정하고 있으나 이에 관련된 자료가 없으므로 어획물의 체장 조성과 von Bertalanffy의 성장식을 이용하여 추정하는 Pauly(1984) 방법을 사용하였다. 최근의 어업에 대한 어획개시연령을 추정하기 위하여 동해구외끌이중형저인망 어업의 2010-2012년간 체장조성자료를 이용하였다.
기름가자미의 자원량을 추정하는데는 성장개념을 포함시킨 생체량을 기초로 한 코호트분석법(Zhang, 1987; Zhang and Sullivan, 1988)이 기존의 다른 모델에 비하여 생물학적 특성을 더 잘 반영한다는 장점을 가지고 있으므로 이 모델을 적용하여 자원량을 추정하였다. 먼저, 마지막 연도의 연령별 자원량은 다음의 모델로 추정하였다,
연도별 최고령의 자원량은
마지막 연도와 연도별 최고령를 제외한 자원량은,
으로 각각 나타낼 수 있다. 여기서,
이 방법에서는 대상연도의 연도별 연령별 어획중량과 연령별 성장률, 순간자연사망계수, 마지막 어획사망계수(terminal
연령별 성장률(
순간자연사망계수는 추정된 성장매개변수를 입력자료로 사용하였으며, 마지막 어획사망계수는 추정된
마지막 어획사망계수의 경우, 동일한 값을 사용하였으므로 1차로 추정된 연도별 연령별 어획사망계수를 사용하여 실제의 자원구조를 반영하는 자원량과 어획사망계수를 구하기 위하여 다음과 같이 반복계산에 의해 보정하였다. 첫째, 어획사망계수는 연령별로 어구의 선택비가 다를 것[
기름가자미의 최대지속적생산량(MSY)과 MSY시의 어획노력량(
Fox 모델은
위 식을 f에 대해 미분하면
가 되며,일 때,
Schaefer 모델은
또한, 적정가입체장(
위의 수식을 사용하여
기름가자미의 생물학적허용어획량은 Zhang and Lee (2001)이 제시한 이용가능 정보의 질적수준에 따른 단계별 ABC 추정시스템을 이용하여 추정하였으며, 본 연구에서는 이용가능한 기름가자미의 정보수준에 따라 3-5단계의 관리수준을 이용하였다.
3단계 추정방법은 최근년도 자원량(
4단계 추정방법은 연도별 어획량과 노력량 자료를 이용하여 최대지속적생산량을 추정하여 적정어획량의 기준으로 사용하였다.
5단계 방법은 가장 간단한 방법으로 단지 대상어종의 연도별 어획량 정보만을 필요로 한다. 적정어획량을 산출함에 있어서 자원량의 변동을 예측할만한 충분한 자료가 없을 경우에 예방적 접근방법이 최근 많이 적용되고 있다. 특히 어느 정도의 자료가 있다 하더라도 자료의 불확실성이나 자원붕괴의 위험성이 존재할 경우 이러한 예방적 접근이 자원붕괴를 방지하는데 효율적일 수 있다.
기름가자미의 순간전사망계수는 2010년부터 2012년까지 3개년에 대한 동해구외끌이중형저인망 어획물의 체장조성 평균치를 사용하여 어획물곡선법(Pauly, 1984)에 의해 1.0542/year로 추정되었고, 이 값을 생잔율(
기름가자미의 순간자연사망계수는 von Bertalanffy 성장매개변수
Estimation of instantaneous coefficient of natural mortality (M) of Glyptocephalus stelleri in the East Sea, Korea
동해구외끌이중형저인망 어업의 2010-2012년(3년간)의 체장조성자료와 von Bertalanffy 성장매개변수(Table 1)를 적용하여 Pauly의 어획물곡선법으로 추정된 동해안 기름가자미의 50% 어구가입연령 혹은 어획개시연령은 1.304 year이었고(Fig. 2), 이 때의 연령을 체장으로 환산한 결과 11.5 cm이었다.
한국 동해산 기름가자미의 자원량을 추정하기 위해 먼저 연령별 순간성장율을 구하였고, Yang et al. (2012)에 의한 연령별 체장조성자료와 체장-체중 관계식을 이용하여 연령별 체중으로 전환한 후 순간성장율을 구하였다(Table 3).
Estimated instantaneous growth rate (Gi) of Glyptocephalus stelleri in the East Sea, Korea
1991-2012년간 연도별 연령별 순간어획사망계수를 추정하였는데, 연령별 성장률, 1차 추정된 어획사망계수로부터 마지막 해의 전(全)연령을 제외한 연령별 평균
Tuned instantaneous rate of fishing mortality by age of Glyptocephalus stelleri in the East Sea, Korea, 1991-2012
마지막 연령의 전(全)연도에 해당되는
Tuned instantaneous rate of fishing mortality by year of Glyptocephalus stelleri in the East Sea, Korea, 1991-2012
생체량을 기초로 한 코호트분석법에 의해 추정된 1991-2012년간 한국 동해안 기름가자미의 연도별 연령별 자원량 변동을 나타내면(Fig. 3, Fig. 4), 1991년에 6,570 mt에서 점차 증가하기 시작하여 1997년 12,800 mt로 최고 수준을 보인 후, 약간 감소하기 시작하여 2004년 10,500 mt까지 감소하였고 최근 2011-2012년 평균 9,800 mt의 수준을 나타내고 있다(Fig. 3). 단위노력당어획량(CPUE) 변동은 2005년 45.7 mt/boat에서 점차 증가하여 2011년 87.6 mt/boat로 가장 높았고 2012년에는 78.4 mt/boat로 약간 감소한 수준을 나타내고 있다(Fig. 1).
연령별 자원량은 3세의 자원량이 전체 자원량의 25%로 가장 높았고 4세가 22%, 2세가 18%, 5세가 13%의 순으로 나타나 2-5세의 자원량이 전체의 78%를 차지하고 있었다. 6세 이후의 자원량은 급격히 감소하여 8-9세의 자원량이 전체의 1.3%로 고령어의 생존율이 매우 낮은 것으로 나타났고(Fig. 4), 연령별 자원량 추정치에서 뚜렷한 탁월연급군은 출현하지 않았다.
2005-2012년의 연도별 어획량 자료와 척당 어획량(CPUE)자료를 Scheafer 모델과 Fox 모델에 적용하여 추정된 동해 안 기름가자미의 최대지속적생산량(MSY)은 각각 3,686 mt (R2=0.54,
Deriso (1987)가 제안한 어획개시체장과
기름가자미의 생물학적허용어획량은 Zhang and Lee (2001)이 제시한 단계별 ABC 추정시스템을 이용하여 추정하였으며, 본 연구에서는 3-5단계의 관리수준을 이용하였다.
3단계 추정방법은 최근년도 자원량(B), 어획사망계수(
4단계 추정방법에서 기름가자미의 생물학적허용어획량은 결정계수(R2) 값이 높은 Fox 모델에 의해 추정된 최대지속적생산량(MSY)의 95% 신뢰구간의 하한값인 3,397 mt를 ABC로 추정하였다.
5단계 방법은 가장 간단한 방법으로 기름가자미의 생물학적 허용어획량을 추정하기 위하여 2010-2012년간의 3개년 어획량자료의 산술평균치(
[Table 6.] ABC estimates of Glyptocephalus stelleri by tiers in the East Sea, Korea
ABC estimates of Glyptocephalus stelleri by tiers in the East Sea, Korea
본 연구는 한국 동해안산 기름가자미의 자원생태학적 특성치, 자원량 및 적정어획량(ABC)를 추정하여 현재 과도어획이 우려되고 있는 기름가자미의 효율적 관리방안 수립에 기초자료로 활용하기 위해 수행되었다.
생잔율을 추정하기 위하여 사용되는 연령조성자료는 원칙적으로 한 코호트(동일발생연급군)에 대해 일생을 통하여 추적한 자료를 사용하여야 하지만, 대부분은 한 코호트에 대한 자료를 얻기가 힘들기 때문에 자원의 평형상태를 가정하여 동해구외끌이중형저인망 어업에 의한 최근 3개년(2010-2012년)의 체장조성자료를 평균하여 Pauly (1984)의 어획물곡선법으로 순간 전사망계수(
순간자연사망계수(
동해구외끌이중형저인망 어업에 의한 기름가자미의 50% 어구가입연령은 대상어구에 대한 체급별 망목선택성 자료가 없는 경우에 체장조성자료를 사용하여 추정가능한 Pauly (1984) 방법을 사용하여 1.304세로 추정되었고 이때의 체장은 11.5 cm로 계산되었다.
기름가자미 어획량은 동해구외끌이중형저인망 47%, 연안자망 36%로 전체 어획량의 83%를 점유하고 있는데(NFRDI, 2012), 동해구외끌이중형저인망 어업이 연안자망 어업에 비해 미성숙개체의 어획비율이 높아 어획강도가 높은 것으로 생각되기 때문에 망목선택성 등 각 어구의 특성이 기름가자미 자원에 미치는 영향은 추후 검토되어야 할 것이다.
본 연구에서 1991년 이전의 기름가자미 어획량(가자미류 어획량 중 기름가자미의 어획비율) 자료가 없어 기름가자미의 자원변동 경향은 파악할 수 없었지만, 기름가자미의 연도별 자원량은 1991년 6,570 mt를 나타낸 이후 점차 증가하기 시작하여 1997년 12,800 mt로 최고값을 보인 후 증감을 반복하면서 점차 감소하여 최근 2011-2012년 평균 10,300 mt의 수준을 보이고 있다. 또한, 기름가자미의 어획량은 1990년대 초반(1991-1993 평균) 약 1,640 mt 수준에서 2012년에 3,560 mt으로 약 2.2배 정도 증가한 것으로 나타났는데, 이는 어획노력량의 증가에 의한 것으로 판단되며, 기름가자미의 자원량에 대한 어획량 수준을 고려해 볼 때 자원이 감소하고 있다는 뚜렷한 징후는 아직까지 보이지 않고 있다.
동해구외끌이중형저인망 어업의 노력량인 척당 어획량과 총어획량 자료를 이용하여 기름가자미의 최대지속적생산량(MSY)을 구한 결과, Schaefer 모델에서 MSY는 3,686 mt,
Deriso (1987)의 적정어획사망계수(
기름가자미의 생물학적허용어획량(ABC)은 최근년도 자원량(
일반적으로, MSY 수준의 자원량이 현재의 자원량과 비교하여 그 차이가 미미하거나 더 클 경우에는 과도개발의 우려가 있으므로 자원보존을 위한 조치가 필요하다(Gulland, 1983). 그러므로 실제 생물학적허용어획량은 Fox 모델에 의한 MSY 값의 95% 신뢰구간의 하한값인 3,397 mt를 ABC로 추정하였다. 따라서, ABC 추정시스템 3-5단계에 의해 ABC는 2,622-3,571 mt로 최종 결정되었다.
이와 같이 본 연구를 통하여 밝혀진 한국 동해안 기름가자미의 자원생태학적 특성치와 어획량 및 자원량의 변동양상, 생물학적허용어획량 등의 자료를 이용하여 현재의 자원상태를 확인해 볼 수 있는 자원평가를 수행하고, 자원평가 결과를 바탕으로 기름가자미 자원보존을 위한 합리적이고 효율적인 자원관리 및 어업관리 정책이 수립되고 시행되어야 할 것이다.