Impacts of Quaternary History on Critical Zone Structure and Processes: Examples and a Conceptual Model From the Intensively Managed Landscapes Critical Zone Observatory

The concept of a critical zone (CZ) supporting terrestrial life has fostered groundbreaking interdisciplinary science addressing complex interactions among water, soil, rock, air and life near Earth’s surface. Pioneering work has focused on the CZ in areas with residual soils and steady-state or erosional topography. CZ evolution in these areas is conceptualized as progressive weathering of local bedrock (e.g. in the flow-through reactor model). However, this model is not applicable to areas in which weathering profiles form in transported materials including the formerly glaciated portion of the Central Lowland of North America. We present a new conceptual model of CZ evolution in landscapes impacted by continental glaciation based on investigations at three study sites in the Intensively Managed Landscapes Critical Zone Observatory (IML-CZO) The IML-CZO is devoted to the study of CZ processes in a region characterized by thick surficial deposits resulting from multiple continental glaciations, with bedrock at depths of up to 150 m. Here the physical (glacial ice, loess, developing soil profiles) and biological (microbes, tundra, forest, prairie) components of the CZ vary significantly in time. Moreover, the spatial relationships between mineral components of the CZ record a history of glacial-interglacial cycles and landscape evolution. We present cross-sections from IML-CZO sites to provide specific examples of how environmental change is recorded by the structure of the mineral components of the CZ. We build on these examples to create an idealized model of CZ evolution through a glacial cycle that represents the IML-CZO sites and other areas of low relief that have experienced continental glaciation. In addition, we identify two main characteristics of CZ structure which should be included in a conceptual model of CZ development in the IML-CZO and similar settings: (1) mineral components have diverse origins and transport trajectories including alteration in past CZs, and, (2) variability in climate, ecosystems, and hydrology during glacial-interglacial cycles profoundly influence the CZ composition, creating a legacy retained in its structure. This legacy is important because the current physical CZ structure influences the occurrence and rates of CZ processes, as well as future CZ responses to land use and climate change.

The concept of a critical zone (CZ) supporting terrestrial life has fostered groundbreaking interdisciplinary science addressing complex interactions among water, soil, rock, air, and life near Earth's surface. Pioneering work has focused on the CZ in areas with residual soils and steady-state or erosional topography. CZ evolution in these areas is conceptualized as progressive weathering of local bedrock (e.g., in the flow-through reactor model). However, this model is not applicable to areas in which weathering profiles form in transported materials including the formerly glaciated portion of the Central Lowland of North America. We present a new conceptual model of CZ evolution in landscapes impacted by continental glaciation based on investigations at three study sites in the Intensively Managed Landscapes Critical Zone Observatory (IML-CZO) The IML-CZO is devoted to the study of CZ processes in a region characterized by thick surficial deposits resulting from multiple continental glaciations, with bedrock at depths of up to 150 m. Here the physical (glacial ice, loess, developing soil profiles) and biological (microbes, tundra, forest, prairie) components of the CZ vary significantly in time. Moreover, the spatial relationships between mineral components of the CZ record a history of glacial-interglacial cycles and landscape evolution. We present cross-sections from IML-CZO sites to provide specific examples of how environmental change is recorded by the structure of the mineral components of the CZ. We build on these examples to create an idealized model of CZ evolution through a glacial cycle that represents the IML-CZO sites and other areas of low relief that have experienced continental glaciation. In addition, we identify two main characteristics of CZ structure which should be included in a conceptual model of CZ development in the IML-CZO and similar settings: (1) mineral components have diverse origins and transport trajectories

INTRODUCTION
The formerly glaciated portions of the Central Lowland of the USA (Fenneman and Johnson, 1946) are among the most intensively managed landscapes and agriculturally productive croplands in the world. The combination of deep fertile soils, temperate climate, and concentrated human interactions to improve growing conditions enables the highest rates of photosynthesis on the Earth (Guanter et al., 2014). Interactions between the atmosphere, hydrosphere, biosphere and lithosphere provide ecosystem services not only in pristine landscapes, but also in this human modified environment. Pioneering research into these interrelated systems at the Earth's surface recognizes a "critical zone" (CZ) that regulates habitats and determines the availability of life-sustaining resources. The CZ was originally defined as stretching from the top of the vegetation canopy to the lower limits of groundwater (e.g., NRC, 2001;Brantley et al., 2007). Yet emerging research suggests the lower CZ limit is better defined as the depth at which meteoric water is in chemical equilibrium with minerals (Riebe et al., 2017). The CZ thus includes the soil, rock, sediment, water, air, organic matter and organisms through the full depth of actively forming soils and extends above to the top of vegetation and below into materials hosting flowing, reactive groundwater. Study of CZ processes requires knowledge of the composition and physical arrangement (structure) of the surficial geologic materials. Furthermore, the structure of the CZ depends on a history stretching back well before the modern anthropogenic land-use changes. Therefore, both ongoing CZ processes and the response of the CZ to future external forcing remain sensitive to the longer history that has shaped the CZ. The Midwest US contains a CZ where current processes are strongly impacted by profound anthropogenic changes, but are ultimately controlled by the events of past global glacial-interglacial climate cycles.
The Intensively Managed Landscapes Critical Zone Observatory (IML-CZO) with field laboratories at Clear Creek Watershed (CCW) in eastern Iowa, Upper Sangamon River Basin (USRB) in central Illinois, and the Minnesota River Valley (MRV) in southern Minnesota is dedicated to studying the CZ processes of land under row crop agriculture. We focus on the interactions in the CZ that control the rates and trajectories of processes responsible for transforming and transporting water, nutrients, sediment, and energy over timescales ranging from minutes to a few million years. Continuous deposits of glacial and glaciofluvial sediment comprise ∼20% of the land surface of the continental United States, including the glaciated portions of the Central Lowland of the Interior Plains (Fenneman and Johnson, 1946;Soller et al., 2009; Figure 1). The IML-CZO is characterized by successive periods of rapid sediment aggradation and erosion associated with continental glaciations punctuated by long (tens to hundreds of thousands of years) interglacial or post-glacial times of surficial weathering and soil development, a fundamentally different landsystem than that of the steady state or erosional topographic conditions of the other eight observatories in the NSF-supported national CZO network (Figure 1). A fundamental challenge for researchers working in the IML-CZO field laboratories has been to understand how the CZ evolved through multiple glacial/interglacial cycles in this region of overall net deposition, and how these landscapes are controlled by past geologic processes. This inheritance of past geologic events in glaciated landscapes of the mid-continent is known as the "glacial legacy, " which still today influences the CZ processes and trajectory of landscape and CZ evolution under intensively managed land-use practices (Phillips and Desloges, 2015).
In an effort to understand the CZ of this region, we confront the fact that the primary CZ conceptual model-progressive weathering of local bedrock moving through a weathering profile in areas of net surface erosion or steady-state topography-, the feed-through reactor of Anderson et al. (2007), has serious limitations in the IML-CZO where thick, continuous glacial sediments blanket the landscape, and in other areas where alluvial, lacustrine, coastal, or eolian sediments predominate. In the feed-through reactor model, as the weathering front moves downward through the subsurface, rock fragments move vertically upward undergoing a progression of chemical and physical weathering reactions. As a result, the CZ weathering record is preserved locally within the soil and weathering profiles. Applying this model across landscapes and regions allows the geological principle of substituting space for time to model the evolution of the CZ in different lithologies under varying environmental conditions. However, in regions with thick glacial deposits, the physical and geochemical properties of geologic materials in the CZ often have little to no affinity with the local bedrock. Here the materials did not form in situ from the local bedrock, but instead are a mix of far-traveled and local sediment, episodically deposited, that underwent varying degrees of weathering prior to and along their transport pathways. Additionally, the significant differences between glacial and interglacial climates, and their associated geomorphological, hydrological, and ecological systems, provide a temporal variability in CZ processes that is of unmistakable FIGURE 1 | Generalized map showing surficial materials (Source of data: U.S Geological Survey, Soller et al., 2009) and the locations of CZO field observatories within the conterminous United States. The glaciated portion of the Central Lowland of the Interior Plains (Source of data: U.S Geological Survey, Fenneman and Johnson, 1946) is indicated by the narrow black line. The US CZO network is indicated by the black stars with the exception of Loquillo, Puerto Rico, which is not within the map area. The bold black oval surrounds the three IML-CZO field sites. Colored areas denote continuous cover of unconsolidated sediment. Green "Glacial" areas are underlain by glacial till, glaciofluvial, glaciolacustrine, and ice contact deposits. Brown "Eolian" are areas where loess or dune sand is at least 20 feet (6 m) thick at the surface. Blue "Coastal" areas contain both fine-grained sediments characteristic of lagoons and estuaries and coarser beach or shoreline deposits. The red "Lacustrine" areas are underlain by lake and playa deposits. Salmon-colored "Alluvial" areas are underlain by alluvial and colluvial deposits along major river valleys and inland basins. The dark gray areas "Organic/Biological" include recent organic sediments (peat) and biological sediments (marl and shelly materials). The white areas include bedrock or residuum developed from weathering of the underlying bedrock.
importance. Our detailed examination motivates us to develop a new conceptual model of CZ evolution in the IML-CZO and the glaciated portions of the Central Lowlands of the Interior Plains that encapsulates two important elements. First, the geologic materials may have diverse origins and distinctive properties acquired during transport and subsequent deposition that do not reflect the composition of the underlying bedrock. Second, the history of change, whether it be in climate, geomorphic processes, or by human-induced modifications is significant and crucial for modeling critical zone structure and processes.
We begin by presenting a series of generalized cross-sections for the three field laboratories of the IML-CZO showing the heterogeneity and distribution of geologic materials in the glaciated central USA. Note that these cross-sections show specific examples of geomorphic settings typical of this CZO and of the Central Lowlands of the Interior Plain more broadly. These cross-sections motivate an idealization of the evolving structure of the CZ over the course of a single glacial cycle. This evolution highlights the key elements that need representation in a conceptual model of the CZ for landscapes characterized by net deposition of glacial sediments (e.g., Kołodynska-Gawrysiak et al., 2017), specifically, the diverse histories of the geologic materials and the dynamic evolution of the CZ through time. Finally we discuss our conceptual model for the CZ of the Central Lowlands of the Interior Plains and assess the applicability of the key elements of this model to CZs more broadly.

GLACIAL HISTORY AND CZ STRUCTURE OF IML-CZO FIELD AREAS
In order to consider the glacial legacy in the IML-CZO, it is necessary to consider the important differences in the history of Quaternary glacial-interglacial cycles at each field laboratory. For example, glaciers last covered CCW almost 500,000 years ago whereas the most recent advance of the Laurentide Ice Sheet in the MRV ended ∼15,000 years ago. The three IML-CZO field observatories encompass the range of glacial histories in the Central Lowlands (Figure 2). During the past 2.6 million years, cycles of physical and chemical weathering, soil formation, landscape development, and ecosystem change in the CZ were by driven by successive glacial advances that reshaped the landscape, altered the surface and subsurface hydrology, and buried or excavated all or part of ancient CZs while depositing "fresh" glacial deposits. The record of landscape development in the stratigraphic record (Willman and Frye, 1970;Hansel and Johnson, 1996;, established the initial geological framework for understanding the subsequent CZ evolution. We illustrate the composition and distribution of Quaternary successions in the IML-CZO using a series of cross sections at a range of scales from the three field laboratories. The major critical zone events in the three areas are summarized in Table 1. The CCW was last glaciated during the early Quaternary (pre-Illinois Episode) and the USRB and MRV were last glaciated during the late Quaternary (Wisconsin Episode). The glacial successions are dominated by stacked, fine-grained till deposited during multiple glaciations. These very poorly sorted sediments (matrix-supported fine-grained diamictons) are intercalated or interspersed with deposits of well-sorted, coarseto fine-grained glaciofluvial and glaciolacustrine sediments that often infill valleys, channels and low areas on the bedrock surface or on surfaces incised in older glacial deposits. Variably thick loess mantles the morainal uplands and outwash plains across landscapes both within and beyond late Wisconsin Episode ice margins. The physical, chemical, and hydrologic properties of the glacial sediments vary due to differences in lithology at the source area, subsequent reworking during transport, and postdepositional alterations that have controlled types and rates of CZ change.

Clear Creek Watershed, Iowa
Glacial landforms from the pre-Illinois Episode glaciation are not preserved, although there are scattered, loess-mantled tabular divides underlain by thick, clayey paleosols, which are likely remains of the latest pre-Illinois Episode "drift" plain (Ruhe, 1969). The bedrock surface in the CCW is crosscut by a palimpsest set of channel segments, up to 60 m deep, which formed during the Tertiary and Early Quaternary. Thick, continuous deposits of sand and gravel do not infill the bedrock channels. The glacial succession generally thins from about 100 to 125m in the headwaters area to 15-30m in the lower part of the watershed where Paleozoic carbonates are often exposed at the land surface. Locally, the thickest glacial successions are associated with the bedrock channels.
The sequence of pre-Illinois Episode (∼2.6-0.5 Ma) glacial deposits includes seven different till sheets, each with an average thickness of 15-20 m, represented by two formations that are differentiated by clay mineral content (Hallberg, 1986;Bettis et al., 2010) (Figure 3). The two lowermost tills, both in the Alburnett Formation, have reversed magnetic polarity (older than 0.78 Ma) and are enriched in kaolinite in the clay fraction. Alburnett Formation tills contain significant amounts of eroded material derived from Tertiary weathering profiles (including kaolinite) eroded from sources north and northwest of CCW during the earliest advances of the Laurentide Ice Sheet (Roy et al., 2004). In addition, the oldest till forms a discontinuous layer of "stony" material on the bedrock surface and contains and abundance of local bedrock clasts (Figure 3). Tills of the overlying Wolf Creek Formation have normal polarity and are enriched in expandable clay minerals. Other than the lowest till, all the tills are similar in their physical properties and have a shared northern to northwestern provenance (Roy et al., 2007). In some places, tills in both formations are buried by deposits of sand and gravel (outwash), or have been completely eroded by meltwater. Following each glacial advance, CZ evolution in the tills produced color changes (primarily reflecting changes in the oxidation state of iron), leaching of primary carbonate minerals, mineral weathering, fracturing, and pedogenicallydeveloped concretions and cementation. During the next glacial advance, erosion by ice or water usually removed all or part of the upland CZ while burying and preserving the CZ formed in valley positions.
The modern landscape the CCW formed as the drainage network expanded and slopes were reshaped during the last two glacial-interglacial cycles occurring over the past half-million years. Although the CCW remained unglaciated during this time period, periglacial processes and distal glaciofluvial processes significantly affected development of the (1) regional drainage system, (2) local topography, and (3) CZ processes. As surface drainage networks, hillslopes, and the shallow groundwater system developed, oxygenated water and organic matter were delivered to the CZ and translocated downward into the unaltered glacial till and deeper part of the glacial succession.
Drainage divides, interfluves, drainage lines, and most hillslope elements on the modern landscape, as well as wellexpressed soils and associated deep weathering profiles were in place at the close of the last interglacial (Sangamon Episode). The most recent (late Wisconsin Episode) glaciation was unique in the US Midwest because of the unprecedented magnitude and areal dispersal of loess Mason et al., 2007). Beyond the late-glacial ice margin, loess composed of carbonaterich silt-loam accumulated in the developing boreal forest and parkland between 22 and 12 ka, which buried the glaciated landscape and associated CZ. In the CCW 3-7 m of loess mantles the uplands and old glacial terraces, completely burying the ancient CZ and its associated clayey paleosols and weathering profiles (Figure 3). This burial altered the chemistry of the preexisting weathering profile and soils through organic matter degradation, additions of calcite and other easily weathered minerals and redox-related changes in iron species. The modern CZ in CCW extends through former CZs that strongly influence the structure and processes operating in the current CZ. For example, the clayey paleosol of the pre-loess CZ acts as an aquitard that "perches" the water table on uplands and focuses Deep weathering profiles and well-integrated drainages developed into Paleozoic and older bedrock.
Pre-Illinois episode glaciation Tills (and outwash) deposited. Erosion of pre-glacial landscape and incision into bedrock. Later deposition masks relief. Pre-Yarmouth soil profiles formed (eroded on uplands, locally preserved in valleys. Deposition and erosion of tills (and outwash) from several glacial advances; bedrock valleys partially filled with outwash and till.
Weakly-developed soil formed on upland and floodplains surfaces (partially preserved).
Tills (and outwash) deposited; varying ice-flow pathways produce distinct compositions. Incision of fluvial drainage occurred during interglacials. Weathering profiles developed in glacial sediments during interglacials are largely not preserved.

Yarmouth episode interglacial
Deep weathering profile and well-integrated drainages developed in glacial sediment.
Deep weathering profile (only locally preserved) and well-integrated drainages developed in glacial sediment.
Weathering during this period is locally contained in lake deposits. Weathering profile rarely preserved in glacial sediments.

Illinois episode glaciation
Likely periglacial reworking of landscape and permafrost during Illinois Episode. Minor addition of loess?
Till, outwash, and glacial lake sediment deposited during multiple glacial advance and retreat cycles.
Possible deposition of tills and outwash.

Sangamon episode interglacial
Elements of modern landscape, including most drainage channels in place. Deep weathering profile continues to develop.
Interglacial weathering profile development (partially preserved) and fluvial network growth.
Interglacial soil development (rarely preserved) in glacial sediments.

Wisconsin episode glaciation
Weathering profile and landscape development continues through middle period. Loess incorporated before ∼27 ka. Periglacial resedimentation before 25 ka. Burial of weathering profile by 3-7 m of loess in a tundra to boreal environment.
Peat and loess deposited prior to glaciation. Partial erosion of Sangamon weathering profiles. Tills (outwash) deposited forming several morainal systems. Periglacial reworking of land surface; later buried by 1-2 m loess in a mixed tundra-boreal environment.
Erosion of older tills and weathering profile by glaciers. Till deposited; uppermost part of unit includes supraglacial material. Sedimentation in Glacial Lake Agassiz and catastrophic drainage as outlet opens. Meltwater carves out Minnesota River (River Warren) through older sediment and bedrock.

Holocene (modern interglacial)
New weathering profile formed into last glacial loess and superimposed on older pre-Illinois Episode weathered sediment in landscape transitioning from hardwood forest to tall grass prairie/savannah.
Weathering of loess (and uppermost till) occurs through a succession of vegetation from forest to prairie. Expansion of fluvial drainage networks.
Knickpoint migration up tributaries. Gully and bluff erosion widespread. Development of modern weathering profile in uppermost glacial till under climatic conditions supporting succession from mixed coniferous/hardwood forest to tall grass prairie.

Human management of landscape
Land use change to row crop agriculture. Accelerated upland soil erosion and accumulation of post-settlement alluvium, burying Holocene soils on floodplains. Humans engineer stream straightening and implement subsurface tile drainage.
Land use change to row crop agriculture. Shallow lakes and wetlands drained. Accelerated upland soil erosion and post-settlement floodplain deposition. Humans engineer stream straightening and implement subsurface tile drainage.
Land use change to row crop agriculture. Shallow lakes and wetlands drained. Intensification of gullying processes. Humans engineer stream straightening and implement subsurface tile drainage.
subsurface flow laterally as the overlying loess thins downslope. The subsurface distribution of the paleosol and its outcrop pattern on the landscape strongly control shallow groundwater flow and the distribution of seasonally saturated intervals in the shallow subsurface that impact weathering and nutrient loss.

Upper Sangamon River Basin, Illinois
The Upper Sangamon River Basin (USRB) was glaciated multiple times during the Pleistocene, including during the most recent (Wisconsin Episode) glaciation (Figure 2). The thickness of unconsolidated deposits, above Paleozoic bedrock, ranges from about 50 to 150 m (Soller et al., 1999;Stumpf and Dey, 2012). Much of the USRB overlies the buried preglacial Mahomet Bedrock Valley (MBV), a prominent feature on the bedrock surface (Horberg, 1945;Kempton et al., 1991). The MBV is part of the Teays-Mahomet Bedrock Valley system that stretches from West Virginia to Illinois (Melhorn and Klempton, 1991) and contained an active river system during the Pliocene and early Pleistocene. Underlying much of the central and southern USRB, the buried MBV ranges from 8 to 32 km wide and is filled with up to 50 m of sediment (Figures 4, 5A), mainly sand to gravelly sand. The coarse-grained fill in this bedrock valley, deposited by glacial meltwaters during the pre-Illinois and Illinois Episodes (Figure 5), contains the regionally important Mahomet Aquifer (Herzog et al., 1995;Roadcap et al., 2011). During the pre-Illinois Episode in the USRB (∼1.2-0.4 Ma), glaciofluvial sediment (or outwash) was repeatedly deposited and excavated from the MBV. Such outwash mainly filled the preglacial bedrock valley and is overlain by till sequences with some glaciolacustrine sediment ( Figure 5C; Kempton et al., 1991;Stumpf and Atkinson, 2015).  Three or more pre-Illinois Episode till sequences, together up to 30 m thick, occur in areas outside the MBV or overlie the pre-Illinois Episode glaciofluvial deposits (Kempton et al., 1991;Stumpf and Atkinson, 2015). Following the pre-Illinois Episode glaciations, a long period of weathering, soil formation and landscape denudation occurred during the Yarmouth interglacial episode (∼0.4-0.2 Ma) (Grimley et al., 2003;. A paleosol is found in some areas of the USRB, but was typically eroded by the initial glacial advance of the Illinois Episode. During the late Illinois Episode (∼0.16-0.13 Ma), glacial ice advanced out of the Lake Michigan Basin across the USRB to southernmost Illinois and westward to eastern Iowa (Figure 2). In Illinois, multiple sequences of glacial deposition, related to ice margin fluctuations, are found statewide (Willman and Frye, 1970;Hansel and McKay, 2010). In the USRB region, multiple depositional sequences contain dense loamy till interstratified with glacial outwash and finer-grained glaciolacustrine sediments ( Figure 5B). Such sediment packages represent one or more glacial advanceretreat sequences (Kempton et al., 1991;Stumpf and Atkinson, 2015;Grimley et al., 2016b). Proglacial outwash below Illinois Episode till, in the MBV, is hydrologically connected with older pre-Illinois Episode outwash (Figures 5A,B, 6) and together comprise the Mahomet Aquifer. During the final glacial retreat of the Illinois Episode, some glacial sequences were incised by meltwaters and created drainageways that were in places reoccupied during the Wisconsin Episode or postglacial period. For example, outwash underlying the present-day Sangamon River Valley ( Figure 5B) may represent a proto-Sangamon River that formed during a deglacial phase of Illinois Episode glaciation (Grimley et al., 2016b). After the Illinois Episode and before the Wisconsin Episode glaciation, a interglacial paleosol formed during the Sangamon Episode (∼0.13-0.05 Ma). The weathering profile of this ancient CZ was partially or completely eroded during the last (Wisconsin Episode) glaciation across most of the USRB, but is preserved in some area in the subsurface ( Figure 5C).
Following the last interglacial, the USRB was covered by glacial ice of the Lake Michigan Lobe of the Laurentide Ice Sheet during the last glacial maximum (Wisconsin Episode), ∼25-20 ka (Grimley, 2016a). This region includes several morainic systems (Figure 4) that mark multiple ice margin positions during the Wisconsin Episode glaciation (Hansel and Johnson, 1996; . At least two till units were deposited as a result of ice margin fluctuations (Wickham, 1979;Grimley et al., 2016b), with peat and proglacial loess preserved in places below the oldest till ( Figure 6C). Glaciofluvial sediment (gravelly sand to fine sand) were deposited in coalescing fans along moraine fronts and in meltwater streams that drained the ice margin. Some of these streams appear to have reoccupied older drainageways or lowlands that existed during the penultimate glaciation. Large proglacial lakes in the USRB during Wisconsin Episode deglaciation were not common, as advancing glaciers flowed down the regional topographic gradient and the ice margin was fairly well-drained by large streams leading to the southwest. However, in the present-day landscape, many small depressions remain poorly drained and stream networks are poorly integrated (although modified by agricultural ditching and drainage tiling) because of drainage disruptions and landscape flattening that resulted from multiple successive glaciations. Immediately following ice retreat, thin deposits of fluvial and lake sediment accumulated in lowlands within the low-relief landscape of the USRB. Thin discontinuous layers of sand and gravel, too minor to show in typical cross-sections (e.g., Figure 5C), are commonly found immediately above Wisconsin Episode till deposits. Loess deposits, 1-2 m thick, accumulated above the till from ∼22 to 16 ka (Nash et al., 2017) and effectively blanketed other glacial sediments. Periglacial climatic conditions for some centuries, immediately after ice margin retreat in the USRB, is evidenced by relict patterned ground and ice-wedge casts (Johnson, 1990). Such periglacial conditions may have led to some re-mobilization of loess and resedimentation into depressional and low-lying areas. The relict periglacial features, along with other aspects of the glacial legacy, today influence the diversity of soil types, soil moisture conditions, and shallow groundwater flow.
The history of Quaternary glacial-interglacial cycles significantly impacts the landscapes and soils found today in the USRB. The CZ here is developed in loess and glacial sediments that accumulated on an older interglacial landscape containing the Sangamon Geosol. Vestiges of this older landscape, including the paleosol and weathered lacustrine or alluvial sediments, are exposed in places along the Sangamon River Valley. Elsewhere, this older landscape is completely buried by deposits of the Wisconsin Episode glaciation and post-glacial streams. In places, late-glacial erosion has incised through modern alluvium into shallowly buried outwash. Additionally, anthropogenic excavations of outwash and alluvium in the Sangamon River Valley for aggregate has penetrated more deeply, perhaps into Illinois Episode outwash. The superposition or connectivity of permeable sedimentary units can provide a hydrologic connection between the surface and groundwater systems and a pathway for cycling nutrients and dissolved ions.
Soil development into loess and underlying glacial sediments has occurred on uplands during the last ca. 16 ka. Over this time period, soil development has occurred under a succession of vegetation including tundra, boreal and deciduous woodland, prairie, and cultivated landscapes. Today, the dominant soil type in the USRB are Mollisols (Mount, 1982), reflecting the dominance of prairie vegetation over much of the landscape during the Holocene. Pre-settlement woodlands were more typical within about 1 km of the Sangamon River valley (INHS, 2003), where Alfisols occur today. Approximately 2-5 m of fine-grained alluvial sediment are found in the valley, of which 0.5-0.9 m was deposited following European settlement in the nineteenth-century (Figure 5C; Grimley et al., 2017). Widespread devegetation and ditching of landscapes to increase agricultural land has resulted in an order of magnitude increase in alluvial sedimentation rates in the Sangamon Valley, averaging ∼5 mm/yr during the last 150 years compared with ∼0.7 mm/yr during pre-settlement times .
Topographic controls on deposition and erosion are evident at several temporal and spatial scales. Over the course of the Pleistocene, the MBV was preferentially filled with sediment so that today's surface topography, dominated by last glacial constructional landforms, no longer reflects this feature (Figures 4, 5A, 6). Stretches of the current Sangamon River occupy a former glacial meltwater channel filled with last glacial outwash. In some areas, the near-surface Wisconsin Episode outwash may be immediately underlain by Illinois Episode (Figure 5B), indicating that glacial meltwater drainage was focused into the same area during two separate glacial events. The modern Sangamon River Valley contains a Holocene alluvial record that indicates increased rates of floodplain sedimentation following the establishment of intensive agriculture in the late 1800s   (Figure 5C). At scales of tens to hundreds of meters, erosion of local topographic highs and sedimentation within swales is evident in soil profiles. Rates of erosion and deposition at this scale were also accelerated by landuse change (e.g., Norton, 1986;Kreznor et al., 1990). Even in low-relief landscapes such as in the USRB, subtle topographic gradients can significantly impact shallow groundwater flow and soil moisture conditions. Glacial moraines, along with finerscale periglacial features and other undulations on the glacial till plain (Figure 4), impose trends in soil drainage, resulting in soil associations (Mount, 1982) tied, in part, to glacial landforms and sediments at the macro-and meso-scales.

Minnesota River Valley, Minnesota
In south-central Minnesota, Quaternary deposits are typically ∼100 m thick (Figure 6). Flat-lying sedimentary Paleozoic bedrock beneath the Quaternary section is exposed in places along the Minnesota River valley and in the lower reaches of some major tributaries of the Minnesota River . Valleys cut into the Paleozoic bedrock section are present, but smaller in scale than the Mahomet Bedrock Valley in the USRB.
The majority of Quaternary sediments are glacial in origin, derived from repeated advances of multiple lobes of the Laurentide Ice Sheet. These distinct lobes were fed by ice from different regions with different bedrock geology and four distinct provenances for glacial tills have been identified based on clast lithology: Riding Mountain (NW), Winnipeg (N), Rainy (NE), and Superior (ENE) (Johnson et al., 2016). In east-central Blue Earth County, the pre-Wisconsin Episode units include unnamed tills of Winnipeg, Superior, and Rainy Provenance   (Figure 6). These tills are stratigraphically below the Browerville Formation, a widespread but discontinuous buried diamicton, that is also pre-Wisconsin Episode and likely pre-Illinois Episode in age (Johnson et al., 2016). Deposits of sorted sand and gravel, interpreted as meltwater deposits of these glaciations are present as discontinuous lenses, especially at the boundaries between till units .
During the Wisconsin Episode glaciation, tills of Riding Mountain provenance and mixed Riding Mountain and Winnipeg provenance were deposited in the east-central Blue Earth County . These tills comprise the New Ulm Formation, which is bounded below by an erosional unconformity where the advancing ice eroded earlier glacial deposits . Till sheets with distinct textures and distributions of clast lithology are recognizable within the New Ulm Formation and were likely deposited by adjacent ice streams draining distinct areas within what is commonly referred to as the Des Moines Lobe . The uppermost portions of the New Ulm Formation include supraglacial tills formed by melting of stagnant ice during the retreat of the Des Moines Lobe .
Glacial Lake Minnesota was a large proglacial lake that formed during the retreat of the Des Moines Lobe as water ponded behind the stagnant ice. Portions of Glacial Lake Minnesota were supraglacial and features such as ice walled lake plains are present where lake sedimentation was constrained by ice during deposition, leaving vaguely circular flat-topped hills of fine-grained sediment standing above subglacial tills . Radiocarbon dates place Glacial Lake Minnesota at a maximum age of between 14.6 and 13.8 ka . Glacial Lake Minnesota drained through a series of outlets with the youngest being to the north along the trend of a bedrock valley (Patterson and Hobbs, 1995). This same spillway was further incised by Glacial River Warren during the catastrophic drainage of Glacial Lake Agassiz.
Glacial Lake Agassiz, the largest proglacial lake in North America during the most recent deglaciation, formed as the Red River Lobe of the Laurentide Ice Sheet retreated north of the drainage divide separating the Mississippi basin from Hudson Bay. During its 5,000 year history, Glacial Lake Agassiz drained through four different outlets including the Minnesota and Mississippi drainages (Teller and Leverington, 2004). Drainage through the southern outlet is referred to as Glacial River Warren and this drainage carved a valley 45-70 m deep, now occupied by the Minnesota River, through pre-Wisconsin tills and Paleozoic saprolite, exposing relatively unweathered bedrock in places (Matsch, 1983). This southern outlet was occupied primarily early in the history of Lake Agassiz from ∼13.4 to 10.9 ka (Matsch, 1983) with a short-lived reoccupation of the southern outlet ending by 9.4 ka (Fisher, 2003). The abrupt downcutting of Glacial River Warren left tributaries including the Blue Earth River and LeSueur River stranded above the trunk stream (Gran et al., 2009). Knickpoints formed in these tributaries and continue to cut back today. In the LeSueur River more than 400 terrace surfaces are present in the basin below the knickpoint which has migrated ∼30 km upstream of the confluence with the Minnesota (Gran et al., 2009) (Figure 7). Within the lower LeSueur basin, bluff erosion and gully formation contribute to high turbidity levels and rapid sedimentation within Lake Peppin, a natural lake on the Mississippi River downstream of the junction with the Minnesota River (Gran et al., 2009(Gran et al., , 2013. Because of its late deglaciation compared with the USRB and CCW, the MRV has very thin or no loess deposited at the surface. Instead, the most recent deposits other than river terraces include modern lake and wetland deposits   (Figure 6). The shallow lakes and wetlands were largely drained to allow for agricultural use and drainage ditches now connect the low-relief uplands to the naturally steep and incised tributaries of the Minnesota River. In the Minnesota River Basin the legacy of glacial history is, therefore, felt in both the subsurface composition of the CZ and in the topography shaped by glacial processes that controls shallow subsurface water flow and the distribution of soils.

IDEALIZED EVOLUTION OF THE CZ THROUGH GLACIAL CYCLES
The CZ structure of each field laboratory in the IML-CZO clearly reflects variability in the climate, geomorphological, biological and geochemical systems that were dominant in the glacialinterglacial cycles. We generalize the histories revealed in the cross sections to produce an idealized model of interactions between the Laurentide Ice Sheet, landscapes of the Central Lowlands of the Interior Plain, and sedimentary processes during the most recent glacial-interglacial cycle. We assume, generally, that the sequence of events was generally similar during each of the Quaternary glacial episodes. One significant difference ais that the flow directions of earlier ice advances were likely controlled by topography on the bedrock surface until deeply incised valleys were filled in. In addition, the preglacial landscapes were covered with thick residuum or saprolite that was incorporated in glacially-transported sediment (Patterson and Boerboom, 1999). In many places, the oldest glacial till contains an abundance of local weathered material derived from preglacial soils and weathering profiles, and fragments of the local bedrock (Willman et al., 1966).
Prior to the last glaciation, the interglacial landscape had an integrated drainage network, variable relief, a hydrologically connected soil cover and associated interglacial weathering profile that was formed in fluvial, lacustrine, colluvial, and eolian is developed in the surficial sediments, including alluvium (black) in the major river valleys. (B) As glaciers (light blue) crossed the region, they actively eroded their beds depositing glacial till (light gray). Sorted fluvial or lacustrine sediments (dark gray) were deposited in front of the ice margin (Continued) FIGURE 7 | partially filling the bedrock valleys. Weathering profiles developed on former land surfaces (light gray squiggles) are preserved where glacial erosion did not entirely remove the older sediments. Vegetation beyond the glacial margin was typical of tundra and boreal forest found in northern latitudes today. (C) Glacial conditions (dry and cold) continue beyond the ice margin. Windblown sand and silt (loess), shown in tan, accumulated, covering soils developed in the bedrock and preglacial sediments. Minimal chemical and physical weathering occurred beneath the glacier. The weathering profile continued to develop in the landscape beyond the glacial margin. In places, soil formation is superimposed on previous weathering profiles. Sparse tundra vegetation likely colonized areas near the ice margin. (D) Following glacier retreat, loess further accumulated on the land surface, completely burying the glacial deposits. The vegetation and soil succession began as tundra was replaced by boreal forest, and later by mixed hardwood/coniferous forest.
(E) As the climatic conditions moderated during the interglacial, tallgrass prairie, and oak-hickory savannas replaced the hardwood forests. Continuing fluvial activity established new drainage networks, and valleys were incised into loess and glacial deposits.
deposits overlying a sequence of older sediments from previous glacial-interglacial cycles ( Figure 7A). The vegetation cover at that time was likelysimilar to the pre-settlement vegetation in the region including oak-hickory forest, savannah, and prairie (Grüger, 1972;Curry et al., 2010). Soils were dominantly fertile Alfisols, and Mollisols, similar in composition to those occurring in the CZO today (Follmer, 1983). As the Laurentide Ice Sheet flowed into the region, the CZ underwent dramatic changes over a relatively short period of geologic time (Figures 7B,C). The regional climate cooled during the onset of glaciation, which led to profound changes in the rates and trajectories of all CZ processes. At the time that the ice sheet margins advanced southward into the region, a parkland or boreal forest ecosystem had colonized the landscape (Garry et al., 1990;Curry and Baker, 2000).
As the ice sheet flowed across the region toward its terminal position, the glaciers actively eroded their beds in many places and compacted or deformed underlying materials to varying extent ( Figure 7C). Subglacial tills and associated meltwater sediments were deposited. Partial or complete weathering profiles (paleosols) developed on former land surfaces were sometimes preserved beneath the ice, primarily in lower parts of former landscape where subglacial erosion was less effective. In front of the ice margin, sorted glaciofluvial and glaciolacustrine sediments accumulated, partially or completely filled pre-existing valleys. Locally, the glaciers hosted diverse microbial life both on the ice surface and beneath (Boetius et al., 2015). During glacial maxima a zone of permafrost extended beyond the ice margin in some areas (Garry et al., 1990) and tundra plants surrounded ice-walled lakes in former supraglacial environments (Curry and Petras, 2011). Cold climate soils (Gelisols) and tundra vegetation occupied landscapes dominated by intensive solifluction, irregular, and fluctuating stream flows and localized eolian deposition (Johnson, 1990;Schwert et al., 1997;Mason, 2015). Beyond the permafrost zone, or in more southern parts of the Central Lowlands, boreal forest and associated cool climate soils (Cryalfs and Spodosols) occupied the landscape.
Critical zone processes were radically different under full glacial conditions than they had been during the interglacial.
Beneath the glacier erosion and deposition significantly modified the CZ structure and added new sediments. The biogeochemical system was altered and shifted from a relatively warm, oxygenrich environment to an oxygen deficient, cool, and biologically limiting condition (Anderson et al., 1997;Skidmore et al., 2000;Boetius et al., 2015). In the tundra zone beyond the ice margin, similarly cool and oxygen deficient conditions dominated all but the uppermost part of the CZ. In the annually thawed part of the CZ (the active zone), components of past CZs were reworked and transformed by mass wasting, eolian and fluvial activity, and enriched in organic matter as a result of significant reduction in the rate of carbon decomposition, although organic carbon inputs were low relative to warmer periods. The CZ was also significantly impacted in the boreal forest zone where dry and cool conditions prevailed, fostering intense seasonal eolian deposition. Deflation of sand and silt from outwash deposits upwind resulted in the accumulation of wind-blown sand and silt (loess). In the CZ, this material was aggraded at varying rates through time, depending upon the distance from the source and multi-centennial fluctuations in climate and/or glacial sediment production Muhs and Bettis, 2003;Nash et al., 2017). CZ processes in this zone were very different than in the periglacial zone and beneath the ice sheet. Here a warmer (but still cool), oxygenated, biologically enriched environment having internal (vertical) connectivity inherited from the interglacial period drove active chemical weathering and maintained the connection between the active weathering zone and deeper parts of the CZ.
Termination of the glacial period, including a shift to warmer climatic conditions, occurred over a much shorter time period than the glacial advance phase (Lisiecki and Raymo, 2005). Glacial retreat was not a continuous process, and often the stagnant ice margins readvanced, overriding previously deglaciated terrain . A series of recessional moraines were formed during these readvances. In the decades and centuries immediately following deglaciation, the dominant geomorphic processes shifted from an active or stagnant ice system to paraglacial, lacustrine, eolian, and fluvial processes. Over millennia, loess blanketed the uplands downwind of the major valleys ( Figure 7D). As the climate moderated and glacial sedimentation ceased, vegetation began to recolonize recently glaciated areasand ecosystemsevolved from tundra ( Figure 7D) to boreal forest, and finally to hardwood forest (e.g., King, 1981;Davis and Shaw, 2001;Curry et al., 2010). These changes in vegetation type were accompanied by an evolution of soils from Gelisols and Cryalfs to Alfisols and Spodosols characteristic of the boreal forest, and later to Alfisols under mixed hardwood/coniferous forest vegetation (Figures 7D,E).
As the climate warmed further during the interglacial, tallgrass prairie, and oak-hickory savannas replaced the hardwood forests ( Figure 7E). As vegetation changed through the period in response to variations in climate and fire regime, soils evolved from the Alfisols of the early postglacial period to Mollisols and transitional soils (Mollic Hapludalfs) in the prairie and savanna (Bartlein et al., 1984;Baker et al., 1992;Nelson et al., 2006). The postglacial fluvial system continued to develop from initial headward extension, through to meander belt development and episodes of alluvial deposition punctuated by channel incision (Knox, 1972;Bettis et al., 2008).
Transformation of ecosystems and climatic conditions during the glacial-interglacial cycle caused long-term shifts in water flux, chemical weathering, and CZ evolution. This idealized progression of CZ evolution illustrates the dynamic nature of the subsurface, and highlights its direct relationship with geomorphic processes, and the complexity of the geologic materials.

EVOLUTION OF THE MODERN CZ
We have illustrated three examples of the structure of the CZ within the formerly glaciated Midwest US (e.g., Figures 3, 5, 6) and shown an idealized evolution of the CZ through a glacial cycle (Figure 7). We argue that the structure of the CZ plays an important role in the current state and future evolution of the CZ within the glaciated portion of the Central Lowlands, now under intensive human management. Management practices vary geographically to accommodate soil and drainage differences associated with glacial, periglacial, and fluvial landforms. Thus, the history of environmental change influences CZ structure in terms of the distribution of materials in the subsurface and the surface topography. Both topography and subsurface sedimentary structure influence ongoing CZ processes and the response of the CZ to anthropogenic forcing. Perhaps most important is the recognition that different areas of the Central Lowlands have different preglacial, glacial, and interglacial histories, creating different CZ structure and influencing future evolution.
Modern soil profiles in the glaciated parts of the Central Lowlands represent the transition between the above and below ground parts of the CZ. Their properties are the cumulative effects of additions, removals, transfers and transformations of energy, mineral, and biological materials in an evolving landscape over the past 15 ka (late glacial and Holocene). These soil profiles (and their horizonation) evolved over changing climate and vegetation systems, but resulted in mainly Alfisol or Mollisol profiles in upland areas of the study area (Mount, 1982;Ruhe, 1984;Baker et al., 1992). In Illinois and Iowa, the base of the A horizon (organic-rich topsoil) occurs at ∼0.1-0.5 m depth and the base of the B horizon (pedogenically altered subsoil) at ∼1-1.5 m depth. The depth to the base of the C horizon (oxidized, slightly altered) depends on the sub-regional geologic history, and typically occurs at ∼3-5 m depth in parts of the region glaciated during the Wisconsin Episode it (Mount, 1982;Ruhe, 1984;Tandarich et al., 1994). In such areas, matrix carbonates (calcite, dolomite) have been leached to 2-3.5 m below the surface. Below the C horizon is a D horizon that consists of essentially unaltered, unoxidized (containing minor pyrite) and calcareous glacial till or other glacial sediment that has generally remained below the long-term water table (Tandarich et al., 1994). Beyond the Wisconsin Episode glacial margin, where the modern soil is formed in thicker loess, the oxidized C horizon is thicker and carbonate leaching ranges from 1.5 to >8m along a west-east precipitation gradient (Ruhe, 1984). In addition to climatic effects, soil profile thicknesses, leaching depths, and oxidation depths also vary on a more local scale with soil drainage, vegetation type, parent material, and soil age (Birkeland and, 1999).
Landscape relief is one way in which history and CZ structure influence ongoing processes and the trajectories for future evolution. Pre-agricultural landscape relief and drainage were strong controls on the pace of agricultural development (e.g., Bogue, 1963) and influenced both land management and system responses to intensive agriculture. This pre-agricultural landscape relief is itself a strong function of the local glacial history. Mature fluvial networks were present in CCW prior to agricultural development because deposition during recent glaciations was limited to loess that blanketed the pre-existing landscape and did not rearrange drainage patterns. These mature fluvial networks create landscape relief that provides energy for sediment transport down hillslopes and, ultimately in to stream networks (Abban et al., 2016). In contrast, fluvial networks in the USRB prior to intensive agriculture, were limited in extent and did not drain large areas of the upland surface (Rhoads et al., 2016). The limited landscape relief in the USRB likely contributes to lower rates of post-settlement alluvial deposition than in CCW  and a fluvial suspended load dominated by near-channel sources rather than agricultural fields (Neal and Anders, 2015). In the MRV, the legacy of drainage of Lake Agassiz strongly controls landscape relief and ongoing erosion and sediment transport that is focused near the entrenched Minnesota River and its tributaries (Gran et al., 2013).
The impact of glacial history and CZ structure on modern CZ processes stretches beyond the topographic influences of landscape relief. For example, differential permeability in the materials within the CZ creates pathways for preferential flow of groundwater in the near surface and in deep aquifers (Atkinson et al., 2014), controlling spatial variability in the flux of reactive water through the CZ. The distribution of sorted sand and gravel and clayey paleosols in the subsurface is thus crucial for predicting subsurface reactions. Sorted sand and gravel deposits occur in a variety of contexts in the IML-CZO. Stacked units of sorted sediment of different ages are present in the USRB at depth within Mahomet Bedrock Valley (preglacial alluvium and pre-Illinois Episode outwash Figure 5B) and on the surface (modern alluvium, Wisconsin Episode and Illinois Episode outwash Figure 5C). In CCW, there is also juxtaposition of sorted sediment of different ages (pre-Illinois Episode outwash and Holocene alluvium Figure 3). Isolated and/or thin bodies of sand and gravel also occur within the glacial units at all IML-CZO sites. The geometry of these deposits and their connectedness is vital for predicting the interactions between surface and subsurface water in the CZ (Stumpf and Dey, 2012;Atkinson et al., 2014). The importance of coarse sediments in the subsurface extends past their role in influencing subsurface weathering and groundwater flow: the water stored within these sediment is a resource. Groundwater within the Mahomet Aquifer is a valuable resource for humans, providing high-quality drinking water for much of east-central Illinois, including the Champaign-Urbana area, and was designated a Sole Source Aquifer (United States Environmental Protection Agency Notice, 2015). In much of the USRB, the aquifer is protected by multiple overlying fine-grained till units from various glaciations; however, a groundwater age of <11,000 years from isotopic data  suggests a hydrologic connection with surface recharge over millennial time scales. Another central Illinois city in the USRB (Decatur, IL) relies on surface water from a reservoir for its drinking water and is thus vulnerable to water quality issues resulting from high nutrient levels in the Sangamon River (Keefer et al., 2010).
As we continue exploring the CZ at the IML-CZO we anticipate finding more examples illustrating how the geologic structure and history of CZ evolution constrain the present function and future evolution of the CZ. The importance of CZ structure to understanding CZ process in the IML-CZO further motivates us to present a new conceptual framework for understanding CZ development in the glaciated Midwest.

The Patchwork CZ: New Metaphor for the IML-CZO and Glaciated Central Lowlands CZ
The CZ of the IML-CZO and the CZs of other depositional landscapes (Figure 1) are not well-described by the feed-through reactor model of the CZ (Anderson et al., 2007). We propose a new metaphor for the CZ of the glaciated portion of the Central Lowlands of the Interior Plain based on the human social systems of this region. Beginning at least 13,000 years ago (Lothrop et al., 2016) a diverse group of people with a variety of histories and origins began to arrive in waves. Acknowledging that the history of people in this region features conflict and colonization, we suggest that the present inhabitants were shaped by their individual and shared experiences of changing local, regional, national, and global human culture. The story of this culture must include both the diversity of its constituents and their ongoing evolution together. The patchwork quilt is a metaphor that has been used to describe American culture (e.g., Chametzky, 1989) which we propose as an image that captures the essential features of the IML-CZO study sites and, more broadly, the CZ of areas that experienced net deposition through cycles of continental glaciation. While the regular division of agricultural fields in the region along evenly spaced section lines easily brings to mind geometrical patchwork quilt designs, the vertical structure of the CZ is more akin to a crazy quilt which features irregularly shaped pieces, exotic fabrics, and embroidery added on top of the fabric (Figure 8).
This image of the patchwork CZ emphasizes two key characteristics of the CZ of the IML-CZO sites: (1) the mineral components of the CZ are spatially variable and have been influenced by their individual histories of transport and weathering, and, (2) CZ processes vary enormously at timescales ranging from glacial cycles (10 5 years) to decades. As is evident from the cross sections presented (Figures 3, 5, 6), the CZs of the IML-CZO are characterized by spatial variability in the mineral and geologic components. The physical and chemical properties of the sediments in the CZ vary due to differences in their source areas, prior weathering histories, processes of erosion and deposition, and postglacial alterations. This multi-sequence evolution of CZ materials before, during and subsequent to transport, is a primary control on CZ structure in the IML-CZO and in other areas with accretionary sedimentary sequences, including areas where the CZ is formed in sedimentary bedrock. The CZ, like the soil profile in its upper part, is a historical body, one whose history must be appreciated to understand its present condition and function. The dramatic and quantifiable impacts of climate change, geomorphic processes, soil development, and vegetation on the CZ structure over the last few million years highlights the need to consider CZ process timescales in reference to variations of climate, vegetation, microbial, and geomorphic processes over those timescales. This point is important for critical zone science even beyond the glaciated Central Lowlands. Climate variability at timescales of centuries to several hundred thousand years is characteristic of the last few million years of Earth history, globally. This variability affected the chemical, physical and biological processes in the CZ across the planet. Similarly, while recent anthropogenic impacts on ecosystem services (Richardson and Kumar, 2017), hydrology, and climate are pervasive in heavily modified landscapes like the IML-CZO, more subtle changes in land use, climate and water are occurring in the CZ elsewhere. The imprint of anthropogenic changes on CZ structure and processes is occurring in tandem with the pre-settlement glacial legacy of the landscape (Kołodynska-Gawrysiak et al., 2017).
We hypothesize that the diversity of materials and variability in climate and ecosystems characteristic of the glaciated Central Lowland influence not only the structure of the CZ but also its current function and its trajectory of future evolution. Future work at the IML-CZO sites will explore the importance of the subsurface permeability contrasts in glacial materials and their geometry in driving current CZ process and the response of the CZ to changing agricultural practices and climate variability. Another target for future research is the importance of smallscale glacial and periglacial topography including hummocky moraines and meter-scale closed depressions in setting the rates of physical, chemical and biological CZ processes. Additionally, soils and weathering profiles provide an integrated record of CZ processes through time and we seek to develop CZ process models capable of hindcasting the observed spatial variability. Moreover, the influences of previous history and temporal variability also impact CZs beyond the areas glaciated by the Laurentide Ice Sheet. Thus, our new conceptual model of the CZ in the glaciated Central Lowlands as a patchwork quilt motivates future research in other regions with thick sedimentary deposits (e.g., Coastal Plain, Great Plains, large basins) for consideration of the mechanisms that shape CZ structure, timescales over which the structure is shaped, and potential for CZ structure to influence near-surface processes and respond to future environmental changes.

AUTHOR CONTRIBUTIONS
All authors contributed to framing of the problem, discussion to create and revise the conceptual model, and the writing and

ACKNOWLEDGMENTS
We thank Harvey Thorleifson and Richard Lively of the Minnesota Geological Survey for providing electronic files of published cross-sections of Blue Earth County, Minnesota as a starting point for our illustration, and Phil Kerr, Iowa Geological Survey for assistance with the CCW cross-section.