{'fraction_non_alphanumeric': 0.026804860614724802, 'fraction_numerical': 0.006433166547533953, 'mean_word_length': 5.361363636363636, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Discrete Element Modeling of Dike-induced Deformation

Introduction: The Tharsis region of Mars is characterized by large volcanic and tectonic centers with distinct sets of graben systems. Many of the radially oriented grabens have been inferred to form in response to intrusion of magmatic dikes. This interpretation is based primarily upon early physical and numerical (boundary element) models that were developed originally to understand surface deformation associated with dike emplace-ment on Earth. In this study, we constructed and analyzed two-dimensional discrete element models to test the hypothesis of shallow dike emplacement and widening as a primary mechanism for the production of grabens on Mars. In particular, our models are designed to explore the extent to which a widening subsurface dike under varying emplacement conditions will induce near-surface and surface deformation. Methodology: The use of discrete element models allows for permanent deformation and material heteroge-neity to be simulated. In contrast, boundary element models are limited by an assumption of homogeneous elastic behavior. We performed discrete element model-ing of dike-related deformation using the code PFC2D [1]. PFC2D is based on the discrete element method that was developed initially for analyzing the mechanical interaction of granular materials [2]. With the addition of particle bonding, the method has been extended to a variety of problems in solid mechanics, including rock fracturing and seismicity [3,4,5,6,7]. The discrete element formulation requires selection of a set of micromechanical properties (e.g., friction, normal and shear bond stiffness and strength) that describe the interaction of the elastic particles with each other and the model boundaries (see [4], for details on the theoretical framework of PFC). The result is a set of micromechanical properties that provides the correct bulk (macroscale) material behavior, which includes both the elastic (recoverable) and inelastic (non-recoverable) components characteristic of natural deformation. For all models, stratigraphic material was modeled by generating particles in an initial model area 10 km wide. The particles have a diameter range of 12.8 – 42.9 m scaled to the kilometer scale model boundaries (smaller particle size ranges generate billions of particles which overwhelm the computations). A half model space was used to conserve computation time. The particles are assigned uniform stiffness and friction values throughout the model area. Blue and red colored layering was added for visualization purposes only and does not affect material behavior. Once the model reached equilibrium (no large out-of-balance forces), the simulated dike was widened. The dike was assumed to have reached a neutral …