bluuebunny/crossref_metadata_2025_split huggingface dataset:
Land use/land cover change prediction in Dak Nong Province based on remote sensing and Markov Chain Model and Cellular Automata
Authors: Nguyen Thi Thanh Huong , Ngo Thi Thuy Phuong 
Published: September 2020

Land use and land cover changes (LULCC) including deforestation for agricultural land and others are elements that contribute on global environmental change. Therefore, understanding a trend of these changes in the past, current, and future is important for making proper decisions to develop in a sustainable way. This study analyzed land use and land cover (LULC) changes over time for Tuy Duc district belonging to Dak Nong province based on LULC maps classified from a set of multi-date satellite images captured in year 2003, 2006, 2009, and 2013 (SPOT 5 satellite images). The LULC spatio-temporal changes in the area were classified as perennial agriculture, cropland, residential area, grassland, natural forest, plantation and water surface. Based on these changes over time, potential LULC in 2023 was predicted using Cellular Automata (CA)–Markov model. The predicted results of the change in LULC in 2023 reveal that the total area of forest will lose 9,031ha accounting of 50% in total area of the changes. This may be mainly caused by converting forest cover to agriculture (account for 28%), grassland (12%) and residential area (9%). The findings suggest that the forest conversion needs to be controlled and well managed, and a reasonable land use plan should be developed in a harmonization way with forest resources conservation.
Thay đổi sử dụng đất và thảm phủ (LULCC) bao gồm cả việc phá rừng để phát triển nông nghiệp và vì các mục đích khác là tác nhân đóng góp vào biến đổi môi trường toàn cầu. Vì vậy hiểu biết về khuynh hướng của sự thay đổi này trong quá khứ, hiện tại và tương lai là quan trọng để đưa ra những quyết định dúng đắn để phát triển bền vững. Nghiên cứu đã phân tích LULCC trong thời gian qua dựa vào các bản đồ sử dụng đất và thảm phủ (LULC) đã được phân loại từ một loạt ảnh vệ tinh đa phổ được thu chụp vào năm 2003, 2006, 2009 (ảnh SPOT 5). Những thay đổi LULC theo thời gian và không gian trong khu vực được phân loại thành đất nông nghiệp với cây dài ngày, cây ngắn ngày, thổ cư, trảng cỏ cây bụi, rừng tự nhiên, rừng trồng và mặt nước. Dựa trên sự thay đổi này theo thời gian, LULC tiềm năng cho năm 2023 đã được dự báo bằng cách sử dụng mô hình CA-Markov. Kết quả dự báo LULCC năm 2023 đã cho thấy tổng diện tích rừng bị mất khoảng 9,031 ha chiếm 50% trong tổng số diện tích thay đổi. Điều này chủ yếu là do chuyển đổi từ rừng tự nhiên sang canh tác nông nghiệp (chiếm 28%), trảng cỏ cây bụi (12%) và khu dân cư (9%). Kết quả cho thấy việc chuyển đổi rừng cần phải được kiểm soát và quản lý tốt và một kế hoạch sử dụng đất hợp lý cần được xây dựng trong sự hài hòa với bảo tồn tài nguyên rừng.
Comparison of sensitivity of three legume species exposed to crude extracts of toxic and non-toxic cyanobacteria
Authors: Pham Thanh Luu 
Published: September 2020

We evaluated the effect of cyanobacterial crude extracts containing microcystin (CCEMC+) from a natural bloom on seed germination and initial development of three economically important legume species: green mung bean Vigna radiata, cowpea Vigna cylindrical and red mung bean Vigna angularis and compared it to crude extracts of cyanobacteria without the toxin (CCEMC–). Results showed that CCEMC+ and CCEMC– caused different effects on seed germination and initial development of the three species. There was a clear inhibition on germination and root growth of the green mung bean exposed to the CCEMC+ (20, 200 and 500 μg/L), indicating that the green mung bean being more sensitive to CCEMC+ when compared to the cowpea and red mung bean. CCEMC+ induced a greater occurrence of abnormal seedlings in the green mung bean, duce to inhibition the germination as well as reduction of fresh weight and root length. The CCEMC– extract caused no harmful effects to germination and seedlings growth of the green mung bean and red mung bean. However, it reduced shoot and root length in cowpea bean, suggesting that the cowpea being more sensitive to both extracts. Our results indicated that the sensitivity in germination and root growth of the green mung bean V. radiata could be used as an indicator to evaluate the toxic effect and monitor the toxin concentration of water contaminated with microcystins.
Nghiên cứu này khảo sát và so sánh các tác động bất lợi của ly trích vi khuẩn lam có chứa và không chứa độc tố lên sự nẩy mầm và sự phát triển ở giai đoạn đầu của ba loại cây họ đậu gồm đậu xanh Vigna radiata, đậu đỏ Vigna angularis và đậu đen Vigna cylindrical. Kết quả cho thấy hai loại ly trích gây ra các tác động khác nhau lên ba loại cây đậu thí nghiệm. Ly trích vi khuẩn lam có chứa độc tố ở nồng độ 20, 200 và 500 μg/L ngăn chặn đáng kể sự nẩy mầm và sự phát triển rể ở đậu xanh. Cả hai loại ly trích có chứa và không chứa độc tố đều ngăn chặn sự nẩy mầm và sự phát triển rể ở đậu đen. Ngược lại ly trích vi khuẩn lam có chứa độc tố ở nồng độ 500 μg/L lại kích thích chiều dài rể, thân mầm và trọng lượng tươi ở đậu đỏ. Kết quả cho thấy đậu đen khá nhạy cảm với cả hai loại ly trích có chứa và không chứa độc tố, trong khi đó đậu xanh nhạy cảm hơn với ly trích có chứa độc tố. Tính nhạy cảm của các loại cây họ đậu khi phơi nhiễm với ly trích vi khuẩn lam có thể được sử dụng để chỉ thị cho sự ô nhiễm và quan trắc độc tố vi khuẩn lam trong môi trường.
Numerical model for estimating greenhouse gas emissions from pulp and paper industrial wastewater treatment systems in Vietnam
Authors: Dang Xuan Hien , Nguyen Thi Van Anh , Nguyen Duc Toan , Dang Thanh Son 
Published: September 2020

At present, it is difficult and costly to measure directly greenhouse gas (GHG) emissions from the wastewater treatment system. Application of model will reduce measurement cost and quickly obtain the forecast data set of GHG emissions. This study developed a mathematical model for both steady and dynamic states to calculate GHG (CO2, CH4, and N2O) emissions from wastewater treatment systems for industrial paper processing. These models are constructed based on mass balance equations of species, including substrate balance equations, biomass balance equations for reactors of treatment systems, stoichiometric coefficiences of species in biochemical reactions and biological processes. The obtained equations were solved based on algorithm of Runge-Kutta and the model was programmed by MATLAB. Results of applying the model to calculate GHG emissions from the paper industrial wastewater treatment system at Bai Bang and Tan Mai plants are as follows: total GHG emissions and emission factor are 3,070.3 kgCO2-eq/day, 0.38 kgCO2-eq/m3, respectively for Bai Bang plant (8,000 m3/day) and 7,413.6 kgCO2-eq/day, 0.74 kgCO2-eq/m3, respectively for Tan Mai plant (10,000 m3/day). The research evaluated a number of influencing factors, such as temperature, flow rate of influent, and substrate concentrations, to GHG emissions at the Tan Mai paper plant.
Hiện nay, việc đo đạc trực tiếp phát thải khí nhà kính (KNK) từ hệ thống xử lý nước thải còn khó khăn và tốn kém. Việc áp dụng mô hình sẽ giảm được chi phí đo đạc và nhanh chóng có được bộ số liệu dự báo một cách tương đối về phát thải KNK. Nghiên cứu đã thiết lập được mô hình toán ở trạng thái ổn định và trạng thái không ổn định để tính toán phát thải khí nhà kính (CO2, CH4, N2O) từ hệ thống xử lý nước thải sản xuất giấy. Các mô hình này dựa trên các phương trình cân bằng chất của các cấu tử bao gồm các phương trình cân bằng cơ chất, các phương trình cân bằng sinh khối trong các bể phản ứng và các hệ số tỷ lượng của các chất tham gia các phản ứng sinh hóa. Các phương trình được giải bằng thuật toán Runge-Kutta và mô hình được lập trình trên ngôn ngữ MATLAB. Mô hình được áp dụng tính toán phát thải khí nhà kính từ hệ thống xử lý nước thải tại nhà máy giấy Bãi Bằng và nhà máy giấy Tân Mai, được kết quả như sau: tổng phát thải khí nhà kính (KNK) và hệ số phát thải là 3.070,3 kg CO2-tđ/ngày, 0,38 kg CO2-tđ/m3 tại Nhà máy giấy Bãi Bằng (8.000 m3/ngày) và 7.413,6 kg CO2-tđ/ngày, 0,74 kg CO2-tđ/m3 nhà máy giấy Tân Mai (10.000 m3/ngày). Nghiên cứu đã đánh giá được một số các yếu tố ảnh hưởng như nhiệt độ, lưu lượng nước thải và nồng độ cơ chất dòng vào đến sự phát thải KNK tại nhà máy giấy Tân Mai.
Study on synthesis of MoS2 modified g-C3N4 materials for treatment of Direct black 38 dye
Authors: Phung Thi Lan , Nguyen Thi Kim Giang 
Published: September 2020

Pure g-C3N4 and MoS2 modified g-C3N4 materials were synthesized using a facile heating method and a low-temperature hydrothermal method, respectively. The obtained samples were characterized by XRD pattern and N2 adsorption-desorption technique at 77K. The adsorption and photocatalytic performance of all obtained samples were investigated by discoloration of direct black 38 dye in the dark and under visible light irradiation. The results showed that all obtained samples exhibited good discoloration efficiency of direct black 38 dye. The two factors including pH values and Mo loading effected mainly on elimination efficiency of direct black 38 dye. MoS2 modified g-C3N4 materials possessed the more enhanced adsorption and photocatalytic performance in comparison to pure g-C3N4 at pH value of 3.5, with adsorbent dosage of 0.1 g/L. Furthermore, it was found that the adsorption process and photo-catalysis simultaneously occurred under visible light irradiation and followed up a pseudo-second-order kinetic reaction of Langmuir - Hinshelwood model.
g-C3N4 và g-C3N4 biến tính bởi MoS2 đã được tổng hợp theo phương pháp nung đơn giản và phương pháp thủy nhiệt ở nhiệt độ thấp tương ứng. Các mẫu tổng hợp đã được đánh giá đặc trưng bởi các phương pháp hiện đại như giản đồ nhiễu xạ tia X, phương pháp hấp phụ-khử hấp phụ N2 ở 77K. Khả năng hấp phụ và quang hóa xúc tác của các vật liệu tổng hợp đã được nghiên cứu bởi quá trình phân hủy màu thuốc nhuộm direct black 38 trong điều kiện bóng tối và chiếu sáng bởi ảnh sáng nhìn thấy của đèn chiếu sáng sợi đốt wolfram (220V-100W). Các kết quả nghiên cứu chỉ ra rằng các mẫu tổng hợp đều có hiệu suất xử lý màu cao đối với thuốc nhuộm direct black 38. Hai yếu tố gồm pH dung dịch và hàm lượng MoS2 ảnh hưởng chính đến hiệu suất xử lý màu direct black 38. g-C3N4 biến tính bởi MoS2 luôn thể hiện hiệu suất hấp phụ và quang hóa cao hơn so với g-C3N4 tinh khiết. Hơn nữa, khi được chiếu sáng bởi ánh sáng nhìn thấy thì quá trình hấp phụ và quá trình quang hóa thuốc nhuộm direct black 38 trên các vật liệu tổng hợp đã xảy ra đồng thời và mô hình Langmuir - Hinshelwood động học bậc 2 đã được đề xuất cho quá trình này.
Removal of arsenic from contaminated groundwater using laterite, sand and ash:  a case study in Son Dong commune, Hoai Duc district, Ha Noi
Authors: Nguyen Quoc Bien , Tran Hai Nam , Nguyen Thi Hoang Ha 
Published: September 2020

Arsenic (As) contaminated groundwater has been a major concern due to the negative impacts to exposed people. This research was conducted to assess and compare the removal efficiency of As from groundwater by laterite, sand, and ash. The experiment was carried out in 14 days in a household scale in Son Dong commune, Hoai Duc district, Ha Noi. Groundwater was pumped directly from a well and flowed through 20 cm (diameter) x 80 cm (length) columns. The initial As concentration in groundwater was 526 µg/L decreasing to an average of 189, 192 and 154 µg/L after being filtrated using sand, ash, and laterite, respectively. Average removal efficiency of sand, ash and laterite during the experiment was 63.3, 63.9, and 70.5%, respectively. Laterite had higher As removal efficiency may be due to higher content of goethite and kaolinite in this sorbent which resulted in better adsorption of As. The concentrations of As in the outflow water were higher than the allowable limit set by the national technical regulation on drinking water quality (QCVN 01:2009/BYT). Therefore, it requires higher volume of sorbents or additional treatment technologies for removal of As from groundwater.
Ô nhiễm asen (As) trong nước ngầm hiện nay là một vấn đề cấp bách đối với xã hội bởi những rủi ro tiềm ẩn với sức khoẻ con người. Nghiên cứu được thực hiện nhằm đánh giá và so sánh khả năng hấp phụ As trong nước ngầm của đá ong, cát và tro. Thí nghiệm được tiến hành trong vòng 14 ngày đặt tại một hộ gia đình có nguồn nước ngầm bị ô nhiễm As thuộc xã Sơn Đồng, huyện Hoài Đức, Hà Nội. Nước ngầm được bơm từ giếng và chảy qua các cột đựng vật liệu có đường kính và chiều dài lần lượt là 20 cm và 80 cm. Nồng độ As ban đầu là 526 µg/L đã giảm xuống còn 189, 192 và 154 µg/L sau khi lọc bằng cát, tro và đá ong. Hiệu suất xử lý As trung bình của cát, tro và đá ong lần lượt là 63,3, 63,9 và 70,5 µg/L. Đá ong xử lý As tốt hơn có thể do hàm lượng goethit và kaolinit cao hơn trong vật liệu hấp phụ này dẫn đến khả năng hấp phụ As tốt hơn. Tuy nhiên hàm lượng As trong nước đầu ra vẫn chưa đạt quy chuẩn cho nước uống (QCVN 01:2009/BYT). Do đó, cần tăng thêm lượng vật liệu hoặc kết hợp với các phương pháp khác để xử lý As hiệu quả hơn.
Allometric relations between biomass and diameter at breast height and height of tree in natural forests at Me Linh Station for Biodiversity, Vinh Phuc Province, Vietnam
Authors: Dang Thi Thu Huong , Do Huu Thu , Trinh Minh Quang , Nguyen Hung Manh , Bui Thi Tuyet Xuan , Nguyen Tien Dung 
Published: September 2020

Stem diameter at breast height (D1.3m) and tree height (H) are commonly used measures of tree growth. Based on correlation analysis between biomass of stem, branches and leaves and stem diameter and height of tree we can identify allometric equation for predicting biomass and carbon sequestration of the vegetation. This study was carried out in the natural forests of Me Linh Station for biodiversity to develop allometric equation between biomass and diameter at breast height and height of tree. The study results indicated that twenty tree species dominate in natural forests in Me Linh Station for Biodiversity and they were selected for sampling. Through the 80 established linear equation models for above and below –ground biomass (AGB and BGB), we found that the biomass of tree species in Me Linh Station for Biodiversity were closely correlated with the diameter factor (R>0.902) and not clearly correlated with the height (correlation coefficient = 0.5498, R2< 0.549). Four regression equations were established, including: Pstem = 25.3051*(D1.3m)0.4627 (R2 : 9.661); Pbranch = 12.1043*(D1.3m)0.5416 (R2 : 9.8); Pleaves = 9.446*(D1.3m)0.5976 (R2 : 0.9363); P total biomass of forest = 25.882*D1.725 with R2: 0.8561) for estimating biomass and carbon sequestration of natural forest at the research site.
Đường kính ngang ngực (D1.3m) và chiều cao (H) cây là hai nhân tố thường được dùng để đánh giá sự phát triển của cây gỗ. Việc xây dựng các phương trình tương quan giữa sinh khối (SK) thân, cành, lá, sinh khối tầng cây gỗ, sinh khối của quần xã thực vật với đường kính và chiều cao cây góp phần rất lớn trong dự báo sinh khối và khả năng hấp thụ khí carbon của thảm thực vật. Kết quả nghiên cứu cho thấy 20 loài cây gỗ chiếm ưu thế trong rừng tự nhiên và chúng được chọn để thu mẫu. Mối tương quan giữa sinh khối với 2 nhân tố điều tra rừng là đường kính ngang ngực và chiều cao cây đã đươc kiểm tra thông qua 80 phương trình tương quan. Nhìn chung, sinh khối có tương quan chặt chẽ với nhân tố đường kính (hệ số tương quan R > 0,902), và không tương quan rõ với nhân tố chiều cao (R < 0,5498). Bốn phương trình tính sinh khối cho thảm rừng tại khu vực nghiên cứu đã được thiết lập: SKthân = 25,3051*(D1,3m)0,4627 (R2: 9,661); SKcành: 12,1043*(D1,3m)0,5416 (R2: 9,8); SKlá: 9,446*(D1,3m)0,5976 (R2: 0,9363) và SKtổng = 25,882*D1,725 with R2: 0,8561).
Die Wahrnehmung der Figur „Adolf Hitler“ in Daniel Levys „Mein Führer- die wirklich wahrste Wahrheit über Adolf Hitler“ durch internationales Filmpublikum. Eine empirische Untersuchung im landeskundlichen DaF-Unterricht.
Authors: Herzner Dominik 
Published: September 2020

Filme beeinflussen die Wahrnehmungen von Figuren. Daniel Levys Film "Mein Führer" (die erste deutschsprachige Komödie über Adolf Hitler) führte zu einer positiven Wahrnehmung der Hauptfigur. In einer empirischen Untersuchung wurde gezeigt, dass sich dieser Effekt bei internationalem Filmpublikum, die einen anderen Zugang zur Thematik Nationalsozialismus haben, verstärkt.
OPTIMALISASI PORTOFOLIO SAHAM MENGGUNAKAN METODE MEAN ABSOLUTE DEVIATION DAN SINGLE INDEX MODEL PADA SAHAM INDEKS LQ-45
Authors: Wulandari Diah , Ispriyanti Dwi , Hoyyi Abdul 
Published: September 2020

Stock investment is the planting of money in a securities that indicates the ownership of a company in order to provide benefits in the future. In obtaining optimal results from stock investments, investors are expected to create a series of portfolios. The portfolio will help investors in allocating some funds in different types of investments in order to achieve optimal profitability. For selection of optimal stocks representing LQ-45 Index, used 2 methods of Mean Absolute Deviation (MAD) method and Single Index Model (SIM) method. In MAD method, 5 best stocks are BBCA with weight 23%, INDF 8%, KLBF 23%, TLKM 23%, and UNVR 23%. While the SIM method of candidate portfolio obtained is AKRA with weight 15,459%, BBCA 48,193%, BBNI 5,028%,KLBF 0,258% and TLKM 31,062%. Portfolio performance meter is used by sharpe ratio. The value of sharpe ratio is 0,36754 for optimal portfolio using MAD method and 0,40782 for optimal portfolio using SIM method, this means that optimal portfolio using SIM method has better performance than MAD. Keywords: Investment, Portfolio, Index LQ-45, Mean Absolute Deviation, Single Index Model, Sharpe Ratio
ANALISIS ANTREAN BUS NONPATAS AKAP DAN AKDP JALUR TIMUR TERMINAL TIRTONADI KOTA SURAKARTA
Authors: Sitomurang Rosalina Aprilda , Sugito Sugito , Mukid Moch. Abdul 
Published: September 2020

The queuing system is a set of customers, services and a set of rules governing the arrival of its customers and services. Queue is a waiting phenomenon that is part of everyday human life. The queue is formed if the number of subscribers to be served exceeds the available service capacity. Queue phenomenon one of them seen in the queue nonpatas buses at Terminal Tirtonadi Surakarta. Nonpatas bus lanes studied include non-purpose buses Surabaya, Karanganyar, Wonogiri, Purwodadi and Pedesaan. The queue displant used is FIFO (First In First Out). For the five nonpatas bus lanes it meets steady state conditions because it has utility value less than 1. The selected model is a model that has the following 4 types of distributions: Erlang, Weibull, Gamma and Lognormal. The queue model generated for the five tracks (ERLA/ERLA/1):(GD/∞/∞) for Surabaya nonpatas buses, (ERLA/WEIB/1):(GD/∞/∞) for Karanganyar nonpatas buses, (GAMM/WEIB/1):(GD/∞/∞) for Wonogiri nonpatas buses, (ERLA/WEIB/1):(GD/∞/∞) for Purwodadi nonpatas buses, (WEIB/LOGN/1):(GD/∞/∞) for Pedesaan nonpatas buses. Based on the value of the system performance measure indicated that the five lines are queue system is good. Keywords: Beta, Erlang, FIFO, Gamma, Steady State Conditions, Lognormal, Queue Model, Queuing Systems, System Performance Measure, Weibull
PEMODELAN DEFORESTASI HUTAN LINDUNG DI INDONESIA MENGGUNAKAN MODEL GEOGRAPHICALLY AND TEMPORALLY WEIGHTED REGRESSION (GTWR)
Authors: Adiningrumh Thea Zulfa , Prahutama Alan , Santoso Rukun 
Published: September 2020

Regression analysis is a statistical analysis method that is used to modeling the relationship between dependent variables and independent variables. In the linear regression model only produced parameter estimators are globally, so it’s often called global regression. While to analyze spatial data can be used Geographically Weighted Regression (GWR) method. Geographically and Temporally Weighted Regression (GTWR) is the development of GWR model to handle the instability of a data both from the spatial and temporal sides simultaneously. In this GWR modeling the weight function used is a Gaussian Kernel, which requires the bandwidth value as a distance parameter. Optimum bandwidth can be obtained by minimizing the CV (cross validation) coefficient value. By comparing the R-square, Mean Square Error (MSE) and Akaike Information Criterion (AIC) values in both methods, it is known that modeling the level of deforestation in protected forest areas in Indonesia in 2013 through 2016 uses the GTWR method better than global regression. With the R-square value the GTWR model is 25.1%, the MSE value is 0.7833 and AIC value is 349,6917. While the global regression model has R-square value of 15.8%, MSE value of 0.861 and AIC value of 361,3328. Keywords : GWR, GTWR, Bandwidth, Kernel Gaussian
KAPABILITAS PROSES DENGAN ESTIMASI FUNGSI DENSITAS KERNEL PADA PRODUKSI DENIM DI PT APAC INTI CORPORA
Authors: Ramadhani Puput , Ispriyanti Dwi , Safitri Diah 
Published: September 2020

The quality of production becomes one of the basic factors of consumer decisions in choosing a product. Quality control is needed to control the production process. Control chart is a tool used in performing statistical quality control. One of the alternatives used when the data obtained is not known distribution is analyzed by nonparametric approach based on estimation of kernel density function. The most important thing in estimating kernel density function is optimal bandwidth selection (h) which minimizes Cross Validation (CV) value. Some of the kernel functions used in this research are Rectangular, Epanechnikov, Triangular, Biweight, and Gaussian. If the process control chart is statistically controlled, a process capability analysis can be calculated using the process conformity index to determine the nature of the process capability. In this research, the kernel control chart and process conformity index were used to analyze the slope shift of Akira-F style fabric and Corvus-SI style on the production of denim fabric at PT Apac Inti Corpora. The results of the analysis show that the production process for Akira-F style is statistically controlled, but Ypk > Yp is 0.889823 > 0,508059 indicating that the process is still not in accordance with the specified limits set by the company, while for Corvus- SI is statistically controlled and Ypk < Yp is 0.637742 < 0.638776 which indicates that the process is in accordance with the specification limits specified by the company. Keywords: kernel density function estimation, Cross Validation, kernel control chart, denim fabric, process capability
ESTIMASI VALUE AT RISK PORTOFOLIO SAHAM MENGGUNAKAN METODE GARCH-COPULA (Studi Kasus : Harga Penutupan Saham Harian Unilever Indonesia dan Kimia Farma Periode 1 Januari 2013- 31 Desember 2016)
Authors: Prasetya Lingga Bayu , Ispriyanti Dwi , Prahutama Alan 
Published: September 2020

Any investment in the stock market will earn returns accompanied by risks. Return and risk has a mutual correlation that equilibrium. The formation of a portfolio is intended to provide a lower risk or with the same risk but provide a higher return. Value at Risk (VaR) is a instrument to analyze risk management. Time series model used in stock return data that it has not normal distribution and heteroscedastisicity is Generalized Autoregressive Conditional Heteroscedasticity (GARCH). GARCH-Copula is a combined method of GARCH and Copula. The Copula method is used in joint distribution modeling because it does not require the assumption of normality of the data and can capture tail dependence between each variable. This research uses return data from stock closing prices of Unilever Indonesia and Kimia Farma period January 1, 2013 until December 31, 2016. Copula model is selected based on the highest likelihood log value is Copula Clayton. Value at Risk estimates of Unilever Indonesia and Kimia Farma's stock portfolio on the same weight were performed using Monte Carlo simulation with backtesting of 30 days period data at 95% confidence level. Keywords : Stock, Risk, Generalized Autoregressive Conditional Heteroscedasticity (GARCH), Copula, Value at Risk