Comparing soil organic carbon, infiltration, and bulk density in various grass communities and management practices on Confederated Salish and Kootena

Progress report for GW24-014

Project Type: Graduate Student
Funds awarded in 2024: $30,000.00
Projected End Date: 03/31/2026
Grant Recipient: Salish Kootenai College
Region: Western
State: Montana
Graduate Student:
Principal Investigator:
Virgil Dupuis
Salish Kootenai College
Principal Investigator:
Maureen McCarthy
Desert Research Institute
Expand All

Project Information

Summary:

The Bison Range is located on the Flathead Indian Reservation in Western Montana, the remaining homelands of the Selis, Qlipse, and Ksanka people. Invasive macrophytes pose significant threat to highly valued cool season Palouse type bunchgrass grasslands dominated by bluebunch wheatgrass (Pseudoroegneria spicata), and rough fescues (Festuca campestris Rydb). Ventenata (Ventenata dubi), a recent annual invasive grass invader has become especially concerning, as it provides poor forage and creates a dense straw cover that smothers native vegetation (Wallace 2015, & Bradley, 2006). Ventenata was first noticed on the Bison Range around 2012, and is now throughout the 16,000-acre Bison Range.

Bison are native to Montana, yet they are relatively recent additions to the Bison Range’s inter-mountain ecosystem. Although they are known for being easier on the land than cattle, little is actually known on how they influence soil health on the Bison Range.

To more fully understand annual invasive grasses and bison grazing effects on soil health, soil organic carbon (SOC), infiltration rates, and bulk density needs to be studied. This study preposes to compare SOC, infiltration rates, and bulk density in two sites with three vegetation functional groups that will be mapped and sampled using the line-point intercept method. One site has a long history of bison grazing while the other nearby site is in an exclosure.

Invasive grasses and overgrazing can reduce organic carbon (Pendell et al., 2018, Koteen et al., 2011, & Harden, 1999). This study will help managers understand the relationship between soil carbon, vegetation, and grazing management on the Bison Range. Proxy measurements estimating total organic carbon, such as soil organic matter, bulk density, and infiltration rates will give insight into microsite variability (Franzluebbers, 2002).

Project Objectives:

Research Objective

  • The first objective is to sample vegetation and soils in three mapped grassland functional groups in one season. We will compare SOC, infiltration, and bulk density and determine which functional group is the most beneficial to soil health.
  • The second objective is to sample soil and vegetation in a bison grazed and ungrazed grassland in one growing season. We will compare the SOC, infiltration, and bulk density to determine whether Bison Range stocking rates are maintaining soil health.

 

Education Objective:

  1. After a field day with the producers and stakeholders, they
    will be able to discuss or describe the link between vegetation
    type and soil health as well as how management such as stock
    numbers and maintaining native grasses can influence soil
    health.  Land managers will be able to utilize data on the
    risks associated with annual invasive grasses to make informed
    management decisions in the following growing seasons.

Following a classroom presentation, SKC students will be able to
discuss the link between vegetation type and soil health as well
as how management such as stock numbers and maintaining native
grasses can influence soil health.

Cooperators

Click linked name(s) to expand/collapse or show everyone's info
  • Peter Bugoni - Producer
  • Virgil Dupis
  • Dr. Adrian Leighton (Researcher)

Research

Materials and methods:

The Bison Range is located in the Mission Valley of Montana and includes a low-rolling mountain, wetlands, a creek, and intermontane Palouse-type grasslands (Gilles & Coffman, 2019). This variety of ecosystems allows for a variety of vegetation and wildlife. Bluebunch wheatgrass and rough fescue are some of the important cool-season grasses on the Bison Range. (Garcia Neto, 2014). There are other grasses, small shrubs, and forb such as Koeleria macrantha (prairie junegrass), Symphoricarpos albus (snowberry), and Balsamorhiza sagittate (arrowleaf balsamroot). (Gilles & Coffman, 2019). Invasive species on the Bison Range include Cirsium arvense (Canada thistle), Bromus tectorum (cheatgrass), Potentilla recta (Sulphur cinquefoil), Dipsacus fullonum (teasel), Ventenata dubiaLepidium draba (whitetop), (Gilles & Coffman, 2019) The main ungulates on the Bison Range putting grazing pressure on the grasslands are bison, elk, white-tailed deer, mule deer, bighorn sheep, and pronghorn, (Garcia Neto, 2014). 

The intermountain valley where the Bison Range is located was created by glacier action and was once a large glacial lake. Soil Valley sediments, lakebed silt, and glacier debris are the main soils in the valley (Smith, 2001). The north side of the mountain tends to be moister than the south side (Gilles & Coffman, 2019). Ventenata can be found in all aspects but it tends to dominate southern aspects. 

The majority of precipitation is in the spring and early summer averaging about 2.5 inches per month in May and June (Gilles & Coffman Gilles & Coffman, 2019). Climate projections for the next several decades predict a rise in temperature of about 2 to 4 degrees globally with a local, 2019).

The grassland environment on the CSKT Bison Range is currently subjected to numerous environmental and climatic variabilities. Since the Bison Range was returned to the Confederated Salish, Kootenai, and Pend d’Oreilles Tribes, the tribes have decided to let the bison choose where to graze within the Bison Range (Peter Bugoni, Personal Communication, 2023). Future changes from climate change are projected to include an increase in annual invasive grasses and loss of cool season grass production (Abatzoglou and Kolden, 2011, & Loka et al., 2018). 

Site selection

The study sites have been selected on the northeastern side of the Bison Range, at the base of a small rolling mountain. The desirable aspects of the site are the presence of ventenata dubia and native grasses, no recent herbicide applications, and no western cedar nearby. The site is relatively flat with few if any gullies and small hills. The soil type is primarily the Ronan series. (USDA, 1978). Each site will be one of two main management groups; bison grazed and a bison exclosure. According to Bison Range staff, the exclosed area has been fenced off from the rest of the Bison Range for the past 50 years. Deer, elk, and possibly pronghorns have been the only ungulate grazers in the fenced area. These sites are in triangles fenced off at the base of Red Sheep Mountain on the Northeast corner of the Bison Range. The grazed area is just northwest of the enclosers and has been grazed by Bison since the Bison Range was established in 1908 (CSKT, Bison Range). 

 

Step 1 Mapping Focal Groups

  The grassland communities were stratified in 2023. The researcher walked through the sites to visually assess patches of vegetation and designate them into general vegetation focal groups. Ventenata dubia, an annual invasive grass prevalent on the Bison Range, will be one focal group (Gilles & Coffman, 2019). Other vegetation focal groups included native grasses such as bluebunch wheatgrass (Pseudoroegneria spicata), rough fescue (Festuca campestris Rydb), prairie junegrass (Koeleria macrantha), western wheatgrass (Pascopyrum smithii), (Gilles & Coffman, 2019). Researchers will use GIS technology to map focal groups. 

Step 2 Vegetation Sampling 

Researchers selected three of the most prominent focal groups; bunch grasses, rhizomatous grasses, and annual invasive grasses. Three plots out of each focal group made up a study site. Researchers estimated vegetation cover and diversity within each plot. They also took five soil samples, and five bulk density samples, and perform three soil infiltration tests from each plot. 

 In June and part of July, vegetation percent cover was estimated using a 15 m line-point intercept (LPI)(De Stephano et al. 2021, and Rangeland Gateways, 2023). A 15 m radius circular plot was placed in the center of the focal groups and split into quadrants. The Researcher walked back and forth, documenting plants not identified with the LPI, such as trees and shrubs, and plants not picked up with the LPI (Sherman and Anderson, 2023 p. 61).

 

Step 3 Soil Sampling

In August and September, the researchers took infiltration rates and soil samples from five random points selected in each plot, using ArcGIS Pro and Google random number generator. The closest microsite (a small patch of vegetation of interest) to the random point was used for sampling. Sarah Holloway and sometimes Jarod Swan or Brad Mitchell took infiltration rates first using the Mini-Disk Infiltrometer. Infiltration rates took from 12 to 80 minutes. The researchers dug micro pits for soil sampling within the same microsite while avoiding the moisture from the infiltration test. Each pit was about 30cmX30cm. The researchers hammered a 4cmX4cm cylinder into the side of the pit at 1 cm and 15 cm below the surface then used a Hori Hori knife to remove the cylinder. The researchers collected small chunks of soil at 30 cm to keep the hole small rather than using the cylinder. They stored soil samples in labeled Ziplock bags. After returning to the office, each soil sample collected that day was weighed and dried in a dehydrator at 160 F for 12 hours. The following day, the soil samples were weighed and re-bagged.

In January 2025, Holloway took the soil to the Desert Research Institute in Reno, Nevada. Holloway worked in the soils lab, with the assistance of Sally Houseman, Kristine Lu, and Alan Villanueva, to test for soil organic carbon (SOC) using the loss on ignition (LOI) method. Holloway split the samples and placed them in a 105-degree Celsius oven the day before LOI testing. The samples remained in the warm oven until each batch was tested. Small batches were weighed and put in a muffle furnace at 400 degrees Celsius for an hour. Holloway cooled the samples for 20 min in a desiccator before weighing them. After the LOI testing, Holloway split several more samples for a laser particle size analysis. She took one surface sample from each plot. Houseman ran the laser particle size analysis and emailed the results.

Data collection

   Soil and vegetation data collection forms were signed, time and date recorded, and original copies stored at the SKC Extension Office.

Schedule of activities

 

  • Fall 2023 -Stratify three grass communities within triangle areas.
  • Spring and summer 2024- Sample Sites
  • Fall 2024- Analyze Data
  • Winter 2025-LOI testing
  • Spring and summer 2025- Data analysis and write up

 

Research Objectives

  • The first objective is to sample vegetation and soils in three mapped grassland functional groups in one season. We will compare SOC, infiltration, and bulk density and determine which functional group is the most beneficial to soil health.
  • The second objective is to sample soil and vegetation in a bison grazed and un-grazed grassland in one growing season. We will compare the SOC, infiltration, and bulk density to determine whether Bison Range stocking rates are maintaining soil health.

Data Management Plan

Data collection

        The collected soil was moved directly from the cylinder into a new bag. Soil was stored in a secure draw in a filing cabinet or in a locked cabinet. Data collected was signed or initialed, timed, and dated. The team was trained to properly identify the common grasses.        

       To protect the data, soil data was written on waterproof paper and researchers used pencil or a waterproof pen. Pictures were taken of data sheets for backup. When the research team returned to the office, they entered the data onto data sheets onto a computer that has a good security system in place including passwords and a firewall. Data sheets are backed up on Google Drive. More than one team member has access to the data so it remains accessible in case someone leaves the team unexpectedly.

To maintain quality control consistent methods were be used at each site. Researchers will measure the completeness if any data is found to be invalid (Holmes et al., 2021). They will ensure it is comparable and is a representative data set by documenting variabilities (Holmes et al., 2021).

Research results and discussion:

The line-point intercept (LPI) data showed that ventenata was the most common plant to get hit with the rod. But even in the plots with the most ventenata, Western wheatgrass and Kentucky bluegrass still grew. The LPI data shows that bunchgrasses supported the most diversity and native forbs. The bunchgrass plots also had the least cover of the desired plants like.

Bunchgrasses Average

   

Rhizomatous grasses Average

   

Annual Invasive Grasses Average

Vegetation

Total

   

Vegetation

Total

   

Vegetation

Total

Ventenata

9.33

   

Ventenata

25.33

   

Ventenata

42

Bluebunch wheatgrass

2.835

   

Bluebunch wheatgrass

1.335

   

Bluebunch wheatgrass

0

Western wheatgrass

5.5

   

Western wheatgrass

19.5

   

Western wheatgrass

9.665

Rough Fescue

19.165

   

Rough Fescue

0

   

Rough Fescue

0

Cheatgrass

0.835

   

Cheatgrass

0.665

   

Cheatgrass

0.665

Canada bluegrass

3.665

   

Canada bluegrass

0.83

   

Canada bluegrass

1

Kentucky Bluegrass

0.835

   

Kentucky Bluegrass

9

   

Kentucky Bluegrass

7.335

Canadian wild rye

0

   

Canadian wild rye

0.165

   

Canadian wild rye

0

Junegrass

0

   

Junegrass

0

   

Junegrass

0

Crested wheatgrass

0.835

   

Crested wheatgrass

0

   

Crested wheatgrass

0.165

Japanese brome

0.166667

   

Japanese brome

2.83

   

Japanese brome

4.83

 

 

   

 

 

   

 

 

Native Forb

4.833334

   

Native Forb

0.835

   

Native Forb

1.335

Non-native Forb

5.665

   

Non-native Forb

2.335

   

Non-native Forb

4.83

Analysis found higher infiltration rates and higher bulk density in Site 2. Infiltration rates were highest in bunchgrasses and the lowest in invasive annuals. Bulk density was higher in rhizomatous grasses than bunch grasses and lowest in annual grasses (Table 1). 

PLOT

Infiltration Rate
MEAN Ml PER MIN

Standard
Deviation

% Soil Moisture

Mean
Bulk Density

Bluebunch1_1

4.391

0.532

9.42

1.09

Bluebunch1_2

1.929

0.597

13.78

1.22

Bluebunch1_3

0.573

0.129

18.43

1.21

Ventenata1_1

1.417

0.507

15.53

1.07

Ventenata1_2

1.162

0.197

20

1.17

Ventenata1_3

1.247

0.530

18.12

0.96

Western wheatgrass1_1

2.201

1.356

12.35

1.1

Western wheatgrass1_2

1.231

0.533

10.46

1.31

Western wheatgrass1_3

2.257

1.411

14.02

1.27

Rough fescue2_1

5.063

0.588

6.05

0.96

Rough fescue2_2

2.503

0.920

10.18

1.34

Rough fescue2_3

1.962

0.548

8.88

1.38

Ventenata2_1

1.881

1.366

7.3

1.35

Ventenata2_2

5.384

1.478

4.8

0.98

Ventenata2_3

2.628

1.516

9.19

0.99

Western wheatgrass2_1

3.152

1.927

18.94

1.21

Western wheatgrass2_2

2.753

0.965

12.56

1.56

Western wheatgrass2_3

4.022

1.001

16.93

1.33

Table 1. Infiltration rates and bulk density by plot

Houseman ran the laser particle size analysis which showed that Site 1 was primarily Sandy Loam, and Site 2 was primarily sandy, though both sites had a mixture of both soil types. Initial LOI data showed that the grazed site (Site 1) had more SOC than Site 2, and that the Invasive annuals had slightly higher SOC in the surface soil and at 15 cm but that the rhizomatous grass had higher SOC at 30 cm than the other sites at the same depth (Table 2).

Bunchgrasses
Bluebunch wheatgrass Site 1
PLOTS 1 2 3 AVERAGE
TOP 2.74 2.68 2.56 2.66
15 CM 1.34 1.47 1.26 1.36
30 CM 1.07 0.9 1.29 1.09
Rough fescue Site 2
PLOTS 1 2 3 AVERAGE
TOP 4.02 2.43 2.3 2.92
15 CM 1.61 1.3 1.07 1.33
30 CM 1.33 0.84 0.55 0.91
Rhizomatous
Western Wheatgrass Site 1
PLOTS 1 2 3 AVERAGE
TOP 2.99 2.23 2.39 2.53
15 CM 2 1.91 1.61 1.84
30 CM 1.41 1.53 1.31 1.42
Western Wheatgrass Site 2
PLOTS 1 2 3 AVERAGE
TOP 2.12 1.66 1.68 1.82
15 CM 1.95 1.2 1.08 1.41
30 CM 1.63 0.97 0.82 1.14
Annual Invasive grass
Ventenata Site 1
PLOTS 1 2 3 AVERAGE
TOP 3 2.58 3.07 2.88
15 CM 1.75 1.5 2.03 1.76
30 CM 1.44 1.04 1.49 1.32
Ventenata Site 2
PLOTS 1 2 3 AVERAGE
TOP 2.38 3.07 2.89 2.78
15 CM 1.52 1.69 1.24 1.48
30 CM 1.02 1.12 0.95 1.03

Table 2. Average soil organic carbon in sites and plots on various depths. 

Further analysis is needed to make comparisons.

 

Participation Summary
1 Producers participating in research

Research Outcomes

Recommendations for sustainable agricultural production and future research:

This research provides baseline data on vegetation cover, soil organic carbon (SOC), infiltration rates, and bulk density that may help in restoration or rehabilitation efforts, potentially improving ecosystem services and increasing forage. Some vegetation analysis showed Western wheatgrass and Kentucky bluegrass were in heavily invested ventenata patches. As both species are rhizomatous grasses, further research could look into using rhizomatous grasses to restore rangelands heavily infested with annual invasive grasses.

Understanding how invasive species and the grazing regime change the soil helps inform management decisions regarding native species, which will benefit the tribes and the public viewing the Bison Range. The data on SOC could inform land managers about management practices and their impact on SOC. Managers could purposely improve soil carbon storage by understanding grazing impacts and plant composition on the soil.  Also, understanding grazing impacts can help inform stocking rates and other grazing management decisions. Further research on stocking rates on cool season grasses and their impact on soil health could be beneficial to rangeland management.

1 Grant received that built upon this project

Education and Outreach

1 Curricula, factsheets or educational tools
1 Tours
1 Webinars / talks / presentations

Participation Summary:

8 Ag professionals participated
Education and outreach methods and analyses:

There have been two opportunities to educate local scientists and the broader community about the research on soil health and grass functional groups on the Bison Range.

While collecting soil samples on the bison range, scientists were touring the Bison Range with Peter Bugoni. They stopped by one of the plots, and Holloway had the opportunity to share information about the research objectives and methods. The conversation with one of the soil scientists led to the decision to collect soil at 30 cm instead of just at 1 cm and 15 cm.

In early April members of Native FEWs Alliance, a group concerned with food, energy, and water availability on tribal lands, visited Salish Kootenai College as part of an annual gathering. In addition to various college faculty and staff presenting the college's involvement in pertinent areas, the SKC Extension Office hosted a poster walk. Holloway presented a poster on microsite variability, comparing soil carbon and other soil health indicators in various grass communities. The poster presentation opened up the opportunity to discuss some of the research findings and the intricacies of dealing with invasive species and carbon storage.

Summer 2025- Tour with scientists

Spring 2025- Poster presentation

Summer/Fall 2025 - Field Trip with producers and stakeholders

Fall 2025 - Classroom presentation

 

Education and outreach results:

Both educational opportunities lead to conversations about how native and non-native species store carbon in the soil. The discussion with a young member of the Native FEWS Alliance, in particular, led to some thought-provoking questions on how we balance the preservation of native species and the desire to store more carbon.

Information Products

    Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and should not be construed to represent any official USDA or U.S. Government determination or policy.