The Heart River Basin is predominantly an agricultural basin in western North Dakota and is approximately 3,350 square miles. The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service and the Grant County Soil Conservation District, completed a study to assess spatial and temporal patterns of water quality in the Heart River Basin. The purpose of this report is to describe the methods and results of a study to evaluate salinity and nutrients in the Heart River Basin in western North Dakota. Water-quality and streamflow data used in the study were compiled from 1970 to 2020 using the National Water Quality Monitoring Council Water Quality Portal and National Water Information System.

Changes in streamflow characteristics were investigated at three sites from 1970 to 2020, and changes in water quality were investigated at four sites from 1974 to 2019. Streamflow analysis indicated decreasing streamflow from 1970 until the late 1990s followed by increasing streamflow through 2020, with the largest increase in the 7-day minimum streamflow or base flow. For the historical water-quality trend period (1974–2019), total dissolved solids, sulfate, sodium, chloride, and sodium adsorption ratio concentrations have increased since the mid-1970s through 2019. Potassium concentrations during the historical period remained mostly constant with some small fluctuations. Calcium and magnesium concentrations increased since the mid-1970s at all sites, except for a decrease at one site between 1974 and 1999. During the recent trend period (1999–2019), increasing concentrations in total dissolved solids, sulfate, sodium, chloride, calcium, magnesium, and sodium adsorption ratios were observed across the Heart River Basin. The magnitude of the increases was smaller at tributary sites compared to main-stem sites. During the recent period, potassium was mostly constant, although small (−0.9 milligram per liter or less) decreases on tributaries and minor (1.3 milligrams per liter) increases on the main-stem sites were detected. Unlike dissolved ion concentrations, significant increases in nutrient concentrations were not detected from 1999 to 2019, but nitrate plus nitrite concentrations most likely decreased upstream from Lake Tschida.

Inverse modeling for period 1 (1974–99) in model zone 1 (Heart River reach from site 5 to site 6) had eight reasonable models that indicated the clay mineral-water interactions and dissolution of evaporites control the geochemistry. Results of the inverse modeling for period 2 (1999–2019) in model zone 1 also had eight reasonable models that indicated that the dissolution of evaporites was the major geochemical control. Results of the geochemical modeling for period 1 (1974–99) in model zone 2 (Heart River and Sweetbriar Creek reach from sites 20 and 21 to site 22) produced seven reasonable models, and the geochemical control of the system was the dissolution of sulfate evaporite minerals. Geochemical modeling results for period 2 (1999–2019) in model zone 2 produced 11 reasonable models and was also controlled by the dissolution of sulfate evaporite minerals. Differences between the two model zones indicated that geology controls some of the water-quality changes in the Heart River Basin.

Loads were estimated for total dissolved solids, sulfate, sodium, and chloride and total phosphorus. Annual loads estimated for the Heart River from 2013 through 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) were generally greatest in 2014 and least in 2016 for total dissolved solids, sulfate, sodium, and chloride. Most of the annual loads of total dissolved solids, sulfate, sodium, and chloride are delivered in March through July in the Heart River at these sites and are likely from snowmelt and spring and summer rains. The mean annual yields of total dissolved solids and sodium from 2013 to 2020 generally were largest in Big Muddy Creek (site 18), whereas yields of sulfate and chloride were largest at Sweetbriar Creek (site 21) compared to the other selected sites in the Heart River Basin. Larger yields of total dissolved solids, sulfate, sodium, and chloride at sites located on Big Muddy Creek and Sweet Briar Creek in the lower Heart River Basin were likely a result of differences in geology and soils upstream from the selected sites.

A mass balance of total dissolved solids, sulfate, sodium, and chloride was estimated for the lower Heart River Basin, specifically the reach below Lake Tschida to Mandan (site 7 to site 22). Intervening flow was the largest contributor to the dissolved ion loads in the lower Heart River Basin and is an important part of understanding the transport of dissolved ions in the basin. The intervening load can include groundwater discharge, irrigation return flow, local runoff, and input from smaller ephemeral tributaries. Tributaries in the lower Heart River Basin contributed portions of the total dissolved solids, sulfate, sodium, and chloride loads at the Heart River near Mandan (site 22) that generally were proportional to the streamflow contributions.

Annual loads for total phosphorus between 2013 and 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) generally were largest in 2019 and smallest in 2016. Most of the total phosphorus loads for main-stem sites 5 and 22 were transported in March, April, and June, likely from snowmelt and early summer rains. The mean annual yields of total phosphorus for 2013–20 were largest on the main-stem site upstream from Lake Tschida (site 5) and Sweetbriar Creek (site 21), whereas the smallest yields were in Big Muddy Creek (site 18). Much of the phosphorus that enters Lake Tschida from the upper basin does not get transported downstream to the lower basin, and much of the phosphorus in the lower basin was attributed to intervening flow.