Kristin Kinic, MS


(505)-718-0193

kristinmkinic@gmail.com


Kristin Kinic's GIS mapping, analysis and cartography portfolio



Cartography: For my thesis, Acequias as Service Ecologies, I created study area maps and Water Flow accumulation graphics.

Flow Accumulation Model: San Augustine


Water Flow Accumulation Analysis


 Study Area Map


Inspiration to study Hydrologic Processes in Las Vegas, New Mexico

The Las Vegas area has several water impoundments, some healthy acequias, remnant acequias. These NAIP images show two different water management tactics: Storrie Lake and the Middle Gallinas Acequias.
Storrie Lake Comparison using NAIP Images

Acequias as Service Ecologies


Watersheds of the Las Vegas NM area with rivers over a 10 meter DEM

Preliminary Research for Master's Degree

NDVI analysis of San Luis, Colorado

JPG files were created by author from NAIP images from www.earthexplorer.usgs.gov with special permission

 


Climate change and water scarcity will bring about circumstances in which voluntary, self-organized irrigation systems will be increasingly important,(Rodriquez, S,et al, 2020). Small scale irrigation projects are often associated with meaningful riparian areas, (Arellano, J.A.2014, Fernald, et al, 2007, Raheem, et al. 2015, Rodriquez, et. al. 2020, SER, 2017).

Globally, restoration best practices include incorporating small scale water systems that are managed locally. Agroecological practices are increasingly part of water infrastructure plans, especially in places affected by drought and seasonal flooding.

The United Nations has an assessment protocol for ecosystems called the Millennium Ecosystem Assessment (MEA). An assessment of acequias(heritage irrigation systems) in the Upper Rio Grande Valley (URGV), using the MEA was completed in 2015, (Raheem,2015).

To wit, acequias are part of communal land use arrangement that includes the watershed uplands, its drainage, and its forested headwaters, (Arellano,2014).

The MEA is established in terms of ecological services in four areas. The acequia system serves all four aspects of the framework. First of four ecological services is provisioning:acequias provide fresh water and food. Second is regulating:acequias recharge aquifers. Third, supporting: acequias enrich soil. Finally, acequias incorporate culturally relevant practices as they are maintained and governed voluntarily, communally, and democratically. In the URGV region, some acequias have been tied to increases ground water recharge throughout the growing season from seepage from lateral ditches.The authors of the MEA assessment called for extensive mapping of acequias as a "next step". Maps created by acequia managers can be enhanced by using geospatial information systems,(GIS), and normalized difference vegetation indices (NDVI).

 

Screen Shot from www.earthexplorer.usgs.gov. Imagery ID m_31755543_13_060_20190912 with user permissions

Drawing by Oliver La Forge, in The Mother Ditch

Screen shot of NAIP image selection process from www.earthexplorer.usgs.gov

Area of Analysis. San Luis Valley, Colorado

Map by Kristin M Kinic


Superscript

What is NDVI (red)

NDVI: Natural Vegetation Difference Index

NDVI has a fixed scale, from 1 to -1. NDVI is calculated by isolating elements of electromagnetic spectrum inherent in aerial and satellite imagery and then layering those elements to filter out "noise" to better visually appreciate healthy vegetation. In NAIP imagery, there are 4 bands readily available for this analysis, and I will be using two of them, in raster form, band 1(red) and band 4 (near infra-red or NIR). The numerical value of each raster will be used in this formula, NDVI(red)= NIR-RED/NIR+RED.

Screen shot of NAIP image selection process from www.earthexplorer.usgs.gov

data

NDVI ANALYSIS OF SAN LUIS VALLEY

NDVI (red)

Statistical Hypotheses:

 H/0 = µNDVI/ace = µNDVI/cpi

H/A = µNDVI/ace > µNDVI/cpi

H/0= There is no difference in the indices of chlorophyll bearing plants in agricultural areas that have been irrigated for maximum density, in paired samples of acequia sites, and overhead pivot irrigation sites in related watersheds.

H/A= Crops that are watered by acequias or central pivot irrigation systems have different NDVI values.

 METHODS: Identified potential agricultural sites for NDVI analysis using NAIP aerial photography of agricultural areas in Southern Colorado and Northern New Mexico. I will use remote sensing to find n= 20 study areas, with paired samples. The pairs will consist of two types of agricultural sites: those with acequias and those with center pivot irrigation. On site assessments of acequia ditches will be conducted to establish appropriate acequia samples, to inventory beneficial, non-cultigen flora. Examples of beneficial flora are red willow, and seges. Water temperatures in ditches will also be taken for future research. To ensure rigor, I will randomize sample points in circular plots with radii of 200m, in Arc Map ™. I will then use these points to calculate NDVI to describe the health and prosperity of cultigens on the ground of both types of irrigation: acequia and pivot irrigation. NDVI has a fixed scale, 1 to -1. NDVI is calculated by isolating elements of electromagnetic spectrum inherent in imagery and then layering those elements to filter out "noise" to better appreciate healthy vegetation. In NAIP imagery, there are 4 bands readily available for this analysis, and I will be using two of them, in raster form, band 1(red) and band 4 (near infra-red or NIR). The numerical value of each raster will be used in this formula, NDVI= NIR-RED/NIR+RED.

 


Steps leading to statistical Analysis using RStudio

  1. create uniform areas that are >4500<5000m^2

  2. separate raster data into bands in Arc Map

  3. compute NDVI in Arc Map

  4. create random samples in the areas in ArcMap

  5. collect numerical data

  6. save numerical data in .csv format

  7. import .csv into R Studio

  8. Use procedural flowchart to select a test

  9. test assumptions of that test

  10. test normality of data and homogeneity of variance

  11. run tests to either a. reject the null hypothesis and b. accept the the alternate hypothesis.




A mixed agricultural area in the San Luis Valley, CO

Results

This is the results of the statistical analysis
A mixed agricultural area in the San Luis Valley, CO

R outputs (below) show some of the process of the statistical analysis. On the left, is one of the 6 theoretical quantile plots I used to visually assess the data. This one is the assessment of the normality of the data from an average of NDVI points from the acequia study areas in the sample part called pixel 1. I sampled each site 3 times, the data sets were named pixel 1, pixel 2, pixel 3. While there are probably better way top sample in NDVI. I used the Pixel Inspector tool that gives the value of pixels in a tabular data frame. I simple entered the first three columns of point (pix 1,2,and 3 in the code) into a csv file, for n=10 sites. In total I had 900 points, n=10(30x3). The points were then averaged. For this analysis I used the type of irrigation (acequia, and center point) as my predictor or independent variable and the NDVI value (a continuous numerical value from 1, to -1) as the response or dependent variable.


Results of the Shapiro Wilkes tests and Fligner tests showed that all of the data was normally distributed with significant p value >0.05. In addition, Fligner tests showed that all of the variances in the data sets were homogeneous. I then conducted a two sample t test.




R OUT PUTS FOR NORMALITY TEST AND VARIANCE HOMOGENEITY, Shapiro.test(ace$pix1)Shapiro-Wilk normality test data: ace$pix1 W = 0.93979, p-value = 0.6644R

> Shapiro.test(ace$pix2) Shapiro-Wilk normality test data: ace$pix2 W = 0.91899, p-value = 0.5234

> Shapiro.test(ace$pix3) Shapiro-Wilk normality test data: ace$pix3 W = 0.90865, p-value = 0.4595

> Shapiro.test(cpi$pix1) Shapiro-Wilk normality test data: cpi$pix1 W = 0.98067, p-value = 0.9382

> Shapiro.test(cpi$pix2) Shapiro-Wilk normality test data: cpi$pix3 W = 0.81485, p-value = 0.1065> Shapiro.test(cpi$pix3) Shapiro-Wilk normality test data: cpi$pix3 W = 0.91112, p-value = 0.4744

> Fligner.test(pix1~Type)Fligner-Killeen test of homogeneity of variances data: pix1 by Type

Fligner-Killeen:med chi-squared = 0.94181, df = 1,

p-value = 0.3318

> fligner.test(pix2~Type) Fligner-Killeen test of homogeneity of variances data: pix2 by Type

Fligner-Killeen:med chi-squared = 0.36934, df = 1,

p-value = 0.5434

> fligner.test(pix3~Type)

Fligner-Killeen test of homogeneity of variances data: pix3 by Type Fligner-Killeen:med chi-squared = 0.50167, df = 1,

p-value = 0.4788

Two Sample t-test

data: ndvi $pix1 by ndvi $Type

t = -0.066394, df = 8, p-value = 0.9487

Null Hypothesis: true difference in means between group acequia and group center point is equal to 0

95 percent confidence interval:

-0.2559594 0.2416327

Map of an Agricultural Area in San Luis, CO. with NDVI valuesNAIP imagery of an Agricultural Area in the San Luis Valley, CONAIP (red band) Imagery of an Agricultural Area in San Luis ,CO

There was no difference in means

The results of the T Test show this with a 95% confidence interval for test area "pix 1"

I can not reject the null hypothesis.

References

Aragon, A. Remote Sensing Evaluation Techniques of Traditional Agriculture and Acequia Systems in the Upper Pecos River Valley, NM, US. Unpub. Thesis NMHU. 2018

Arellano, J.E. Enduring Acequias: Wisdom of the Land, Knowledge of the Water. UNM. 2014

Benavides, D. A Retrospective on NMAA's Policy Achievements. Noticias de las Acequias. Santa Fe. New Mexico Acequia Association. 2019

Clark, I. Water in New Mexico: A history of its management and use. Albuquerque: University of New Mexico Press.1987

Fernald, Baker, et.al. 2007. Hydrologic, Riparian and Agroecosystem Functions of Traditional Acequia Irrigation Systems. Journal of Sustainable Agriculture, Vol30(2). https://jsa.hawthornepress.com.2007. Accessed on 10/07/2021

Hall, G. E. Steve Reynolds -Portrait of a State Engineer as a Young Artist, 38 Nat. Resources J. 537, 538–561. 1998

Huang, Tang, et.al. Commentary Revue NDVI. For. Res32(1)[10]. 2021

R Core Team. R: A Language and Environment for Statistical Computing.R Foundation for Statistical Computing. Vienna, Austria. 2021. https://www.R-project.org.

Raheem, Gonzales et.al. Framework for Assessing Ecosystem Service in Acequia Irrigation. Wiley. 2015

Rodriguez, ,S; Rivera, J, et.al Key Concepts for Multi Dis. Approach to Acequias. NMSU.2021

Society for Restoration. International Primer on Restoration Methods. 2017

Scarborough, V. The Flow of Power. SAR. 2011

Tuwani-Smith. 2014. Decolonizing Research Methods. Routledge.

USDA. 2021.NAIP. https://www.fsa.usda.gov/programs-and-services/aerial-photography/imagery-programs/naip-imagery/index [9]

.

Superscript

Kristin Kinic

DBA: Urbana and the Art of Bricolage.

(505)-718-0193

kristinmkinic@gmail.com