Nghiên cứu một số đặc điểm sol khí khu vực việt nam qua số liệu aeronet và tư liệu viễn thám

1 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Kết luận văn trung thực chưa cơng bố cơng trình khác Tác giả luận văn Đỗ Ngọc Thúy MỤC LỤC LỜI CAM ĐOAN MỤC LỤC DANH MỤC CÁC KÝ HIỆđU, CÁC CHỮ VIẾT TẮT DANH MỤC CÁC BẢNG DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ SOL KHÍ (AEROSOL) 12 1.1 Định nghĩa số đặc điểm chung sol khí 12 1.2 Một số đặc điểm quang học sol khí 16 1.2.1 Độ dày quang học sol khí 16 1.2.2 Hệ số Angstrom sol khí 17 1.2.3 Phân bố kích thước hạt 18 1.3 Tổng quan tình hình nghiên cứu sol khí giới nước 19 1.3.1 Tổng quan tình hình nghiên cứu giới 19 1.3.2 Tổng quan tình hình nghiên cứu nước 29 CHƯƠNG PHƯƠNG PHÁP NGHIÊN CỨU SOL KHÍ BẰNG CƠNG NGHỆ VIỄN THÁM 30 2.1 Phương pháp nghiên cứu sol khí áp dụng mạng trạm AERONET 30 2.1.1 Mô hình hóa q trình truyền xạ 30 2.1.2 Mơ hình vi vật lý đặc tính quang học sol khí 32 2.1.3 Thuật toán nghịch đảo xạ Mặt trời xạ bầu trời 34 2.1.3.1 Thống mơ hình thuận (forward model) với phép nghịch đảo 34 2.1.3.2 Thực phép nghịch đảo 36 2.2 Phương pháp nghiên cứu sol khí qua vệ tinh MODIS 38 2.2.1 Giới thiệu vệ tinh MODIS/TERRA MODIS/AQUA 38 2.2.2 Thuật tốn tính tốn sản phẩm đặc điểm sol khí MODIS 40 CHƯƠNG MỘT SỐ ĐẶC ĐIỂM CỦA SOL KHÍ TRÊN KHU VỰC VIỆT NAM 56 3.1 PHÂN BỐ THEO KHÔNG GIAN VÀ THỜI GIAN CỦA AOD THEO SỐ LIỆU MODIS 56 3.1.1 Phân bố theo không gian AOD khu vực Việt Nam 56 3.1.2 Phân bố theo thời gian AOD khu vực Việt Nam 60 3.2 ĐẶC ĐIỂM SOL KHÍ KHU VỰC VIỆT NAM TỪ SỐ LIỆU AERONET 63 3.2.1 Đặc điểm sol khí Bắc Giang 63 3.2.2 Đặc điểm sol khí Bạc Liêu 71 KẾT LUẬN 77 TÀI LIỆU THAM KHẢO 79 PHỤ LỤC 83 DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT ACE : Aerosol Characterizatiom Experiment AERONET : AErosol RObotic NETwork AOD : Aerosol Optical Depth BOBMEX-99 : Bay of Bengal Monsoon Experiment in 1999 EOS : the Earth Observing System IGAC : International Global Atmospheric Chemistry INDOEX : INDian Ocean EXperiment LUT : Look Up Table MISR : Multi-angle Imaging Spectroradiometer MODIS : MODerate resolution Imaging Spectroradiometer MSC : Meteorological Service of Canada NCAR : National Center for Atmospheric Research NDVI : Normalized Difference Vegetation Index POLDER :Polarization and Directionality of the Earth’s Reflectances PWC : Precipitable Water Content QC : Quality Control TCO : Total Column Ozone TOA : Top Of Atmospheric DANH MỤC CÁC BẢNG Bảng 2.1a Đặc điểm vệ tinh MODIS 38 Bảng 2.1b Khả ứng dụng vệ tinh MODIS 39 Bảng 2.2 Phân bố kích thước albedo tán xạ đơn sử dụng bảng tra MODIS thuật toán áp dụng khu vực đất liền 46 Bảng 3.1 Độ dày quang học trung bình năm ba khu vực lãnh thổ Việt Nam giai đoạn 2001 – 2010…………………………………………………58 Bảng 3.2 Độ dày quang học trung bình tháng ba khu vực lãnh thổ Việt Nam giai đoạn 2001-2010 61 Bảng 3.3 Phân bố kích thước hạt lớn bán kính hạt tương ứng trạm Bắc Giang giai đoạn 2003-2009 70 Bảng 3.4 Phân bố kích thước hạt lớn bán kính hạt tương ứng trạm Bạc Liêu giai đoạn 2003-2009 76 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1 Sol khí có nguồn gốc từ núi lửa 13 Hình 1.2 Ảnh vệ tinh MODIS vùng nhiệt đới Bắc Đại Tây Dương 3/2004 13 Hình 1.3 Hình dạng sol khí 15 Hình 1.4 Kích thước hạt sol khí 16 Hình 1.5 Hệ thống quan trắc thí nghiệm INDOEX 21 Hình 1.6 Sơ đồ mạng trạm AERONET 23 Hình 1.7 Sơ đồ mạng lưới AEROCAN 24 Hình 1.8 Sơ đồ mạng lưới GAWPFR 25 Hình 1.9 Bản đồ phân bố trạm lidar thuộc mạng lưới EARLINET lãnh thổ Châu ÂU 26 Hình 2.1 Sơ đồ thuật tốn MODIS áp dụng khu vực đất liền 43 Hình 2.2 Bản đồ tỉ lệ độ dày quang học tổng cộng gán “khơng bụi” sol khí dạng mịn khu vực đất liền tính trung bình tháng 45 Hình 2.3 Sơ đồ thuật toán MODIS áp dụng đại dương 51 Hình 3.1 Độ dày quang học trung bình năm khu vực Việt Nam giai đoạn 2001-2010………………………………………………………………… 52 Hình 3.2 Độ dày quang học trung bình tháng khu vực Việt Nam giai đoạn 2001-2010 62 Hình 3.3 Số lần quan trắc độ dày quang học sol khí trạm Bắc Giang giai đoạn 2003 - 2009 64 Hình 3.4 Số ngày quan sát AOD trạm Bắc Giang giai đoạn 20032009 64 Hình 3.5 (a) Biến trình ngày, (b) Biến trình năm độ dày quang học sol khí trạm Bắc Giang giai đoạn 2003-2009 65 Hình 3.6 (a) Biến trình ngày, (b) Biến trình năm thơng số Angstrom trạm Bắc Giang giai đoạn 2003-2009 67 Hình 3.7 Phân bố kích thước hạt sol khí trung bình tháng trạm Bắc Giang giai đoạn 2003-2009 69 Hình 3.8 Số lần quan trắc độ dày quang học sol khí trạm Bạc Liêu giai đoạn 2003-2009 71 Hình 3.9 Số ngày quan sát AOD bước sóng 500nm trạm Bạc Liêu giai đoạn 2003-2009 72 Hình 3.10 (a) Biến trình ngày, (b) Biến trình năm độ dày quang học sol khí trạm Bạc Liêu giai đoạn 2003 – 2009 73 Hình 3.11 (a) Biến trình ngày, (b) Biến trình năm thơng số Angstrom trạm Bạc Liêu giai đoạn 2003- 2009 73 Hình 3.12 Phân bố kích thước trung bình tháng trạm Bạc Liêu giai đoạn 2003-2009 76 MỞ ĐẦU Tính cấp thiết đề tài Những thập kỷ gần đây, mức độ phát thải sol khí (các hạt thể rắn lỏng tồn lơ lửng không khí) vào khí ngày tăng liên quan đến q trình phát triển cơng nghiệp quốc gia Nồng độ sol khí khơng khí tăng lên tác động trực tiếp tới sức khoẻ đời sống người giảm chất lượng khơng khí, ngồi cịn tác động gián tiếp thông qua ảnh hưởng tới thời tiết, khí hậu [6] Các hạt sol khí tán xạ hấp thụ xạ làm cho lớp (bề mặt) khí ấm lên làm cho bề mặt Trái Đất lạnh Các hạt sol khí làm tăng số hạt nhân ngưng kết hình thành hạt nước nhỏ hơn, dẫn đến tăng tán xạ phản xạ mây Các hạt nước nhỏ làm hạn chế va chạm liên kết, kéo dài thời gian tồn mây ngăn cản lớn lên hạt nước mây tạo hạt mưa Ảnh hưởng trực tiếp sol khí làm giảm lượng xạ trung bình tồn cầu giới hạn khí 0.35W/m2, tăng lớp khí 3.0W/m2, giảm bề mặt Trái đất 3.4W/m2 [9] Trong vùng châu Á gió mùa, trung bình năm, lượng xạ khí (mặt đất) tăng (giảm) 10-20W/m2 Màn mây nâu ABCs (Atmospheric Brown Clouds) cấu thành từ chất ô nhiễm cácbon đen, cácbon hữu cơ, tro, bụi chất hấp thụ sunfat, ngăn cản xạ mặt trời tới mặt đất làm giảm 50% nóng lên tồn cầu tăng khí nhà kính [17] Việt Nam nằm khu vực quốc gia phát triển, trình phát thải chất vào khí đa dạng thành phần, phong phú số lượng Thêm vào chế độ hoàn lưu khu vực Việt Nam phức tạp, nên vận chuyển sol khí từ vùng khác Thế giới đến Việt Nam khuếch tán sol khí từ Việt Nam vào khí đa dạng Nhận biết vấn đề trên, học viên chọn hướng nghiên cứu số đặc điểm sol khí khu vực Việt Nam làm chủ đề cho luận văn tốt nghiệp Một số đặc điểm sol khí sở cho nghiên cứu sức khỏe môi trường, biến đổi khí hậu, đánh giá chất lượng khơng khí, xác định mơ hình truyền xạ nghiên cứu nguồn xạ từ Trái Đất Đề tài đưa phù hợp cần thiết cho trình nghiên cứu vấn đề liên quan tới khí Nghiên cứu tiến hành dựa chuỗi tư liệu viễn thám từ vệ tinh MODIS (MODerate Resolution Imaging Spectroradiometer – quang phổ kế tạo ảnh có độ phân giải trung bình) số liệu từ mạng trạm quan trắc sol khí tồn cầu AERONET (AErosol RObotic NETwork) NASA Bắc Giang Bạc Liêu thuộc quyền quản lý Viện Vật lý địa cầu Mục tiêu đề tài Đưa quy luật biến đổi độ dày quang học sol khí không gian thời gian khu vực Việt Nam qua phân tích số liệu MODIS AERONET Nghiên cứu làm rõ số đặc điểm quang học sol khí qua số liệu mạng trạm AERONET đặt Bắc Giang Bạc Liêu Đối tượng phạm vi nghiên cứu Đề tài nghiên cứu số đặc điểm sol khí độ dày quang học sol khí, hệ số Angstrom, phân bố kích thước hạt sol khí khu vực Việt Nam Nội dung nghiên cứu Nghiên cứu sở lý thuyết phương pháp viễn thám sử dụng để quan trắc sol khí 10 Nghiên cứu quy luật phân bố theo thời gian không gian độ dày quang học sol khí khu vực Việt Nam dựa số liệu MODIS Nghiên cứu phân tích đặc điểm quang học sol khí từ số liệu ARONET Bắc Giang Bạc Liêu bao gồm độ dày quang học dải phổ, hệ số Angstrom, phân bố kích thước hạt sol khí Phương pháp nghiên cứu Phương pháp nghiên cứu sol khí qua thiết bị quang phổ kế đa dải phổ từ mặt đất Trong đề tài sử dụng thiết bị đo hai trạm Bắc Giang Bạc Liêu để xác định thông số quang học sol khí Lợi phương pháp cho phép xác định chi tiết điểm Phương pháp viễn thám sử dụng số liệu vệ tinh MODIS/TERRA để nghiên cứu sol khí Phương pháp cho phép nghiên cứu khu vực rộng, bổ trợ cho phương pháp quang phổ kế đặt mặt đất Trong luận văn, thực việc nghiên cứu thuật toán nghịch đảo áp dụng mạng trạm AERONET thuật tốn áp dụng xử lí số liệu MODIS Các phương pháp thống kê, phân tích khơng gian, thời gian sử dụng để xác định quy luật số đặc điểm sol khí Ý nghĩa khoa học thực tiễn đề tài nghiên cứu Ý nghĩa khoa học: đề tài khai thác sử dụng hiệu số liệu vệ tinh số liệu quan sát mặt đất để phục vụ mục đích nghiên cứu sol khí cho khu vực Việt Nam Các kết nghiên cứu quy luật biến đổi sol khí có ý nghĩa quan trọng làm tiền đề cho nhiều nghiên cứu lĩnh vực khí tượng, khí hậu mơi trường Ý nghĩa thực tiễn: kết nghiên cứu đề tài sử dụng để hiệu chỉnh ảnh hưởng khí tới đặc điểm bề mặt xác định 84 clc; clear all; %%%!!!!!! %%% SO LIEU MODIS KHU VUC VIET NAM !!!!!!!!! %%%%%%%% %% Doc so lieu ngay, tinh trung binh thang, trung binh nam, ve thi clc; clear all; addpath E:\fonction aodmy=ones(2160,10); for my=1:10; yy=2000+my; aody=ones(180,10); for mm=1:12; if yy==2000 | yy==2004 | yy==2008; n1=[31 29 31 30 31 30 31 31 30 31 30 31 else n1=[31 28 31 30 31 30 31 31 30 31 30 31]; end n=n1(mm); aodm=ones(180,n); for dd=1:n dm1=[0 31 59 90 120 151 181 212 243 273 304 334]; dm2=[0 31 60 91 121 152 182 213 244 274 305 335]; if yy==2004 | yy==2008; day=dm2(mm)+dd; else day=dm1(mm)+dd; end yyst=int2str(yy); dayst=int2str(day); if day=10&day0 85 aod = fread(f1,'float32'); aod=reshape(aod,360,180); aod(find(aod','LineWidth',0.8,'MarkerSize',3.5); xlim([0 13]); set(gca,'XTick',1:1:13) set(gca,'XTickLabel',{'1','2','3','4','5','6','7','8','9','10','11','12','13'},'FontSize',10) ylim([0 2.2]); set(gca,'YTick',0:0.2:2.2) set(gca,'YTickLabel',{'0','0.2','0.4','0.6','0.8','1.0','1.2','1.4','1.6','1.8','2.0','2.2'},'FontSize',10) legend('1020nm','870nm','675nm','500nm','440nm','380nm') %%%% BIEN TRINH NGAY %%%% for i=1:24 landoNgay1020(i)=0; tongNgay1020(i)=0; TBNgay1020(i)=0; landoNgay870(i)=0; tongNgay870(i)=0; TBNgay870(i)=0; landoNgay675(i)=0; tongNgay675(i)=0; TBNgay675(i)=0; landoNgay500(i)=0; tongNgay500(i)=0; TBNgay500(i)=0; landoNgay440(i)=0; tongNgay440(i)=0; TBNgay440(i)=0; landoNgay380(i)=0; tongNgay380(i)=0; TBNgay380(i)=0; end 92 for j=1:r if AOD1020(j)~=999; landoNgay1020(hh1(j)+1)=landoNgay1020(hh1(j)+1)+1; tongNgay1020(hh1(j)+1)=tongNgay1020(hh1(j)+1)+AOD1020(j); end if AOD870(j)~=999; landoNgay870(hh1(j)+1)=landoNgay870(hh1(j)+1)+1; tongNgay870(hh1(j)+1)=tongNgay870(hh1(j)+1)+AOD870(j); end if AOD675(j)~=999; landoNgay675(hh1(j)+1)=landoNgay675(hh1(j)+1)+1; tongNgay675(hh1(j)+1)=tongNgay675(hh1(j)+1)+AOD675(j); end if AOD500(j)~=999; landoNgay500(hh1(j)+1)=landoNgay500(hh1(j)+1)+1; tongNgay500(hh1(j)+1)=tongNgay500(hh1(j)+1)+AOD500(j); end if AOD440(j)~=999; landoNgay440(hh1(j)+1)=landoNgay440(hh1(j)+1)+1; tongNgay440(hh1(j)+1)=tongNgay440(hh1(j)+1)+AOD440(j); end if AOD380(j)~=999; landoNgay380(hh1(j)+1)=landoNgay380(hh1(j)+1)+1; tongNgay380(hh1(j)+1)=tongNgay380(hh1(j)+1)+AOD380(j); end end for i=1:24 if landoNgay1020(i)~=0; TBNgay1020(i)=tongNgay1020(i)/landoNgay1020(i); end if landoNgay870(i)~=0; TBNgay870(i)=tongNgay870(i)/landoNgay870(i); end 93 if landoNgay675(i)~=0; TBNgay675(i)=tongNgay675(i)/landoNgay675(i); end if landoNgay500(i)~=0; TBNgay500(i)=tongNgay500(i)/landoNgay500(i); end if landoNgay440(i)~=0; TBNgay440(i)=tongNgay440(i)/landoNgay440(i); end if landoNgay380(i)~=0; TBNgay380(i)=tongNgay380(i)/landoNgay380(i); end end for i=1:24 if (TBNgay1020(i)==0); TBNgay1020(i)=NaN; end if (TBNgay870(i)==0); TBNgay870(i)=NaN; end if (TBNgay675(i)==0); TBNgay675(i)=NaN; end if (TBNgay500(i)==0); TBNgay500(i)=NaN; end if (TBNgay440(i)==0); TBNgay440(i)=NaN; end if (TBNgay380(i)==0); TBNgay380(i)=NaN; end end figure(2); x=0:1:23 plot(x,TBNgay1020,'-bp','LineWidth',0.8,'MarkerSize',3.5);hold on; grid on plot(x,TBNgay870,'-ro','LineWidth',0.8,'MarkerSize',3.5); plot(x,TBNgay675,'-gs','LineWidth',0.8,'MarkerSize',3.5); plot(x,TBNgay500,'-yd','LineWidth',0.8,'MarkerSize',3.5); plot(x,TBNgay440,'-c^','LineWidth',0.8,'MarkerSize',3.5); plot(x,TBNgay380,'-m>','LineWidth',0.8,'MarkerSize',3.5); xlim([5 18]) set(gca,'XTick',6:1:18) set(gca,'XTickLabel',{'6','7','8','9','10','11','12','13','14','15','16','17','18'},'FontSize',10) ylim([0 2]) set(gca,'YTick',0:0.2:2.0) set(gca,'YTickLabel',{'0','0.2','0.4','0.6','0.8','1.0','1.2','1.4','1.6','1.8','2.0'},'FontSize',10) legend('1020nm','870nm','675nm','500nm','440nm','380nm') %% SO LAN QUAN TRAC AOD for i=1:r if (DATA(i,9)==999); DATA(i,9)=NaN; end if (DATA(i,10)==999); DATA(i,10)=NaN; end 94 if (DATA(i,11)==999); DATA(i,11)=NaN; end if (DATA(i,17)==999); DATA(i,17)=NaN; end if (DATA(i,20)==999); DATA(i,20)=NaN; end if (DATA(i,22)==999); DATA(i,22)=NaN; end end aod1020=DATA(1:r,9); aod870=DATA(1:r,10); aod675=DATA(1:r,11); aod500=DATA(1:r,17); aod440=DATA(1:r,20); aod380=DATA(1:r,22); plot(aod870,'-r'); plot(aod675,'-g'); plot(aod500,'-y'); plot(aod440,'-c'); plot(aod380,'-m'); figure(333);clf; plot(aod1020,'-b'); hold on; xlabel('Ngay Julian','FontSize',10) ylabel('Do day quang hoc cua sol khi','FontSize',10) title('DO DAY QUANG HOC CUA SOL KHI TRAM BAC GIANG') legend('1020nm','870nm','675nm','500nm','440nm','380nm') %%% ANGSTROM TRONG DAI PHO 440-870 load E:\LuanVanThSy\SoLieu\BacGiangAeronet\AeronetBG.mat angstrom2003=a2003(:,42); angstrom2004=a2004(:,42); angstrom2005=a2005(:,42); angstrom2006=a2006(:,42); angstrom2007=a2007(:,42); angstrom2009=a2009(:,42); angstrom=DATA(:,42); [r,c]=size(angstrom) %!! BIEN TRINH NAM CUA ANGSTROM (2003-2009) for i=1:12 %! tao mang so lieu lando(i)=0; tong(i)=0; TB(i)=0; end for j=1:r if angstrom(j)~=999; lando(mm1(j))=lando(mm1(j))+1; tong(mm1(j))=tong(mm1(j))+angstrom(j);end end 95 for i=1:12 if lando(i)~=0; TB(i)=tong(i)/lando(i); end end figure(1); x=1:1:12 plot(x,TB,'-ok','LineWidth',1.0,'MarkerSize',3.5);hold on; grid on; xlim([0 13]) set(gca,'XTick',1:1:13) set(gca,'XTickLabel',{'1','2','3','4','5','6','7','8','9','10','11','12','13'},'FontSize',10) ylim([0.6 2.2]) set(gca,'YTick',0.6:0.2:2.2) set(gca,'YTickLabel',{'0.6','0.8','1.0','1.2','1.4','1.6','1.8','2.0','2.2'},'FontSize',10) %!BIEN TRINH NGAY for i=1:24 lando(i)=0; tong(i)=0; TB(i)=0; end for j=1:r if angstrom(j)~=999; lando(hh1(j)+1)=lando(hh1(j)+1)+1; tong(hh1(j)+1)=tong(hh1(j)+1)+angstrom(j); end end for i=1:24 if lando(i)~=0; TB(i)=tong(i)/lando(i); end end for i=1:24 if (TB(i)==0); end TB(i)=NaN; end 96 figure(2); x=0:1:23 plot(x,TB,'-ko','LineWidth',1,'MarkerSize',3.5);hold on; grid on xlim([5 18]) set(gca,'XTick',6:1:18) set(gca,'XTickLabel',{'6','7','8','9','10','11','12','13','14','15','16','17','18'},'FontSize',10) ylim([0.4 2.4]) set(gca,'YTick',0.4:0.2:2.4) set(gca,'YTickLabel',{'0.4','0.6','0.8','1.0','1.2','1.4','1.6','1.8','2.0','2.2','2.4'},'FontSize',10) %%%% PHAN BO KICH THUOC %%%% clc; clear all; chdir E:\LuanVanThSy\SoLieu\BacGiangAeronet\PhanboKT a=load('BGPhanBo.txt'); b=load('BanKinh.txt'); [r,c]=size(a) mm(:,:)=a(:,2); yy(:,:)=a(:,1); data(:,:)=a(:,3:c); [r1,c1]=size(data) %!! TRUNG BINH NAM for k=1:12 for j=1:22 lando(k,j)=0; tong(k,j)=0; tb(k,j)=0; end end for i=1:r1 for j=1:22 if data(i,j)~=99999; lando(mm(i),j)=lando(mm(i),j)+1; tong(mm(i),j)=tong(mm(i),j)+data(i,j);end end end 97 for k=1:12 for j=1:22 if lando(k,j)~=0; tb(k,j)=tong(k,j)/lando(k,j); end end end figure(100); for k=1:12 if k= =1; semilogx1=semilogx(b,tb(k,:),'-bo','LineWidth',0.5,'MarkerSize',3.5); grid on; elseif k= =2; semilogx1=semilogx(b,tb(k,:),'-rs','LineWidth',0.5,'MarkerSize',3.5); elseif k= =3; semilogx1=semilogx(b,tb(k,:),'-go','LineWidth',0.5,'MarkerSize',3.5); elseif k= =4; semilogx1=semilogx(b,tb(k,:),'-ms','LineWidth',0.5,'MarkerSize',3.5); elseif k= =5; semilogx1=semilogx(b,tb(k,:),'-co','LineWidth',0.5,'MarkerSize',3.5); elseif k= =6; semilogx1=semilogx(b,tb(k,:),':bs','LineWidth',1,'MarkerSize',3.5); elseif k= =7; semilogx1=semilogx(b,tb(k,:),':ro','LineWidth',1,'MarkerSize',3.5); elseif k= =8; semilogx1=semilogx(b,tb(k,:),':gs','LineWidth',1,'MarkerSize',3.5); elseif k= =9; semilogx1=semilogx(b,tb(k,:),':mo','LineWidth',1,'MarkerSize',3.5); elseif k= =10; semilogx1=semilogx(b,tb(k,:),':cs','LineWidth',1,'MarkerSize',3.5); elseif k= =11; semilogx1=semilogx(b,tb(k,:),'-bs','LineWidth',0.5,'MarkerSize',3.5); else semilogx1=semilogx(b,tb(k,:),'-ro','LineWidth',0.5,'MarkerSize',3.5); end end hold on; 98 xlim([0.01 100]); xlabel(' r ylim([0 0.18]) (micromet)','FontSize',10) ylabel(' dV(r)/dln r (micromet^3/micromet^2)','FontSize',10) title('PHAN BO KICH THUOC PHAN TU THEO THE TICH TRAM BAC GIANG') legend('thang 1','thang 2','thang ','thang 4','thang 5','thang 6','thang 7','thang 8','thang 9','thang 10','thang 11','thang 12') %!!TRUNG BINH THANG for k=1:12 figure(k) semilogx1=semilogx(b,tb(k,:),'-ks','LineWidth',1,'MarkerSize',3.5);grid on; xlim([0.01 100]) ylim([0 0.18]) xlabel(' r (micromet)','FontSize',12) ylabel(' dV(r)/dln r (micromet^3/micromet^2)','FontSize',12) k1= num2str(k) tit = strcat(' THANG ',k1) title(tit) end ... theo thời gian AOD khu vực Việt Nam 60 3.2 ĐẶC ĐIỂM SOL KHÍ KHU VỰC VIỆT NAM TỪ SỐ LIỆU AERONET 63 3.2.1 Đặc điểm sol khí Bắc Giang 63 3.2.2 Đặc điểm sol khí Bạc Liêu ... học sol khí khơng gian thời gian khu vực Việt Nam qua phân tích số liệu MODIS AERONET Nghiên cứu làm rõ số đặc điểm quang học sol khí qua số liệu mạng trạm AERONET đặt Bắc Giang Bạc Liêu Đối tư? ??ng... tư? ??ng phạm vi nghiên cứu Đề tài nghiên cứu số đặc điểm sol khí độ dày quang học sol khí, hệ số Angstrom, phân bố kích thước hạt sol khí khu vực Việt Nam Nội dung nghiên cứu Nghiên cứu sở lý thuyết