Erosions on the southern Tibetan Plateau: Evidence from in-situ cosmogenic nuclides 10Be and 26Al in fluvial sediments

Investigating topographic and climatic controls on erosion at variable spatial and temporal scales is essential to our understanding of the topographic evolution of the orogen.In this work,we quantified millennial-scale erosion rates deduced from cosmogenic 10Be and 26Al concentrations in 15 fluvial sediments from the mainstream and major tributaries of the Yarlung Zangbo River draining the southern Tibetan Plateau (TP).The measured ratios of 26Al/10Be range from 6.33 ± 0.29 to 8.96 ± 0.37,suggesting steady-state erosion processes.The resulted erosion rates vary from 20.60 ± 1.79 to 154.00 ± 13.60 m Myr-1,being spatially low in the upstream areas of the Gyaca knickpoint and high in the downstream areas.By examining the relationships between the erosion rate and topographic or climatic indices,we found that both topography and climate play significant roles in the erosion process for basins in the upstream areas of the Gyaca knickpoint.However,topography dominantly controls the erosion processes in the downstream areas of the Gyaca knickpoint,whereas variations in precipitation have only a second-order control.The marginal Himalayas and the Yarlung Zangbo River Basin (YZRB) yielded significantly higher erosion rates than the central plateau,which indicated that the landscape of the central plateau surface is remarkably stable and is being intensively consumed at its boundaries through river headward erosion.In addition,our 10Be erosion rates are comparable to present-day hydrologic erosion rates in most cases,suggesting either weak human activities or long-term steady-state erosion in this area.     

风沙流的输沙率沿程变化规律

将风力作用下的地表侵蚀分解为气流直接侵蚀和跃移沙粒冲击侵蚀两个部分,在冲击侵蚀率中考虑颗粒与床面之间以及颗粒之间碰撞影响的前提下,结合地表蚀积平衡方程从理论上导出了风沙流的输沙率沿程变化公式。公式的预测值与实测资料吻合较好。以该公式为基础,进一步得到了流体侵蚀率、冲击侵蚀率和总侵蚀率各自沿程的理论变化表达式以及饱和路径长度公式,并且发现流体侵蚀率沿程逐渐减小,而冲击侵蚀率和总侵蚀率沿程变化规律都是开始时增大、然后减小,存在极大值;流体侵蚀累积量在总累积侵蚀量中所占的比例并不是固定不变,而是随气流速度的增大而增大。     

用风洞实验方法计算土壤风蚀量的时距问题

通过对延津县沙质壤土原状土进行风蚀风洞模拟实验，探讨了土壤风蚀量计算中的时距问题。结果表明，前人使用“风蚀模数”概念来计算土壤风蚀量会导致较大误差，原因在于同一风蚀事件中单位时间内的风蚀量随吹蚀时间的延续而递减，在风洞实验计算风蚀量时使用“累积风蚀量”概念将会有更精确的结果。     

不同砾石覆盖度戈壁床面风蚀速率定量模拟

通过风洞实验,利用称重传感器自动记录风蚀观测样方重量变化过程,对供沙条件下不同砾石覆盖度戈壁床面风蚀速率进行了定量模拟研究。结果表明:砾石覆盖度是影响戈壁风蚀速率的关键因子,戈壁床面风蚀速率随砾石覆盖度增加按指数规律递减。各实验风速下,砾石覆盖度>50%时,戈壁床面风蚀速率随砾石覆盖度增加而减小量有限,甚至无风蚀发生;而盖度从10%到50%时,风蚀速率显著减小。因此,两种实验粒径砾石(3 cm与4 cm)至少在50%盖度时才能达到较好的风蚀防治效果。戈壁风蚀防护机理主要是砾石覆盖度的增加增大了砾石间沙粒的临界起动剪切风速,而且减少了作用在砾石间可蚀地表的剪切压。与沙质对照床面相比,10%~90%砾石盖度戈壁床面沙粒临界起动剪切风速增大了0.8~3.4倍,只有0.5％~28%的剪切压作用在砾石间可蚀地表。     

地球化学示踪在现代土壤侵蚀研究中的应用

土壤侵蚀导致土地退化、农产品产量和品质下降、水环境污染等诸多问题 ,引起各方面的广泛关注 ,而且这个问题随着持续增长的人口压力和农业生产方式深刻变革的影响而日益严重。由于目前广泛使用的高差法、遥感研究法、RUSLE、野外调查法、水土流失监测点法等土壤侵蚀研究方法存在着一些不足 ,因此 ,核素地球化学示踪法、稀土元素示踪法、土壤地球化学指标法等地球化学示踪方法应运而生 ,在土壤侵蚀研究中发挥其独特作用。在比较了各种传统的土壤侵蚀研究方法的优势和不足的基础上 ,作者认为 :1 137Cs、 2 10 Pb可用于研究 30、 40年左右平均土壤侵蚀速率 ,7Be可用于示踪季节性土壤侵蚀和一次降雨事件的土粒运移 ,并说明这些核素的示踪原理和侵蚀速率的定量依据 ;2人为施放稀土元素 ( REE)示踪法适合定量一次或多次降雨事件的侵蚀速率 ,可用于土壤侵蚀理论的研究 ,尤其适合于室内模拟实验的研究 ;3土壤磁化率法和土壤理化性质指标法研究土壤侵蚀是一种处于起步阶段的侵蚀研究方法 ,适用于土壤质量退化研究 ,但要用于土壤侵蚀研究特别是侵蚀速率的定量研究还很不成熟。     

Soil erosion in the Falkland Islands: an assessment

Observations and data concerning soil erosion on clay-rich, sand-rich and organic soil materials are reported for the Falkland Islands, South Atlantic, where soil erosion is both widespread and currently active, though rather uneven in its distribution. The most extensively eroded areas are located on coastal sand deposits and where sandy soils occur at inland sites. There is evidence that some erosion is natural and was initiated before the Islands were settled; other erosion is of more recent origin and can be associated with land management practices, principally grazing management, stocking rates, and the frequency and intensity of pasture burning. The impact of soil erosion on the agricultural economy of the Islands appears to be negligible and it is not perceived to be an important environmental issue. This lack of concern probably arises because the on-farm and off-farm costs of erosion have never been determined and there are no data relating to rates of erosion. It is suggested that a research programme be established aimed at monitoring erosion rates, assessing land use practices that initiate erosion, and assessing the costs of erosion. Furthermore, guidelines for soil conservation should be established. If pasture burning is to continue, a recommended code of practice should be drawn up. Given the importance of the sheep industry to the Falkland Islands, it is desirable that the grazing potential be maintained.     

Denudation rates and a topography-driven rainfall threshold in northern Chile: Multiple cosmogenic nuclide data and sediment yield budgets

The quantification of geomorphic process rates on the outcrop- and the orogen-scale is important to describe accurately the interaction between the relative effects of erosion, tectonics and climate on landscape evolution. We report single and paired cosmogenic nuclide (¹⁰Be, ²⁶Al and ²¹Ne) derived erosion rates and exposure ages on hillslope interfluves from the tectonically active western central Andes that show a distinct spatial variation. A positive correlation of erosion rates with elevation and present-day rainfall rates is observed. Erosion rates at lower altitudes–the hyperarid Coastal Cordillera and the Western Escarpment with the northern part of the Atacama Desert–are extremely low and of the order of 10–100 cm/My (nominal exposure ages 1–6 My). In contrast, erosion rates at higher altitudes–the semiarid Western Cordillera–range up to 4600 cm/My (nominal exposure ages 0.02–0.1 My). This latter average long-term bedrock erosion rate record, suggested to be coupled to an orographically controlled pattern of rainfall, is also reflected in the pattern of denudation rates derived from a short-term decadal record of limited sediment yield data. Specifically, denudation rates calculated from sediment flux data are of a similar order of magnitude as erosion rates deduced from long-lived cosmogenic nuclides from bedrock hillslope interfluves of the Western Cordillera. Nevertheless, the production and the supply of sediment from the western Andean slope are very limited.Analysis of multiple cosmogenic nuclides allows simultaneous determination of erosion rates and exposure ages but also reveals complex exposure histories of non-bedrock samples, such as boulders or amalgamated clast samples. Notably, this study shows that saturation of nuclides, usually assumed in studies where only a single nuclide is analyzed, is rather the exception than the rule, as revealed by erosion island plots. Constant erosion that started much later than the formation age of the rocks or episodic erosion by spalling can partially explain non-steady-state concentrations and more complicated exposure scenarios. Furthermore, the use of multiple nuclides with different half-lives allowed us to infer that no significant variations in long-term erosion rates have occurred and that at the Western Escarpment erosion rates have been low and constant for most of the late Neogene. Nevertheless, the time intervals necessary to reach steady-state concentrations for cosmogenic nuclides can be quite different from those needed for landscapes to reach steady state.     

The use of multilevel modeling in evaluating controls on erosion rates on inter-tidal shore platforms,Kaikoura Peninsula,South Island,New Zealand

Repeated measurements of erosion rates on seashore platforms on Kaikoura Peninsula, South Island, New Zealand were gathered over a 28-month period using a micro-erosion meter (MEM). The resultant hierarchical dataset was analysed using multilevel regression modeling. Results indicate that there is persistent variation of rates across measurement occasions within individual points and significant variation is also found across sites (or local ‘contexts’). There is little evidence for significant variation across individual points. The explanatory variables of platform/rock type, distance along platform and seasonality all have a significant effect on the erosion rates and aid the explanation of some of the variation in erosion rates across sites.     

CARBONATE EROSION RATES IN SOUTHWESTERN WISCONSIN

Between 1977 and 1979 three sets of erosional weight loss tablets were em-placed in the karst of southwestern Wisconsin to measure both potential and actual rates of carbonate erosion. Tablets of Jamaican White Limestone, which has been used in many other karst areas, were used to gauge relative erosion potential while tablets of the local bedrock, Prairie du Chien dolomite, were used to assess absolute rates. Tablets of Yugoslavian limestone were emplaced as part of a world-wide study of comparative rates. Tablets were recovered in 1983 and 1984. Although potential erosion rates are similar to those in other carbonate areas, the sluggish dissolution of the indigenous dolomite retards karst development. Soil cover increases erosion rates and rates at a given depth are consistent over sites of up to 10 m². The highest erosion rates were recorded by the Yugoslavian tablets but betweensite variations caution that results based on a single site may be unrepresentative. Over 37% of all tablets were affected by mechanical weathering processes, probably freeze-thaw action.     

SHORELINE CHANGE AT THOMPSON ISLAND,BOSTON HARBOR,MASSACHUSETTS, 1938-1977

A photogrammetric and sediment analysis is presented to illustrate the relationship between beach erosion and seacliff recession on Thompson Island, Boston Harbor, Massachusetts. Aerial photographs taken in 1938, 1952, 1963, and 1977 were measured to determine rates of shoreline change around the island. The 39-year average rate of beach erosion is 0.3 m/yr ± 2% with an average rate of cliff recession at 0.2 m/yr ± 2%. Rates of beach erosion between six orientations that reflect principal wave approach to the island were not found to be significantly different (0.05) for the 39-year period. The rank order correlation between beach erosion and cliff recession for these six orientations was moderately well correlated (0.63). To determine possible controls of cliff erosion, 58 sediment samples were collected from glacial cliffs along the shoreline. The textural composition was determined, and then tested with discriminant function analysis. Partial correlation analysis between beach erosion and cliff recession holding a surrogate for sediment size constant improved the rank order from 0.63 to 0.84. The results indicate that coarser-grained cliffs recede at faster rates, but with less erosion occurring on adjacent beaches. Alternatively, finer-grained cliffs recede at slower rates, but with greater erosion occurring on adjacent beaches.     

Relief-rejuvenation and topographic length scales in a fluvial drainage basin, Napf area, Central Switzerland

This paper explores how, and to what extent, a phase of relief-rejuvenation modifies the mode of surface erosion in an approximately 63 km² drainage basin located at the northern border of the Swiss Alps (Luzern area). In the study area, the retreat of the Alpine glaciers at the end of the Last Glacial Maximum (LGM) caused base level to lower by approximately 80 m. The fluvial system adapted to the lowered base level by headward erosion. This is indicated by knickzones in the longitudinal stream profiles and by the continuous upstream narrowing of the width of the valley floor towards these knickzones. In the headwaters above these knickzones, processes are still to a significant extent controlled by the higher base level of the LGM. There, frequent exposure of bedrock in channels and especially on hillslopes implies that sediment flux is to a large extent limited by weathering rates. In the knickzones, however, exposure of bedrock in channels implies that sediment flux is supply-limited, and that erosion rates are controlled by stream power.The morphometric analysis reveals the existence of length scales in the topography that result from distinct geomorphic processes. Along the tributaries where the upstream sizes of the drainage basins exceed 100,000–200,000 m², the mode of sediment transport and erosion changes from predominantly hillslope processes (i.e., landsliding, creep of regolith, rock avalanches and to some extent debris flows) to processes in channels (fluvial processes and debris flows). This length scale reflects the minimum size of the contributing area for channelized processes to take over in the geomorphic development (i.e., threshold size of drainage basin). This threshold size depends on the ratio between production rates of sediment on hillslopes, and export rates of sediment by processes in channels. Consequently, in the headwaters, erosion rates and sediment flux, and hence landscape evolution rates, are to a large extent limited by weathering processes. In contrast, in the lower portion of the drainage basin that adjusts to the lowered base-level, rates of channelized erosion and relief formation are controlled mainly by stream power. Hence, this paper shows that base-level lowering, headward erosion and establishment of knickzones separate drainage basins in two segments with different controls on rates of surface erosion, sediment flux and relief formation.     

耕作土壤表面的空气动力学粗糙度及其对土壤风蚀的影响

土壤表面粗糙度是影响耕作土壤抗风蚀能力的一个重要因素。根据风速廓线计算得到的空气动力学粗糙度,可以简捷而有效地刻画土壤表面的空气动力学性质。风洞模拟实验表明,耕作土壤表面的空气动力学粗糙度主要取决于暴露地表的土块直径,在土块大致均匀分布的条件下,直径愈大,空气动力学粗糙度愈大。土壤风蚀速率则随空气动力学粗糙度的增大而迅速减小,二者具有良好的相关性。     

地表岩石侵蚀速率对宇生核素暴露测年影响的研究

利用原地生宇生核素测定暴露年代时,通常会假设地貌体侵蚀速率为0。研究表明,该假设会低估地貌体的真实暴露年代。搜集2009~2012年全球不同区域56个岩石样品的宇生核素¹⁰Be测年数据,探讨侵蚀速率为0对于侵蚀速率为0.5、1以及2 mm/ka的样品,在不同暴露尺度上对暴露年代计算的影响幅度。结果表明,对于1×10⁴a尺度的样品暴露年代可能低估约0.5%,1%,2%;对于10×10⁴a尺度的样品可能低估约5%,7%,20%;对于50×10⁴a尺度的样品可能低估约40%,70%甚至100%以上。     

Relations between denudation,glaciation, and sediment deposition: implications from the Plio‐Pleistocene Central Alps

Despite abundant data on the early evolution of the Central Alps, the latest stage exhumation history, potentially related to relief formation, is still poorly constrained. We aim for a better understanding of the relation between glaciation, erosion and sediment deposition. Addressing both topics, we analysed late Pliocene to recent deposits from the Upper Rhine Graben and two modern river sands by apatite fission‐track and (U‐Th‐Sm)/He thermochronology. From the observed age patterns we extracted the sediment provenance and paleo‐erosion history of the Alpine‐derived detritus. Due to their pollen and fossil record, the Rhine Graben deposits also provide information on climatic evolution, so that the erosion history can be related to glacial evolution during the Plio‐Pleistocene. Our data show that Rhine Graben deposits were derived from Variscan basement, Hegau volcanics, Swiss Molasse Basin, and the Central Alps. The relations between glaciation, Alpine erosion, and thermochronological age signals in sedimentary rocks are more complex than assumed. The first Alpine glaciation during the early Pleistocene did not disturb the long‐term exhumational equilibrium of the Alps. Recent findings indicate that main Alpine glaciation occurred at ca. 1 Ma. If true, then main Alpine glaciation was coeval with an apparent decrease of hinterland erosion rates, contrary to the expected trend. We suggest that glaciers effectively sealed the landscape, thus reducing the surface exposed to erosion and shifting the area of main erosion north toward the Molasse basin, causing sediment recycling. At around 0.4 Ma, erosion rates increased again, which seems to be a delayed response to main glaciation. The present‐day erosion regime seems to be dominated by mass‐wasting processes. Generally, glacial erosion rates did not exceed the pre‐glacial long‐term erosion rates of the Central Alps.     

Modelling the impact of land-use change and farm dam construction on hillslope sediment delivery to rivers at the regional scale

Landscapes in southeastern Australia have changed dramatically since the spread of European colonisation in the 19th century. Due to widespread forest clearance for cultivation and grazing, erosion and sediment yields have increased by a factor of more than 150. In the 20th century, erosion and sediment yield were reduced again due to an increasing vegetative cover. Furthermore, during the last decades, thousands of small farm dams were constructed to provide drinking water for cattle. These dams trap a lot of sediment, thereby further reducing sediment delivery from hillslopes to river channels. Changes in sediment delivery since European colonisation are documented in sediment archives. Within this study, these changing rates in hillslope erosion and sediment delivery were modelled using a spatially distributed erosion and sediment delivery model (WATEM/SEDEM) that was calibrated for Australian ecosystems using sediment yield data derived from sedimentation rates in 26 small farm dams. The model was applied to the Murrumbidgee river basin (30,000 km²) under different land-use scenarios. First, the erosion and sediment yield under pre-European land-use was modelled. Secondly, recent land-use patterns were used in the model. Finally, recent land-use including the impact of farm dams and large reservoirs was simulated. The results show that the WATEM/SEDEM model is capable of predicting the intensity of the geomorphic response to changes in land-use through time. Changes in hillslope erosion and hillslope sediment delivery rates are not equal, illustrating the non-linear response of the catchment. Current hillslope sediment supply to the river channel network is predicted to be 370% higher compared to the pre-European settlement period, yet farm dams have reduced this back to 2.5 times the pre-19th century values. The role of larger reservoirs is even more important as they have reduced the current sediment supply downstream to their pre-European values, thus completely masking the increased hillslope erosion rates from land-use change. However, the model does so far not include valley widening and sediment storage in river systems. Therefore, modelled rates of sediment delivery are lower than observed values.     

Tree roots — Methodological review and new development in dating and quantifying erosive processes

Exposed roots have been used in Dendrogeomorphology since the 1960s to determine erosion rates. This was done by using the root axis as the relative position of the former soil layer and then relating this measure to the age of the root in order to quantify the amount of soil erosion over time. Recently, wood anatomical research revealed the possibility of determining the first year of exposure of roots due to specific anatomical changes within the respective annual ring. As a consequence, using exposed roots that are still in contact to the soil surface, it is possible to reconstruct the size and position of the root at the time of exposure. Consequently, a new equation is presented to calculate the thickness of the eroded soil layer allowing a detailed reconstruction of erosion rates. The technique also helps to determine if the root was exposed by erosional processes or just by its ongoing secondary growth without any erosion involved.     

Erosion rates in badland areas recorded by collectors,erosion pins and profilometer techniques (Ebro Basin,NE-Spain)

In badland areas of the Ebro Basin, in a semiarid climate, two erosion plots (257 m²; 5° slope and 128 m²; 23° slope) on exposed Tertiary clays were monitored over two years (Nov. 1991–Nov. 1993). This material is characterized by high sodium absorption ratios which lead to high soil dispersivity. The dominant erosion processes in both plots are rilling and sheet erosion. Rainfall intensity was recorded at a weather station, connected to a data-logger, sediment production for single events was collected in tanks, and ground lowering was measured every six months by erosion pins and microtopographic profile gauge techniques. Significant runoff was produced only by rainfall events above 5 mm. Another threshold at 20 mm rain was noted. For rainfalls higher than 20 mm, the 23° slope plot shows a greater runoff response than the 5° one. Rainfall events exceeding this threshold showed a higher sediment production for the steeper slope. In the relationship between precipitation and sediment concentration, an envelope curve can be drawn indicating that any rainfall event of a given amount and intensity has a maximum sediment concentration which we speculate to be a function of the runoff sediment transport capacity. Runoff response and sediment yield in the studied plots are controlled by the rainfall and soil characteristics and their seasonal variations. In both plots, the erosion pins show that erosion rates in rill areas are 25–50% higher than in the interrill areas. Sediment yield recorded by collector devices was higher than the rates measured by erosion pins. The erosion rates based on rill cross-sections by profilometers were higher than the ones recorded by collectors.     

From water to tillage erosion dominated landform evolution

While water and wind erosion are still considered to be the dominant soil erosion processes on agricultural land, there is growing recognition that tillage erosion plays an important role in the redistribution of soil on agricultural land. In this study, we examined soil redistribution rates and patterns for an agricultural field in the Belgian loess belt. ¹³⁷Cs derived soil erosion rates have been confronted with historical patterns of soil erosion based on soil profile truncation. This allowed an assessment of historical and contemporary landform evolution on agricultural land and its interpretation in relation to the dominant geomorphic process. The results clearly show that an important shift in the relative contribution of tillage and water erosion to total soil redistribution on agricultural land has occurred during recent decades. Historical soil redistribution is dominated by high losses on steep midslope positions and concavities as a result of water erosion, leading to landscape incision and steepening of the topography. In contrast, contemporary soil redistribution is dominated by high losses on convex upperslopes and infilling of slope and valley concavities as a result of tillage, resulting in topographic flattening. This shift must be attributed to the increased mechanization of agriculture during recent decades. This study shows that the typical topographical dependency of soil redistribution processes and their spatial interactions must be accounted for when assessing landform and soil profile evolution.     

Effectiveness of grade-control structures in reducing erosion along incised river channels: the case of Hotophia Creek, Mississippi

Herein, we undertake a geomorphological analysis in which spatial and temporal trends of bed and bank erosion along an 18-km length of Hotophia Creek, Mississippi, are estimated for the period between 1961 and 2050. The evaluation was undertaken for two scenarios of channel response to channelization during 1961–1963. One scenario represents the ‘actual’ response of the channel and includes the effects of installing a series of grade-control structures (GCS) between 1980 and 1996, while the other represents a hypothetical scenario in which the channel is left to adjust naturally. This allows the effectiveness of GCS in reducing in-channel erosion to be assessed. The analysis relies on the availability of channel survey data to develop empirical bed and bank response models for each adjustment scenario, supplemented by bank stability modelling to predict future rates of bank erosion. Results indicate that channel erosion rates decline nonlinearly with respect to time since 1961, for both adjustment scenarios. However, by the year 2050, the “with” GCS adjustment scenario results in the cumulative removal of some 663,000 (9%) extra tonnes of sediment relative to the “without” GCS scenario. Most (63%) of this excess is derived from enhanced bed erosion during 1976–1985 and 1985–1992, with the remainder derived from increased bank erosion during 1985–1992. Detailed analysis of the patterns of erosion and deposition, and their association with the GCS, provides evidence to support the view that GCS installed along Hotophia Creek have, for the most part, been ineffective in reducing channel erosion rates. This is because the GCS were installed too late to prevent bed degradation, caused by the 1961–1963 channelization, migrating upstream. In addition, some structures have disrupted the downstream transmission of bed material from eroded reaches upstream, exacerbating bed degradation and bank erosion in incised reaches downstream.     

Modeling watershed-scale sequestration of soil organic carbon for carbon credit programs

Impending risks associated with climate change have forced the global community to devise tradable pollution permit or “cap and trade” approaches to control the release of greenhouse gases. In the U.S, soils have the potential to offset about 10 percent of annual CO₂ emissions; however, if carbon credits are to be included in greenhouse gas control programs, soil organic carbon (SOC) sequestration rates associated with agricultural land uses must be computed at a watershed scale. The Soil Water Assessment Tool (SWAT) water quality model, the Water Erosion Prediction Project (WEPP) erosion model, and the CENTURY 4.0 a soil carbon model were used to simulate carbon sequestration rates for 160 crop-tillage rotations in 272 sub-basins of the Big Creek watershed (12,300 hectares). Under annual crops, only no-till in a corn-soybean rotation, on low to moderate slopes results in net gains in SOC. Substantial annual rates of SOC sequestration occur only under perennial crops such as Conservation Reserve Program (CRP; 0.14 t/ha without erosion; 0.08 with erosion), pasture (0.67 t/ha without erosion; 0.58 with erosion), hay (0.88 t/ha without erosion; 0.52 with erosion), and forest (2.66 t/ha without erosion; 2.49 with erosion). Erosion thus has a large effect on the spatial distribution of field-measured SOC by moving it down slope and increasing its spatial variability. Because of this, carbon credit programs should be based on field practices, thus targeting the locations where the sequestration of atmospheric carbon actually occurs and minimizing monitoring costs. Developing model-based estimates of SOC sequestration rates of field practices at many locations would thus greatly serve the needs of carbon crediting programs.