1
ĐẶT VẤN ĐỀ
Rừng ngập mặn là một hệ sinh thái rừng đặc biệt mà Việt Nam là một
trong những Quốc gia được thiên nhiên ban tặng. Rừng ngập mặn chỉ có ở các
nước nhiệt đới, á nhiệt đới và có vai trò bảo vệ môi trường, con người, đặc
biệt trong bảo vệ bờ biển các vùng duyên hải.
Việt Nam, với bờ biển dài hơn 3620 km, phù sa bồi đắp nhiều tạo điều
kiện thuận lợi cho sự phát triển rừng ngập mặn kinh tế biển. Rừng ngập mặn
có vai trò quan trọng trong việc bảo vệ môi trường sống cho các loài động
thực vật và cung cấp nguồn thức ăn cho con người, bên cạnh đó rừng ngập
mặn có vai trò chắn sóng, tăng lượng bồi đắp phù sa, điều hòa không khí – là
nhân tố góp phần chống lại biến đổi khí hậu.
Hiện nay, quá trình phát triển kinh tế xã hội với tốc độ đô thị hóa diễn
ra ngày càng mạnh, cùng với tốc độ gia tăng dân số nên con người đã khai
thác và sử dụng rừng ngập mặn vào nhiều mục đích khác nhau làm cho diện
tích rừng ngập mặn ngày càng bị thu hẹp, thành phần loài thực vật và chất
lượng rừng ngày càng bị suy giảm.
Gần đây, vấn đề nuôi trồng thủy sản của nhân dân địa phương đang gặp
nhiều khó khăn do tình trạng ô nhiễm môi trường sống tại các khu đầm và
vùng nước ven biển. Tình trạng bệnh dịch của các loài thủy sản đang ngày
càng xuất hiện với tần suất lớn hơn, các bãi bồi ngày càng bị thu hẹp, nguyên
nhân chủ yếu là do rừng ngập mặn đã bị khai thác quá mức làm giảm chức
năng bảo vệ môi trường sống, hơn nữa vấn đề nuôi trồng thủy sản đang có sự
mâu thuẫn rất lớn giữa các cộng đồng nên hiệu quả đạt thấp.
2
Với vị trí địa lý là một bán đảo, thị xã Quảng Yên là một trong những
khu vực có diện tích rừng ngập mặn lớn của tỉnh Quảng Ninh. Rừng ngập
nghiệp, các cá nhân... đánh giá được hiện trạng của các quá trình, các thực thể
tự nhiên, kinh tế - xã hội thông qua các chức năng thu thập, quản lý, truy vấn,
phân tích và tích hợp các thông tin được gắn với một nền hình học (bản đồ)
nhất quán trên cơ sở toạ độ của các dữ liệu đầu vào.
Có nhiều cách tiếp cận khác nhau khi định nghĩa GIS. Nếu xét dưới góc
độ hệ thống, thì GIS có thể được hiểu như một hệ thống gồm các thành phần:
con người, phần cứng, phần mềm, cơ sở dữ liệu và quy trình-kiến thức
chuyên gia nơi tập hợp các quy định, quy phạm, tiêu chuẩn, định hướng, chủ
trương ứng dụng của nhà quản lý, các kiến thức chuyên ngành và các kiến
thức về công nghệ thông tin.
Khi xây dựng một hệ thống GIS ta phải quyết định xem GIS sẽ được
xây dựng theo mô hình ứng dụng nào, lộ trình và phương thức tổ chức thực
hiện nào. Chỉ trên cơ sở đó người ta mới quyết định xem GIS định xây dựng
sẽ phải đảm đương các chức năng trợ giúp quyết định gì và cũng mới có thể
có các quyết định về nội dung, cấu trúc các hợp phần còn lại của hệ thống
4
cũng như cơ cấu tài chính cần đầu tư cho việc hình thành và phát triển hệ
thống GIS. Với một xã hội có sự tham gia của người dân và quá trình quản lý
thì sự đóng góp tri thức từ phía cộng đồng đang ngày càng trở nên quan trọng
và càng ngày càng có vai trò không thể thiếu.
1.1.1. Khái niệm về GPS:
Hệ thống Định vị Toàn cầu (Global Positioning System - GPS) là hệ
thống xác định vị trí dựa trên vị trí của các vệ tinh nhân tạo, do Bộ Quốc
phòng Hoa Kỳ thiết kế, xây dựng, vận hành và quản lý. Trong cùng một thời
điểm, tọa độ của một điểm trên mặt đất sẽ được xác định nếu xác định được
khoảng cách từ điểm đó đến ít nhất ba vệ tinh.
1.1.2. Khái niệm về viễn thám
nhau tạo thành Hệ thống thông tin địa lý. Nếu thiếu một trong hợp phần trên
thì GIS sẽ ngừng hoạt động hoặc chỉ là một hệ thống chết. Theo đà phát triển
của khoa học kỹ thuật bốn hợp phần trên cũng được được phát triển mạnh mẽ,
tạo lên một Hệ thông tin địa lý cũng được phát triển hơn, thực hiện được các
chức năng ưu việt hơn, nhanh hơn, mạnh hơn, tiện lợi và dễ dàng hơn trong
vận hành và sử dụng.
Như vậy ta có thể hiểu một cách khái quát về GIS như sau: Hệ thống
thông tin địa lý là tập hợp có tổ chức của phần cứng, phần mềm máy tính, dữ
liệu địa lý và các thủ tục của người sử dụng nhằm trợ giúp việc thu nhận, lưu
trữ, quản lý, xử lý, phân tích, hiện thị các thông tin không gian từ thế giới
6
thực để giải quyết các vấn đề tổng hợp thông tin cho các mục đích của con
người đặt ra.
Ngoài ra còn một số định nghĩa về GIS của một số tác giả như sau:
- Theo Burrough (1986) thì GIS là “Tập hợp các công cụ để thu nhập,
lưu trữ, tra cứu, chuyển đổi và biểu thị các dữ liệu không gian từ thế giới
thực”.
- Parker (1988) định nghĩa GIS như một “kĩ nghệ thông tin nhằm lưu
trữ, phân tích và biểu thị dữ liệu không gian và phi không gian”.
- Aronoff (1989) quan niệm GIS là “bất kỳ một phương thức trên sách
tra khảo hoặc máy tính dùng để lưu trữ, thao tác các dữ liệu tham chiếu địa
lý”
1.2.1. Ứng dụng viễn thám và GIS trên thế giới
Các nước tiên tiến ứng dụng viễn thấm đầu tiên là: Mỹ, Nga, Ấn Độ,
Canada, Nhât, và mới đây có thêm Trung Quốc … việc ứng dụng kết hợp
công nghệ viễn thám và GIS đã trở thành công nghệ hoàn chỉnh, được sử
dụng rộng rãi không chỉ để theo dõi, kiểm kê, dự báo, quản lý các tài
Định hay nghiên cứu sự sạt lở bờ sông Tiền, sông Hậu thuộc tỉnh Vĩnh Long.
Các công trình này đã thể hiện được tính ưu việt của tư liệu ảnh viễn thám kết
hợp với GIS trong khả năng ứng dụng của nó. Các công trình nghiên cứu biến
động tập trung vào việc sử dụng các tư liệu ảnh đa thời gian cũng như ảnh của
nhiều vệ tinh khác nhau để nghiên cứu biến động lớp phủ thực vật như biến
động tài nguyên rừng, biến động sử dụng đất, biến động đường bờ biển, biến
động rừng ngập mặn, nghiên cứu biến động về hiện trạng sử dụng đất....
8
Công cụ viễn thám và GIS với các phần mềm có các chức năng phân
tích, xử lý dữ liệu mạnh ngày càng được ứng dụng rộng rãi trong nhiều lĩnh
vực khác nhau. Trong những năm gần đây, GIS được sử dụng nh một công cụ
đắc lực trong công tác nghiên cứu, quản lý, khai thác cũng nh bảo vệ môi
trường. Các bản đồ quy hoạch, bản đồ xói mòn đất tiềm năng, bản đồ hiện
trang rừng, hiện trạng giao thông.v.v. có thể được thành lập bằng công cụ
GIS.
Đối với việc nghiên cứu vùng duyên hải nói chung và nghiên cứu đất
ngập nước nói riêng, công cụ GIS cũng được áp dụng, tuy chưa nhiều nhưng
đã mang lại những hiệu quả đáng kể. Từ năm 1990, GIS đã được một số tác
giả thuộc Viện Địa lý ứng dụng để đánh giá tài nguyên ven đồng bằng sông
Hồng. Năm 1996 tác giả Nguyễn Hoàn, Vũ Văn Phái và những người khác đã
tiến hành nghiên cứu biến đổi địa hình và quá trình hình thành các cồn bãi ở
khu vực cửa sông Hồng, trong đó sử dụng GIS như một trong các phương
pháp nghiên cứu chính để phân tích ảnh vệ tinh, phân tích số liệu địa hình
thuộc nhiều giai đoạn khác nhau để xác định biến động sử dụng đất và tài
nguyên theo không gian và thời gian.
Có thể nói việc kết hợp nghiên cứu tư liêu lịch sử, nghiên cứu, phỏng
vấn thực địa với ứng dụng viễn thám và GIS trong đánh giá, phân tích và dự
MỤC TIÊU, ĐỐI TƯỢNG, NỘI DUNG
VÀ PHƯƠNG PHÁP NGHIÊN CỨU
2.1. Mục tiêu nghiên cứu
2.1.1. Mục tiêu chung
Góp phần nâng cao hiệu quả quản lý sử dụng hợp lý tài nguyên rừng
ngập mặn.
2.1.2. Mục tiêu cụ thể
Ứng dụng công nghệ không gian địa lý GIS vào việc giải đoán ảnh viễn
thám thành lập bản đồ hiện trạng tài nguyên rừng tại thị xã Quảng Yên – Tỉnh
Quảng Ninh giai đoạn 1990 – 2015.
Dự báo xu thế biến động tài nguyên rừng và đất rừng ngập mặn trong
10 năm tiếp theo.
Để xuất giải pháp quản lý sử dụng hợp lý rừng và đất rừng ngập mặn
tại khu vực nghiên cứu.
2.2. Đối tượng điều tra khảo sát
Đối tượng điều tra khảo sát của đề tài là dải đất ngập nước ven biển
thuộc khu vực thị xã Quảng Yên – tỉnh Quảng Ninh bao gồm: Điều tra diện
tích rừng, đất ngập mặn của các năm từ 1990 – 2015 từ đó xác định sự biến
động tài nguyên rừng và đất ngập mặn. Cùng nằm trong tình trạng chung của
các vùng đất ngập nước trên thế giới cũng như ở Việt Nam, đây là khu hệ sinh
thái rất nhạy cảm, dễ bị biến đổi bởi các tác động nhân sinh cũng như tự
11
nhiên, đặc biệt là các tác động nhân sinh. Đất ngập nước ở khu vực này có thể
được sử dụng cho các mục đích khác nhau như rừng ngập mặn (mọc tự nhiên
hoặc trồng), ao đầm nuôi trồng thuỷ sản, hệ thống thuỷ văn, đường giao
thông, các bãi bùn hay các bãi bồi còn để trống.
2.3. Phạm vi nghiên cứu
Khai thác và sử dựng hợp lý tài nguyên rừng ngập mặn.
Trồng và phục hồi lại diện tích rừng ngập mặn.
Thay đổi phương pháp nuôi trồng thủy sản trong rừng ngập mặn.
Ngăn chặn các hành vi phá rừng, đất ngập mặn làm mục đích khác.
2.4.4.2. Giải pháp quản lý
Thành lập ban quản lý rừng ngập mặn chuyên trách.
Quy hoạch quản lý sử dụng tổng thể khu vực rừng và đất rừng ngập
mặn.
2.4.4.3 Giải pháp về mặt xã hội
Nâng cao nhận thức cho người dân địa phương trong công tác quản lý
bảo vệ rừng ngập mặn.
Đảm bảo kế sinh nhai cho người dân địa phương.
2.5. Phương pháp nghiên cứu
13
Một số phương pháp đã được sử dụng trong quá trình thực hiện đề tài
như: Thu thập và nghiên cứu tài liệu có liên quan đến vùng nghiên cứu,
phương pháp GPS, viễn thám và hệ thông tin địa lý (GIS), phương pháp điều
tra khảo sát thực địa. Trong đó GPS, Viễn thám và Hệ thông tin địa lý là công
cụ chính để thực hiện các công việc trong suất quá trình nghiên cứu. Các kiến
thức và kỹ thuật, kinh nghiệm về GPS, giải đoán ảnh viễn thám được sử dụng
trong quá trình liên kết dữ liệu (số hoá các đối tượng từ ảnh vệ tinh vùng
nghiên cứu) đầu vào với Hệ thông tin địa lý. Thực địa là bước quan trọng
nhằm kiểm chứng lại kết quả của công việc giải đoán để có thể đưa ra các bổ
sung, chỉnh sửa cần thiết.
2.5.1. Phương pháp kế thừa tài liệu
Thu thập và nghiên cứu các tài liệu các liên quan đến vùng nghiên cứu.
Tài liệu được thu thập chủ yếu được lưu trữ tại phòng kinh tế và hạt kiểm lâm
hiện cùng một đối tượng. Đó chính là quá trình phân lớp tự động.
2.5.2.1. Phương pháp xây dựng khóa giải đoán
Khóa giải đoán được thực hiện qua việc xác định phân tích hình ảnh theo
cấp phổ màu tại các khu vực nghiên cứu trong mỗi khoảng thời gian khác
nhau đối với mỗi trạng thái khác nhau. Để thành lập khóa giải đoán chúng ta
cần áp dụng nghiêm nghặt phương pháp lựa chọn khóa đặc trưng cho từng đối
tượng theo các nguyên tắc về sự khác nhau trong “tone” màu. Mỗi một đối
tượng khác nhau sẽ có sự phản xạ các bước sóng ánh sáng khác nhau, do vậy
chúng tạo ra những độ sắc nét, độ đậm nhạt khác nhau trong mỗi loại màu
sắc. Đây là cơ sở quan trọng để thành lập các khóa giải đoán.
15
Khóa giải đoán là tập hợp các chuẩn dùng để giải đoán một đối tượng
nhất định. Kết quả giải đoán phụ thuộc vào khóa giải đoán. Mục đích của việc
sử dụng khóa giải đoán là làm chuẩn hóa các kết quả giải đoán của nhiều
người khác nhau. Thông thường khóa giải đoán do những người có kinh
nghiệm và hiểu biết thành lập dựa trên những vùng nghiên cứu thử nghiệm đã
được điều tra kỹ lưỡng. Tất cả 8 chuẩn giải đoán cùng với các thông tin về
thời gian chụp, tỷ lệ ảnh, mùa chụp đều phải đưa vào khóa giải đoán.
Tư liệu ảnh dùng để giải đoán tốt nhất là ảnh tổng hợp màu. Đặc điểm cơ
bản của ảnh tổng hợp màu là sự mã hóa bằng màu sắc các khác biệt về phổ
của các đối tượng. Ở đây chuẩn giải đoán chính là sự tương phản màu được
nhấn mạnh nhờ sự lựa chọn một cách có ý thức phương án tổng hợp màu.
Trong trường hợp tư liệu gốc thỏa mãn các điều kiện kỹ thuật nếu sử dụng
phương án tổng hợp màu chuẩn và điều kiện xử lý hóa ảnh chặt chẽ thì màu là
một chuẩn giải đoán tương đối ổn định.
Nhờ khả năng phân biệt cao của màu sắc mà nó có thể truyền đạt các
khác biệt về phổ của đối tượng, ảnh tổng hợp màu có tính trực quan sinh động
thông tin địa lý cùng các dữ liệu phụ trợ khác. Kết quả giải đoán từng thời
điểm sẽ được phân tích bằng phương pháp tính bảng chéo (crossing table) để
tính ra biến động. Tích hợp các thông tin viễn thám với các thông tin kinh tế
xã hội khác sẽ giúp cho việc đánh giá mối liên quan giữa biến động tài
nguyên rừng, đất ngập mặn với các yếu tố thời gian và kinh tế xã hội mà ta
quan tâm.
2.5.2.4. Các bước tiến hành
17
2.5.2.4.1.. Thu thập số liệu
Các loại dữ liệu khác nhau được thu thập trong quá trình thực hiện đề
tài được phân chia như sau:
a. Số liệu không gian:
- Ảnh vệ tinh Landsat TM tổ hợp màu tự nhiên khu vực thị xã Quảng
Yên chụp từ năm 1990 - 2015 độ phân giải không gian 30x30 m
b. Số liệu phi không gian:
Đây là loại số liệu thuộc tính của khu vực cũng nh của đối tượng
nghiên cứu. Các báo cáo, văn bản, luận văn, tạp chí khoa học… có liên quan
tới vùng nghiên cứu được tham khảo để hình thành lên cái nhìn tổng quan về
khu vực nghiên cứu. Đó là các loại số liệu sau:
- Đặc điểm tự nhiên vùng nghiên cứu
- Điều kiện kinh tế xã hội của khu vực
- Số liệu về tình hình khai thác và sử dụng tài nguyên môi trường khu
vực cửa sông Hồng
Các loại số liệu trên là dữ liệu đầu vào cho Hệ thống thông tin địa lý.
Những dữ liệu này được nhập vào máy tính, trở thành cơ sở dữ liệu để các
phần mềm GIS có thể thực hiện được các chức năng phân tích, tính toán.
2.5.2.4.2. Nhập dữ liệu
định, vị trí thường ở trên các cồn cát hay các
bãi bồi phía ngoài đê
4
5
6
Khu
nuôi
trồng Màu xanh nhạt, hoa văn mịn, cấu trúc dạng
thuỷ sản (KNTTS)
mảnh hay thửa có bờ bao quanh
KNTTS + RNM
Sù xen kẽ của hai đối tượng, hoa văn chấm thô
Đất trống
Màu trắng (cát bồi) hoặc xám (cát và phù sa).
Tông ảnh sáng hoặc sẫm, hoa văn mịn.
Bãi bùn, cát ngập Màu trắng đục hơi đậm,hoặc đen xám hoa văn
7
triều
đồ hiện trạng.
b. Số hoá các đối tượng nghiên cứu
Công việc này được thực hiện chủ yếu bằng phần mềm AcrGis,
Mapinfo. Trước hết là gán cho ảnh vùng nghiên cứu các điểm toạ độ khống
chế, sau đó dựa trên chìa khóa giải đoán đã được xây dựng để số hoá các đối
tượng không gian trong ảnh, lập lên bản đồ vector hiện trạng sử dụng đất ở
khu vực nghiên cứu trong thời kỳ 1990 và 2015.
Trong quá trình số hoá, các dữ liệu thuộc tính liên quan như diện tích,
chu vi của đối tượng vùng, độ dài của đối tượng đường... được tính toán tự
động bằng máy và liên kết với các đối tượng. Đồng thời các đối tượng không
gian được phân loại và gán cho mã số riêng (ID). Chỉ số ID là một trường đặc
biệt của lớp (layer) thông tin trong AcrGis, Mapinfo cũng như trong các phần
mềm GIS khác. Nó được dùng để liên kết dữ liệu thuộc tính và dữ liệu không
gian cùng loại với nhau.
2.5.2.4.3. Chồng lớp và phân tích dữ liệu
Các đối tượng không gian trong ảnh của vùng nghiên cứu ở các năm
khác nhau sau khi được số hoá và gán thông tin thuộc tính thì được chuyển
sang phần mềm AcrGis để xử lý. Phần mềm GIS AcrGis với khả năng phân
tích không gian (spatial analys) mạnh và chức năng chồng lớp (overlayers) sẽ
đưa ra được kết quả biến động sử dụng đất giữa hai thời điểm 1990 và 2015:
20
- Các lớp thông tin không gian về các đối tượng trong vùng nghiên cứu
ở hai thời kỳ được chồng lên nhau bằng chức năng overlay để tìm ra quy luật
biến đổi các đối tượng không gian này sang đối tượng không gian khác giữa
hai thời kỳ.
- So sánh các số liệu thuộc tính nh diện tích, tên, ID của các đối tượng
không gian tại hai thời điểm 1990 và 2015. Thành lập bảng số liệu biến đổi và
Ranh giới:
+ Phía Bắc giáp thành phố Uông Bí và huyện Hoành Bồ.
+ Phía Nam giáp đảo Cát Hải và cửa Nam Triệu.
+ Phía Đông giáp thành phố Hạ Long và Vịnh Hạ Long.
+ Phía Tây giáp huyện Thủy Nguyên, thành phố Hải Phòng.
3.1.2. Địa hình, địa mạo
Riêng phần phía Bắc nằm trong khu vực chuyển tiếp của vùng núi
Đông Bắc nên mang tính chất của nhóm đất đồi núi. Theo đặc tính phân loại,
Quảng Yên có các nhóm đất chính sau:
+ Đất đồi núi có diện tích 6100 ha chiếm 18,3% diện tích, phân bố ở
khu vực phía Bắc thị xã tập trung ở các địa phương: phường Minh Thành,
phường Đông Mai và một phần ở xã Sông Khoai, phường Cộng Hòa, phường
Tân An, xã Tiền An, xã Hoàng Tân. Đất bao gồm chủ yếu là các loại đất
feralit vàng đỏ trên đá macma axit và đất feralit nâu vàng, xám vàng trên các
đá trầm tích phiến thạch, sa thạch, đá vôi. Đất có tầng dày trung bình 60 80cm, thành phần cơ giới từ thịt nhẹ đến trung bình, độ pH từ 4 - 4,5; hiện
chủ yếu là đất rừng và đất trồng cây ăn quả.
23
+ Đất đồng bằng có diện tích gần 13.528 ha chiếm 44 % diện tích đất
đai, gồm chủ yếu là đất phù sa cổ và đất phù sa cũ nằm trong đê, phân bố ở
hầu hết các xã, phường trong thị xã nhưng tập trung ở khu vực Hà Nam. Đất
có thành phần cơ giới trung bình đến nặng, độ PH dưới 4,5 hàm lượng mùn
trung bình. Một số điểm nội đồng đất trũng bị ngập nước mùa mưa nên lầy
mạnh, đất chua hàm lượng mùn thấp. Hiện đất được sử dụng chủ yếu để trồng
cây lương thực thực phẩm, trồng lúa hai vụ cho năng suất khá cao.
+ Đất bãi bồi cửa sông, ven biển gồm các loại đất mặn và đất cát có
diện tích gần 12.300 ha chiếm 37,1% diện tích, phân bố ở các khu vực ven
biển và cửa sông, tập trung nhiều các khu đầm Nhà Mạc, đầm Soài, Cái Tráp,
hợp đối với sản xuất nông nghiệp do nước bị nhiễm mặn.
Quan trọng nhất là dòng chính của sông Bạch Đằng chảy ở phía Tây
ngăn cách Thị xã Quảng Yên với Thành Phố Hải Phòng và các chi lưu chảy
vào thị xã là sông Chanh, sông Nam, các sông này đều đổ ra biển ở khu cửa
Nam Triệu - Lạch Huyện. Phần phía Đông thị xã còn có một số sông nhỏ
khác như sông Hốt, sông Bến Giang và sông Bình Hương nhưng các sông này
đều ngắn, diện tích lưu vực nhỏ, chủ yếu trong phạm vi thị xã.
Nguồn nước ngầm khá phong phú, mạch nước ngầm thường nằm ở độ
sâu 5 - 6 mét, khu vực Hà Nam và ven biển nước bị nhiễm mặn ít sử dụng
được, khu vực Hà Bắc nước ngọt đủ để khai thác sử dụng cho sinh hoạt.
25
Thuận lợi lớn nhất về nguồn nước của Thị xã là có hồ Yên Lập, là hồ
thuỷ lợi lớn của tỉnh có dung tích thường xuyên 127, 5 triệu m3, dung tích
hữu ích 113, 2 triệu m3 với kênh chính dẫn nước cho Thị xã dài 28,4 km.
Nguồn cấp nước từ hồ Yên Lập dồi dào và hiện là nguồn nước chủ yếu cho
sản xuất và sinh hoạt của dân cư trong thị xã, đặc biết là khu vực Hà Nam
luôn được chống xâm nhập mặn.
3.1.5. Các nguồn tài nguyên khác
3.1.5.1. Tài nguyên du lịch
- Về tuyến du lịch có 3 tuyến bao gồm: Tuyến Trung tâm TX Quảng
Yên (Bảo tàng Bạch Đằng, Hai cây lim Giếng Rừng, Bãi cọc Bạch Đằng, Đền
thờ Trần Hưng Đạo, miếu Vua Bà, Bến đò Cổ, chợ Rừng); Tuyến Trung tâm
TX Quảng Yên - Đảo Hà Nam (điểm di tích lịch sử Bạch Đằng, làng nghề
đan ngư cụ truyền thống thôn Hưng Học thuộc phường Nam Hòa, đi thuyền
nan trên sông cửa đình Cốc, đình Phong Cốc, nhà thờ họ Lê); Tuyến Trung
tâm TX Quảng Yên- phường Đông Mai, Minh Thành (điểm Bác Hồ dừng
chân, Rừng Thông, Thác Mơ).