The bearing capacity of a foundation determines its ability to support loads without experiencing shear failure or excessive settlement. Traditionally, geotechnical analysis has treated this capacity as a relatively constant parameter, influenced primarily by instantaneous material properties and load conditions. However, the reality beneath our feet is far less static. Time-dependent degradation, creep, weathering, and stress redistribution can all significantly alter the mechanical properties of weak rock masses over time, thus compromising foundational stability.

This latest research utilizes advanced numerical modeling techniques to simulate the behavior of weak rock masses under varying load and environmental conditions. The study’s simulations incorporate parameters such as creep deformation, progressive failure mechanisms, and time-dependent reductions in shear strength. By integrating these complex processes into a cohesive numerical framework, the authors have achieved a more precise and predictive model of foundation bearing capacity that evolves over the lifespan of a structure.

One of the standout features of this investigation is its focus on weak rock masses — geological formations characterized by reduced intact strength, pervasive fracturing, and an innate susceptibility to environmental influences. Despite their prevalence, these materials have often been sidelined in engineering analyses due to the challenges in characterizing their behavior accurately. The insights provided by this study, therefore, mark a significant advancement in geotechnical engineering.

At the heart of the research is a sophisticated finite element model that realistically captures the mechanical responses of weak rock under sustained loading. The model’s architecture considers time-dependent phenomena such as creep, wherein rock masses under stress gradually deform, and subcritical crack growth, a process by which microfractures propagate slowly, reducing structural integrity over extended periods. These processes, often overlooked in traditional design standards, are crucial for anticipating long-term risks.

Moreover, the authors highlight how environmental factors, including moisture fluctuations and temperature variations, modulate the time-dependent behavior. Water infiltration can accelerate weathering and chemical alteration, softening rock interfaces and encouraging delayed failure. This environmental coupling necessitates integrating hydromechanical factors alongside the purely mechanical ones in any limit-state analysis.

Quantitatively, the simulations suggest that the effective bearing capacity of foundations on weak rock may degrade by significant margins—sometimes up to 30%—within decades post-construction under moderate loading conditions. Such reductions underscore the potential for catastrophic failures if temporal effects are not rigorously accounted for during the design phase. The researchers advocate for a recalibration of safety factors to incorporate these delayed strength losses explicitly.

The study’s implications extend beyond academic curiosity into real-world applications, as infrastructure increasingly encroaches on geologies riddled with weak rock masses. Urban development in mountainous or sedimentary basin regions, mining operations, and even renewable energy installations such as wind turbines all demand a granular understanding of foundation performance over time. Misestimating bearing capacity could lead to costly repairs, structural collapses, or even loss of life.

Crucially, the numerical approach presented proposes a pathway toward more resilient design codes. By parametrizing time-dependent effects, engineers and decision-makers can generate predictive maintenance schedules, implement real-time monitoring systems, and improve early warning mechanisms. This digital foresight tool allows better resource allocation and proactive mitigation strategies, ultimately enhancing societal safety margins.

The multidisciplinary nature of this work also signals a shift in how geotechnical engineering integrates with environmental sciences. The coupling of mechanical modeling with hydrogeological and chemical processes represents a new frontier in foundation research, reflective of the complex reality governing subsurface conditions. Such integrative frameworks pave the way for holistic infrastructure risk assessments.

While the study is numerical in essence, its authors emphasize the necessity of experimental validation. Field studies, in situ testing, and long-term monitoring campaigns are vital to calibrate and verify the modeled parameters accurately. Future research is likely to focus on refining these models with empirical datasets sourced from diverse weak rock environments worldwide.

In conclusion, the findings by Mamaghani, Zaheri, and Ranjbarnia substantially enrich our comprehension of time-dependent effects on foundation bearing capacity within weak rock masses. Their work offers a vital cautionary tale against static design assumptions while providing engineers with the tools to anticipate and counteract time-induced degradation effectively. The influence of their approach on both academic research and practical engineering practices is poised to be profound and enduring.

As infrastructure resilience takes center stage globally, harnessing these predictive insights becomes a linchpin for sustainable construction in geologically challenging contexts. This study not only bridges a critical knowledge gap but also empowers engineers with an enhanced predictive shield against the hidden risks lurking beneath our built environment, encouraging a smarter, safer future.

Subject of Research: Numerical investigation of time-dependent effects on foundation bearing capacity on weak rock masses.

Article Title: Numerical investigation of time-dependent effects on the bearing capacity of foundations on weak rock masses.

Article References:
Motamedi Mamaghani, F., Zaheri, M. & Ranjbarnia, M. Numerical investigation of time-dependent effects on the bearing capacity of foundations on weak rock masses. Environ Earth Sci 84, 610 (2025). https://doi.org/10.1007/s12665-025-12625-2

Image Credits: AI Generated

The seismic tomography model underpinning this research was assembled from an extensive dataset collected between 2007 and 2013, involving 387 broadband land stations spread throughout the Ibero-Maghrebian region. Importantly, the study incorporated data from 24 ocean-bottom seismometers deployed offshore Southwest Iberia during the NEAREST experiment, alongside instruments from the TOPOMED project, enhancing offshore ray coverage. Analyzing over 25,000 arrival-time residuals from 451 teleseismic events with magnitudes exceeding 5.5, researchers meticulously constructed a detailed velocity model extending from the crust-mantle boundary (the Moho) to depths of 800 kilometers. This comprehensive dataset was foundational in resolving fine-scale structures in the upper mantle, critical for identifying delamination signatures.

Seismic wave travel times were initially aligned and filtered to isolate relevant signals, followed by an adaptive stacking procedure that refined the seismic phase arrival estimates. The innovative use of the FMTOMO package allowed for the iterative inversion of these time residuals, solving the forward travel-time problem through a grid-eikonal method known as the Fast Marching Method. This approach uncovers three-dimensional variations in seismic wave velocity that correspond to temperature, compositional, and structural heterogeneities deep within the Earth. Crucially, the inversion accounted for crustal effects using a prior three-dimensional model (PRISM3D), which corrected for topographic and velocity variations in the crust that might otherwise distort mantle imaging.

To evaluate the robustness of the tomography results, a synthetic spike resolution test was performed. This numerical experiment introduced pairs of velocity anomalies with known properties into the starting model to test whether these features could be reliably recovered by the inversion process. The successful identification of these synthetic anomalies in the output model confirmed the high accuracy and resolution of the seismic imaging, particularly for uppermost mantle structures near 90 kilometers depth. This accomplishment lends strong support to the interpretation that the imaged high-velocity anomalies offshore Southwest Iberia represent genuine lithospheric structures consistent with delaminated oceanic material.

Complementing the seismic imaging, the research team employed advanced numerical modeling to explore the mechanics of oceanic plate delamination. Using the computational platform Underworld, they simulated the coupled processes of momentum, mass conservation, and thermal evolution under realistic boundary conditions. The models solve governing equations of fluid dynamics and heat transfer, incorporating nonlinear rheologies governed by temperature, pressure, and strain-rate-dependent viscosity. Importantly, the mechanical behavior includes viscoplastic deformation, with yielding determined by a Drucker–Prager criterion that accounts for frictional failure and plastic strain weakening. This detailed formulation allows the model to capture complex interactions between brittle fracture, ductile flow, and thermal weakening—all critical for realistic simulation of lithospheric peeling.

The numerical experiments were conducted within a large two-dimensional domain measuring 2,800 kilometers in length and 660 kilometers deep, discretized with thousands of finite elements to ensure fine spatial resolution—down to 1.25 kilometers within the lithosphere. The modeled geometry includes two contrasting oceanic plates: a thicker, older Africa-like plate beneath the southern part of the model, and a younger, thinner Eurasia-like plate to the north. Notably, the younger plate incorporates a serpentinized mantle layer, a low-viscosity zone prone to weakening, reflecting real geological observations from seismic refraction profiles. These contrasting lithospheric features create conditions conducive to delamination under tectonic compression.

To simulate natural convergence, a slow northward velocity of 8 millimeters per year was imposed on the African-like plate, with the Eurasian-like plate fixed in place, replicating the Cenozoic Africa–Eurasia plate motions. Multiple scenarios were tested, varying the presence and thickness of serpentinized layers and vertical weak zones that represent inherited faults or fractures. Models with two vertical weak zones evenly spaced and a 10-kilometer-thick serpentinized weak layer best matched observed seismic data and geological constraints, faithfully reproducing the delamination process. These results underscore the critical role of preexisting lithospheric weaknesses and compositional heterogeneities in facilitating such complex tectonic behavior.

The simulations illuminate the dynamic progression of delamination, showing that gravitational forces and induced stresses cause the dense oceanic lower lithosphere to detach and sink into the mantle. This peeling away disrupts isostatic equilibrium and modifies mantle flow patterns, potentially triggering volcanism and seismicity. Notably, the study also explored the influence of stopping convergence after 18 million years, finding that delaminated blocks may continue sinking under gravity alone, highlighting the interplay between tectonic forcing and buoyancy-driven dynamics. This insight refines previous conceptions of delamination duration and its feedbacks with surface tectonic processes.

The decision to pursue a two-dimensional modeling approach was strategic. The elongated geometry of the delaminating structure, oriented perpendicular to the convergence direction, supports the assumption of plane-strain symmetry. Moreover, focusing on a simplified framework enabled systematic parametric studies of key controlling mechanisms without the computational burden and complexity of full three-dimensional modeling. While three-dimensional effects are expected in nature, this minimalistic modeling provided essential physical understanding, serving as a proof-of-concept to test hypotheses derived from seismic observations.

Advanced rheological formulations underpin the simulations, with effective viscosity calculated via experimentally derived flow laws that incorporate activation energy and volume, stress exponent, and temperature dependence. The models capture the transition from ductile creep at high temperatures and pressures to brittle failure at shallower depths. Incorporation of strain weakening mimics the progressive loss of strength as deformation accumulates, reproducing realistic lithospheric weakening that fosters delamination initiation. These physically based constitutive laws enhance model fidelity and predictive power, bridging laboratory rheology and geodynamic processes.

Thermomechanical coupling is central to the model, with temperature evolution governed by an advection-diffusion equation incorporating shear heating and adiabatic heating terms. Shear heating arises from viscous deformation work, while adiabatic heating relates to compression under mantle conditions. These thermal effects modify viscosity and density distributions, feeding back into deformation patterns and delamination progression. This coupling mirrors natural conditions where thermal and mechanical processes are intertwined, adding another layer of realism to the model outcomes.

The integration of seismic tomography and numerical modeling in this study represents a pioneering approach in geosciences, shedding light on complex lithosphere-mantle interactions offshore Southwest Iberia. The high-resolution seismic images confirm the presence of a dense, high-velocity anomaly interpreted as a delaminated oceanic slab fragment descending into the mantle, while the sophisticated simulations reveal governing physical mechanisms and key parameters controlling the process. This dual methodology sets a benchmark for future multidisciplinary investigations of plate dynamics in regions where direct observation is impossible.

Findings from this research have profound implications beyond Southwest Iberia. Understanding oceanic plate delamination is fundamental to deciphering tectonic regime changes, intraplate volcanism, seismic hazard, and mantle-driving forces globally. The observed link between inherited lithospheric fabrics, serpentinization, and delamination initiation offers new perspectives on how plate weakening modulates large-scale Earth dynamics. Moreover, the study emphasizes the importance of integrating seismological, geological, and numerical evidence to unravel deep Earth processes, inspiring a holistic paradigm in geodynamics.

Looking ahead, the research team envisions extending their framework to incorporate three-dimensional geometries, anisotropic material properties, and coupling with surface processes such as erosion and sedimentation. Such enhancements will enable even more detailed reconstructions of lithospheric evolution and its surface manifestations. Furthermore, applying similar approaches to other convergent margins worldwide can test the ubiquity and variability of delamination phenomena, providing a richer understanding of the Earth’s tectonic mosaic.

The multidisciplinary nature of this breakthrough underscores the synergy emerging in Earth sciences as computational power increases and data acquisition evolves. Combining dense seismic observations with cutting-edge modeling tools allows scientists to “see” and simulate hidden processes shaping our planet’s lithosphere. Studies like this herald a new era where theory, observation, and computation converge to solve longstanding geodynamic puzzles, with promising impacts on hazard assessment, resource exploration, and fundamental Earth science.

In conclusion, the seismic evidence for oceanic plate delamination revealed offshore Southwest Iberia not only solves an important regional geodynamic mystery but also opens pathways for novel explorations into plate tectonics and mantle convection. The coupling of meticulous seismic imaging and robust numerical simulations demonstrates the power of integrative science in unraveling the deep Earth’s secrets. As we push the limits of resolution and computational sophistication, our planetary understanding becomes ever clearer, revealing the dynamic tapestry beneath our feet.

Subject of Research:

Article Title:

Article References:
Duarte, J.C., Riel, N., Civiero, C. et al. Seismic evidence for oceanic plate delamination offshore Southwest Iberia. Nat. Geosci. (2025). https://doi.org/10.1038/s41561-025-01781-6

Image Credits: AI Generated

DOI:

Keywords:

At the core of this pioneering research lies the nuanced understanding that real-world slopes rarely exhibit uniform mechanical properties. Instead, anisotropy—the variation of material strength and deformation characteristics with direction—plays a crucial role in influencing how slopes fail and how landslides propagate downhill. Chang, Lu, Yeh, and their team have harnessed advanced numerical modeling techniques to capture these directional effects, providing a window into landslide dynamics that could transform risk assessment and mitigation strategies worldwide.

The team’s approach integrates anisotropic material behavior directly into dynamic simulation models that replicate the initiation, acceleration, and eventual deposition phases of landslide movement. By encoding inherent directional dependencies of slope materials into the governing equations, their framework surpasses traditional isotropic models that assume uniformity in every direction. This distinction is vital as anisotropic conditions frequently arise from geological layering, foliation, or sedimentary structures that predispose slopes to fail preferentially along certain planes, resulting in distinctive runout patterns.

One of the most striking insights from the study is the influence of anisotropy on the velocity profile and distance traveled by landslide masses. Unlike isotropic slopes, where runout distance tends to follow more predictable and symmetrical distributions, anisotropic slopes display marked asymmetries—exhibiting directional biases that can either accelerate or decelerate different portions of a moving slide mass. This heterogeneity in flow behavior highlights the importance of incorporating directional mechanical properties when simulating potential landslide scenarios, especially for predictive hazard zoning.

Moreover, the researchers identified that anisotropy could amplify or attenuate the energy dissipation mechanisms intrinsic to granular flow during runout. This has profound implications for the calculation of impact forces exerted on infrastructure and populated areas situated in landslide-prone regions. Accurate prediction of these forces is essential for designing effective protective barriers and early warning systems, making the modeling framework introduced by Chang and colleagues a valuable tool for engineers and policymakers alike.

To validate their numerical models, the team conducted a series of simulations comparing landslide runout on slopes with varying degrees and orientations of anisotropy. These simulations revealed that even subtle variations in mechanical properties can drastically alter runout behavior, emphasizing the sensitivity of landslide dynamics to underlying geological heterogeneities. This sensitivity underlines the need for detailed site investigations that characterize anisotropic properties, supporting the integration of field data into simulation workflows.

The methodology employed deviates markedly from traditional empirical and isotropic computational approaches by exploiting anisotropic constitutive models capable of reflecting real-world conditions more faithfully. These constitutive models describe the stress-strain relationships that govern material deformation and failure, accommodating phenomena such as directional shear strength and variable permeability. By embedding these models into finite element or finite difference frameworks, the researchers produced simulations that closely mirror observed landslide behaviors.

Beyond theoretical refinement, the practical applications of this research are extensive. In mountainous regions where anisotropic structures are prevalent due to complex geological histories, the ability to accurately predict landslide runout becomes a critical component of disaster preparedness. Urban planners and civil engineers can harness such precise models to optimize land use planning, ensuring that residential and infrastructural developments avoid zones of exacerbated landslide risk.

Furthermore, climate change-driven increases in intense precipitation events are likely to trigger more frequent and larger landslides. Understanding anisotropy’s role in modulating runout is therefore timely. Incorporating anisotropic effects into early warning systems can enhance their predictive accuracy, potentially saving lives by providing better estimates of both timing and spatial extent of landslide occurrences. This represents a significant leap forward in integrating material science with environmental hazard modeling.

The study also opens avenues for cross-disciplinary collaboration, integrating geological fieldwork, remote sensing data, and advanced computational mechanics. The authors highlight the importance of combining micro-scale material characterization with macro-scale numerical simulations, bridging gaps between laboratory analysis and field-scale hazard prediction. This multi-scale methodology is poised to become a standard approach in the future of landslide research.

As computational power continues to grow, so too does the potential for increasingly sophisticated and high-resolution models. The numerical framework pioneered in this study is readily adaptable to incorporate further complexities, such as pore-water pressure effects, vegetation influence, and evolving topography during landslide progression. Each of these factors can interplay with anisotropy to shape landslide dynamics, underscoring the ongoing need for integrated modeling efforts.

In conclusion, the groundbreaking research by Chang and his colleagues ushers in a new paradigm in landslide science. By illuminating the intricate effects of anisotropic material properties on runout behavior, they have advanced our ability to model, predict, and ultimately mitigate one of nature’s most destructive forces. These findings not only deepen scientific understanding but also provide critical tools for safeguarding vulnerable communities in an era of environmental uncertainty.

As landslide hazards continue to exact heavy tolls worldwide, innovations such as these hold the promise of more resilient future landscapes. The integration of anisotropic considerations into numerical modeling marks a transformative milestone, elevating both the accuracy and applicability of landslide risk assessments. It is a vivid example of how scientific insight paired with technological prowess can pave the way for safer coexistence with Earth’s dynamic and sometimes perilous terrain.

Subject of Research: Landslide runout behavior on anisotropic slopes through numerical modeling.

Article Title: Numerical modeling of landslide runout behavior for an anisotropic slope.

Article References:
Chang, KT., Lu, CA., Yeh, PT. et al. Numerical modeling of landslide runout behavior for an anisotropic slope. Environ Earth Sci 84, 407 (2025). https://doi.org/10.1007/s12665-025-12403-0

Image Credits: AI Generated