Refraction Tomography Research Articles

Risk Factors Associated with Keratoconus in an Iranian Population.

To determine associated factors for keratoconus (KCN) in the Iranian population. In this retrospective case-control study, 100 KCN patients and 200 age- and sex-matched individuals, who were either candidates for photorefractive keratectomy or healthy referrals from the Torfeh Eye Hospital, were included as the case and control groups, respectively. KCN patients were all registered at the Iranian National Registry of Keratoconus (KCNRegⓇ). Demographic characteristics, patients' symptoms and their habits, as well as systemic and ocular disorders were documented. Clinical examinations included best corrected visual acuity (BCVA) and refractive error measurements, biomicroscopic examination, and corneal imaging. In this case group, the frequency of mild, moderate, and severe KCN was 38%, 28%, and 34%, respectively. Parental consanguinity (odds ratio [OR] = 1.758, P = 0.029), positive familial history in patients' first degree (OR = 12.533, P 0.001) and second degree (OR = 7.52, P 0.001) relatives, vernal keratoconjunctivitis (OR = 7.510, P = 0.003), severe eye rubbing (OR = 10.625, P 0.001), and systemic diseases including migraine, hypertension, and thyroid disease (OR = 6.828, P = 0.021) were found as associated factors for KCN. Lesser frequency of KCN was observed in patients with Fars ethnicity (OR = 0.583, P = 0.042), with higher levels of wealth indices (OR = 0.31, P 0.001) and higher levels of education (OR = 0.18, P = 0.024). Severe eye rubbing, vernal keratoconjunctivitis, parental consanguinity and positive familial history of KCN, low socioeconomic status, and low levels of education were significantly associated with KCN in our study population.

Integrated geophysical methods for boulder delineation to improve mining

AbstractNear‐surface boulders can pose serious challenges to opencast mining. They often introduce complexities, delays in drilling, blasting and excavation programmes, which subsequently decrease mining efficiency, increase mining risks and costs. The location of subsurface boulders and the identification of other geological features that may impact mining activities (e.g. fractures, the presence of iron‐rich ultramafic pegmatites and the variation in weathering across a mining region) are necessary to reduce the challenges posed by these geological features, therefore optimizing mining efficiency. In this study, magnetics, electrical resistivity tomography, seismic refraction tomography, ground penetrating radar and borehole data are integrated for boulder delineation and mapping of other geological features that may impact mining using an unmined section at Tharisa Mine, Bushveld Complex (South Africa), as a test site. The results obtained from the different geophysical techniques are found to complement each other and successfully delineate boulders, fractures, iron‐rich ultramafic pegmatites and the variation in weathering and layering across the area. The incorporation of geophysical results can thus improve mining efficiency, while reducing mining risks and costs.

Quantification of structural heterogeneity in H&E stained clear cell renal cell carcinoma using refractive index tomography.

Clear cell renal cell carcinoma (ccRCC) is a common histopathological subtype of renal cancer and is notorious for its poor prognosis. Its accurate diagnosis by histopathology, which relies on manual microscopic inspection of stained slides, is challenging. Here, we present a correlative approach to utilize stained images and refractive index (RI) tomography and demonstrate quantitative assessments of the structural heterogeneities of ccRCC slides obtained from human patients. Machine-learning-assisted segmentation of nuclei and cytoplasm enabled the quantification at the subcellular level. Compared to benign regions, malignant regions exhibited a considerable increase in structural heterogeneities. The results demonstrate that RI tomography provides quantitative information in synergy with stained images on the structural heterogeneities in ccRCC.

Toward high-fidelity imaging: Dynamic matching FWI and its applications

Full-waveform inversion (FWI) is firmly established within our industry as a powerful velocity model building tool. FWI carries significant theoretical advantages over conventional velocity model building methods such as refraction and reflection tomography. Specifically, by solving a nonlinear inverse problem through the wave equation, FWI is able to recover a broadband velocity model containing both high and low spatial wavenumbers, thus extending the approximation of residual moveout correction inherent in traditional velocity model building approaches. Moreover, FWI is capable of inverting information from the entire wavefield (i.e., early arrivals, reflections, refractions, and multiple energy) rather than from a subset as in conventional approaches (i.e., first break and primary reflections), thereby availing itself of more information to better constrain its model estimate. However, these theoretical benefits cannot be realized easily in practice because various complexities of real seismic data often conspire to violate algorithmic assumptions, leading to unsatisfactory results. Dynamic matching FWI (DMFWI) is a newly developed algorithm that solves an inversion problem that maximizes the cross correlation of two dynamically matched data sets — one recorded and the other synthetic. Dynamic matching of the two data sets de-emphasizes the amplitude impact, which allows the algorithm to focus on minimizing their kinematic differences rather than amplitude in the data-fitting process. The multichannel correlation makes the algorithm robust for data with low signal-to-noise ratio. Applications of DMFWI across different types of acquisition and geologic settings demonstrate that this novel FWI approach can resolve complex velocity errors and provide high-quality migrated images that exhibit a high degree of geologic plausibility. Additionally, reflectivity images can be obtained in a straightforward manner as natural byproducts through computation of the directional derivative of the inverted FWI velocity models.

Geophysical Imaging of the Shallow Geyser and Hydrothermal Reservoir Structures of Spouter Geyser, Yellowstone National Park: Geyser Dynamics I

AbstractYellowstone National Park (YNP) is home to roughly 500 geysers, making it the most concentrated geyser field in the world. Recent studies exploring the mechanics of geyser eruptions utilizing laboratory models were limited by a lack of knowledge of the subsurface geometry of the geyser features. This study presents results from active hydrogeophysical surveys, including ground penetrating radar, nuclear magnetic resonance (NMR), seismic refraction, electrical resistivity tomography, and transient electromagnetics to image subsurface geyser structure and constrain the geophysical response of the reservoir structure of Spouter Geyser, Yellowstone National Park. Previous geophysical studies on similar geyser systems in Yellowstone characterize the hydrothermal reservoir that supplies the geyser with hydrothermal fluids and vapors as high porosity, low density structures. Whereas our imaging highlights the hydrothermal conduits and reservoir as high resistivity and high velocity structures, interpreted as silica precipitate decreasing the porosity and increasing the bulk modulus (i.e., from unconsolidated media to semi/fully consolidated). Ground penetrating radar identifies sinter thickness around the geyser, NMR provides porosity measurements of the geyser reservoir, and transient electromagnetics maps the depth to the Biscuit Basin Rhyolite bedrock throughout the field site. Electrical resistivity tomography and seismic refraction image the shallow geyser conduit structures to ∼15 m depth and establish the extent of the geyser hydrothermal reservoir structure. Furthermore, a spatial correlation of the electrical resistivity and seismic velocity results establishes the interpreted hydrothermal geyser reservoir as a high resistivity, high velocity structure located to the northeast of Spouter Geyser at depths greater than 15 m.

Pre‐stack depth imaging techniques for the delineation of the Carosue Dam gold deposit, Western Australia

AbstractIn this study, we explore the latest generation of seismic imaging algorithms (migration) that have been successful in the oil and gas exploration industry and apply them in the more challenging hard rock environment. Seismic migration is a crucial processing step required to build an accurate image of the Earth's subsurface. The seismic method applied in a hard rock environment has a specific set of challenges: complex geological settings often comprised of steeply dipping interfaces and heterogeneities (faults, fracture zones, thrusts, etc.); spatially variable zones of alteration, low intrinsic signal‐to‐noise ratio; complex near‐surface conditions (the weathered overburden has a very high contrast in seismic properties with base formations). Here, we present how a dense source–receiver 3D grid in combination with the latest pre‐stack depth imaging techniques can image the geology with a remarkable level of detail and to depths as shallow as the top of fresh rock. We also propose a comprehensive velocity model building strategy applied specifically to the Carosue Dam deposit, Western Australia, which is mainly characterized by volcaniclastic and volcanic rocks within the Carosue Basin disrupted by a complex system of gold‐bearing faults. To arrive at an optimized velocity model for pre‐stack depth imaging, we combine drillhole data with tomographic refinement. Excellent correlation between resultant seismic images and refraction tomography, vertical seismic profiling and drillholes can be attributed to the integrated velocity model building workflow used for depth migration.

Seismic refraction tomography using underwater and land based seismic data for evaluation of foundation of civil structures

Inadequate subsurface information around the proposed structure often leads to disruptions to the plan of construction, resulting in delays in the implementation of design, and finally hampering the progress of the work. Therefore, geophysical and geotechnical investigations prior to any civil construction should be adopted to characterize the rock mass conditions at a given site. Quick results on engineering evaluation of sites are imminent due to tight timelines and when site conditions do not favor readily deployable testing methods. Suitable tailor-made acquisition techniques need to be adopted to arrive at the best evaluation possibility. Two such cases are discussed here, wherein one site had proposed barrage alignment, abandoned after laying foundations for four pillars and in the second site, foundation of three pillars of a road bridge was unexpectedly met with poor rock conditions. In order to investigate and evaluate the foundation rock conditions, Seismic Refraction (SR) survey was carried out in both the cases. SR data were acquired using tailor-made array across the river using hydrophones and geophones on water and land respectively. This paper discusses the pragmatic approach in data acquisition for the given site conditions along with the pertinent data processing using Seismic Refraction Tomography (SRT) technique, analysis, interpretation of the results correlated with the geotechnical data at the same site. The outcome of the investigations and the benefits of the pragmatic approach for evaluation of foundation in such cases are aided in design modifications of the proposed civil structures.

Unified Modulation Pattern Analysis (UMPA) algorithm for 1D sensitive X-ray phase contrast imaging techniques

X-ray phase contrast imaging (XPCI) techniques are sensitive to refraction (differential-phase) and small-angle X-ray scattering (dark-field) signals, not measurable with conventional absorption imaging techniques. Among XPCI techniques, edge illumination (EI), grating interferometry (GI), and speckle-based imaging (SBI) make use of wavefront markers, such as absorbing masks with periodical apertures or random diffusers, to encode refraction and dark-field signals induced by the sample. The Unified Modulated Pattern Analysis (UMPA) provides an algorithmic solution to extract the transmission, refraction, and dark-field images from EI, GI, and SBI datasets where the wavefront marker is directly resolved by the employed detection system. In its original implementation, UMPA has been designed for XPCI techniques sensitive to refractions along two axes. This work presents a modified version of the algorithm to extend its applicability to all the existing XPCI techniques that use wavefront markers with sensitivity to refraction limited along one direction (UMPA-1D). The algorithm, written in C++ and Cython and parallelized with OpenMP, enables fast reconstruction times that are particularly convenient for large tomographic datasets. The validity of the UMPA-1D has been demonstrated using both simulated images and real acquisitions with an EI setup in beam-tracking mode.

Modeling of refractive images during the propagation of laser beams in optically inhomogeneous liquid media

Modeling of refractive images during the propagation of laser beams in optically inhomogeneous liquid media

Optimized Artificial Intelligence for Enhanced Ectasia Detection Using Scheimpflug-Based Corneal Tomography and Biomechanical Data

To optimize artificial intelligence (AI) algorithms to integrate Scheimpflug-based corneal tomography and biomechanics to enhance ectasia detection. Multicenter cross-sectional case-control retrospective study. A total of 3886 unoperated eyes from 3412 patients had Pentacam and Corvis ST (Oculus Optikgeräte GmbH) examinations. The database included 1 eye randomly selected from 1680 normal patients (N) and from 1181 "bilateral" keratoconus (KC) patients, along with 551 normal topography eyes from patients with very asymmetric ectasia (VAE-NT), and their 474 unoperated ectatic (VAE-E) eyes. The current TBIv1 (tomographic-biomechanical index) was tested, and an optimized AI algorithm was developed for augmenting accuracy. The area under the receiver operating characteristic curve (AUC) of the TBIv1 for discriminating clinical ectasia (KC and VAE-E) was 0.999 (98.5% sensitivity; 98.6% specificity [cutoff: 0.5]), and for VAE-NT, 0.899 (76% sensitivity; 89.1% specificity [cutoff: 0.29]). A novel random forest algorithm (TBIv2), developed with 18 features in 156 trees using 10-fold cross-validation, had a significantly higher AUC (0.945; DeLong, P < .0001) for detecting VAE-NT (84.4% sensitivity and 90.1% specificity; cutoff: 0.43; DeLong, P < .0001) and a similar AUC for clinical ectasia (0.999; DeLong, P=.818; 98.7% sensitivity; 99.2% specificity [cutoff: 0.8]). Considering all cases, the TBIv2 had a higher AUC (0.985) than TBIv1 (0.974; DeLong, P < .0001). AI optimization to integrate Scheimpflug-based corneal tomography and biomechanical assessments augments accuracy for ectasia detection, characterizing ectasia susceptibility in the diverse VAE-NT group. Some patients with VAE may have true unilateral ectasia. Machine learning considering additional data, including epithelial thickness or other parameters from multimodal refractive imaging, will continuously enhance accuracy. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Application of Resistivity and Seismic Refraction Tomography for Landslide Stability Assessment in Vallcebre, Spanish Pyrenees

Geophysical surveys are a noninvasive reliable tool to improve geological models without requiring extensive in situ borehole campaigns. The usage of seismic refraction tomography (SRT), electrical resistivity tomography (ERT) and borehole data for calibrating is very appropriate to define landslide body geometries; however, it is still only used occasionally. We present here the case of a Spanish Pyrenees slow-moving landslide, where ERT, SRT and lithological log data were integrated to obtain a geological three-dimensional model. The high contrasts of P-wave velocity and electrical resistivity values of the upper materials (colluvial debris and clayey siltstone) provided accurate information on the geometry of the materials involved in the landslide body, as well as the sliding surface. Geophysical prospecting allowed us to identify the critical sliding surface over a large area and at a reduced cost and, therefore, gives the geophysical method an advantage over borehole data. The three-dimensional model was used to carry out stability analyses of a landslide in 2D and 3D, which, coherently with previous studies, reveal that the lower part is more unstable than the upper units.

Experimental demonstration of ultrahigh sensitivity Talbot-Lau interferometer for low dose mammography

Objective. Even though the techniques used for breast cancer identification have advanced over the years, current mammography based on x-rays absorption, the ‘gold standard’ screening test at present, still has some shortcomings as concerns sensitivity and specificity to early-stage cancers, due to poor differentiation between tumor and normal tissues, especially in the case of the dense breasts. We investigate a possible additional technique for breast cancer detection with higher sensitivity and low dose, x-ray phase-contrast or refraction-based imaging with ultrahigh angular sensitivity grating interferometers, having several meters length. Approach. Towards this goal, we built and tested on a mammography phantom, a table-top laboratory setup based on a 5.7 m long Talbot-Lau interferometer with angular sensitivity better than 1 μrad. We used a high-power x-ray tungsten anode tube with a 400 μm focal spot, operated at 40 kVp and 15 mA with a 2 mm aluminum filter. Main results. The results reported in our paper confirm the ultrahigh sensitivity and dose economy possible with our setup. The visibility of objects simulating cancerous formations is strongly increased in the refraction images over the attenuation ones, even at a low dose of 0.32 mGy. Notably, the smallest fiber of 400 μm diameter and calcifications specs of 160 μm in diameter are detected, even though the spatial resolution at the object of our magnification M ∼ 2 setup with a 400 μm source spot is only ∼250 μm. Significance. Our experiments on a mammography phantom illustrate the capabilities of the proposed technique and can open the way toward low-dose interferometric mammography.

Near-Surface Geophysical Characterization of Lithologies in Corfu and Lefkada Towns (Ionian Islands, Greece)

Lefkada and Corfu old towns are located in the western part of Greece, in the Ionian Sea. Their proximity to the Hellenic subduction zone (HSZ) is the reason for their intense seismicity. The main goal of this study was the estimation of the geotechnical characteristics of the subsurface, with the contribution of applied geophysical techniques. Therefore, seismic refraction tomography (SRT) and multichannel analysis of surface waves (MASW) were applied. A total of thirty-three (33) seismic and geoelectrical profiles were performed in both towns in order to evaluate the geotechnical characteristics of the subsurface formations. Additionally, subsurface resistivity distributions were investigated with the application of electrical resistivity tomography (ERT). Some important elastic moduli were calculated through the combination of estimated seismic wave velocities and laboratory density measurements. The horizontal distribution of seismic velocities and mechanical properties (σ, E, K, G) of Corfu town was illustrated in maps, for the depth of 5 m. The geophysical interpretation also revealed that Lefkada’s subsurface consists of only one compact geological formation, with little or no variation of its geophysical-geotechnical characteristics. Beyond that, the ground type classifications for the two towns were determined according to the European Committee for Standardization Eurocode 8, based on VS30 values.

Enhanced Diagnostics for Corneal Ectatic Diseases: The Whats, the Whys, and the Hows.

There are different fundamental diagnostic strategies for patients with ectatic corneal diseases (ECDs): screening, confirmation of the diagnosis, classification of the type of ECD, severity staging, prognostic assessment, and clinical follow-up. The conscious application of such strategies enables individualized treatments. The need for improved diagnostics of ECD is related to the advent of therapeutic refractive procedures that are considered prior to keratoplasty. Among such less invasive procedures, we include corneal crosslinking, customized ablations, and intracorneal ring segment implantation. Besides the paradigm shift in managing patients with ECD, enhancing the sensitivity to detect very mild forms of disease, and characterizing the inherent susceptibility for ectasia progression, became relevant for identifying patients at higher risk for progressive iatrogenic ectasia after laser vision correction (LVC). Moreover, the hypothesis that mild keratoconus is a risk factor for delivering a baby with Down's syndrome potentially augments the relevance of the diagnostics of ECD. Multimodal refractive imaging involves different technologies, including Placido-disk corneal topography, Scheimpflug 3-D tomography, segmental or layered tomography with layered epithelial thickness using OCT (optical coherence tomography), and digital very high-frequency ultrasound (VHF-US), and ocular wavefront. Corneal biomechanical assessments and genetic and molecular biology tests have translated to clinical measurements. Artificial intelligence allows for the integration of a plethora of clinical data and has proven its relevance in facilitating clinical decisions, allowing personalized or individualized treatments.

Integrated Agrogeophysical Approach for Investigating Soil Pipes in Agricultural Fields

Soil erosion is one of the most significant challenges for soil management and agri-food production threatening human habitat and livelihood. Although soil erosion due to surficial processes is well-studied, erosion due to subsurface processes such as internal soil pipes has often been overlooked. Internal soil pipes directly contribute to the total soil loss in agricultural fields and impede agricultural sustainability. Locating, measuring, and mapping internal soil pipes and their networks are vital to assessing the total soil loss in agricultural fields. Their hidden and uncorrelated nature of subsurface occurrences constricts the applicability of manual and remote sensing-based detection techniques. Non-invasive agrogeophysical methods can overcome these limitations with detailed subsurface pictures and high spatial resolution. In this study, the applicability of three agrogeophysical methods including seismic refraction tomography (SRT), electrical resistivity tomography (ERT), and ground-penetrating radar (GPR) was tested at Goodwin Creek, an experimental field site with established internal soil pipes. SRT showed low P and S wave velocities anomalies in soil pipe-affected zones. ERT results indicated the location of soil pipes with high resistivity anomalies. However, both SRT and ERT lack resolution to identify individual soil pipes. GPR diffraction hyperbolas and their apexes however effectively identified individual soil pipes. The agrogeophysical anomalies for soil pipes were compared with the low penetration resistance of the cone penetrologger (CPL) results. Correspondence between low PR in CPL and agrogeophysical anomalies verify the locations of internal soil pipe-affected zones. Moreover, the fragipan layer is identified below the soil pipe-affected zone by all three methods.

Dual Optical and Acoustic Negative Refraction in Phoxonic Crystals

We report dual optical and acoustic negative refraction based on a defect-free phoxonic crystal within a triangular lattice. The phoxonic negative refraction is achieved based on abnormal dispersion effect, by intentionally creating convex equal-frequency contours for both photonic and phononic modes. As a potential application, negative refraction imaging for both photonic and phononic modes is also achieved. Numerical simulations based on the finite element method demonstrate the coexistence of negative refraction and the resultant imaging for electromagnetic and acoustic waves. Compared with the defect-based bandgap effects that need low fault tolerance, phoxonic negative refraction relying on passbands has considerable advantages in realizing controllable propagation of waves. The new scheme for the simultaneous control of electromagnetic and acoustic waves provides a potential platform for designing novel phoxonic devices.

Characterization of an earth dam in Brazil from seismic refraction tomography and multichannel analysis of surface waves

Effective methodologies to obtain information regarding the internal features in an earth dam are fundamental for rapid and technically correct decision making. Although several geophysical methods have potential applicability for investigations on dams, the contribution of seismic methods is well known, since the P and S-waves velocities can be associated with elastic modules of geotechnical interest. The study area was the Paranoá dam, the largest dam in Brazil's federal capital. We analyzed the internal structure of the dam and identified anomalies that could be associated with water saturation. We used traveltimes tomography to obtain Vp models and the MASW method to obtain Vs models. Four profiles of seismic data were acquired downstream of the dam. Borehole's information was used to corroborate with the geophysical interpretation of the profiles. The dam's abutments and the foundation ground could be identified within most models, and the velocity values obtained were mainly related to the clay material that forms the dam's massif, with higher values associated with the quartzite rocks. The rockfill material was well marked as a low velocity zone. Saturated zones were interpreted as local anomalies of high Vp/Vs ratio. An analysis of the depth range obtained within velocity models suggests that the depth of investigation from both methodologies are likely to be site specific, rather than exclusive matter of source and receiver instrumentation.

Underwater acoustic metamaterials.

Acoustic metamaterials have been widely investigated over the past few decades and have realized acoustic parameters that are not achievable using conventional materials. After demonstrating that locally resonant acoustic metamaterials are capable of acting as subwavelength unit cells, researchers have evaluated the possibility of breaking the classical limitations of the material mass density and bulk modulus. Combined with theoretical analysis, additive manufacturing and engineering applications, acoustic metamaterials have demonstrated extraordinary capabilities, including negative refraction, cloaking, beam formation and super-resolution imaging. Owing to the complexity of impedance boundaries and mode transitions, there are still challenges in freely manipulating acoustic propagation in an underwater environment. This review summarizes the developments in underwater acoustic metamaterials over the past 20years, which include underwater acoustic invisibility cloaking, underwater beam formation, underwater metasurfaces and phase engineering, underwater topological acoustics and underwater acoustic metamaterial absorbers. With the evolution of underwater metamaterials and the timeline of scientific advances, underwater acoustic metamaterials have demonstrated exciting applications in underwater resource development, target recognition, imaging, noise reduction, navigation and communication.

Detecting sub-basalt Mesozoic sediments and active tectonics in the Narmada-Tapti Rift Zone, central India

To obtain structural information below basalt cover using conventional seismic reflection has always been a challenge due to the high impedance of the basaltic layer. This study enumerates P-wave refraction tomography to image the shallow structure along three long offset profiles in the Deccan Volcanic Province, one of the world's largest flood basalt volcanic regions. The travel time inversion is performed to obtain a starting model for the wavefield inversion, and synthetic tests are used to estimate the resolution. Our results show an abrupt increase of crystalline basement depth from 1.5 to 3 km across the Narmada South Fault. The shallow basement depth coincides with a 70-km-wide sub-basalt (>1.5 km deep) seismic shadow zone between the Narmada and Tapti rivers, associated with the sub-basalt graben caused by a crustal extension that contains sediments or a combination of basalt and sediments. The basaltic layer appears to be highly heterogeneous in velocity and thickness, and the degree of heterogeneity increases laterally and vertically with depth within the layer. It is quite likely that these sub-basalt sediments could be the potential zone for hydrocarbon, similar to the exposed Vindhyan in the north and Gondwana sediments in the south of our study area. Another intriguing feature is the Barwani Sukta Fault, which controls the tectonics of the region by dividing the basement configuration. We also compare our tomographic model with those determined from previously obtained seismic, gravity, and surface mapping. The tomographic image clearly shows that numerous faults are reaching up to the basement and extending laterally. We suggest that these neo-tectonic faults may be capable of generating large-magnitude earthquakes, including the 21 May 1997 (mb 6.0) Jabalpur earthquake.

High-resolution spatiotemporal analysis of single serotonergic axons in an in vitro system.

Vertebrate brains have a dual structure, composed of (i) axons that can be well-captured with graph-theoretical methods and (ii) axons that form a dense matrix in which neurons with precise connections operate. A core part of this matrix is formed by axons (fibers) that store and release 5-hydroxytryptamine (5-HT, serotonin), an ancient neurotransmitter that supports neuroplasticity and has profound implications for mental health. The self-organization of the serotonergic matrix is not well understood, despite recent advances in experimental and theoretical approaches. In particular, individual serotonergic axons produce highly stochastic trajectories, fundamental to the construction of regional fiber densities, but further advances in predictive computer simulations require more accurate experimental information. This study examined single serotonergic axons in culture systems (co-cultures and monolayers), by using a set of complementary high-resolution methods: confocal microscopy, holotomography (refractive index-based live imaging), and super-resolution (STED) microscopy. It shows that serotonergic axon walks in neural tissue may strongly reflect the stochastic geometry of this tissue and it also provides new insights into the morphology and branching properties of serotonergic axons. The proposed experimental platform can support next-generation analyses of the serotonergic matrix, including seamless integration with supercomputing approaches.