Soil provenance using ATR-FTIR spectroscopy augmented with chemometrics–A pilot study

기술 요약: 화학계량학을 결합한 ATR-FTIR 분광법을 이용한 토양 기원 판별

문제 제식
법의학적 토양 검사는 범죄 현장, 피해자, 용의자를 연결하고 지리적 기원을 결정하며 범죄 현장을 재구성하는 데 있어 중요한 구성 요소이다. 토양 분석은 시각적 및 미세 구조적 검사에서 정교한 유기 및 무기 화학 분석으로 진화해 왔으나, 가스 크로마토그래피-질량 분석법(GC-MS), 고성능 액체 크로마토그래피(HPLC), 원자 흡수 분광법과 같은 전통적인 기술은 비용이 많이 들고 시간이 오래 걸리며 광범위한 시료 준비가 필요하다. 또한, 토양은 다양한 지리적, 생물학적, 인위적 요인의 영향을 받는 유기 및 무기 성분의 복잡한 혼합물이기 때문에, 서로 다른 지역의 샘플을 구별하는 것은 매우 까다로운 작업이다. 따라서 샘 samples의 조성 프로파일을 기반으로 샘플을 구별할 수 있는 신속하고 비파괴적이며 비용 효율적인 방법이 필요하다.

방법론
본 파일럿 연구에서는 다양한 지리적 위치의 토양 샘플을 구별하기 위해 화학계량학 도구와 결합된 감쇄 전반사 푸리에 변환 적외선(ATR-FTIR) 분광법을 사용하였다.

• 시료 채취 및 준비: 인도 전역의 15개 서로 다른 지리적 지역(잠무 카슈미르, 히마찰 프라데시, 펀자브, 델리 및 케랄라의 여러 지역 포함)에서 총 90개의 토양 샘플을 수집하였다. 각 지역에서 6개의 대표 샘플을 수집하였다. 샘플은 공기 건조 후 이물질을 제거하기 위해 2mm 메쉬로 체질하였으며, 분석을 위해 보관하였다.
• 기기 분석: 스펙트럼 데이터는 다이아몬드 결정이 장착된 Nicolet iS20 ATR-FTIR 분광광도계를 사용하여 획득하였다. 측정은 중간 적외선 영역(4000–600 cm⁻¹)에서 4 cm⁻¹의 해상도로 수행되었으며, 샘플당 평균 32회의 스캔을 실시하였다. 균질한 압력을 보장하기 위해 샘플을 결정에 밀착시켰으며, 스캔 사이에는 분광 등급의 아세톤으로 결정을 세척하였다.
• 화학계량학적 분석: 스펙트럼 데이터는 베이스라인 보정, 평활화(Savitzky–Golay), 평균 중심화(mean centering)를 통해 전처리되었다. 두 가지 다변량 통계 기법이 적용되었다:
  • 주성분 분석 (PCA): 스펙트럼 공간 분포를 시각화하고 자연스러운 클러스터링을 식별하기 위해 사용된 비지도 차원 축소 기법이다.
  • 부분 최소 제곱 판별 분석 (PLS-DA): 스펙트럼 특성을 기반으로 샘플 분류 모델을 개발하기 위해 사용된 지도 학습 분류 방법이다.
• 검증: 모델 성능은 교차 검증(무작위 20개 세그먼트), 평균 제곱근 오차(RMSE), 분류 정확도를 사용하여 평가되었다.

주요 결과

• 스펙트럼 특성: 대표 스펙트럼은 1010–1045 cm⁻¹ 및 870–954 cm⁻¹ 범위에서 주요 피크를 나타냈다. 1033 cm⁻¹ 피크는 Al-OH 변형(카올리나이트) 및 방향족 C-H 굽힘에 의한 것으로 간주되었으며, 870–945 cm⁻¹ 범위는 벤조산, 피라노스 고리 및 셀룰로오스에 해당했다. 본 연구는 높은 반복성과 재현성을 확인하였으며, 표준 편차는 각각 ±0.18 및 ±0.19였다.
• PCA 성능: PCA 모델은 제1주성분(PC1)이 전체 분산의 99%를 설명함을 밝혔다. 점수 플롯(score plot)은 북부 지역(Class 1–6) 샘플에 대해 뚜렷한 클러스터링을 보여주었으며, 케랄라 지역(Class 7–15) 샘플은 지리적 근접성과 유사한 지질학적 형성으로 인해 조밀한 클러스터를 형성했다. 케랄라 하위 집합에 대한 2차 PCA 분석은 구별 능력을 향상시켰다. PCA 모델의 전체 분류 정확도는 **86.11%**였다. 오분류는 주로 지리적으로 가까운 케랄라 지역 사이에서 관찰되었다(예: Class 8이 Class 12로 오분류됨).
• PLS-DA 성능: PLS-DA 모델은 변별을 위한 세 가지 유의미한 요인을 식별하였으며, 제1요인이 분산의 99%를 설명했다. PCA와 유사하게, 모델은 북부 샘플을 효과적으로 분리했으나 케랄라 샘플은 밀집된 클러스터를 형성했다. 케랄라 하위 집합(샘플 7–15)을 위해 타겟팅된 PLS-DA 모델이 개발되었으며, 이는 캘리브레이션 과정에서 높은 R² 값(0.9838)과 낮은 RMSE(0.308)를 나타냈다. 전체 PLS-DA 모델의 분류 정확도는 **85.56%**를 달agt했다.
• 오분류 패턴: 오분류는 주로 동일한 주 내의 인접한 위치 간에 발생했다. PLS-DA 모델에서 Class 1과 Class 7은 정확도가 0%였으며(모든 샘플이 각각 Class 4와 Class 8로 오분류됨), Class 12는 83.33%의 정확도를 보였다(한 개의 샘플이 Class 10으로 오분류됨). Class 2, 3, 4, 5, 6, 8, 9, 10, 11, 13, 14, 15는 100%의 정확도를 달성했다.

주요 기여

• 방법론적 검증: 본 연구는 ATR-FTIR 분광법이 화학계량학을 통해 보완될 때, 전통적인 법의학적 토양 분석 기술을 대체할 수 있는 유효한, 비파괴적이며 신속한 대안이 될 수 있음을 입증한다.
• 차별화 능력: 본 연구는 15개의 뚜렷한 지리적 지역으로부터의 토양 샘플을 성공적으로 차별화하였으며, 높은 분류 정확도(PCA 86.11%, PLS-DA 85.56%)를 달성하였다.
• 타겟 모델링: 본 연구는 초기 광범위 모델이 놓칠 수 있는 미세한 스펙트럼 차이를 해결하기 위해 지리적으로 밀집된 샘플(케랄라와 같은)에 대한 2차 타겟 모델을 개발하는 것의 유용성을 강조한다.
• 스펙트럼 마커: 특정 진동 피크(예: 1000 cm⁻¹ 부근의 Si-O 격자 신축)의 식별은 샘플 구별을 담당하는 화학적 구성 성분을 해석하는 기초를 제공한다.

의의 및 주장
저자들은 본 파일럿 연구가 ATR-FTIR 분광법과 화학계량학적 분석의 결합이 법의학적 맥락에서 토양 샘플의 지리적 구별을 위한 신뢰할 수 있는 도구로서의 잠재력을 효과적으로 보여준다고 주장한다. 연구는 개발된 모델이 높은 R-제곱 값과 낮은 RMSE를 바탕으로 견고함을 입증하여 실제 사례 적용에 적합하다고 단언한다. 지리적으로 인접한 지역의 샘플들이 유사한 스펙트럼 프로파일을 공유하여 미세한 오분류가 발생할 수 있음을 인정하면서도, 저자들은 이 접근 방식이 더 저렴하고 빠르며 시료 준비가 적게 든다는 점에서 기존 방법보다 상당한 이점을 제공한다고 결론짓는다. 이 연구는 범죄 현장, 용의자, 피해자를 연결하고 출처를 규명하는 데 있어 이러한 기술의 활용을 뒷받침한다.