Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Congo Basin Water Balance and Terrestrial Fluxes Inferred From Satellite Observations of the Isotopic Composition of Water Vapor

Large spatio-temporal gradients in the Congo basin vegetation and rainfall are observed. However, its water-balance (evapotranspiration minus precipitation, or ET - P) is typically measured at basin-scales, limited primarily by river-discharge data, spatial resolution of terrestrial water storage measurements, and poorly constrained ET. We use observations of the isotopic composition of water vapor to quantify the spatio-temporal variability of net surface water fluxes across the Congo Basin between 2003 and 2018. These data are calibrated at basin scale using satellite gravity and total Congo river discharge measurements and then used to estimate time-varying ET - P over four quadrants representing the Congo Basin, providing first estimates of this kind for the region. We find that the multi-year record, seasonality, and interannual variability of ET - P from both the isotopes and the gravity/river discharge based estimates are consistent. Additionally, we use precipitation and gravity-based estimates with our water vapor isotope-based ET - P to calculate time and space averaged ET and net river discharge within the Congo Basin. These quadrant-scale moisture flux estimates indicate (a) substantial recycling of moisture in the Congo Basin (temporally and spatially averaged ET/P > 70%), consistent with models and visible light-based ET estimates, and (b) net river outflow is largest in the Western Congo where there are more rivers and higher flow rates. Our results confirm the importance of ET in modulating the Congo water cycle relative to other water sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations

Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Shortwave Cloud-SST Feedback Amplifies Multi-Decadal Pacific Sea Surface Temperature Trends: Implications for Observed Cooling

Climate models struggle to produce sea surface temperature (SST) gradient trends in the tropical Pacific comparable to those seen recently in nature. Here, we find that the magnitude of the cloud-SST feedback in the subtropical Southeast Pacific is correlated across models with the magnitude of Eastern Pacific multi-decadal SST variability. A heat-budget analysis reveals coupling between cloud-radiative effects, circulation, and SST gradients in driving multi-decadal variability in the Eastern Pacific. Using this relationship and observed feedback estimates, we find that internal Eastern Pacific SST variability is underestimated in most models. Adjusting for model bias increases the likelihood of generating a cooling trend at least as large as observations in preindustrial control simulations by ~56% on average. If models underestimate climate “noise,” as our results suggest, this bias should be accounted for when attributing the relative importance of forced versus unforced changes in the climate.

58 GEOSCIENCES↗
DOE OSTI2024

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimating the Electric Fields Driving Lightning Dart Leader Development With BIMAP‐3D Observations

Abstract In this paper, a numerical dart leader model has been implemented to understand the leader's development and the corresponding electric field changes observed by the 3D broadband mapping and polarization (BIMAP‐3D) system. The model assumes the extending leader channel is equipotential and has a linear charge distribution induced by an ambient electric field. The charge distribution induced by the ambient field can be used to model the electric field change at the ground. We then find the ambient electric field which best fits the field change measurements at the two BIMAP stations. The estimated ambient electric field decreases in the direction of dart leader propagation. Our observations and modeling results are consistent with our earlier hypothesis that dart leader speed is proportional to the electric field at the leader tip. The model also supports our earlier analysis that leader speed variations near branch junctions were due to previous charge deposits near the junctions. The modeled tip electric field is generally lower than the breakdown field unless the pre‐dart‐leader channel has a significant temperature of 3,000 K. This is consistent with the fact that dart leaders typically do not form new branches into the virgin air. Furthermore, the tip field is generally close to the negative streamer stability field at ambient temperatures, explaining the nature of the narrow and well‐defined channel structure. In addition to the charge distribution and the ambient and tip electric field, the development of the channel potential and current distribution are also presented.

3D lightning mapping↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimating Interplanetary Magnetic Field Conditions at Mercury's Orbit From MESSENGER Magnetosheath Observations Using a Feedforward Neural Network

Abstract Mercury's small magnetosphere is embedded in the dynamic and intense solar wind environment characteristic of the inner heliosphere. Both the magnitude and orientation of the interplanetary magnetic field (IMF) significantly influence the solar wind‐magnetospheric interaction at Mercury, driving phenomena such as magnetic reconnection. The MErcury Surface, Space Environment, Geochemistry and Ranging (MESSENGER) spacecraft provided in‐situ magnetic field measurements of the solar wind, the magnetosheath, and the magnetosphere along each orbit. However, it is a challenge to directly assess the IMF's impact on Mercury's plasma environment due to the temporal separation between observations within the solar wind and the magnetosphere, especially in the absence of an upstream monitor. Here, we present a feedforward neural network (FNN) trained on a subset of magnetosheath observations to estimate the strength and orientation of the IMF upstream of the bow shock. Utilizing magnetosheath magnetic field, cylindrical spatial coordinates, and heliocentric distance measurements, the FNN predicts upstream IMF conditions with an score of 0.70 and mean averaged error of 5.3 nT, thereby greatly decreasing the temporal separation between IMF estimates and magnetospheric measurements throughout the MESSENGER mission. This approach yields IMF estimates for all magnetosheath data measured by MESSENGER, providing a useful tool for future investigations of the IMF impact on Mercury's magnetosphere. This method will be integrable with the dual‐spacecraft BepiColombo magnetosheath measurements, providing useful estimates of upstream IMF conditions particularly during the extended periods in which neither spacecraft sample the solar wind. Our results demonstrate the utility of machine learning techniques on advancing space science research.

Bowers, Charles F.↗
DOE OSTI2025

Direct Observations of Solute Dispersion in Rocks With Distinct Degree of Sub‐Micron Porosity

Abstract The transport of chemical species in rocks is affected by their structural heterogeneity to yield a wide spectrum of local solute concentrations. To quantify such imperfect mixing, advanced methodologies are needed that augment the traditional breakthrough curve analysis by probing solute concentration within the fluids locally. Here, we demonstrate the application of asynchronous, multimodality imaging by X‐ray computed tomography (XCT) and positron emission tomography (PET) to the study of passive tracer experiments in laboratory rock cores. The four‐dimensional concentration maps measured by PET reveal specific signatures of the transport process, which we have quantified using fundamental measures of mixing and spreading. We observe that the extent of solute spreading correlate strongly with the strength of subcore‐scale porosity heterogeneity measured by XCT, while dilution is enhanced in rocks containing substantial sub‐micron porosity. We observe that the analysis of different metrics is necessary, as they can differ in their sensitivity to the strength and forms of heterogeneity. The multimodality imaging approach is uniquely suited to probe the fundamental difference between spreading and mixing in heterogeneous media. We propose that when multi‐dimensional data is available, mixing and spreading can be independently quantified using the same metric. We also demonstrate that one‐dimensional transport models have limited predictive ability toward the internal evolution of the solute concentration, when the model is solely calibrated against the effluent breakthrough curves. The data set generated in this study can be used to build realistic digital rock models and to benchmark transport simulations that account deterministically for rock property heterogeneity.

Kurotori, Takeshi [Department of Chemical Engineer↗
DOE OSTI2025

Opportunistic Short‐Term Water Uptake Dynamics by Subalpine Trees Observed via In Situ Water Isotope Measurements

Abstract Variations in tree water sources are important to understand in semi‐arid ecosystems because climatic shifts towards lower snowpack and increased drought affect water availability in subalpine forests of the western US. Here, we use daily in situ measurements of stable isotopes ( 2 H & 18 O) in soil and tree stem water, soil matric potential and sap flow to study tree water uptake dynamics. We instrumented three soil profiles down to 90 cm, as well as three aspen and engelmann spruce trees near Gothic, Colorado, in the East River watershed. We observed the fate of natural isotopic variations in rainfall, soil, and plants from June to October 2022, and in August 2023 we conducted a 2 H labeled irrigation experiment. Our observations showed that all studied aspen trees compensated for water scarcity in the shallow soil by shifting the dominant water source at 60(±20) cm to ⅔ of uptake from 90 cm within a few days of a dry period. Both species relied on snowmelt stored in the subsoil to sustain transpiration. Intense rainfall caused the plant water uptake to shift partially to top soil layers within 2 days. Spruce transpiration was lower and relied more on snowmelt, because rainfall infiltration was low in the spruce stand due to high canopy interception. Our findings highlight the important role of snowmelt stored in the deep soil layers for subalpine forest drought response and the dominant fate of monsoonal rainfall to become transpiration rather than recharging groundwater and streams in the Upper Colorado River. Plain Language Summary There is a need to understand how trees in mountainous regions respond to dry conditions that lead to water scarcity, because climate projections suggest that such conditions will become more frequent in the future. Here we present a novel data set of measurements of daily stable isotopes of water across soil profiles and in tree stems of aspen and spruce. Our data show that when the upper soil dried out, aspen trees shifted to using water from deeper layers (beneath 60 cm) to keep transpiring. For spruce trees the uptake pattern is less clear, but both types of trees mainly used snowmelt stored in the deeper soil layers to survive the dry summer. After heavy rain, aspen and spruce trees switched to using water from the top 20 cm of soil. However, for spruce, only some rain reached the soil because the dense tree canopy intercepted it, so spruce trees stayed more dependent on snowmelt and used less water overall. This study shows how important deep snowmelt water is for helping forests survive dry periods and suggests that most summer rain is quickly used by trees rather than replenishing streams and groundwater in the headwaters of the Colorado River. Key Points Tree water resources changed within a few days from snow dominated to higher share of rainfall as soils wetted up after a dry period Compensatory plant water uptake by aspen from the deep layer (90 cm), while uptake from soil depths that became drier (60 cm) declined Strong differences between water sources and availability beneath aspen and spruce, respectively

Sprenger, Matthias↗
DOE OSTI2025

A U.S. Scientific Community Vision for Sustained Earth Observations of Greenhouse Gases to Support Local to Global Action

Managing carbon stocks in the land, ocean, and atmosphere under changing climate requires a globally‐integrated view of carbon cycle processes at local and regional scales. The growing Earth Observation (EO) record is the backbone of this multi‐scale system, providing local information with discrete coverage from surface measurements and regional information at global scale from satellites. Carbon flux information, anchored by inverse estimates from spaceborne Greenhouse Gas (GHG) concentrations, provides an important top‐down view of carbon emissions and sinks, but currently lacks global continuity at assessment and management scales (<100 km). Partial‐column data can help separate signals in the boundary layer from the overlying atmosphere, providing an opportunity to enhance surface sensitivity and bring flux resolution down from that of column‐integrated data (100–500 km). Based on a workshop held in September 2024, the carbon cycle community envisions a carbon observation system leveraging GHG partial columns in the lower and upper troposphere to weave together information across scales from surface and satellite EO data, and integration of top‐down/bottom‐up analyses to link process understanding to global assessment.

Parazoo, Nicholas C. [California Institute of Tech↗
DOE OSTI2025

Multidecadal Fluctuations in the Observed ENSO‐Tropical Cyclone Teleconnection

Abstract El Niño‐Southern Oscillation (ENSO) is a skillful predictor for seasonal tropical cyclone (TC) activity in most TC basins. This study examines recent changes in the observed ENSO‐TC teleconnection strength, as measured by ENSO modulation of hurricane frequency. We find that the ENSO‐North Atlantic TC teleconnection fluctuated over time, with the strongest relationship occurring from the 1980s to the mid‐2000s. In the western and eastern North Pacific, the ENSO‐TC teleconnection has strengthened in recent decades. Periods with a strong ENSO‐TC teleconnection are associated with more favorable environmental conditions for TCs, with higher values of genesis potential indices. Positive phases of the Atlantic Multidecadal Oscillation coincided with periods of strong ENSO‐TC teleconnections in the Atlantic and North Pacific basins. A weaker Atlantic ENSO‐TC relationship was associated with negative phases of the Pacific Decadal Oscillation and the North Atlantic Oscillation. This research reveals climate conditions that modulate ENSO's utility for seasonal TC prediction. Plain Language Summary El Niño‐Southern Oscillation (ENSO) is a useful predictor for seasonal tropical cyclone (TC) activity in many basins. Here we found that the strength of the ENSO‐TC teleconnection, represented as the correlation between ENSO and the number of hurricanes and accumulated cyclone energy, has changed in the historical record. The ENSO‐TC teleconnection in the North Atlantic fluctuated over time, with a weak relationship during the 1960s and 1970s and a strong relationship during the 1980s to mid‐2000s. Meanwhile, the ENSO‐TC teleconnection strengthened in the North Pacific in recent decades, with strong teleconnections after the 1980s in the western North Pacific and after the 2000s in the eastern North Pacific. Periods of strong ENSO‐TC teleconnections are associated with more favorable environmental conditions for TCs, including higher values of genesis potential indices and higher mid‐tropospheric humidity, as well as positive phases of the Atlantic Multidecadal Oscillation. Additionally, the negative phase of the Pacific Decadal Oscillation leads to strong/weak ENSO‐TC teleconnections in the eastern North Pacific and North Atlantic, respectively. Furthermore, a negative North Atlantic Oscillation is associated with a weak ENSO‐North Atlantic TC teleconnection. This research highlights variations in ENSO's effectiveness for seasonal TC prediction. Key Points The observed impact of ENSO on tropical cyclone (TC) activity exhibits multidecadal fluctuations The ENSO‐TC teleconnection was strong in the Atlantic from the 1980s to mid‐2000s and strengthened over the North Pacific in recent decades The ENSO‐TC teleconnection is stronger in the Atlantic and North Pacific basins during a positive Atlantic Multidecadal Oscillation

ENSO↗
DOE OSTI2026

Long-term observations reveal non-negligible temperature contributions to moist static energy during the MJO transition from suppressed to active phases

The evolution of moist static energy (MSE) is widely used to understand the organization and propagation of the Madden-Julian oscillation (MJO). Past studies, largely based on reanalysis or short-term observations, have highlighted humidity as the dominant driver of MJO evolution. Using 14 years of bias-corrected radiosonde observations from the Department of Energy – Atmospheric Radiation Measurement (DOE-ARM) Facility site at Manus in the western Pacific, we show that temperature has a significant influence on the vertically-integrated MSE during MJO’s transition from suppressed to active phase. Specifically, during the transition phase in boreal winter (spring), ~48% (~31%) of the vertically-integrated MSE anomalies are attributed to temperature anomalies. This temperature contribution is confined to the mid- and upper troposphere, where temperature and MSE anomalies exhibit a strong correlation. These findings have important implications for MJO theory and provide new constraints for evaluating MJO transition processes in models.

Alex, Nirmal↗
DOE OSTI2026

Surface Observations From Atmospheric Radiation Measurement Sites Constrain the Anthropogenic Contribution to Cloud Droplet Number

Uncertainty in anthropogenic forcing driven by aerosol-cloud interactions (aci) limits our ability to infer the sensitivity of the Earth system to forcing from historical records. The driver of aci is the change in cloud droplet number concentration (N d ) due to changes in aerosol serving as cloud condensation nuclei (CCN). Here, we combine a perturbed parameter ensemble run in a global Earth system model with observations of CCN and single-layer-cloud N d at surface sites in the Azores, the Southern Great Plains, and Ascension Island to provide a constraint on the anthropogenic contribution to present-day N d . These observational lines of evidence constrain the preindustrial to present-day change in N d to be between 11 and 43 cm −3 . This is consistent with the upper end of some previous estimates but has a higher minimum perturbation, pointing to a stronger historical aerosol cooling.

ARM↗
DOE OSTI2026

First EarthCARE CPR Observations of Dynamics–Microphysics Coupling in Marine Cold‐Air Outbreak Clouds

Marine cold-air outbreak (MCAO) clouds begin as shallow cloud streets near the ice edge and evolve into broken open-cell convection downstream. Their rapid transitions in cloud structure and microphysics are closely linked to vertical motions, but this interplay has not previously been observed from space. The EarthCARE satellite, carrying the first spaceborne Doppler radar (Cloud Profiling Radar, CPR), now measures hydrometeor sedimentation velocity and vertical air motion. Using CPR observations over the Norwegian and Barents Seas (December 2024–May 2025), we analyze MCAO and non-MCAO clouds. We find that strong updrafts in MCAO clouds coincide with larger supercooled liquid water (SLW) path (LWP). These conditions are accompanied by increases in sedimentation velocity and its vertical gradient, reflecting rapid riming-driven ice growth. In contrast, non-MCAO clouds show lower riming efficiency at similar LWP due to relatively weaker dynamical support. This emphasizes the importance of dynamics in shaping polar cloud structure and microphysics.

Kim, Jiseob [McGill Univ., Montreal, QC (Canada)] ↗
DOE OSTI2020

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using as a model system ZnO, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Direct observation of elemental fluctuation and oxygen octahedral distortion-dependent charge distribution in high entropy oxides

The enhanced compositional flexibility to incorporate multiple-principal cations in high entropy oxides (HEOs) offers the opportunity to expand boundaries for accessible compositions and unconventional properties in oxides. Attractive functionalities have been reported in some bulk HEOs, which are attributed to the long-range compositional homogeneity, lattice distortion, and local chemical bonding characteristics in materials. However, the intricate details of local composition fluctuation, metal-oxygen bond distortion and covalency are difficult to visualize experimentally, especially on the atomic scale. Here, we study the atomic structure-chemical bonding-property correlations in a series of perovskite-HEOs utilizing the recently developed four-dimensional scanning transmission electron microscopy techniques which enables to determine the structure, chemical bonding, electric field, and charge density on the atomic scale. The existence of compositional fluctuations along with significant composition-dependent distortion of metal-oxygen bonds is observed. Consequently, distinct variations of metal-oxygen bonding covalency are shown by the real-space charge-density distribution maps with sub-ångström resolution. The observed atomic features not only provide a realistic picture of the local physico-chemistry of chemically complex HEOs but can also be directly correlated to their distinctive magneto-electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Annual cycle observations of aerosols capable of ice formation in central Arctic clouds

The Arctic is warming faster than anywhere else on Earth, prompting glacial melt, permafrost thaw, and sea ice decline. These severe consequences induce feedbacks that contribute to amplified warming, affecting weather and climate globally. Aerosols and clouds play a critical role in regulating radiation reaching the Arctic surface. However, the magnitude of their effects is not adequately quantified, especially in the central Arctic where they impact the energy balance over the sea ice. Specifically, aerosols called ice nucleating particles (INPs) remain understudied yet are necessary for cloud ice production and subsequent changes in cloud lifetime, radiative effects, and precipitation. Here, we report observations of INPs in the central Arctic over a full year, spanning the entire sea ice growth and decline cycle. Further, these observations are size-resolved, affording valuable information on INP sources. Our results reveal a strong seasonality of INPs, with lower concentrations in the winter and spring controlled by transport from lower latitudes, to enhanced concentrations of INPs during the summer melt, likely from marine biological production in local open waters. This comprehensive characterization of INPs will ultimately help inform cloud parameterizations in models of all scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observation of electron orbital signatures of single atoms within metal-phthalocyanines using atomic force microscopy

Resolving the electronic structure of a single atom within a molecule is of fundamental importance for understanding and predicting chemical and physical properties of functional molecules such as molecular catalysts. However, the observation of the orbital signature of an individual atom is challenging. We report here the direct identification of two adjacent transition-metal atoms, Fe and Co, within phthalocyanine molecules using high-resolution noncontact atomic force microscopy (HR-AFM). HR-AFM imaging reveals that the Co atom is brighter and presents four distinct lobes on the horizontal plane whereas the Fe atom displays a “square” morphology. Pico-force spectroscopy measurements show a larger repulsion force of about 5 pN on the tip exerted by Co in comparison to Fe. Our combined experimental and theoretical results demonstrate that both the distinguishable features in AFM images and the variation in the measured forces arise from Co’s higher electron orbital occupation above the molecular plane. The ability to directly observe orbital signatures using HR-AFM should provide a promising approach to characterizing the electronic structure of an individual atom in a molecular species and to understand mechanisms of certain chemical reactions.

36 MATERIALS SCIENCE↗