Surface Altitude Research Articles

The Effects of Bedrock Topography and Soil Permeability on Saturated Zone Distribution in a Mountainous Steep‐Slope Area

ABSTRACTSaturation development and distribution at the soil–bedrock interface are important for predicting shallow landslide occurrence. Previous studies have indicated that saturation is generated in bedrock depressions and valleys and that bedrock groundwater seepage generates locally saturated areas. However, the effects of soil permeability, which is known to be heterogeneously distributed, on saturation development and distribution are poorly understood. In this study, we performed unprecedented high‐resolution (approximately 50 cm grid) soil pore water pressure and soil temperature monitoring using 141 tensiometer–thermocouple sets in a plot measuring approximately 5 × 4 m to investigate the effects of topography and bedrock groundwater seepage on saturation development and distribution. We then measured permeability distribution of two soil profiles, including at the soil–bedrock interface, using the Guelph Permeameter method (GP method) for comparison with saturated zone distribution and saturation duration. The results indicated that a perennial saturated area was formed by bedrock groundwater seepage and was distributed downstream from a certain bedrock surface altitude in the lower region of the study plot. After a peak of rainfall, the perennial saturated area expanded upslope owing to the increased seepage. In areas without the influence of bedrock groundwater, saturation was observed to retreat rapidly at high permeability points and persist over long periods at low permeability points; however, the saturation duration was inconsistent with the bedrock surface topography. Therefore, it is suggested that the bedrock altitude controls the saturation distribution generated by bedrock groundwater, whereas the distribution of saturation that is associated with direct rainwater infiltration may be controlled by the permeability distribution during recession periods. Although the plot size was small, the unprecedented high‐resolution observations suggest that the permeability distribution, rather than the bedrock topography, may control the saturated zone distribution following rainfall.

Suspended Sediment Concentration Estimation along Turbid Water Outflow Using a Multispectral Camera on an Unmanned Aerial Vehicle

Optical remote sensing using unmanned aerial vehicles (UAVs) is proposed to monitor changes in marine environments effectively. Optical measurements were performed using a UAV multispectral camera (RedEdge, five spectral wavelengths of 475, 560, 668, 717, and 842 nm) with high spatial (5 cm) and temporal (1 s) resolutions to monitor the rapidly changing suspended sediment concentration (SSC) in the Saemangeum coastal area on the western coast of Korea. To develop the SSC algorithm, optical field, and water sample measurements were obtained from outside (11 stations) and inside (three stations) regions separated by a seawall, accounting for 100 measurements from 2018 to 2020. Accordingly, the remote sensing reflectance (Rrs) was estimated at each sampling station and used to develop the SSC algorithm based on multiple linear regression. The algorithm reasonably estimated the SSC with an R2 and root mean square error of 0.83 and 4.27 (mg L−1), respectively. Continuous individual UAV measurements over the coastal area of Saemangeum were combined to generate a wider SSC map. For the UAV observational data, the atmospheric influence at each altitude was reduced to the surface altitude level using a relative atmospheric correction technique. The SSC map enabled front monitoring of SSC fluctuations caused by discharge water due to the sluice gate opening. These results demonstrated the usability of the UAV-based SSC algorithm and confirmed the possibility of monitoring rapid SSC fluctuations.

A Possible Cause for Preference of Super Bolt Lightning Over the Mediterranean Sea and the Altiplano

AbstractExceptionally high‐energy lightning strokes >106 J (X1000 stronger than average) in the very low‐frequency band between 5 and 18 kHz, also known as superbolts (SB), occur mostly during winter over the North‐East Atlantic, the Mediterranean Sea, and over the Altiplano in South America. Here we compare the World‐Wide Lightning Location Network database with meteorological and aerosol data to examine the causes of lightning stroke high energies. Our results show that the energy per stroke increases sharply as the distance between the cloud's charging zone (where the cloud electrification occurs) and the surface decreases. Since the charging zone occurs above the 0°C isotherm, this distance is shorter when the 0°C isotherm is closer to the surface. This occurs either due to cold air mass over the ocean during winter or high surface altitude in the Altiplano during summer thunderstorms. Stroke energy decreases with the warm phase of the cloud, as proxied by the cloud base temperature, and increases with a more developed cloud, as proxied by the cloud top temperature, but to a much lesser extent than the distance between the surface and 0°C isotherm. Aerosols play no significant role. It is hypothesized that a shorter distance between the charging zone and the ground represents less electrical resistance that allows stronger discharge currents.

The Mass Balance of the Kozelsky Glacier in Kamchatka for 1977–2022

Received May 26, 2023; revised June 6, 2023; accepted June 27, 2023The change in the volume of the Kozelsky Glacier in Kamchatka for the period 1977–2022 (1977–2015 and 2015–2022) was estimated using historical data and modern DEM. During this period, the area of the glacier did not change much. At the same time, its length increased by about 0.7 km, while the width decreased over its almost whole extent. The volume of the glacier decreased by 34.15 ± 6.74 million m3, and its surface became lower by 17.30 m, on the average. The cumulative mass balance amounted 14.70 ± 3.94 m w.e., and the mean annual value –0.33 m w.e. yr–1. In the last 45 years, the ice loss and redistribution to lower hypsometric levels took place on the Kozelsky Glacier. In 1977–2015, the average area change in the altitude of the glacier surface was equal to –17.84 m, the volume decreased by 35.21 ± 7.20 million m3, the cumulative mass balance amounted –15.16 ± 4.17 m w.e., and the mean annual balance –0.40 m w.e. yr–1. In the period 2015–2022, an elevation of the glacier surface was recorded by 0.59 ± 1.55 m on the average, the volume increased by 1.01 ± 2.65 million m3, the cumulative mass balance amounted to 0.50 ± 1.35 m w.e., and the mean annual balance – to 0.07 m w.e. yr–1. During the last decade, a slowdown in the movement of the glacier front down the valley was recorded. In 2012–2022, the glacier front advanced with a velocity of about 5.2 m/year, while it was 17.9 m/year in 1977–2007, and 20.0 m/year in 2007–2012. The current climatic conditions are not favorable for development of glaciers. In 1977–2022, a trend of the summer air temperature rise was observed with a relatively stable amount of precipitation falling during the cold period. The almost continuous (except 1978–1981) advance of the glacier in 1977–2022 can be explained by the influence of the volcanic factor. A thick surface moraine covers more than 2/3 of the glacier area and, thus, prevents the surface ablation. Increased seismic activity associated with active volcanism promotes the ice movement.

Effects of urbanization on the structure of plant-flower visitor network at the local and landscape levels in the northern Argentinian Yungas forest

Human population and cities are growing fast, with a concomitant modification of the land surface. Urbanization is driving biodiversity loss and biological homogenization, which impacts human wellbeing. In this study, we evaluated the influence of urbanization on flower visitor assemblage using an interaction network approach. We assessed the effect of different variables at the local and landscape scales on community parameters and network metrics along a gradient of urbanization located in a subtropical montane Yungas forest. We found that local variables affected the richness of flower visitors, which increased with greater flower coverage, high stability of floral resources, and the proportion of exotic plants. Moreover, local variables affected the diversity, nestedness (NODF), and specialization (H2) of the interaction network. Landscape variables, such as altitude and proportion of impervious surface (a proxy of urbanization), affected both the richness of flower visitors and specialization. The effect of urbanization on the richness of flower visitors differed across the altitudinal gradient, with higher impact at higher altitudes. In conclusion, our results indicate that local and landscape variables affect community parameters and the structure of plant-flower visitor networks to different extents and strengths.

An Improved Intersection Method for Positioning and Its Application to Dalk Glacier, East Antarctica

Stake measurement is a significant ground-observation method used to monitor glacier dynamics in the cryosphere and validate the reliability of remotely sensed investigations. However, some stake observation records obtained by total stations have not been efficiently utilized via the classic forward intersection method. The forward intersection method requires a pair of cognominal observation records and cannot function if stakes are only observed from one control point at a stake measurement, which commonly occurs in measurements. Therefore, we propose an improved intersection method that combines stake trajectories and sightlines, the trajectory-sightline intersection method, which requires only a single observation record. In the processing of stake identification, location and elevation from the trajectory-sightline intersection method are constrained by the velocity and surface altitude of the glacier. Then, the trajectory-sightline intersection method is applied to the stake measurement of Dalk Glacier in East Antarctica. The results show that the trajectory-sightline intersection method extracted an additional 44 intersections, thus increasing the data utilization ratio by 20.35%. Moreover, leave-p-out cross-validation indicated that the trajectory-sightline intersection method, with an average offset of 0.88±0.54 m in the horizontal position and 0.26±0.12 m at the elevation, has similar accuracy to the differential global positioning system (DGPS) applied to Livingston Island. The ice flow velocities of Dalk Glacier calculated from trajectory-sightline intersections agree with those from forward intersections (the mean velocity difference is 0.47 m·year−1). Therefore, the trajectory-sightline intersection method is an effective method for improving data utilization and obtaining stake locations and glacier velocities.

Enhanced Interannual Variability in Temperature during the Last Glacial Maximum

Abstract Using all relevant climate experiments archived in phases 3 and 4 of the Paleoclimate Modeling Intercomparison Project (PMIP3/4), we examine the interannual variability change in global-scale surface air temperature and associated physical mechanisms during the Last Glacial Maximum (LGM). The results show that relative to the preindustrial period, the LGM interannual temperature variability increased by 20% globally, which was mainly attributed to the large-scale increase in the meridional temperature gradient, especially at midlatitudes. Larger magnitudes of change occurred in areas where the underlying surface properties, such as the surface altitude, land–sea distribution, and ice sheet extent, differed from those in the preindustrial period. In addition, the relationship between changes in the meridional temperature gradient and the interannual temperature variability became closer in the winter hemisphere. In the tropical land regions, changes in interannual temperature variability are mainly related to the adjustment of latent and sensible heat fluxes during the LGM.

Spatio-temporal variations and ecological risks of organochlorine pesticides in surface waters of a plateau lake in China.

Qilu Lake is one of the 9 plateau lakes in Yunnan, China, with a lake surface altitude of 1796.62m above sea level. In spite of the importance and agriculturally-intensive phenomenon in Qilu Lake Basin, few studies have provided a modern evaluation of pesticide residues and potential effects to local aquatic organisms. The primary goal of this study was to determine the spatio-temporal variations of organochlorine pesticides (OCPs) in this area, and to further assess the related ecological risks. Of the 25 OCPs analyzed, 14 were detected, and the concentrations of ∑25OCPs were highest in the upstream of rivers, followed by regions close to the lake shore, and the lowest concentrations were found in Qilu Lake in every season except winter. The concentrations of ∑25OCPs were the highest in summer, and the lowest in winter. OCP concentrations in spring and in autumn were similar. 4,4'-DDD, γ-HCH, HCB, trans-chlordane, and cis-chlordane were 5 OCPs with relatively high risk in Qilu Lake Basin. Interestingly, higher OCP concentrations do not necessarily correspond to higher ecological risk levels. Low predicted no-effect concentration (PNEC) values and relatively high toxicity of these OCPs led to their high risk quotient (RQ) values. This work further illustrated that although OCPs have been banned for many years, they were still frequently detected in surface waters, and caused risks to aquatic animals.

Rivers' Water Level Assessment Using UAV Photogrammetry and RANSAC Method and the Analysis of Sensitivity to Uncertainty Sources.

Water-level monitoring systems are fundamental for flood warnings, disaster risk assessment and the periodical analysis of the state of reservoirs. Many advantages can be obtained by performing such investigations without the need for field measurements. In this paper, a specific method for the evaluation of the water level was developed using photogrammetry that is derived from images that were recorded by unmanned aerial vehicles (UAVs). A dense point cloud was retrieved and the plane that better fits the river water surface was found by the use of the random sample consensus (RANSAC) method. A reference point of a known altitude within the image was then exploited in order to compute the distance between it and the fitted plane, in order to monitor the altitude of the free surface of the river. This paper further aims to perform a critical analysis of the sensitivity of these photogrammetric techniques for river water level determination, starting from the effects that are highlighted by the state of the art, such as random noise that is related to the image data quality, reflections and process parameters. In this work, the influences of the plane depth and number of iterations have been investigated, showing that in correspondence to the optimal plane depth (0.5 m) the error is not affected by the number of iterations.

Flood Detection Using Backpropagation Neural Network Method

Lack of river and watershed management will cause problems and disasters. One of it is the flood that can cause physical, social and economic loss. So countermeasures or flood anticipation are needed by using the Early Warning System (EWS) to provide early information if a flood is going to occur. This study uses five input indicators: temperature, humidity, water discharge, water surface altitude and rainfall data that will produce output in the form of notifications and alarms for the Early Warning System (EWS). Then the input and output data configuration will be processed using a Backpropagation Neural Network. Data used is data recorded in real-time on the research object for two weeks with the composition of training and testing data with a percentage of 80% and 20%. The best backpropagation neural network model used has the input of 5 neurons layer architecture, 15 neurons as the hidden layer and three neurons as the output layer. The flood prediction result uses the Backpropagation Neural Network method, has an RMSE score performance of 2.16e-21 and a percentage success testing system of 91.33%. It shows that the model has an excellent accuracy level.

Effects of Charged Martian Dust on Martian Atmosphere Remote Sensing

In this paper, the extinction property and optical depth of charged Martian dust at infrared band 3 THz–300 THz are studied using the Mie scattering theory. It is found that the extinction coefficients of Martian atmospheric dust and the dust optical depth (DOD) of the Martian atmosphere can be amplified significantly as the dust particles are charged. This extinction amplification has a peak, called amplification resonance, which shifts toward the upper left of the r-q parameter plane with increasing frequency. Here, r denotes the particle radius and q denotes the particle’s total net charge. The amplification of the Martian DOD is more significant at high altitudes than at low altitudes because the particles at high altitudes are smaller. For example, at an altitude of 30–50 km, the dust optical depth at 30 THz can be increased by 60–200%. However, at 3 THz–10 THz, the DOD at the near surface altitude (0–10 km) can still be enhanced by ~80%. This implies that by treating the Martian dust as uncharged particles, the dust density constructed from the Martian DOD data might be overestimated. The estimation error of the dust density of the Martin atmosphere may be reduced by counting the enhancement of the DOD that is caused by charged dust.

Global validation of clear-sky models for retrieving land-surface downward longwave radiation from MODIS data

Surface longwave radiation (SLR) plays an important role in the energy budget of the Earth's climate system. Remote sensing provides various data sources to retrieve SLR on a large scale and with high spatial resolution (e.g., 1 km). Multiple retrieval methods of surface downward longwave radiation (SDLR) based on satellite thermal infrared data produce different retrieval results for the same scenario. Therefore, validating these models is necessary to understand their characteristics and limitations. To this end, ground-based measurements were used to provide independent validation of six widely used SDLR models for clear sky conditions over 41 Baseline Surface Radiation Network (BSRN) stations worldwide. The Wang2020 model had the best overall performance (bias of −5.480 W/m2, root-mean-square errors [RMSE] of 23.226 W/m2, R2 of 0.879), and Tang2008 model had similar retrieval capability. The errors of LST had limited influence on the retrieval accuracy of SDLR models. When using the near-surface air temperature, the retrieval accuracy of the Zhou2007 model was significantly improved with a range of ~9.5 W/m2 for the RMSE. The uncertainty of TCWV had significant effect on all model performances, wherein the Zhou2007 model had stronger error resilience of TCWV. Moreover, MODIS TIR TCWV data provided better performance than NIR TCWV in most situations, and thereby are preferred to use in the SDLR retrieval. Surface altitude had a lesser impact on SDLR retrieval than terrain effects. All models overestimated SDLR for peak stations in mountainous areas, with biases reaching 56.614 W/m2 and RMSE reaching 63.909 W/m2. Land cover type also had a significant effect on retrieval accuracy; model performances were poorer in the desert and barren where atmospheric conditions are extremely dry and hot. Remote sensing SDLR data with high accuracy are needed for hydrological, agricultural, and climate change applications. The results of this study provide a reference for the SDLR retrieval accuracy based on clear-sky models.

Mercury in the sediments of selected peatlands in Małopolska region

The mercury content in peat layers of seven peat bogs in the Małopolska region (Wolbrom in the Kraków-Częstochowa Upland and Otrębowskie Brzegi, Puścizna Wysoka, Puścizna Mała, Przybojec, Puścizna Długopole and Bór na Czerwonem in the Orawa-Nowy Targ Basin) was determined to a depth of 100 cm at a resolution of 1–2 cm. Geochemical background and level of anthropogenic mercury enrichment were determined for each peat bog. To establish the interdependence between mercury concentration, mineral matter content and selected metals, geochemical diagrams and models of the system were drawn up representing the correlational relationships between the selected geochemical variables for three of the peat bogs in the Orawa--Nowy Targ Basin (Otrębowskie Brzegi, Przybojec and Bór na Czerwonem). It was found that, comparing against other European raised bogs, the range of mercury concentrations in the peat layer down to 50 cm in the peatlands of Małopolska is similar to those in Norway, Sweden, Denmark and Spain, and decidedly lower than in the peatlands in Scotland, Belgium and Czechia. Furthermore, the sites in the Orawa-Nowy Targ Basin were clearly lower, the lower the altitude of the peat bog surface, which is probably associated with local topoclimatic conditions.

An Early Miocene Lowland on the Northeastern Tibetan Plateau

The northeastern Tibetan Plateau (NE TP) has long been thought to be the last part of the Plateau to be raised, but this assumption has been challenged by recent analyses of fossil leaf energy, which have pointed to the possibility that the present surface altitude of ∼3,000 m above sea level (asl) in the Qaidam Basin (QB) was attained during the Oligocene. Here, for the first time, we present a record of glycerol dialkyl glycerol tetraethers (GDGTs) from a well-dated Cenozoic section in the QB. This record appears to demonstrate that the mean annual average paleotemperature of the QB was 28.4 ± 2.9°C at ∼18.0 Ma. This would suggest that the paleoelevation of the QB was only ∼1,488 m asl at that time and that a ∼1,500 m uplift was attained afterwards, in agreement with the massive shortening of the QB and the rapid drying of inland Asia since the late Miocene.

Bionic Nonsmooth Drag Reduction Mathematical Model Construction and Subsoiling Verification.

In this study, a bionic nonsmooth drag-reducing surface design method was proposed; a mathematical model was developed to obtain the relationship between the altitude of the nonsmooth drag-reducing surface bulges and the spacing of two bulges, as well as the speed of movement, based on which two subsoiler shovel tips were designed and verified on field experiments. The mechanism of nonsmooth surface drag reduction in soil was analyzed, inspired by the efficient digging patterns of antlions. The nonsmooth surface morphology of the antlion was acquired by scanning electron microscopy, and a movement model of the nonsmooth surface in soil was developed, deriving that the altitude of the nonsmooth drag-reducing surface bulge is proportional to the square of the distance between two bulges and inversely proportional to the square of the movement speed. A flat subsoiler shovel tip and a curved tip were designed by applying this model, and the smooth subsoiler shovel tips and the pangolin scale bionic tips were used as controls, respectively. The effect of the model-designed subsoilers on drag reduction was verified by subsoiling experiments in the field. The results showed that the resistance of the model-designed curved subsoiler was the lowest, the resistance of the pangolin scale bionic subsoiler was moderate, and the resistance of the smooth surface subsoiler was the highest; the resistance of the curved subsoiler was less than the flat subsoilers; the resistance reduction rate of the model-designed curved subsoiler was 24.6% to 33.7% at different depths. The nonsmooth drag reduction model established in this study can be applied not only to the design of subsoilers but also to the design of nonsmooth drag reduction surfaces of other soil contacting parts.

Atmospheric Correction of True-Color RGB Imagery with Limb Area-Blending Based on 6S and Satellite Image Enhancement Techniques Using Geo-Kompsat-2A Advanced Meteorological Imager Data

This study aims for producing high-quality true-color red-green-blue (RGB) imagery that is useful for interpreting various environmental phenomena, particularly for GK2A. Here we deal with an issue that general atmospheric correction methods for RGB imagery might be breakdown at high solar/viewing zenith angle of GK2A due to erroneous atmospheric path lengths. Additionally, there is another issue about the green band of GK2A of which centroid wavelength (510 nm) is different from that of natural green band (555 nm), resulting in the unrealistic RGB imagery. To overcome those weakness of the RGB imagery for GK2A, we apply the second simulation of the satellite signal in the solar spectrum radiative transfer model look-up table with improved information considering altitude of the reflective surface to reduce the exaggerated atmospheric correction, and a blending technique that mixed the true-color imagery before and after atmospheric correction which produced a naturally expressed true-color image. Consequently, the root mean square error decreased by 0.1–0.5 in accordance with the solar and view zenith angles. The green band signal was modified by combining it with a veggie band to form hybrid green which adjust centroid wavelength of approximately 550 nm. The original composite of true-color RGB imagery is dark; therefore, to brighten the imagery, histogram equalization is conducted to flatten the color distribution. High-temporal-resolution true-color imagery from the GK2A AMI have significant potential to provide scientists and forecasters as a tools to visualize the changing Earth and also expected to intuitively understand the atmospheric phenomenon to the general public.

Response of terrace deposit thickness to climate change and tectonic deformation: An example of the Liyuan River in the Northeast Tibetan Plateau

AbstractThe relationships between tectonics, climate and sedimentation are important for understanding landscape evolution. In order to elucidate these relationships, we conducted a detailed survey of the thickness of terrace deposits and the altitude of terrace surfaces along the Liyuan River valley, in the Northeast (NE) Tibetan Plateau. We have dated the formation ages of terraces T3 and T4 to 42 ± 3.2 and 93.7 ± 7.5 ka, respectively, using optically stimulated luminescence (OSL) method. By compiling the deposit thickness of terraces, palaeoclimate records and local tectonism, we found that variations in precipitation account for the temporal differences in thickness/rates of sediment deposition in terraces T4 and T3. However, the spatial distribution of terrace deposit thickness was determined by tectonism, including faulting and folding, which changed the local base level. In order to maintain a foreland‐ward‐dipping gradient, the Liyuan River adjusts its sedimentary processes via spatially differentiated aggradation.

3D MODELING AND GIS ANALYSIS FOR AERODROME FOREST OBSTACLE MONITORING

Abstract. The article discusses methods for constructing and using digital photogrammetric and cartographic models as a basis for growing tree height control and plantation planning in aerodrome areas. Forests or gardens in the take-off and landing flight areas, exceeding special limitation surfaces, are dangerous obstacles and intended to cut down. Tree and bush vegetation should be under periodic monitoring because of their growth. The research was aimed to determine the maximum permissible obstacle height and tree age when it reaches the obstacle limitation surface altitude. For these purposes, it is proposed to use geospatial modeling and geoinformation analysis methods. As a basis for geospatial models, remote sensing optical stereo images were used. The allowable height is calculated as a difference between 3D obstacle limitation surface and the earth surface altitude values. The article presents the study results for a Belarus climatic zone, where the tree species predictive age in reaching the maximum permissible height is calculated. The main goal of the technology is to manage the aerodrome forest plantation growth without further labor-intensive monitoring, while ensuring the safety of aircraft flights.

POSITIVE AND NEGATIVE ROUGHNESS ACCORDING TO LOCAL DIFFERENCES BETWEEN DEM SURFACE AND 3D REFERENCE PLANES

Abstract. The irregularities of the earth’s surface are quantified by means of roughness measurements using Digital Elevation Models (DEM’s). This article presents a roughness measurement method that is based on the calculation of the difference of altitude existing between a plane passing through the centre of a moving window and the altitude of the DEM surface inside this window. This method differs from the measure of the standard deviation and best fit plane, in the sense that it considers all difference values, positives or negatives. The measurement is done in a 3 × 3 or a 5 × 5 moving window and contemplates inside this window the plane which passes through the centre of the window and the highest pixel located in the border or perimeter of this window. According to the 3D configuration of the DEM surface inside the moving window, the sum of all the differences is positive or negative, allowing to discriminate the local morphology independently of the global roughness. The roughness variable which distinguishes negative and positive values allows to classify accurately landscape units such as watersheds, riverbeds, volcanic assemblages as well as landforms associated with tectonic structures.

Tools for Comprehensive Assessment of Fluid-Mediated and Solid-State Alteration of Carbonates Used to Reconstruct Ancient Elevation and Environments

Carbonates are ubiquitous in the rock record and provide a broad array of stable isotope-based paleoclimatic proxies (i.e., δ18O, δ13C, ∆17O, ∆47, ∆48) that provide information on stratigraphy, carbon cycling, temperature, hydrology, and the altitude of ancient land surfaces. Thus, carbonates are an essential archive of environmental and topographic histories of continental terranes. However, carbonate minerals are highly susceptible to post-depositional alteration of primary isotopic values via fluid-mediated and solid-state reactions. We propose a hierarchical suite of techniques to comprehensively assess alteration in carbonates, from essential and readily accessible tools to novel, high-resolution techniques. This framework provides a means of identifying preserved textures in differentially altered samples that contain high-value environmental information. To illustrate this progressive approach, we present a case study of Tethyan nearshore carbonates from the Paleocene Tso Jianding Group (Tibet). We demonstrate the utility of each technique in identifying chemical and crystallographic indicators of post-depositional alteration at progressively finer spatial scales. For example, secondary ionization mass spectrometry (SIMS) oxygen isotope maps of micrite and bioclasts reveal significant isotopic heterogeneity due to grain-scale water-rock exchange in textures that were labeled “primary” by optical inspection at coarser spatial resolution. Optical and cathodoluminescence microscopy should be the minimum required assessment of carbonate samples used in stable isotope analyses, but supplemented when necessary by SIMS, PIC mapping, and other yet untapped technologies that may allow distinction of primary and altered fabrics at finer spatial resolutions.