검색 전체 메뉴
PDF
맨 위로
OA 학술지
Fisheries resources management of crucian carp based on assessment of fish stock and potential yield in the mid-upper system of Seomjin River
  • 비영리 CC BY-NC
  • 비영리 CC BY-NC
ABSTRACT
Fisheries resources management of crucian carp based on assessment of fish stock and potential yield in the mid-upper system of Seomjin River
KEYWORD
allowable biological catch , Carassius auratus , fisheries resources management , potential yield , stock assessment , yield per recruit
  • INTRODUCTION

    Fish population as a top predator of the aquatic food web plays a major role in aquatic systems and is also of great economic importance (Kalff 2002). Their dynamics have a large influence on the others aquatic organisms; especially, commercial fish species have a direct bearing on the changes of fisheries resources. Therefore, it is essential to manage the fish population systematically and scientifically for efficient utilization of fisheries and eco-system sustainability.

    The fisheries resources management in freshwater of Korea had relied on technical measures. These technical measures included restriction of mesh or setting size, boat licenses, closed fishing areas/seasons, and gear limitations (Kim 2008). Recently, to manage their fisheries, some governments, including Korea, Japan, USA, Australia, and Mexico, adopted a newer fisheries management measure (e.g., output control) based on potential yield and stock assessment in seawater (Lopez-Martinez et al. 2002, Ryu et al. 2005, Lee et al. 2009, Griffiths 2010, Lee et al. 2011). This management measure requires stock assessment methodologies for estimating potential yield (Zhang and Lee 2001). The Potential Yield (PY) will represent an appropriate harvest rate that is applied to a current fish stock size estimated. Ideally, the PY will produce a biologically average yield over the long term and avoid overfishing during periods when the stock is low. Therefore, PY is essential for sustainable utilization and effective management of the stock (Fazli et al. 2009). To properly estimate PY, it is essential to obtain accurate information of target species such as the ecological population characteristics, optimal levels of catch and fishing effort, biomass and age at first capture (Jang et al. 2010, 2011b).

    Crucian carp (Carassius auratus) are distributed throughout North Korea, China, Japan, Africa, and Eurasia (Nakamura 1969, Choi et al. 1990, Nelson 2006). In Korea, crucian carps are widely distributed and caught in freshwater, comprising a large proportion of fisheries resources production in inland water (Seo and Byun 2004). Crucian carps are deemed an iconic species and are greatly targeted by commercial and recreational fishers as they grow to large sizes. Notwithstanding the ecological and commercial significance of crucian carp, scientific information on population ecological characteristics and stock is inadequate (Han et al. 2001, Jang et al. 2011a). Therefore, in order to suggest fisheries resource management implications in freshwater ecosystem and study applicability, crucian carps were selected in this study as one of fish that have commercial value and are distributed widely.

    The specific aims of this study are as follows: (i) to assess annual the stock of crucian carp in the mid-upper system of the Seomjin River, (ii) to estimate the potential yield that is considered the Allowable Biological Catch (ABC), and (iii) to investigate and evaluate the fishing intensity at the current level using yield-per-recruit analysis.

    MATERIALS AND METHODS

      >  Stock Assessment

    Sampling areas were located in the mid-upper system of the Seomjin River (Fig. 1). This area was relatively not disturbed by the artificial environmental factors, and commercial catch activities were strictly prohibited. It has been reported that there are many various species in this area (Jang et al. 2009). Sampling was carried out at 9 sites from August 2008 to April 2009 seasonally, considering fish habitat, water depth, and other factors (Fig. 1). All Sampling type were selected the random block sampling.

    Study on potential yield and biomass of crucian carp could facilitate sustainable utilization of fisheries resources in the way that stock diagnosis enable the quantitative assessment and variation prediction of biomass. The method for estimating stock followed the “swept area method” (Pezzuto et al. 2008, Shin 2009). This restricts the sampling area using seine net (mesh of 10 × 10 mm) and specimens were collected by cast net (mesh of 7 × 7 mm; π × 1.2 m2 = 4.5 m2), fyke net (mesh of 30 × 30 mm; 500 m2), and divers. The divers were sent in to search the ground-fish of deepwater. A total of 402 fish were collected and the catch was 75,233 g.

    In the field, collected specimens were measured up to the nearest 1 mm (total length, TL; body length, BL) and up to the nearest 0.1 g (body weight). Date of collection, length and weight of each specimen were recorded. All specimens were identified according to Choi et al. (1990) and Kim and Park (2002). After identification of species, specimens were randomly selected and transferred to the laboratory for analysis. To analyze age determination, the scales were removed from 93 species at the interval of 10 mm body length. The catch data were applied as input data for estimating stock (Table 1).

    [Table 1.] Catch in number at total length (TL) of crucian carp in the midupper system of the Seomjin River

    label

    Catch in number at total length (TL) of crucian carp in the midupper system of the Seomjin River

    The method for estimating stock followed the “swept area method”, the following equation was used (Pezzuto et al. 2008, Shin 2009).

    image

    where B is the biomass, is an average density (individuals/m2) , is the target species of average weight (g) , A is the area of investigated waters (m2), and q is a catchability. The estimated values of average density and average weight were 0.018 ± 0.001 indviduals/m2 and 187.1 ± 11.2 g, respectively. The catchability was calculated by the assumption that import and export of fish in restricted area by the seine net were equal each to each and all specimens should be caught in sampling area. The area of investigated waters and fishing ground were calculated by Rural Agriculture Water Resource Information System (https://rawris.ekr.or.kr). They were calculated at approximately 22,347 m2 and 2,259,900 m2, respectively.

    In this study, parameters of population ecological characteristics that were estimated from the prior research were applied as input data. In the previous study, age structure of crucian carp has been assessed by otolith. Biological parameters for crucian carp are summarized in Table 2 (Jang et al. 2011a).

    [Table 2.] Parameters estimates used in yield-per recruit and Potential Yield for crucian carp in the mid-upper system of the Seomjin River

    label

    Parameters estimates used in yield-per recruit and Potential Yield for crucian carp in the mid-upper system of the Seomjin River

      >  Age at first capture (tc)

    In general, tc of a stock was estimated directly from fishing experiments either by attaching a small-meshed cover over the cod-end or from the size composition of the catches of nets of smaller meshes caught at the same time and place (Gulland 1983). Pauly method, based on a length-converted catch curve, was used because of the lack of information on fishing experiment. The method involves the following assumption: instantaneous coefficient of total mortality (Z) is same in all age group. When Z was estimated by linearized catch curve, data under age at full recruitment was not used. Virtual catch in number under age at the full recruitment was estimated by extrapolation method. Therefore, it calculates the gap between actual catch in number and virtual catch in number, at all ages. The value of gap means compositive effect between recruitment and setting mesh size (Zhang 2010). In this study, age at first capture was estimated by Pauly method (Pauly 1983, 1984).

      >  Optimal fishing intensity and optimal age at first capture

    To determine the optimal fishing mortality and the optimal age at first capture, the following equation was used (Beverton and Holt 1957).

    image

    where asymptotic weight (W), growth coefficient (K), and age at zero length (t0) are parameters of von Bertalanffy growth function; R is the number of fish alive at time tr, the age at recruitment, which was the youngest age in the catch, the minimum length in the catch was identified as 8.0 cm. converted tr was 0.88 age. Y is the yield and M is the instantaneous coefficient of natural mortality. Un is a summation parameter equal to +1, -3, +3, and -1 for n = 0, 1, 2, and 3, respectively.

      >  Potential Yield

    The Potential Yield (PY) that was defined as the allowable biological catch (ABC) provides an optimal catch level about fisheries resources exploited or unexploited in specified areas. To estimate the PY, this study was undertaken using the improved Fisheries Resources Potential Estimation System (FRPES) (Zhang and Na 2006). It was composed of a five-tier-classification system, which was modified from the six-tier-system used under the US Fisheries Management Plan for the North Pacific groundfish fisheries (NFMS/AFSC 1998). The estimation of FABC (instantaneous rate of fishing mortality that provide for Allowable Biological Catch) of stock in Fisheries Resources Potential Estimation System must take into consideration the quantity and quality of information available and the exploitation history of fisheries resources. For tier 1-3, FABC was determined, and the following equation was used to determine PY.

    image

    where FABC is the instantaneous coefficient of fishing mortality for ABC determined by the available information and the stock status, r is a recruit age, and tL is a maximum fishing age. PYr was computed with the following equation:

    image

    where B is the estimate of biomass. The study was adapted to the 3rd tier classification (Zhang and Lee 2001).

    image

    An estimate of F0.1, which means the fishing mortality corresponding to 10% of the slope at the origin when no fishing occurs, was based on the following equation differentiated from the yield per recruit (YPR) equation of the Beverton and Holt (1957).

    image

    RESULTS AND DISCUSSION

      >  Age at first capture (tc)

    The age at first capture (tc) of crucian carp was estimated to be 1.468 year from the length converted catch curve of the Pauly method (Pauly 1984). Body length that was converted by estimating tc was calcurated as 10.8 cm (Total length 13.1 cm). Crucian carp was known to grow to 14-16 cm in total length during whole 1 year (Son and Song 2006), but, in this study, the results were somewhat differ (Fig. 2).

      >  Optimal f??ishing intensity and optimal age at first capture

    Fishery impact has traditionally been assessed using yield-per-recruit (YPR) analysis (Baranov 1914). YPR is estimated for different fishing mortality rates and different ages at entry to fishery (Allen 2001). Yield isopleth is developed to afford a visual representation of the influence of fishing mortality and age at entry on YPR. This analysis, which represents the tradeoff between growth and natural mortality, can be used to identify whether fish have been caught before their full growth potential (i.e., while their growth rate is still higher than natural mortality) (Wang et al. 2009). The management implications can apply the yield isopleth to help determine the relevant mesh or setting size and fishing effort.

    The estimates of YPR against fishing mortalities (F) and tc for crucian carp are summarized in Fig. 3. Yield isopleths of crucian carp were constructed to examine the response of yield-per-recruit estimates to changes in tc and F (Fig. 4). The current average yield-per-recruit (with F = 0.067 year-1, and tc = 1.468 year) was 15.999 g per recruit, indicating that the fishery was operating below the maximum-yield-per-recruit at 75.617 g, when tc was 2.0 year. Fixing tc at the current level, maximum-yield-per-recruit was 70.903 g when F increased more than 1.890 year-1, but this F was out of range of the maximum-yield-per-recruit line. Fixing F at the current level when tc decreased to 0.359 year, maximum-yield-per-recruit was 17.876 g which resulted in an increase of 11%, i.e., 1.88 in yield-per-recruit.

    Reference points such as Fmax and F0.1 have often applied to evolve fisheries management implications. Table 3 delineated Fmax and F0.1 values from age 1 to 5 of crucian carp, which were calculated from the Beverton and Holt model. The current instantaneous coefficient of fishing mortality (F) did not exceed biological reference points. When tc was 2 year, the yield-per-recruit was the highest at Fmax and F0.1 of 75.067 g and 59.346 g, respectively. Fmax and F0.1 decreased until 5. The yield-per-recruit corresponding to each Fmax value increased to age 2 but the yield-per-recruit gradually declined after the age 2. Several researchers have supported designating Fmax as target reference points and F0.1 as a threshold reference point in order to obtain near optimal yields while guarding against stock collapse (Fazli et al. 2009). The result of the base YPR analysis in the study with tc of 1.468 year suggests that the annual fishing intensity of crucian carp in the mid-upper system of the Seomjin River is lower than the target reference point Fmax as well.

    [Table 3.] Yield per recruit (Y/R) of crucian carp, in the mid-upper system of the Seomjin River under harvest strategies corresponding to Fmax and F0.1

    label

    Yield per recruit (Y/R) of crucian carp, in the mid-upper system of the Seomjin River under harvest strategies corresponding to Fmax and F0.1

      >  Stock assessment and potential yield

    Study on potential yield and biomass of crucian carp could facilitate sustainable utilization of fisheries resources in the way that stock diagnosis enabled the quantitative assessment and variation prediction of biomass. This study calculated biomass using the Swept area method. The annual biomass of cruician carp was estimated as 7,608 kg and biomass per unit area was 3.37 g m-2. The biomass per unit area was not relatively more plentiful than the result of preceded researchthat was performed in lake Namyangand estimated the result as 4.68 g m-2 (Jang et al. 2008). For estimating PY, information on dependable annual biomass, F0.1, and M for crucian carp plus environmental factors affecting recruitment is required. If instantaneous coefficient of fishing mortalities (F) was set at current level, PY was estimated as 343 kg. However, in order to manage fisheries resources of crucian carp population efficiently, PY would be predicted to increase to 2,531 kg, provided that the FABC was set at 0.618 year-1 as derived from the estimated F0.1 for crucian carp (Fig. 5). Also, yield per recruit would increase to 62 g. Considering the current stock status of crucian carp, we suggest the following method based on our result in order to sustain the maximum sustainable yield of crucian carp population in the mid-upper system of the Seomjin River.

    First, make a restriction on minimum fishing body length at 11 cm and over.

    Second, maintain a fishing intensity of 10 times level at the current.

    To establish the ecological sustainability of fish harvest, the PY estimated must be compared with actual catch figure. If total catch exceeds the PY level, fisheries resources obviously are not sustainable. Thus, Estimation of PY value must use to the threshold point, not the target reference point (Caddy and Mahom 1995). If PY values were applied to target reference level, it will occur for overfishing and collapse of fisheries.

    A limited knowledge of the fisheries resources ecology of some freshwater fishes has led to overexploitation and inadequate management in freshwater (Dankel et al. 2008). Varied knowledge of the fisheries resources ecology may help establish appropriate management implication for freshwater fishes, and increasing fishery and biological information may enable a higher level evaluation to be carried out (i.e., quantitative stock assessment) (Zischke 2012). In order to suggest an effective fisheries resources management system, this study was undertaken using stock assessment and potential yield analyses of crucian carp population in the freshwater ecosystem of Korea for the first time.

    CONCLUSIONS

    Study on potential yield and biomass of crucian carp could facilitate sustainable utilization of fisheries resources in the way that stock diagnosis enables the quantitative assessment and variation prediction of biomass. The annual biomass of cruician carp was estimated as 7.6 t (1t = 1,000 kg) in the mid-upper system of the Seomjin River. PY would be predicted to 2.5 t, and this value can be used to establish the allowable biological catch for fisheries resources management. In order to manage fisheries resources of crucian carp population efficiently, if age at first capture (tc) and instantaneous coefficient of fishing mortalities (F) were set at 2 year and 0.618 year-1, the yield per recruit and total allowable catch would be predicted to increase to 62 g and 2,531 kg by about 3.9 times and 7.3 times, respectively.

    Crucian carps (C. auratus) are a significant recreational and commercial species nationwide, and there will be escalating need to evaluate the stock of crucian carp more regularly to guarantee that they can be utilized sustainably. The quality of assessments will considerably depend upon accurate biological and fishery information, as well as an acceptable grasp of stock structure, if stock assessments are required. Future study should base on this information (Griffiths et al. 2010) (Fig. 6).

참고문헌
  • 1. Allen RL 2001 Organization, Functions, and Achievement of the Inter-American Tropical Tuna commission google
  • 2. Baranov FI 1914 The capture of fish by gillnets [Mater Poznaniyu Russ Rybolov] Vol.3 P.56-99 google
  • 3. Beverton RJH, Holt SJ 1957 On the Dynamics of Explotied Fish Populations. Fishery Investigations. Series II google
  • 4. Caddy JF, Mahom R 1995 Reference points for fisheries management google
  • 5. Choi KC, Jeon SR, Kim IS, Son YM 1990 Colored Illustrations of the Freshwater Fishes of Korea google
  • 6. Dankel DJ, Skagen DW, Ulltang O 2008 Fisheries management in practice: review of 13 commercially important sifh stocks [Rev Fish Biol Fish] Vol.8 P.201-233 google
  • 7. Fazli H, Zhang CI, Hay DE, Lee CW 2009 Stock assessment and management implications of anchovy kilka (Clupeonella engrauliformis) in Iranian waters of the Caspian Sea [Fish Res] Vol.100 P.103-108 google cross ref
  • 8. Griffiths SP 2010 Stock assessment and efficacy of size limits on longtail tuna (Thunnus tonggol) caught in Australian waters [Fish Res] Vol.102 P.248-257 google cross ref
  • 9. Griffiths SP, Pepperell JG, Tonks ML, Sawynok W, Olyott L, Tickell S, Zischke MT, Burgess J, Jones E, Joyner D, Lynne J, Makepeace C, Moyle K 2010 Biology, fisheries and status longtail tuna (Thunnus tonggol), with special reference to recreational fisheries in Australian waters google
  • 10. Gulland JA 1983 Fish Stock Assessment: a Manual of Basic Methods google
  • 11. Han KH, Jin DS, Yoo DJ, Baek SR, Hwang DS 2001 The early life history of the crucian carp, Carassius auratus in Korea [Bull Yosu Natl Univ] Vol.16 P.379-387 google
  • 12. Jang SH, Ryu HS, Lee JH 2009 Fish fauna and community structure in the mid-upper region of the Seomjin River [Korean J Limnol] Vol.42 P.394-403 google
  • 13. Jang SH, Ryu HS, Lee JH 2010 Stock assessment and mangement implications of the Korean aucha perch(Coreopera herzi) in freshwater: (1) estimation of population ecological characteristics of Coreperca herzi in the midupper system of the Seomjin River [Korean J Limnol] Vol.43 P.82-90 google
  • 14. Jang SH, Ryu HS, Lee JH 2011 Estimation on population ecological characteristics of crucian carp, Carassius auratus in the mid-upper system of the Seomjin River [Kor J Env Eco] Vol.25 P.318-326 google
  • 15. Jang SH, Ryu HS, Lee JH 2011 Stock assessment and mangement implications of the Korean aucha perch (Coreopera herzi) in freshwater: (2) estimation of potential yield assessment and stock of Coreoperca herzi in the mid-upper system of the Seomjim River [Korean J Limnol] Vol.44 P.172-177 google
  • 16. Jang SH, Zhang CI, Na JH, Kim SW, An KG, Lee JJ, Lee JH 2008 A analysis of trophic structure in lake Namyang using the ecopath modeling [Korean J Limnol] Vol.41 P.144-154 google
  • 17. Kalff J 2002 Limnology google
  • 18. Kim DY 2008 A study on the reinforcement of the competitive power of Korean Inland fisheries [J Fish Bus Adm] Vol.39 P.111-137 google
  • 19. Kim IS, Park JY 2002 Freshwater Fish of Korea google
  • 20. Lee HH, Maunder MN, Piner KR, Methot RD 2011 Estimating natural mortality within a fisheries stock assessment model: an evaluation using simulation analysis based on twelve stock assessments [Fish Res] Vol.109 P.89-94 google cross ref
  • 21. Lee SI, Yang JH, Yoon SC, Chun YY, Kim JB, Cha HK, Choi YM 2009 Biomass estimation of sailfin sandfish, Arctoscopus japonicus in Korean waters [Kor J Fish Aquat Sci] Vol.42 P.487-493 google cross ref
  • 22. Lopez-Martinez J, Arreguın-Sanchez F, Morales-Azpeitia R, Salinas-Zavala C 2002 Stock assessment and potential yield of the rock shrimp, Sicyonia penicillata, fishery of Bahia Kino, Sonora, Mexico [Fish Res] Vol.59 P.71-81 google cross ref
  • 23. Nakamura M 1969 Cyprinid fishes of Japan. Studies on the life history of Cyprinid fishes of Japan [Reg Inst Nat Resour Spec Pub] Vol.4 P.1-455 google
  • 24. Nelson JS 2006 Fishes of the World google
  • 25. 1998 To redefine acceptable biological catch and overfishing. Environmental Assessment/Regulatory Impact Review for Amendment 44 to the Fishery Management Plan google
  • 26. Pauly D 1983 Length-converted catch curves. A powerful tool for fisheries research in the torpics (Part II) [ICLARM Fishbyte] Vol.2 P.17-19 google
  • 27. Pauly D 1984 Fish Population Dynamics in Tropical Waters. A Manual for Use with Programmable Calculators. ICLARM Stud Rev google
  • 28. Pezzuto PR, Alvarez-Perez JA, Wahrlich R 2008 The use of the swept area method for assessing the seabob shrimp Ziphopenaeus kroyeri(Heller, 1862) biomass and removal rates based on artisanal fishery-derived data in southern Brazil: using depletion models to reduce uncertainty [Lat Am J Aquat Res] Vol.36 P.245-257 google cross ref
  • 29. Ryu JG, Kim DY, Lee JS, Kim SJ 2005 The Study on Utilization of Bioeconomic Model for Effectiveness Analysis of Fisheries Management Means google
  • 30. Seo BI, Byun HB 2004 A Philological study on Carassius auratus (L.) [J Appl Orient Med] Vol.4 P.15-28 google
  • 31. Shin YS 2009 A method for estimating potential fishery yield in coastal waters. MS Thesis google
  • 32. Son YM, Song HB 2006 Freshwater fishes of Geum river, Korea google
  • 33. Wang SP, Maunder MN, Aires-da-Silva A, Bayliff WH 2009 Evaulating fishery impacts: application to bigeye tuna (Thunnus obesus) in the eastern Pacific Ocean [Fish Res] Vol.99 P.106-111 google cross ref
  • 34. 2010 Marine Fisheries Resources Ecology google
  • 35. Zhang CI, Lee JB 2001 Stock assessment and management implications of horse mackerel(Trachurus japonicus) in Korean waters, based on the relationships between recruitment and the ocean environment [Prog Oceanogr] Vol.49 P.513-537 google cross ref
  • 36. Zhang CI, Na JH 2006 A studies on estimating potential yield in inland fisheries [Proceedings of 2006 Symposium for Fishery & Resource] P.243-244 google
  • 37. Zischke MT 2012 A review of the biology, stock structure, fisheries and status of wahoo (Acanthocybium solandri), with reference to the Pacific Ocean [Fish Res] Vol.119/120 P.13-22 google cross ref
OAK XML 통계
이미지 / 테이블
  • [ Fig. 1. ]  Map showing the sampling sites (from St. 1 to St. 9, marked as closed circles) in the upper system of the Seomjin River.
    Map showing the sampling sites (from St. 1 to St. 9, marked as closed circles) in the upper system of the Seomjin River.
  • [ Table 1. ]  Catch in number at total length (TL) of crucian carp in the midupper system of the Seomjin River
    Catch in number at total length (TL) of crucian carp in the midupper system of the Seomjin River
  • [ ] 
  • [ Table 2. ]  Parameters estimates used in yield-per recruit and Potential Yield for crucian carp in the mid-upper system of the Seomjin River
    Parameters estimates used in yield-per recruit and Potential Yield for crucian carp in the mid-upper system of the Seomjin River
  • [ ] 
  • [ ] 
  • [ ] 
  • [ ] 
  • [ ] 
  • [ Fig. 2. ]  E Estimation of the selection ogive of crucian carp from length-converted catch curve analysis using the method of Pauly (1984).
    E Estimation of the selection ogive of crucian carp from length-converted catch curve analysis using the method of Pauly (1984).
  • [ Fig. 3. ]  Yield per recruit (Y/R) curve of the crucian carp. (a) Y/R against the age at first capture (tc) for various fishing mortalities (F); (b) Y/R against the fishing mortalities (F) for various the age at first capture (tc).
    Yield per recruit (Y/R) curve of the crucian carp. (a) Y/R against the age at first capture (tc) for various fishing mortalities (F); (b) Y/R against the fishing mortalities (F) for various the age at first capture (tc).
  • [ Fig. 4. ]  Contour plot of yield per recruit of crucian carp in the midupper system of the Seomjin River. P indicates the current state of fishing mortality (F) and age at first capture (tc).
    Contour plot of yield per recruit of crucian carp in the midupper system of the Seomjin River. P indicates the current state of fishing mortality (F) and age at first capture (tc).
  • [ Table 3. ]  Yield per recruit (Y/R) of crucian carp, in the mid-upper system of the Seomjin River under harvest strategies corresponding to Fmax and F0.1
    Yield per recruit (Y/R) of crucian carp, in the mid-upper system of the Seomjin River under harvest strategies corresponding to Fmax and F0.1
  • [ Fig. 5. ]  Variation of potential yield for crucian carp, Carassius auratus, sampled in the mid-upper system of the Seomjin River.
    Variation of potential yield for crucian carp, Carassius auratus, sampled in the mid-upper system of the Seomjin River.
  • [ Fig. 6. ]  A flow diagram showing key information required for stock assessment and management that is known from this research (open box), and still unknown (dashed box) for crucian carp (Griffiths et al., 2010).
    A flow diagram showing key information required for stock assessment and management that is known from this research (open box), and still unknown (dashed box) for crucian carp (Griffiths et al., 2010).
(우)06579 서울시 서초구 반포대로 201(반포동)
Tel. 02-537-6389 | Fax. 02-590-0571 | 문의 : oak2014@korea.kr
Copyright(c) National Library of Korea. All rights reserved.