The integration of bacteria-based self-healing mortars has emerged as a promising solution to address repair due to recurring cracks and preserving masonry durability. Building upon a recent pilot study demonstrating the efficacy of a self-healing agent in the repair of masonry made with cement-based mortar, this follow-up study explores the potential of integrating the added-in healing agent in a pre-bagged cement-lime mortar - more commonly used in masonry applications. Through bond wrench tests and a 30-day healing period involving wet-dry cycles, the study evaluates aesthetic and flexural bond strength recovery of couplets built with solid clay bricks. Results showed that the addition of the agent altered the initial flexural bond strength, with bacteria-based masonry couplets four times stronger than the plain reference ones - without containing the agent. The mortar’s colorwas also affected. Additionally, bacteria-based specimens demonstrated automatic repair, restoring up to 33% of the original flexural bond strength, while referencemasonry couplets showed no evidence of autonomous healing. However, instances of leaching, possibly attributed to the agent’s substrate, prompted a revision of the strategy employed for the healing environment. Further research will specifically target the observed leaching issue by exploring the effects of multiple healing environments.

Subsidence caused by natural or human-induced factors can occur unevenly, resulting in differential settlements. Existing unreinforced masonry (URM) buildings are susceptible to damage from differential settlements. However, the extent of the damage varies between structures, depending on factors such as the magnitude and pattern of the settlements, along with the features of the building and the properties of the underlying soil. Non-linear finite element analyses (NLFEA) are often used for studying the damage response, accounting for variability in soil and structural features. This study uses 6912 NLFEA, including 8 variations in façade geometry, 3 masonry materials, 2 soils, 2 shallow foundation systems, and 72 settlement patterns, to develop fragility curves for URM buildings undergoing subsidence. Old Dutch URM buildings with strip foundations are modelled using 2D plane-stress façade models, accounting for non-linear smeared shearing, cracking and crushing of masonry and 3D effects of transverse walls. Settlement troughs are applied at a non-linear soil-foundation interface, with angular distortion (β) progressively increasing to quantify settlement intensity and building deformation. As β increases, the NLFE models exhibit progressive cracking damage, with severity objectively assessed through the parameter Ψ considering crack width, length, and number. Then, the distortion β is used as the demand parameter to develop the fragility curves. The analysis shows that long façades are twice as likely to experience 5 mm cracks from settlement damage compared to short façades under an applied β of 2 ‰ (1/500). For this applied β, proposed as an acceptable limit for many structures in the Eurocode, half of the models exhibit cracks up to 5 mm wide. Therefore, while 1/500 may be considered safe for structural integrity, it can still lead to noticeable damage. Light damage occurs even at angular distortion values below 0.5 ‰ (1/2000), with 10 % of models showing cracks up to 1 mm wide.

The structural response of masonry walls during flood events is a critical concern for the flood resilience of (Dutch) buildings, as they typically constitute part of the load-bearing structure. This study investigates the out-of-plane behaviour of a full-scale single-wythe fired-clay-brick masonry wall under out-of-plane hydrostatic pressure and debris impact loads. Experimental tests were conducted on a 2.7 × 2.7 m masonry wall subjected to a vertical pre-compression and simultaneously varying water levels and debris impacts at the Flood Proof Holland facility in Delft, the Netherlands. Results demonstrated that the wall remained within the linear-elastic regime up to a water depth of approximately 90 cm when the interior side was dry. Beyond this threshold, crack initiation and stress redistribution occurred, leading to significant deformation. On the basis of calibrated models, failure was predicted at approximately 150 cm water depth for a fully restrained wall. Debris impact tests showed that soft debris, represented by a floating log, caused negligible additional damage, whereas repeated impacts with a steel cube (hard debris) resulted in progressive cracking and local failure, particularly at higher water levels. Numerical models, including analytical, linear-elastic finite element method (FEM), and non-linear FE approaches, were calibrated using the experimental data. While one-way bending models predicted conservative failure thresholds, two-way, non-linear models accurately captured the wall’s deformation and cracking behaviour, demonstrating the importance of lateral boundary constraints in determining wall capacity and stability. The findings emphasise that traditional masonry walls in Dutch buildings can safely withstand water depths up to 90 cm without significant damage. However, higher water levels or hard debris impacts pose substantial risks, highlighting the need for improved flood resilience strategies. Future work should focus on cavity wall systems, leakage effects, and the behaviour of walls with openings.

Cycles of settlement and uplift beneath existing masonry structures can lead to visible cracks, which not only affect the aesthetic appearance and functionality of the building but can also compromise its structural integrity and undermine the occupants' sense of safety. These cyclic ground movements can be triggered by seasonal actions, such as fluctuation in the groundwater table. In the Netherlands, many existing masonry structures on shallow foundations rest directly on the subsurface, making them vulnerable to cyclic ground movements. Settlement and uplift cycles cause “breathing” masonry cracks, which open and close over time without fully sealing. This study uses finite element analyses to investigate and assess the damage of structures subjected to cyclic quasi-static ground movements. A case study is presented for the analysis, featuring the geometry of an existing low-rise masonry structure with an age exceeding 50 years. A 3D non-linear shell-element model is used to evaluate the structural response, featuring an unreinforced strip foundation and including the non-linear tensile softening and cracking behaviour of masonry. Heaving and sinking displacements are applied to a non-linear interface simulating the soil-foundation interaction at the bottom of the strip foundation. The intensity of the ground displacements is quantified by their angular distortion. A damage parameter objectively assesses the severity of damage by considering the number, length, and width of cracks. Results indicate that repeated cycles of settlement (and uplift) have been observed to cause irreversible cracking damage in the model, with crack widths ranging from 1 to 5 mm, progressively increasing over time. Damage occurring during settlement is, on average, twice as severe as that during uplift. Overall, cycles of settlement and uplift may induce cracking damage up to twice as high as that caused by cycles of settlement alone, depending on the magnitude and shape of the ground movements.

Masonry buildings in the Netherlands are especially prone to damage in the form of small cracks. This is because the masonry is unreinforced, the foundations are shallow and often also unreinforced, the bedding is composed of soft soils like peat or clay, dilation joints are missing in older or historical structures, and current loading conditions, such as earthquake vibrations, were never considered in the design of the buildings. The latter includes mining operations for salt and gas that have led to subsidence and induced seismicity. Moreover, farming policy and water management, in combination with regional subsidence, have led to varying groundwater table levels which, in turn, cause wetting and drying of sensitive soils. This process is exacerbated by more extreme seasons of precipitation and drought because of climate change, leading to swelling and compaction of the ground underneath buildings. To understand building damage in this context, it is necessary to evaluate the combined effects of these various hazards. Their actions can be decomposed into vibrations caused by earthquakes and ground deformations. The former can be characterized by the PGV or PGA of the vibrations, and the latter by the induced curvature of the soil surface and/or by the horizontal strains at the surface because of deformations deep in the underground. Moreover, repeated earthquake events and seasonal soil subsidence or heave lead to cyclic actions. The contribution and interaction of these loads causing progressive damage to masonry buildings have been the focus of an extensive modelling study with detailed non-linear models of the buildings and the soil. The slow soil deformations were analyzed first and served as the starting point for subsequent, repeated vibrations. For example, a horizontal strain of 0.1 mm/m caused by mining, in combination with an angular distortion of 1/2000 due to local soil compaction, can produce cracks of about 1 to 2 mm wide in a particular masonry façade. The damage is then aggravated by an earthquake vibration in the order of 5 mm/s, which is further increased by about 10% with a repeated event. The expected final damage may include multiple cracks of up to 3 mm. In this manner, the combination of all actions can lead to the establishment of conservative thresholds to prevent or limit damage to existing structures.

Within the context of light damage to unreinforced masonry structures, recent tests have shown that the cracking behaviour of calcium-silicate brick masonry walls makes them more vulnerable to in-plane loads when compared against fired-clay brick walls. To further explore this observation, four nominallyidentical walls have been tested. Two of the specimens (3m wide and 2.7m tall) were built with calcium-silicate bricks and twowith fired-clay bricks. Additionally, two boundaries were compared: a top cantilever boundary, and a doubled-clamped configuration. The quasi-static, in-plane, two-way cyclic tests imposed small (0.03 to 0.1%), repeated drifts on the walls to investigate the initiation and propagation of small cracks. To monitor the cracking behaviour, high resolution Digital Image Correlation was applied. At the end of the tests, large drifts up to 2% were exerted to compare the near-collapse behaviour of the walls. The tests revealed that thewalls with the more restrictive boundary, deforming mostly in shear, behaved the stiffest and also developed cracks earlier than the cantilever walls. Additionally, this constraint also led to more vertical cracks that split bricks, while the cantilever walls saw more horizontal and diagonal cracks along mortar joints and at mortar-brick interfaces. While the calcium-silicate-proved to be more brittle than the fired-clay masonry for the cantilever test, the claymasonry exhibited similar brick-splitting cracks in the double-clamped configuration. In general, there were fewer but wider cracks in the calcium-silicate specimens, while the clay brick samples showed less localisation and more smeared-crack behaviour. In terms of stiffness, the calcium-silicate walls were initially stiffer and achieved a higher capacity. Moreover, these walls also presented a higher hysteresis associated with more frictional failures. In sum, while cracks on the calcium-silicate walls were confirmed to be more serious, their increased stiffness could lead to smaller drifts during dynamic loading, and the walls would develop less damage; this requires further study.

Terraced buildings with cavity walls are among the most common types of construction in the northern part of the Netherlands. Since 1980, the inner walls of these buildings have been constructed using either calcium silicate bricks (214 × 102 × 75 mm) with thick mortar joints (10 mm) or, more recently, calcium silicate blocks (437 × 198 × 100 mm) with thin mortar joints (3 mm). The shear properties of these units play a crucial role in the seismic response of buildings, particularly in regions like Groningen, which is prone to seismic activity due to artificial extraction. This study investigates the shear interface behavior of these two types of masonry units by testing multiple triplet samples under varying levels of normal stress at the interface. The results provide detailed insights into the shear properties of both brick and block masonry, offering valuable data for enhancing the accuracy of numerical simulations and predicting the structural capacity of these types of masonry buildings.

Soil heterogeneity, due to variations in the subsurface stratigraphy or properties within a layer, can trigger or amplify differential settlements that affect buildings and infrastructure and can thus lead to (increase in) damage. The state-of-the-art mainly focuses on the effect of heterogeneous properties within a layer on engineering problems. From this, it is known that the variation in properties can increase the vulnerability of a structure. However, nearly always variations in the soil lithological conditions are disregarded, while they can influence subsidence potentially even more. Lithological variations are relevant both at the scale of individual buildings as well as different scales (city, regional, country), for which often detailed soil information is not available. Thus, for a better prediction of potential building damage related to subsidence, knowledge about the scale and influence of lithological variations is needed. This paper describes an approach to quantify and investigate the influence of lithological heterogeneity at the scale of a single building. Moreover, this exploratory study evaluates the influence of lithological heterogeneity on the spatial variability of settlements, intending to upscale the approach to regional application. Two independent datasets at high resolution (site-specific) and low resolution (national level) are used to retrieve the stratigraphic conditions for the area selected for the analyses. One-, Two- and Three-dimensional numerical models, based on the collected information are used to simulate the consolidation process and settlement due to a uniform load imposed on the surface level of the study area. Additional analyses investigate the influence of loading conditions and groundwater table. The parameter “correlation length” is used to quantify the spatial variability of the soil layer thickness and then of the computed settlements. The analyses reveal that the spatial variability of the soil strata thickness matches that of the computed settlements, ranging from 2 to 10 meters. In other words, the lithological variability of the soil leads to differential settlements occurring at the scale of man-made structures such as houses, roads, and embankments. Thus, the results encourage including the contribution of lithological heterogeneity in models and predictions of differential settlement at the scale of individual structures. Moreover, the statistical properties, in terms of mean, spread and distribution shape, of the settlement computed through in-situ specific models, match with those derived at the national scale. These results are expected to support the identification of areas potentially influenced by lithological soil heterogeneity, thus showing potential for upscaling to regional or national levels.

This study presents a semi-automated, data-informed framework for selecting parameter-consistent numerical models to approximate the in-plane behaviour of clay and calcium silicate masonry walls. A comprehensive experimental campaign has been executed on full-scale unreinforced calcium silicate and clay masonry walls at Delft University of Technology. The in-plane response of these walls was evaluated based on stiffness, strength, damage intensity at equivalent drift levels, and the overall impact of the damage. The findings indicate that unreinforced calcium silicate masonry walls are more prone to damage through the brick units, while cracks in unreinforced clay masonry walls predominantly align with mortar joints. Calcium silicate walls tend to develop larger and more prominent cracks, often requiring the replacement of individual bricks for a complete repair. In contrast, the damage in clay walls is typically easier to address through repointing of mortar joints. A parametric finite element analysis was performed to investigate these failure mechanisms, systematically varying input parameters to generate 3,456 numerical simulations. Each model permutation was evaluated to select models that closely approximate observed experimental responses. Unlike conventional calibration methods, this framework systematically explores possible combinations of input and output parameters to identify numerical models that replicate key structural behaviours. The preliminary results demonstrate that multiple parameter combinations can yield numerical responses closely matching experimental observations, providing a structured approach for improving masonry modelling practices.

Damage assessment for masonry structures subjected to settlement is crucial for ensuring structural safety, guiding repairs, and preserving the built environment. Non-linear finite element modelling offers an effective approach for this purpose, though balancing model complexity, computational cost, and predictive reliability remains a key challenge. This study addresses the absence of a systematic comparison between macro- and simplified micro-modelling strategies for such analyses, clarifying their respective strengths, limitations, and sensitivity to key parameters. The performance and accuracy of semi-coupled NLFEM models are compared in simulating the response of a 1/10th scaled masonry façade under settlement, available from prior research. The two approaches considered are: simplified micro-modelling, where bricks are represented as expanded blocks with non-linear interfaces for mortar joints and their contact edges, and macro-modelling, where masonry is homogenised into an equivalent orthotropic composite material. The macro-models employ two well-established constitutive models, the Total Strain Rotating Crack Model (TSRCM) and the Engineering Masonry Model (EMM), to capture the non-linear cracking behaviour of masonry. Sensitivity analyses assess the influence of base interface models and the interface’s tangential stiffness. The results show how the selection of the modelling approach depends on the analysis objective: The macro-model with the Engineering Masonry Model best predicts damage severity, deviating by only 10% from the experiment, further improved by calibrating the minimum head-joint tensile strength. While all models yield similar predictions for vertical displacements of the façade, the TSRCM better captures overall and horizontal displacements, whereas the simplified micro-model more accurately represents the crack pattern. The EMM-based macro-models are the most computationally efficient, with TSRCM requiring 1.5 times the CPU time of EMM, and the micro-model requiring twice as much. The analysis also shows that the TSRCM-based macro-model is more sensitive to variations in the type of base interface models and base interface tangential stiffness, convergence criteria, incremental-iterative procedure, and analysis settings, whereas the EMM macro-model and the simplified micro-model are less affected. By identifying the strengths and limitations of each modelling approach, this study supports informed modelling choices for a more reliable assessment of settlement damage, contributing to the effective protection of existing masonry structures.

Cracks represent a prevalent form of damage in masonry structures, posing not only aesthetic concerns but also compromising structural durability; therefore, they are undesirable and need to be repaired. The repointing technique is traditionally implemented in this context, especially in historical masonry. However, this method fails to provide a long-term solution, leaving structures vulnerable to future damage. The paper investigates the applicability of a bio-based self-healing mortar to enable autonomous repair of masonry. This innovative mortar, developed to repair concrete structures, was implemented to explore the capacity of couplets to recover their original bond capacity and aesthetic aspect after multiple damaging events. Specimens built with calciumsilicate and clay bricks were subjected to subsequent cracking cycles using a crack-mouth-opening -displacement controlled bond-wrench test. Experimental results showed that self-repair, in terms of bond restoration and aesthetic filling of cracks, occurs even after multiple cracking cycles when the bio-based mortar is used with both types of bricks, optimizing the autogenous healing (intrinsic) of cement-based mortars. The effectiveness varied also according to the types of brick and healing environment used, e.g. under humid conditions (RH ~ 95%), 50% vs 80% of the original capacity was regained in fully separated couplets made respectively with clay and calcium-silicate bricks.

In this study, 2D and 3D modelling strategies are used to represent the behaviour of historical masonry buildings on strip foundations undergoing settlements. The application focuses on a two-story building, typical of the Dutch architectural heritage. An improved 2D modelling is presented: It includes the effect of the lateral walls to replicate the response of the detailed 3D models. The masonry strip foundation is modelled and supported by a no-tension interface, which represents the soil-foundation interaction. Two settlement configurations, hogging and sagging, are applied to the models, and their intensity is characterized using their angular distortion. The improved 2D model that includes the stiffness and weight of the lateral walls agrees in terms of displacements, stress and damage with the detailed 3D models. Conversely, the simplified 2D façade models without lateral walls exhibit different cracking, and damage from 2 to 7 times lower at an applied angular distortion of 2‰ (1/500). The improved 2D model requires less computational and modelling burden, resulting in analyses from 9 to 40 times faster than the 3D models. The results prove the importance of identifying and including the 3D effects that affect the response of structures subjected to settlements.

The July 2021 flood heavily affected many inhabitants, buildings and critical infrastructure throughout Germany, Belgium and the Netherlands. Specifically, the Ahr Valley (Germany) showcased the destructive power associated with these extreme events. Hence, this region was the focus of a field survey, aiming at describing the flood-induced damage to buildings and assessing the possible underlying processes that led to structural failures. The field assessment revealed a close connection between building failures and (1) local flow depths and velocities, (2) building location, (3) distance from the riverbank and (4) construction type. Although it is difficult to identify the exact causes that induced failures, the detailed assessment revealed that damages mainly originated from local scour and hydraulic loads, often unevenly distributed around buildings. Importantly, many buildings were significantly affected by (large) floating debris impacts and damming, both responsible for additional loads, highlighting their importance in flood-resistant building design. Furthermore, data showed that buildings near the riverbanks and in the upstream part of villages were more severely damaged. Altogether, data provide a better understanding of the flood processes that lead to building failures, fostering future research towards the development of safer protection measures and more effective flood risk management strategies.

This report details the data verification process for the information utilized in the publica- tion ”Empirical fragility and ROC curves for masonry buildings subjected to settlements” [3]. The data are stored in the 4TU.ResearchData repository [4]. Data collection and anal- ysis took place between 2020 and 2021, with the final version being published in 2023. A thorough check was conducted before and after the publication of the manuscript to ensure data quality. While the first check was conducted by the Authors of [3] during the initial data collection, this report focuses on the second check, carried out after the pub- lication. In this document, the adopted methodology for the data verification is detailed and the results of the procedure are presented.

This report begins with Section 2, which outlines the methodology used for the data verification, whereas Section 3 summarises the results. Section 4 offers a discussion and conclusions. Finally, appendix A contains the calculations that support the verification checks.

In the Netherlands, subsidence due to different causes is linked to damage to the ubiquitous masonry structures. Finite element (FE) analyses can be used to assess the response of the structures subjected to settlements. This paper presents the comparison between three-dimensional FE modelling strategies to investigate the response of an unreinforced masonry building on a strip foundation. The aim is to investigate whether different modelling approaches demonstrate consistent results. The soil-structure system is modelled employing two strategies: a coupled model, in which the structure is tied to the soil volume, and an uncoupled approach that divides the soil and structure into two sub-systems. Two displacement fields, imposed at the bottom of the soil volume, idealize various shapes of the subsidence troughs, with increasing intensity measured by their distortion. Non-linear interfaces are used to simulate the soil-foundation interaction, and their stiffness values vary based on the type of model. The displacements, interface stresses and crack patterns of the selected modelling strategies are consistent. The interface types do not influence the response of the façade, whereas the shape of the settlement does play a key role. The uncoupled models exhibit, on average, slightly higher values of damage than coupled models for a given imposed distortion. The two modelling strategies require almost the same computational time and show similar convergence. Because of the limited contribution of small soil volumes in uncoupled models, the superstructure sub-system can be directly utilized to assess the response of structures undergoing vertical displacements, thereby reducing the modelling burden.

Drift limits are useful thresholds; during design or retrofitting analyses, engineers can compare the expected behaviour of a structure to drift limits that predict when the structure will reach a certain condition. This helps ensure that structures satisfy specified performance goals when exposed to certain hazards. Masonry walls are susceptible to damage from lateral in-plane actions such as wind or earthquake loading; ensuring that in-plane drift remains sufficiently small will help limit this damage. Drift limits based on crack-based damage are scarce, however, with DS1 limits being extrapolated from higher damage grades based on structural strength capacity or ductility. In this work, crack-based damage is evaluated on a multitude of full-scale experimental walls surveyed with digital image correlation. This method observes the initiation and propagation of cracking. Cyclically incremental in-plane tests provide a range of drift-damage relationships. These are explored with machine learning to determine influential predictors and ultimately establish drift limits for light damage. Two types of brick masonry are explored: fired-clay and calcium-silicate. For the latter, light damage begins at an in-plane drift of 0.5 mm/m and can extend to 4.8 mm/m (or 0.48%) for the former before the masonry surpasses light damage and reaches structural damage grades. In comparison to drift limits set by other authors and (international) guidelines to characterise light damage, significant damage, or the ultimate capacity of masonry walls, the resulting drift limits for light damage from this work are set directly on the basis of experiments and are in good agreement with other authors. Most importantly, all the consulted values for ultimate capacity are much larger than the upper threshold for light damage determined herein, with limits for significant damage in the same order of magnitude. This result verifies the accuracy of the experimental crack-based characterisation used to establish the drift thresholds.

To quantify damage to unreinforced masonry structures, sensitive to tensile stresses and thus to vibrations caused by induced earthquakes in the north of the Netherlands, a bottom-up or physical approach has been followed. First, the crack-based damage expected on these common structures, was defined on a continuous scale (Ψ) based on width, length and number of cracks, so that it could also be quantified. Then, the initiation of cracks was investigated on full-scale masonry walls, single-wythe and of fired-clay bricks, using high resolution Digital Image Correlation; at 0.1 mm in width, cracks become visible. Further, the propagation of cracks, as they widen and lengthen, was monitored during repeated and cyclic testing. At repeated, identical in-plane drift, cracks still propagate and damage increases, though the increase is minor.

The experiments on walls and spandrels, displaying horizontal, diagonal, and vertical cracks, were used to calibrate numerical, finite-element-method models. These reproduced the behaviour of the tests also in terms of crack patterns and propagation, besides stiffness, strength, and hysteresis. Then, the models were adapted to explore the effect of earthquake vibrations, also in combination with existing damage caused by settlement-like actions. In this manner, the effect of initial damage could be quantified. Several other parameters were varied, such as the material strength, the geometry of the masonry walls, the soil properties in a soil-structure interaction interface, the record, PGV, and repetition of the vibrations, and the intensity of initial damage.

Relationships between these parameters and damage were captured into a surrogate model, which was used in a MonteCarlo simulation to determine the probability of damage. The fragility curves reveal, for instance, that fired-clay brick walls with no visible pre-damage (Ψ0=0) have a 5% chance of visible damage (Ψ≥1) at a PGV of, for example, 10 mm/s, a probability that rose to 20% if the walls had undetectable pre-existing damage (Ψ0=0.5). The probability of exceeding aesthetic damage (Ψ≈2.5) at this PGV is less than 1%. A lower PGV is associated with a smaller probability of damage.

Furthermore, it was concluded that repeated events lead to an increase in damage of about 10% for five similar events. This increase is not noticeable. In a sequence of events, similar events accumulate little damage and increases appear when larger events are experienced by a masonry structure.

In the pursuit of introducing bacteria-based self-healing mortar for masonry repair, this study examined the potential of incorporating a poly-lactic acid (PLA) agent—already established in concrete repair—into cement-lime mortars, typical of historical constructions. Testing prisms constructed with varying lime/cement ratios revealed decreased flexural and compressive strength in high-cement-concentration mortars upon the addition of the agent; for mortars with high lime concentration, however, the agent led to an increase in both strengths. Furthermore, the agent's potential to self-repair was confirmed by allowing the remaining portions of tested samples to heal under humid conditions. Irrespective of mortar composition, cracks were resealed thus confirming the aesthetic, and potentially watertightness, restoration.

This paper demonstrates the use of non-linear finite element modelling to investigate the response of structures subjected to different shapes of subsidence-related ground settlements. The approach is presented with reference to a two-storey unreinforced masonry façade resting on a shallow foundation. Eight realistic settlement shapes, based on field and literature data, are applied in the model with increasing intensity. The intensity of the subsidence profiles is characterized using their (angular) distortion. The extent of the induced damage on the façade is objectively and directly quantified by a damage parameter, based on the number of cracks, their length and opening. The performance of different settlement indicators and corresponding limiting values, typically employed in the state of the art, is in this paper discussed in relation to the damage modelling strategy; these are observed to be dependent on the shape of the settlement profiles. The aim of this paper is thus to provide insight into the extent to which the vulnerability of masonry buildings depends on the shape of the subsidence pattern and may serve as a warning not to use (deterministic) damage indicators such as angular distortion without considering the settlement shape.