10. Lời cảm ơn
3.2.2.Đánh giá diện tích biến ñộ ng sử dụng ñấ t theo không gian
• Biến ñộng sử dụng ñất giai ñoạn 2001-2013 theo ñường giao thông chính Trong bảng 3.10 cho thấy diện tích ñất xây dựng và diện tích ñất ao hồ tăng trên toàn bộ các biên khoảng cách tới ñường giao thông chính. Đối với ñất xây dựng, càng gần ñường giao thông thì diện tích càng biến ñộng mạnh theo chiều hướng tăng, mức tăng cao nhất là 1121.88ha nằm trong biên khoảng cách tới ñường giao thông
chính nh h n 1km. t chuyên lúa, hoa màu và màu lúa giảm t t cả các biên khoảng cách t i ng giao thông, trong biên khoảng cách n ng giao thông d i 1km t lúa giảm m nh nh t lên t i 898.22ha. t màu lúa giảm nhi u nh t trong biên khoảng cách n ng giao thông t 1-2.5km v i 500.16ha di n tích.
Bảng 3. 10. Biến ñộng diện tích (ha) sử dụng ñất theo ñường giao thông chính
Khoảng cách S d ng t (ha) <=1 1_2.5 >2.5 Xây d ng 1121.88 924.60 414.51 Lúa -898.22 -489.35 -422.70 Hoa màu -46.38 -74.48 -83.56 Màu lúa -396.36 -500.16 -76.84 Ao h 224.14 95.29 146.56 t tr ng 4.48 75.18 19.72 Sông -12.66 -32.84 7.18
• Bi n ng s d ng giai o n 2001-2013 t theo khoảng cách t i Hà N i
Bảng 3. 11. Biến ñộng diện tích (ha) theo khoảng cách tới Hà Nội
Khoảng cách (km) S d ng t (ha) <=5 5- 10 >10 Xây d ng 653.54 1160.11 647.35 Chuyên lúa -553.25 -848.38 -408.64 Hoa màu -63.31 -83.99 -57.12 Màu lúa -140.28 -461.20 -371.87 Ao, h 91.40 189.67 184.33 t tr ng 56.38 42.85 0.14 Sông -46.29 4.93 3.04
Hình 3. 13. Biểu ñồ diện tích biến ñộng theo khoảng cách tới Hà Nội
Theo bảng 3.11 và hình 3.13 cho th y di n tích t xây d ng t ng lên m nh nh t trong ai khoảng cách t 5 n 10 km t i Hà N i v i di n tích là 1160.11 ha, ng th i trong ai khoảng cách này thì di n tích t chuyên lúa là di n tích t b m t nhi u nh t v i con s là 848.38ha. Trong ai khoảng cách t i Hà N i nh h n 5km, di n tích t lúa c ng b m t nhi u nh t lên t i 553.25 ha. t màu lúa b giảm m nh nh t nh ng khu v c cách Hà N i khoảng cách h n 10 km v i di n tích là 371.87 ha.
• Bi n ng s d ng t giai o n 2001 - 2013 theo m c ô th hóa n m 2011 T bảng 3.12 và hình 3.14 cho th y di n tích t xây d ng t ng theo t t cả
các m c ô th hóa và di n tích t xây d ng t ng nhi u nh t khi m c ô th hóa trong khoảng t 0.4÷1.2, nh ng theo bi u ph n tr m di n tích bi n ng th hi n thì m c ô th hóa càng m nh, t xây d ng càng bi n i nhi u và t ng m nh nh t là khu v c có m c ô th l n h n 3.3 v i m c t ng lên t i h n 30% so v i t ng di n tích c a các lo i t. Ng c l i v i t xây d ng, t chuyên lúa liên t c giảm. Theo bảng 3.12, di n tích t lúa c ng giảm m nh nh t m c ô th hóa trong khoảng 0.4÷1.2, nh ng theo bi u ph n tr m di n tích bi n ng thì di n tích t lúa l i giảm m nh nh t khi m c ô th hóa cao nh t.
D i n tí ch ( h a) Khoảng cách -1000.00 -500.00 0.00 500.00 1000.00 1500.00 <=5km 5- 10km >10km
Bảng 3. 12. Diện tích biến ñộng sử dụng ñất theo mức ñộñô thị hóa năm 2011 M c ô th hóa S d ng t (ha) <=-0.2 -0.2÷0 0÷0.4 0.4÷1.2 1.2÷3.3 >3.3 Xây d ng 296.64 268.20 669.81 766.93 124.27 339.41 Chuyên lúa -224.89 -149.46 -558.66 -572.12 -48.63 -259.18 Hoa màu -45.32 -42.60 -69.73 -56.73 2.11 0.00 Màu lúa -113.07 -162.14 -171.90 -367.57 -65.72 -86.66 Ao, h 61.78 62.18 121.87 194.06 18.44 6.42 t tr ng 19.91 10.60 25.27 42.06 4.60 0.00 Sông 4.95 13.21 -16.66 -6.64 -35.07 0.00 Hình 3. 14. Biểu ñồ phần trăm diện tích biến ñộng theo mức ñộñô thị hóa
3.3. Xác ñịnh các yếu tốảnh hưởng ñến biến ñộng sử dụng ñất nông nghiệp tại
Đông Anh, Hà Nội
Khi nghiên c u bi n ng s d ng t ngoài vi c phát hi n ra không gian và di n tích bi n ng s d ng t thì vi c nghiên các y u t tác ng n bi n ng s d ng t c ng là nhi m v quan tr ng c r t nhi u nhà khoa h c quan tâm.
S d n g t (%) M c ô th hóa -40.0 -30.0 -20.0 -10.0 0.0 10.0 20.0 30.0 40.0 <=-0.2 -0.2 ÷ 0 0 ÷ 0.4 -0.4 ÷ 1.2 -1.2 ÷ 3.3 >3.3
Nghiên c u y u t ảnh h ng n bi n ng s d ng t nh m m c ích l ng hóa c các y u t ảnh h ng n bi n ng s d ng t và giải thích c ch t ng tác trong m i quan h gi a chúng v i bi n ng s d ng t.
3.3.1. Hồi quy logistic
K t quả phân tích ch ng 1 cho th y có r t nhi u ph ng pháp ánh giá các y u t tác ng n bi n ng s d ng t, m i lo i hình chuy n i s d ng t s có nh ng y u t tác ng riêng r , nh ng các lo i hình chuy n i s d ng t l i t ng tác v i nhau và xảy ra ng th i. Vì vậy, việc lựa chọn phương pháp cho phép phân tích quan tâm ñến tính phù hợp của dữ liệu và ñộ chính xác của mô hình. Dựa vào những yêu cầu ñó, luận án ñã lựa chọn mô hình hồi quy logistic ña thức ñánh giá tác ñộng của các yếu tố ñô thị hóa ñến biến ñộng sử dụng ñất. Mô hình hồi quy ña thức là mô hình hồi quy mở rộng của mô hình hồi ñơn thức [70].
Mô hình ñơn thức với biến phụ thuộc là biến nhị phân với hai trạng thái 0 và 1 có phương trình hồi quy sau:
Ε Φ Γ Η Ι ϑ (3. 1)
trong ñó: là hệ số chặn, Κ là hệ số của phương trình hồi quy, xi là các biến ñộc lập
Xác suất ñể Y = 1 nghĩa là xác suất sự kiện xảy ra ñược xác ñịnh theo công thức:
(ΛΜΝΟΠΘΡΣΤΥςΩ,ΞΨΖ[∴) (]⊥_αβχδεφγη,ιϕκλµ)
ν(ο)
π(θ) (3. 2) Mô hình hồi quy ña thức tương tự như mô hình hồi quy ñơn thức nhưng biến phụ thuộc là biến có nhiều hơn hai trạng thái. Giả sử biến Y có k trạng thái là 0, 1, 2,..., k, ñòi hỏi sẽ có k-1 phương trình hồi quy, mỗi phương trình của từng trạng thái so với trạng thái tham chiếu 0. Các phương trình của hồi quy ña thức so với trạng thái tham chiếu 0 như sau:
} ∼ ! ∀ #∃ % &∋ ( (3. 4) Và
) ∗ +, − ./ 0 12 3 (3. 5)
Xác su t xảy ra c a t ng tr ng thái c xác nh theo công th c:
45(6) 78(9) :;(<) (3. 6) =>(?) ≅Α(Β) Χ∆(Ε) ΦΓ(Η) (3. 7) Và Ιϑ(Κ) ΛΜ(Ν) ΟΠ(Θ) ΡΣ(Τ) (3. 8) Trong h i quy logistic s d ng ph ng pháp xác su t c c i (maximum likelihood) c tính các tham s c a ph ng trình h i quy nh m cho giá tr t i a hóa khả n ng c a các bi n quan sát. (adsbygoogle = window.adsbygoogle || []).push({});
phù h p c a mô hình c xác nh thông qua pseudo R- squared (kí hi u là R2) có giá tr trong khoảng t 0 - 1. T s ó c xác nh theo công th c sau:
Υ (3. 9)
trong ó:
- ] là giá tr c a mô hình có bi n c lập,
- [C] là giá trị của mô hình không có biến ñộc lập.
Giá trị R2 trong khoảng từ 0.2 ñến 0.4 cho thấy mức ñộ phù hợp cao của dữ liệu với mô hình.Trong việc ñánh giá các yếu tố tác ñộng ngoài tính phù hợp của mô hình thì còn ñòi hỏi xem xét ñến hệ số hồi quy có ý nghĩa khi p-value<0.05 với mức ñộ tin cậy 95%, p-value<0.1 với mức tin cậy 90% [71, 139].
T s Odds ratio c s d ng giải thích mô hình. Giá tr t s này c tính theo công th c sau:
ςΩ
(3. 10) trong ó: Ξ là h s c a ph ng trình h i quy
T s này là m t phép o m i quan h xác su t bi n c c a t ng tr ng thái v i xác su t c a bi n c so sánh và là k t quả óng góp c a t ng bi n c lập ñến từng trạng thái.
Trong nghiên cứu biến ñộng sử dụng ñất, mô hình hồi quy ñược sử dụng ñể ño ñạc mối quan hệ giữa các yếu tố tự nhiên và yêu tố xã hội với biến ñộng sử dụng ñất, ñồng thời xác ñịnh ñược những yếu tố nào là tác nhân thúc ñẩy và những yếu tố nào làm hạn chế tới quá trình chuyển ñổi sử dụng ñất trong khu vực. Mô hình hồi quy ñược sử dụng với biến phụ thuộc Y là biến ñộng sử dụng ñất ñược ñặc trưng bởi 5 trạng thái thay ñổi sử dụng ñất chính của khu vực nghiên cứu so với trạng thái tham chiếu là không có thay ñổi sử dụng ñất có giá trị bằng 0. Số lượng mẫu phụ thuộc vào số lượng biến ñộc lập, số lượng mẫu tối thiểu ñược xác ñịnh là n≥ 10(p+1) trong ñó p là số lượng biến ñộc lập. Trong luận án, số ñiểm mẫu ñược lấy là 4984 ñiểm ngẫu nhiên và trải ñều trên khu vực nghiên cứu. Các ñiểm nằm trong thôn nào thì sẽ mang ñặc ñiểm tự nhiên và kinh tế xã hội của thôn ñó.
Hình 3. 15. Xác ñịnh các yếu tố của ảnh hưởng ñến biến ñộng sử dụng ñất nông nghiệp
Các yếu tố tác ñộng ñến biến ñộng sử dụngñất nông nghiệp
Yếu tốñịa lý, tự nhiên
Mức ñộñô thị hóa
Biến ñộng
sử dụng ñất
nông nghiệp Mức ñộ hoạt ñộng nông nghiệp
Điều kiện kinh tế - xã hội
3.3.2. Cơ sở dữ liệu GIS của mô hình
Ph n t ng quan ch ng 1 cho th y có r t nhi u nhân t tác ng n thay i s d ng t. Trong tr ng h p nghiên c u cho khu v c Đông Anh, d a trên c tr ng ô th hóa c a khu v c ven ô l a ch n các y u t ánh giá.
• Bi n ph thu c:
Các lo i hình bi n ng s d ng t chính c xác nh t trong m c 3.2 theo hai giai o n 2001-2006 và 2006-2013.
• Bi n c lập:
Khi ñánh giá tác ñộng của ñô thị hóa ñến sử dụng ñất nông nghiệp trong khu vực luận án ñánh giá từ tổng thể của ñô thị hóa ñến chi tiết của từng hợp phần của ñô thị hóa ñến biến ñộng sử dụng ñất. Do vậy, ngoài các biến tổng hợp của ñô thị hóa thì các biến hợp phần bao gồm các biến về ñiều kiện kinh tế xã hội và các biến về ñiều kiện tự nhiên cũng ñược xem xét trong nghiên cứu. Các biến ñộc lập ñược xây dựng thành các lớp trong cơ sở dữ liệu GIS gồm có:
- Mức ñộ ñô thị hóa (F1) và mức ñộ hoạt ñộng nông nghiệp (F2) của từng thôn ñược xác ñịnh thông qua phân tích thành phần chính trong phụ lục 5.
- Khoảng cách ñến ñường giao thông chính (GT) ñược tính toán thông qua một tập hợp các vùng ñệm liên tiếp có khoảng cách 300 m từ hệ thống ñường giao thông trong khu vực nghiên cứu. Khoảng cách này là khoảng cách ñược ño trung bình từ các khu dân cư bám trục giao thông chính.
- Khoảng cách ñến Hà Nội (HN) ñược tính toán thông qua tập hợp vùng ñệm liên tiếp có khoảng cách 3 km tính từ ñiểm khởi ñầu nối kết giữa Hà Nội với Đông Anh qua trục giao thông Nội Bài - Hà Nội. Khoảng cách này ñược xác ñịnh khi ño từ ñiểm khởi nối ñến khu công nghiệp Nam Thăng Long, là một trong những khu công nghiệp ñầu tiên và lớn nhất của huyện.
- Khoảng cách ñến khu dân cư (DC) ñược xác ñịnh thông qua tập hợp vùng ñệm liên tiếp có khoảng cách 200m tính từ trọng tâm của khu vực dân cư mỗi thôn [129].
Hình 3. 16. Yếu tố tự nhiên của khu vực nghiên cứu (adsbygoogle = window.adsbygoogle || []).push({});
- Thay i c c u lao ng theo ngành c xác nh t hi u s lao ng theo ngành nông nghi p (P_NN), theo ngành th ng m i d ch v (P_DV) và theo ngành công nghi p xây d ng(P_CN) giai o n 2001 n 2006 và 2006 n 2011.
- Thay i c c u ngu n thu nhập theo ngành ñược xác ñịnh từ hiệu số hộ thu nhập theo ngành nông nghiệp (I_NN), theo ngành thương mại dịch vụ (I_DV) và theo ngành công nghiệp xây dựng(I_CN) ở hai giai ñoạn từ năm 2001 ñến 2006 và từ năm 2006 ñến 2011.
- Thay ñổi mật ñộ dân số (DPP) ñược xác ñịnh từ hiệu mật ñộ dân số giai ñoạn 2001 ñến 2006 và 2006 ñến 2011.
ñường giao thông chính 0 2 Km 0 - 300m 300 - 600m 600 - 900m 900 - 1200m 1200 - 1500m > 1500m
Khoảng cách tới ñường giao thông
ñường giao thông chính 0 2 Km Khoảng cách tới Hà Nội 0 - 3km 3 - 6km 6 - 9km 9 - 12km 12 - 15km > 15km
Khoảng cách tới khu dân cư
ñường giao thông chính 0 2 Km 0 - 200m 200 - 400m 400 - 600 600 - 800m > 800m
Bảng 3. 13. Biến phụ thuộc và biến ñộc lập
Loại biến Kiểu
Biến phụ thuộc
Không chuy n i 0
Chuy n i t t chuyên lúa thành t xây d ng 1
Chuy n i t t chuyên lúa thành t màu lúa 2
Chuy n i t t chuyên lúa thành t ao h 3
Chuy n i t t màu lúa thành t xây d ng 4
Chuy n i t t màu lúa thành t lúa 5
Biến ñộc lập
M c ô th hóa (F1) R i r c
M c ho t ng trong l nh v c nông nghi p (F2) R i r c
Đi u ki n t nhiên
Khoảng cách n Hà N i Liên t c
Khoảng cách n ng giao thông chính Liên t c
Khoảng cách n khu dân c Liên t c
Đi u ki n kinh t xã h i
C c u ngành lao ng R i r c
C c u ngu n thu nhập Rời rạc
Mật ñộ dân số Rời rạc
3.3.3. Các yếu tố tác ñộng chuyển ñổi sử dụng ñất nông nghiệp tại Đông Anh, Hà Nội
Kết quả chạy mô hình hồi quy cho hai giai ñoạn 2001-2006 và 2006-2013 xác ñịnh các tác nhân ảnh hưởng tới biến ñổi sử dụng ñất cho khu vực. Kết quả chạy mô hình ñược trình bày trong phụ lục 6. (adsbygoogle = window.adsbygoogle || []).push({});
• Giai ñoạn 2001-2006
T k t quả c a mô hình v i 6 bi n ph thu c v i giá tr R2 = 0.240, cho th y d li u hoàn toàn phù h p v i mô hình. Giá tr p-value là giá tr mang ý ngh a minh ch ng th ng kê cho vi c có ho c không có s ảnh h ng c a các y u t n chuy n i c a các lo i hình s d ng t trong giai o n này. V i giá tr P-value 0.1 i u này ch ng minh có m i quan h mật thiết giữa biến ñộc lập và biến phụ thuộc với mức ý nghĩa 90%.
Bảng 3. 14. Tổng hợp hệ số của phương trình hồi quy từ kết quả chạy mô hình giai ñoạn 2001 -2006 với giá trị P-value ≤0.1
Y X Lúa – xây dựng Lúa – Màu lúa Lúa – Ao, hồ Màu lúa - Xây dựng Màu lúa – Lúa F1_06 0.0976 0.1456 F2_06 0.0471 -0.3076 DPP 0.0003 GT -0.0030 -0.0006 DC -0.0037 -0.0003 0.0002 -0.0046 HN -0.0002 0.0001 0.0001 0.00001
Bảng tổng hợp kết quả chạy mô hình hồi quy cho thấy mức ñộ ñô thị hóa (F1_06) liên quan ñến hai loại hình chuyển ñổi ñó là chuyển ñổi từ ñất chuyên lúa sang ñất xây dựng và chuyển ñổi từ ñất màu lúa sang ñất xây dựng. Đối với cả hai loại hình chuyển ñổi này ñều có hệ số của phương trình hồi quy mang dấu dương thể hiện quan hệ thuận giữa mức ñộ ñô thị hóa với chuyển ñổi từ ñất chuyên lúa sang ñất xây dựng và chuyển ñổi từ ñất màu lúa sang ñất xây dựng. Nghĩa là khi mức ñộ ñô thị hóa tăng lên thì khả năng chuyển ñổi của hai loại hình này cũng tăng lên. Ví dụ: nếu mức ñộ ñô thị hóa tăng lên 1 ñơn vị thì giá trị odds ratio của loại
hình chuy n i t t chuyên lúa sang xây d ng6 là 1.1025, có ngh a k v ng khả
n ng chuy n i t ng 1.1l n.