Formation of todorokite from vernadite in Ni-rich hemipelagic sediments

Todorokite is considered to form from vernadite in nature and commonly concentrates nickel. However, this mineralogical transformation has never been imaged nor explained mechanistically, and its effect on the uptake of nickel has never been quantified at the molecular-level. We have characterized these reactions at the macroscopic, microscopic, nanoscopic and atomic scales in a marine manganese concretion by combining transmission electron microscopy, electron and X-ray microprobe analysis, powder and micro X-ray diffraction, and Mn and Ni K-edge EXAFS spectroscopy. The concretion was collected during the Ticoflux II expedition near the Nicoya Peninsula, Costa Rica, and is representative of Mn deposits in hemipelagic sediments. It consists of 5 to 25 μm aggregates, shaped like sea-urchins, with a core of 7Å-vernadite (1.0 wt% Ni), a rim of 10Å-vernadite (3.8 wt% Ni), and an outermost region of todorokite fibers (1.9 wt% Ni) that extend outwards. The crystallites of 7Å-vernadite are single- to bi-layered, with hexagonal layer symmetry (a = b = 2.83 Å), and an average structural formula of . The crystallites of 10Å-vernadite contain 10 to 20 layers semi-coherently stacked in the ab plane and uniformly separated in the [0 0 1] direction by ∼9 Å due to the intercalation of hydrated Mg²⁺ cations. The average structural formula of 10Å-vernadite is if the layers contain vacancy sites, or alternately , if they contain Mn³⁺. The average formula of todorokite is .A genetic model is proposed based on combining these new data with previously published results. The thermodynamically unstable 7Å-vernadite transforms via dissolution-recrystallization to semi-ordered Mg-rich 10Å-vernadite. Nickel is released from dissolved biogenic silica or reduced organic matter, and taken up mainly in the Mn layer of 10Å-vernadite. Interlayer magnesium serves as a template to the further topotactic transformation of 10Å-vernadite to todorokite. The dimension of the todorokite tunnels in the [0 0 1] direction is uniform and determined by the size of the hydrated Mg²⁺ ion (8.6 Å). The tunnel dimension in the [1 0 0] direction depends on the density of Mg²⁺ in the interlayer and the superstructure of the phyllomanganate layer. If the parent phyllomanganate contains high amounts of Mg²⁺ (i.e., high layer charge), or Mn³⁺ and Mn⁴⁺ cations ordered following the Mn³⁺-Mn⁴⁺-Mn⁴⁺ sequence as in synthetic triclinic birnessite, then the tunnel dimension is ideally 3 × 3 octahedral chain widths in both crystallographic directions. Otherwise, the tunnel dimension is incoherent in the [1 0 0] direction (i.e., T(3,n) tunnel structure), as has been observed in all natural todorokites. Natural todorokite is defective because the precursor natural phyllomanganates either have a layer charge deficit below 0.33e per octahedral site, or rarely are triclinic birnessite. The abundance of Mg in seawater and its key role in converting phyllomanganate to tectomanganate with T(3,n) tunnel structure explain why todorokite is common in marine ferromanganese oxides, and seldom present in terrestrial environments. The topotactic phase transformation described here is the only known route to todorokite crystallization. This implies that all natural todorokites may be authigenic because they are formed in situ from a phyllomanganate precursor.     

Confirmation of mass-independent Ni isotopic variability in iron meteorites

We report high-precision analyses of internally-normalised Ni isotope ratios in 12 bulk iron meteorites. Our measurements of ⁶⁰Ni/⁶¹Ni, ⁶²Ni/⁶¹Ni and ⁶⁴Ni/⁶¹Ni normalised to ⁵⁸Ni/⁶¹Ni and expressed in parts per ten thousand (?) relative to NIST SRM 986 as and , vary by 0.146, 0.228 and 0.687, respectively. The precision on a typical analysis is 0.03?, 0.05? and 0.08? for , and , respectively, which is comparable to our sample reproducibility. We show that this ‘mass-independent’ Ni isotope variability cannot be ascribed to interferences, inaccurate correction of instrumental or natural mass-dependent fractionation, fractionation controlled by nuclear field shift effects, nor the influence of cosmic ray spallation. These results thus document the presence of mass-independent Ni isotopic heterogeneity in bulk meteoritic samples, as previously proposed by Regelous et al. (2008) (EPSL 272, 330-338), but our new analyses are more precise and include determination of ⁶⁴Ni. Intriguingly, we find that terrestrial materials do not yield homogenous internally-normalised Ni isotope compositions, which, as pointed out by Young et al. (2002) (GCA 66, 1095-1104), may be the expected result of using the exponential (kinetic) law and atomic masses to normalise all fractionation processes. The certified Ni isotope reference material NIST SRM 986 defines zero in this study, while appropriate ratios for the bulk silicate Earth are given by the peridotites JP-1 and DTS-2 and, relative to NIST SRM 986, yield deviations in , and of −0.006?, 0.036? and 0.119?, respectively. There is a strong positive correlation between and in iron meteorites analyses, with a slope of 3.03 ± 0.71. The variations of Ni isotope anomalies in iron meteorites are consistent with heterogeneous distribution of a nucleosynthetic component from a type Ia supernova into the proto-solar nebula.     

Coupled mobilization of dissolved organic matter and metals (Cu and Zn) in soil columns

Dissolved organic carbon (DOC) is a key component involved in metal displacement in soils. In this study, we investigated the concentration profiles of soil-borne DOC, Cu and Zn at various irrigation rates with synthetic rain water under quasi steady-state conditions, using repacked soil columns with a metal-polluted topsoil and two unpolluted subsoils. Soil solution was collected using suction cups installed at centimeter intervals over depth. In the topsoil the concentrations of DOC, dissolved metals (Zn and Cu), major cations (Ca²⁺ and Mg²⁺) and anions ( and ) increased with depth. In the subsoil, the Cu and Zn concentrations dropped to background levels within 2 cm. All compounds were much faster mobilized in the first 4 cm than in the rest of the topsoil. DOC and Cu concentrations were higher at higher flow rates for a given depth, whereas the concentrations of the other ions decreased with increasing flow rate. The decomposition of soil organic matter resulted in the formation of DOC, , and and was the main driver of the system. Regression analysis indicated that Cu mobilization was governed by DOC, whereas Zn mobilization was primarily determined by Ca and to a lesser extent by DOC. Labile Zn and Cu²⁺ concentrations were well predicted by the NICA-Donnan model. The results highlight the value of high-resolution in-situ measurements of DOC and metal mobilization in soil profiles.     

First-principles estimates of equilibrium magnesium isotope fractionation in silicate, oxide, carbonate and hexaaquamagnesium(2+) crystals

Equilibrium mass-dependent magnesium isotope fractionation factors are estimated for a range of crystalline compounds including oxides, silicates, carbonates, and salts containing the solvation complex. Fractionation factors for the gas-phase species Mg and MgO are also presented. Fractionation factors are calculated with density functional perturbation theory (DFPT), using norm-conserving pseudopotentials. The results suggest that there will be substantial inter-mineral fractionation, particularly between tetrahedrally coordinated Mg²⁺ in spinel (MgAl₂O₄) and the more common octahedrally coordinated Mg²⁺-sites in silicate and carbonate minerals. Isotope fractionations calculated for Mg²⁺ in hexaaquamagnesium(2+) salts are in good agreement with previous fractionation models of based on large molecular clusters (Black et al., 2007), but show possibly more significant disagreement with a more recent study (Rustad et al., 2010). These models further suggest that solvated , in the form of , will have higher ²⁶Mg/²⁴Mg than coexisting magnesite and dolomite. Calculated fractionations are consistent with Mg-isotope fractionations observed in peridotite mineral separates and inorganic carbonate precipitates. Predicted large, temperature-sensitive spinel-silicate fractionations, in particular, may find use in determining equilibration temperatures of peridotites and other high-temperature rock types.     

Flow and nutrient dynamics in a subterranean estuary (Waquoit Bay, MA, USA): Field data and reactive transport modeling

A two-dimensional (2D) reactive transport model is used to investigate the controls on nutrient (, , PO₄) dynamics in a coastal aquifer. The model couples density-dependent flow to a reaction network which includes oxic degradation of organic matter, denitrification, iron oxide reduction, nitrification, Fe²⁺ oxidation and sorption of PO₄ onto iron oxides. Porewater measurements from a well transect at Waquoit Bay, MA, USA indicate the presence of a reducing plume with high Fe²⁺, , DOC (dissolved organic carbon) and PO₄ concentrations overlying a more oxidizing -rich plume. These two plumes travel nearly conservatively until they start to overlap in the intertidal coastal sediments prior to discharge into the bay. In this zone, the aeration of the surface beach sediments drives nitrification and allows the precipitation of iron oxide, which leads to the removal of PO₄ through sorption. Model simulations suggest that removal of through denitrification is inhibited by the limited overlap between the two freshwater plumes, as well as by the refractory nature of terrestrial DOC. Submarine groundwater discharge is a significant source of to the bay.     

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

The effect of sulfur on the partitioning of Ni and other first-row transition elements between olivine and silicate melt

The effect of sulfur dissolved as sulfide (S²⁻) in silicate melts on the activity coefficients of NiO and some other oxides of divalent cations (Ca, Cr, Mn, Fe and Co) has been determined from olivine/melt partitioning experiments at 1400 °C in six melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), and in derivatives of these compositions at 1370 °C, obtained from the six CMAS compositions by substituting Fe for Mg (FeCMAS). Amounts of S²⁻ were varied from zero to sulfide saturation, reaching 4100 μg g⁻¹ S in the most sulfur-rich silicate melt. The sulfide solubilities compare reasonably well with those predicted from the parameterization of the sulfide capacity of silicate melts at 1400 °C of O’Neill and Mavrogenes (2002), although in detail systematic deviations indicate that a more sophisticated model may improve the prediction of sulfide capacities.The results show a barely discernible effect of S²⁻ in the silicate melt on Fe, Co and Ni partition coefficients, and also surprisingly, a tiny but resolvable effect on Ca partitioning, but no detectable effect on Cr, Mn or some other lithophile incompatible elements (Sc, Ti, V, Y, Zr and Hf). Decreasing Mg# of olivine (reflecting increasing FeO in the system) has a significant influence on the partitioning of several of the divalent cations, particularly Ca and Ni. We find a remarkably systematic correlation between and the ionic radius of M²⁺, where M = Ca, Cr, Mn, Fe, Co or Ni, which is attributable to a simple relationship between size mismatch and excess free energies of mixing in Mg-rich olivine solid solutions.Neither the effect of S²⁻ nor of Mg#^ol is large enough by an order of magnitude to account for the reported variations of obtained from electron microprobe analyses of olivine/glass pairs from mid-ocean ridge basalts (MORBs). Comparing these MORB glass analyses with the Ni-MgO systematics of MORB from other studies in the literature, which were obtained using a variety of analytical techniques, shows that these electron microprobe analyses are anomalous. We suggest that the reported variation of with S content in MORB is an analytical artifact.Mass balance of melt and olivine compositions with the starting compositions shows that dissolved S²⁻ depresses the olivine liquidus of haplobasaltic silicate melts by 5.8 × 10⁻³ (±1.3 × 10⁻³) K per μg g⁻¹ of S²⁻, which is negligible in most contexts. We also present data for the partitioning of some incompatible trace elements (Sc, Ti, Y, Zr and Hf) between olivine and melt. The data for Sc and Y confirm previous results showing that and decrease with increasing SiO₂ content of the melt. Values of average 0.01 with most falling in the range 0.005-0.015. Zr and Hf are considerably more incompatible than Ti in olivine, with and about 10⁻³. The ratio / is well constrained at 0.611 ± 0.016.     

Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Influence of Al for Si substitution in the tetrahedral network of silicate melts

The influence on olivine/melt transition metal (Mn, Co, Ni) partitioning of substitution in the tetrahedral network of silicate melt structure has been examined at ambient pressure in the 1450-1550 °C temperature range. Experiments were conducted in the systems NaAlSiO₄-Mg₂SiO₄- SiO₂ and CaAl₂Si₂O₈-Mg₂SiO₄-SiO₂ with about 1 wt% each of MnO, CoO, and NiO added. These compositions were used to evaluate how, in silicate melts, substitution and ionization potential of charge-balancing cations affect activity-composition relations in silicate melts and mineral/melt partitioning.The exchange equilibrium coefficient, , is a positive and linear function of melt Al/(Al + Si) at constant degree of melt polymerization, NBO/T. The is negatively correlated with the ionic radius, r, of the M-cation and also with the ionization potential (Z/r², Z = electrical charge) of the cation that serves to charge-balance Al³⁺ in tetrahedral coordination in the melts. The activity coefficient ratio, (γ_M/γ_Mg)^melt, is therefore similarly correlated.These melt composition relationships are governed by the distribution of Al³⁺ among coexisting Q-species in the peralkaline (depolymerized) melts coexisting with olivine. This distribution controls Q-speciation abundance, which, in turn, controls (γ_M/γ_Mg)^melt and . The relations between melt structure and olivine/melt partitioning behavior lead to the suggestion that in natural magmatic systems mineral/melt partition coefficients are more dependent on melt composition and, therefore, melt structure the more alkali-rich and the more felsic the melt. Moreover, mineral/melt partition coefficients are more sensitive to melt composition the more highly charged or the smaller the ionic radius of the cation of interest.     

Ammonium loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

Comparative carbonate and hydroxide complexation of cations in seawater

Comparative concentrations of carbonate and hydroxide complexes in natural solutions can be expressed in terms of reactions with bicarbonate that have no explicit pH dependence (). Stability constants for this reaction with n = 1 were determined using conventional formation constant data expressed in terms of hydroxide and carbonate. Available data indicate that stability constants appropriate to seawater at 25 °C expressed in the form are on the order of 10^4.2 for a wide range of cations (M^z+) with z = +1, +2 and +3. Φ₁ is sufficiently large that species appear to substantially dominate MOH^z−1 species in seawater. Evaluations of comparative stepwise carbonate and hydroxide stability constant behavior leading to the formation of n = 2 and n = 3 complexes suggest that carbonate complexes generally dominate hydroxide complexes in seawater, even for cations whose inorganic speciation schemes in seawater are currently presumed to be strongly dominated by hydrolyzed forms (). Calculated stability constants, and , indicate that the importance of carbonate complexation is sufficiently large that carbonate and hydroxide complexes would be generally comparable even if calculated Φ₂ and Φ₃ values are overestimated by two or more orders of magnitude. Inclusion of mixed ligand species in carbonate-hydroxide speciation models allows cation complexation intensities (M_T/[M^z+]) to be expressed in the following form:

     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

Coprecipitation of chromate with calcite: Batch experiments and X-ray absorption spectroscopy

Batch experiments, combined with in situ spectroscopic methods, are used to examine the coprecipitation of Cr(VI) with calcite, including partitioning behavior, site-specific distribution of Cr on the surface of calcite single crystals, and local coordination of Cr(VI) in the calcite structure. It is found that the concentration of Cr incorporated in calcite increases with increasing Cr concentration in solution. The calculated apparent partition coefficient, , is highest at low Cr solution concentration, and decreases to a constant value with increasing Cr solution concentration. DIC images of the surface of calcite single crystals grown in the presence of exhibit well-defined growth hillocks composed of two pairs of symmetrically nonequivalent vicinal faces, denoted as + and −, which reflect the orientation of structurally nonequivalent growth steps. Micro-XRF mapping of the Cr distribution over a growth hillock shows preferential incorporation of Cr into the—steps, which is considered to result from differences in surface structure geometry. XANES spectra confirm that incorporated Cr is hexavalent, and no reduction of Cr(VI) in the X-ray beam was observed up to 2 days at room temperature. EXAFS fit results show the incorporated Cr(VI) has the expected first shell of 4 O at ∼1.64 ± 0.01 Å, consistent with . Best fit results show that the second shell is split with ∼2.5 Ca at ∼3.33 ± 0.05 and ∼2.2 Ca at ∼3.55 ± 0.05 Å, which confirms the incorporation of chromate into calcite. Consideration of possible local coordination indicates that significant distortion or disruption is required to accommodate in the calcite structure.