To send content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about sending content to .
To send content items to your Kindle, first ensure firstname.lastname@example.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about sending to your Kindle.
Note you can select to send to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Shallow granular avalanches on slopes close to repose exhibit hysteretic behaviour. For instance, when a steady-uniform granular flow is brought to rest it leaves a deposit of thickness
on a rough slope inclined at an angle
to the horizontal. However, this layer will not spontaneously start to flow again until it is inclined to a higher angle
, or the thickness is increased to
. This simple phenomenology leads to a rich variety of flows with co-existing regions of solid-like and fluid-like granular behaviour that evolve in space and time. In particular, frictional hysteresis is directly responsible for the spontaneous formation of self-channelized flows with static levees, retrogressive failures as well as erosion–deposition waves that travel through the material. This paper is motivated by the experimental observation that a travelling-wave develops, when a steady uniform flow of carborundum particles on a bed of larger glass beads, runs out to leave a deposit that is approximately equal to
. Numerical simulations using the friction law originally proposed by Edwards et al. (J. Fluid Mech., vol. 823, 2017, pp. 278–315) and modified here, demonstrate that there are in fact two travelling waves. One that marks the trailing edge of the steady-uniform flow and another that rapidly deposits the particles, directly connecting the point of minimum dynamic friction (at thickness
) with the deposited layer. The first wave moves slightly faster than the second wave, and so there is a slowly expanding region between them in which the flow thins and the particles slow down. An exact inviscid solution for the second travelling wave is derived and it is shown that for a steady-uniform flow of thickness
it produces a deposit close to
for all inclination angles. Numerical simulations show that the two-wave structure deposits layers that are approximately equal to
for all initial thicknesses. This insensitivity to the initial conditions implies that
is a universal quantity, at least for carborundum particles on a bed of larger glass beads. Numerical simulations are therefore able to capture the complete experimental staircase procedure, which is commonly used to determine the
curves by progressively increasing the inclination of the chute. In general, however, the deposit thickness may depend on the depth of the flowing layer that generated it, so the most robust way to determine
is to measure the deposit thickness from a flow that was moving at the minimum steady-uniform velocity. Finally, some of the pathologies in earlier non-monotonic friction laws are discussed and it is explicitly shown that with these models either steadily travelling deposition waves do not form or they do not leave the correct deposit depth
When a layer of static grains on a sufficiently steep slope is disturbed, an upslope-propagating erosion wave, or retrogressive failure, may form that separates the initially static material from a downslope region of flowing grains. This paper shows that a relatively simple depth-averaged avalanche model with frictional hysteresis is sufficient to capture a planar retrogressive failure that is independent of the cross-slope coordinate. The hysteresis is modelled with a non-monotonic effective basal friction law that has static, intermediate (velocity decreasing) and dynamic (velocity increasing) regimes. Both experiments and time-dependent numerical simulations show that steadily travelling retrogressive waves rapidly form in this system and a travelling wave ansatz is therefore used to derive a one-dimensional depth-averaged exact solution. The speed of the wave is determined by a critical point in the ordinary differential equation for the thickness. The critical point lies in the intermediate frictional regime, at the point where the friction exactly balances the downslope component of gravity. The retrogressive wave is therefore a sensitive test of the functional form of the friction law in this regime, where steady uniform flows are unstable and so cannot be used to determine the friction law directly. Upper and lower bounds for the existence of retrogressive waves in terms of the initial layer depth and the slope inclination are found and shown to be in good agreement with the experimentally determined phase diagram. For the friction law proposed by Edwards et al. (J. Fluid. Mech., vol. 823, 2017, pp. 278–315, J. Fluid. Mech., 2019, (submitted)) the magnitude of the wave speed is slightly under-predicted, but, for a given initial layer thickness, the exact solution accurately predicts an increase in the wave speed with higher inclinations. The model also captures the finite wave speed at the onset of retrogressive failure observed in experiments.
Filamentary structures can form within the beam of protons accelerated during the interaction of an intense laser pulse with an ultrathin foil target. Such behaviour is shown to be dependent upon the formation time of quasi-static magnetic field structures throughout the target volume and the extent of the rear surface proton expansion over the same period. This is observed via both numerical and experimental investigations. By controlling the intensity profile of the laser drive, via the use of two temporally separated pulses, both the initial rear surface proton expansion and magnetic field formation time can be varied, resulting in modification to the degree of filamentary structure present within the laser-driven proton beam.
Rapid shallow granular flows over inclined planes are often seen in nature in the form of avalanches, landslides and pyroclastic flows. In these situations the flow develops an inversely graded (large at the top) particle-size distribution perpendicular to the plane. As the surface velocity of such flows is larger than the mean velocity, the larger material is transported to the flow front. This causes size segregation in the downstream direction, resulting in a flow front composed of large particles. Since the large particles are often more frictional than the small, the mobility of the flow front is reduced, resulting in a so-called bulbous head. This study focuses on the formation and evolution of this bulbous head, which we show to emerge in both a depth-averaged continuum framework and discrete particle simulations. Furthermore, our numerical solutions of the continuum model converge to a travelling wave solution, which allows for a very efficient computation of the long-time behaviour of the flow. We use small-scale periodic discrete particle simulations to calibrate (close) our continuum framework, and validate the simple one-dimensional (1-D) model with full-scale 3-D discrete particle simulations. The comparison shows that there are conditions under which the model works surprisingly well given the strong approximations made; for example, instantaneous vertical segregation.
Small perturbations to a steady uniform granular chute flow can grow as the material moves downslope and develop into a series of surface waves that travel faster than the bulk flow. This roll wave instability has important implications for the mitigation of hazards due to geophysical mass flows, such as snow avalanches, debris flows and landslides, because the resulting waves tend to merge and become much deeper and more destructive than the uniform flow from which they form. Natural flows are usually highly polydisperse and their dynamics is significantly complicated by the particle size segregation that occurs within them. This study investigates the kinematics of such flows theoretically and through small-scale experiments that use a mixture of large and small glass spheres. It is shown that large particles, which segregate to the surface of the flow, are always concentrated near the crests of roll waves. There are different mechanisms for this depending on the relative speed of the waves, compared to the speed of particles at the free surface, as well as on the particle concentration. If all particles at the surface travel more slowly than the waves, the large particles become concentrated as the shock-like wavefronts pass them. This is due to a concertina-like effect in the frame of the moving wave, in which large particles move slowly backwards through the crest, but travel quickly in the troughs between the crests. If, instead, some particles on the surface travel more quickly than the wave and some move slower, then, at low concentrations, large particles can move towards the wave crest from both the forward and rearward sides. This results in isolated regions of large particles that are trapped at the crest of each wave, separated by regions where the flow is thinner and free of large particles. There is also a third regime arising when all surface particles travel faster than the waves, which has large particles present everywhere but with a sharp increase in their concentration towards the wave fronts. In all cases, the significantly enhanced large particle concentration at wave crests means that such flows in nature can be especially destructive and thus particularly hazardous.
Giant electromagnetic pulses (EMP) generated during the interaction of high-power lasers with solid targets can seriously degrade electrical measurements and equipment. EMP emission is caused by the acceleration of hot electrons inside the target, which produce radiation across a wide band from DC to terahertz frequencies. Improved understanding and control of EMP is vital as we enter a new era of high repetition rate, high intensity lasers (e.g. the Extreme Light Infrastructure). We present recent data from the VULCAN laser facility that demonstrates how EMP can be readily and effectively reduced. Characterization of the EMP was achieved using B-dot and D-dot probes that took measurements for a range of different target and laser parameters. We demonstrate that target stalk geometry, material composition, geodesic path length and foil surface area can all play a significant role in the reduction of EMP. A combination of electromagnetic wave and 3D particle-in-cell simulations is used to inform our conclusions about the effects of stalk geometry on EMP, providing an opportunity for comparison with existing charge separation models.
In this study, we report the characterization of a 304L stainless steel cylindrical projectile produced by additive manufacturing. The projectile was compressively deformed using a Taylor Anvil Gas Gun, leading to a huge strain gradient along the axis of the deformed cylinder. Spatially resolved neutron diffraction measurements on the HIgh Pressure Preferred Orientation time-of-flight diffractometer (HIPPO) and Spectrometer for Materials Research at Temperature and Stress diffractometer (SMARTS) beamlines at the Los Alamos Neutron Science CEnter (LANSCE) with Rietveld and single-peak analysis were used to quantitatively evaluate the volume fractions of the α, γ, and ε phases as well as residual strain and texture. The texture of the γ phase is consistent with uniaxial compression, while the α texture can be explained by the Kurdjumov–Sachs relationship from the γ texture after deformation. This indicates that the material first deformed in the γ phase and subsequently transformed at larger strains. The ε phase was only found in volumes close to the undeformed material with a texture connected to the γ texture by the Shoji–Nishiyama orientation relationship. This allows us to conclude that the ε phase occurs as an intermediate phase at lower strain, and is superseded by the α phase when strain increases further. We found a proportionality between the root-mean-squared microstrain of the γ phase, dominated by the dislocation density, with the α volume fraction, consistent with strain-induced martensite α formation. Knowledge of the sample volume with the ε phase from the neutron diffraction analysis allowed us to identify the ε phase by electron back scatter diffraction analysis, complementing the neutron diffraction analysis with characterization on the grain level.
The Binary Population and Spectral Synthesis suite of binary stellar evolution models and synthetic stellar populations provides a framework for the physically motivated analysis of both the integrated light from distant stellar populations and the detailed properties of those nearby. We present a new version 2.1 data release of these models, detailing the methodology by which Binary Population and Spectral Synthesis incorporates binary mass transfer and its effect on stellar evolution pathways, as well as the construction of simple stellar populations. We demonstrate key tests of the latest Binary Population and Spectral Synthesis model suite demonstrating its ability to reproduce the colours and derived properties of resolved stellar populations, including well-constrained eclipsing binaries. We consider observational constraints on the ratio of massive star types and the distribution of stellar remnant masses. We describe the identification of supernova progenitors in our models, and demonstrate a good agreement to the properties of observed progenitors. We also test our models against photometric and spectroscopic observations of unresolved stellar populations, both in the local and distant Universe, finding that binary models provide a self-consistent explanation for observed galaxy properties across a broad redshift range. Finally, we carefully describe the limitations of our models, and areas where we expect to see significant improvement in future versions.
A comparative study is presented of the chemistry and crystallography of zinc-bearing strunzites from Hagendorf Süd, Bavaria, Germany and the Sitio do Castelo mine, Folgosinho, Portugal. Electron microprobe analyses of samples from the two localities show quite different cation substitutions. The Hagendorf Süd mineral is a Zn-bearing ferristrunzite, with compositional zoning due to Zn2+ replacing predominantly Fe3+ as well as minor Mn2+, whereas the Portugese mineral is a Zn-bearing strunzite, in which Zn2+ replaces Mn2+, with minor replacement of Fe3+ by Mn3+. Zincostrunzite, with dominant Zn in the interlayer octahedrally coordinated site, is a new strunzite-group mineral that has been characterized at both locations. Analysis of single-crystal synchrotron data for zinc-bearing ferristrunzite and zincostrunzite crystals from Hagendorf Süd show that the structures of both minerals contain zeolitic water in the interlayer region. The formula for strunzite-group minerals containing the zeolitic water is MFe23+(PO4)2(OH)2·6.5H2O, M=Fe, Mn, Zn. This formulation agrees with that found for zincostrunzite from the Sitio do Castelo mine, but differs from that reported previously for strunzite, MFe2+(PO4)2(OH)2·6H2O, which has no interlayer water. Interestingly, the zincostrunzites from the two localities differ in the location of the interlayer water molecule, with a corresponding difference in the H bonding.
The Savage-Hutter equations of granular avalanche flows are a hyperbolic system of equations for the distribution of depth and depth-averaged velocity components tangential to the sliding bed. They involve two phenomenological parameters, the internal and the bed friction angles, which together define the earth pressure coefficient which assumes different values depending upon whether the flow is either diverging or contracting. Because of the hyperbolicity of the equations, since velocities may be supercritical, shock waves are often formed in avalanche flows. Numerical schemes solving these free surface flows must cope with smooth as well as non-smooth solutions. In this paper the Savage-Hutter equations in conservative form are solved with a shock-capturing technique, including a front-tracking method. This method can perform for parabolic similarity solutions for which the Lagrangian scheme is excellent, and it is even better in other situations when the latter fails.
One means of preventing areas from being hit by avalanches is to divert the flow by straight or curved walls or tetrahedral or cylindrical-type structures. Thus, there arises the question how a given avalanche flow is changed regarding the diverted-flow depth and flow direction. In this paper a report is given on laboratory experiments performed for gravity-driven dense granular flows down an inclined plane obstructed by plane wall and tetrahedral wedge. It was observed that these flows are accompanied by shocks induced by the presence of the obstacles. These give rise to a transition from super-to subcritical flow of the granular avalanche, associated with depth and velocity changes. It is demonstrated that with an appropriate shock-capturing integration technique for the Savage-Hutter theory, the shock formation for a finite-mass granular flow sliding from an inclined plane into a horizontal run-out zone is well described, as is the shock formation of the granular flow on either side of a tetrahedral protection structure.
Introduction: Effective trauma resuscitation requires a coordinated team approach, yet there is a significant risk for error. These errors can manifest from sequential system-, team- and knowledge based failures, defined as latent safety threats (LSTs). In situ simulation (ISS), a point-of-care training strategy, provides a novel prospective approach to identify factors that impact patient safety. This study quantified and formulated a hierarchy of LSTs during risk-informed ISS trauma resuscitations. Methods: At a Level 1 trauma centre, we conducted 12 multi-disciplinary, unannounced ISSs to prospectively identify trauma-related LSTs. Four, risk-informed scenarios were developed based on 5 recurring themes found within the trauma program’s morbidity and mortality process. The actual, on-call trauma team participated in the study. Simulations were video recorded with 4 cameras, each positioned at a different angle. Using a framework analysis methodology, human factors experts transcribed and coded the videos. Thematic structure was established deductively based on existing literature and inductively based on observed ISS events. All LSTs were prioritized for future patient safety, systems and ergonomic interventions using the Healthcare Failure Mode and Effect Analysis (HFMEA) matrix. Results: We identified 893 LSTs from 12 simulations. LST analysis resulted in 8 themes subcategorized into 43 codes. Themes were associated with team-, knowledge- or system-related issues. The following themes emerged: situational awareness, provider safety, mental model alignment, team/individual responsibility, team resources, equipment considerations, workplace environment and clinical protocols. The HFMEA hazard scoring process identified 13 high priority codes that required urgent attention and intervention to mitigate negative patient outcomes. Conclusion: A prospective, video-based framework analysis represents a novel and robust approach to LST identification within trauma care. Patterns of LSTs within and between simulations provide a high degree of transparency and traceability for an inter-professional trauma program review. Hazard matrix scoring facilitates the classification and prioritization of human factors interventions intended to improve patient safety.
Geophysical granular flows, such as avalanches, debris flows, lahars and pyroclastic flows, are always strongly influenced by the basal topography that they flow over. In particular, localised bumps or obstacles can generate rapid changes in the flow thickness and velocity, or shock waves, which dissipate significant amounts of energy. Understanding how a granular material is affected by the underlying topography is therefore crucial for hazard mitigation purposes, for example to improve the design of deflecting or catching dams for snow avalanches. Moreover, the interactions with solid boundaries can also have important applications in industrial processes. In this paper, small-scale experiments are performed to investigate the flow of a granular avalanche over a two-dimensional smooth symmetrical bump. The experiments show that, depending on the initial conditions, two different steady-state regimes can be observed: either the formation of a detached jet downstream of the bump, or a shock upstream of it. The transition between the two cases can be controlled by adding varying amounts of erodible particles in front of the obstacle. A depth-averaged terrain-following avalanche theory that is formulated in curvilinear coordinates is used to model the system. The results show good agreement with the experiments for both regimes. For the case of a shock, time-dependent numerical simulations of the full system show the evolution to the equilibrium state, as well as the deposition of particles upstream of the bump when the inflow ceases. The terrain-following theory is compared to a standard depth-averaged avalanche model in an aligned Cartesian coordinate system. For this very sensitive problem, it is shown that the steady-shock regime is captured significantly better by the terrain-following avalanche model, and that the standard theory is unable to predict the take-off point of the jet. To retain the practical simplicity of using Cartesian coordinates, but have the improved predictive power of the terrain-following model, a coordinate mapping is used to transform the terrain-following equations from curvilinear to Cartesian coordinates. The terrain-following model, in Cartesian coordinates, makes identical predictions to the original curvilinear formulation, but is much simpler to implement.
Observations of slab-avalanche releases in alpine terrain have led to the hypothesis that rocky outcrops can influence the spatial distributions of temperature and heal flow in dry alpine snow covers and Unis control the local distribution of depth hoar. We investigate the effects of terrain on crystal growth by using a two-dimensional finite-clemcnl model of heal How coupled with a model of crystal growth from vapor. We used the model to examine the influence of snow properties, terrain geometry and snow depth on this phenomenon. The efleel is stronger in the early winter than in the late winter, because the rock has then had time to cool. In all cases, we found that depth-hoar growth occurs preferentially over the rock. This suggests that snow-pit investigations made over soil can he misleading if rocky outcrops are present.
Kummerite, ideally Mn2+Fe3+A1(PO4)2(OH)2.8H2O, is a new secondary phosphate mineral belonging to the laueite group, from the Hagendorf-Süd pegmatite, Hagendorf, Oberpfalz, Bavaria, Germany. Kummerite occurs as sprays or rounded aggregates of very thin, typically deformed, amber yellow laths. Cleavage is good parallel to ﹛010﹜. The mineral is associated closely with green Zn- and Al-bearing beraunite needles. Other associated minerals are jahnsite-(CaMnMn) and Al-bearing frondelite. The calculated density of kummerite is 2.34 g cm 3. It is optically biaxial (-), α= 1.565(5), β = 1.600(5) and y = 1.630(5), with weak dispersion. Pleochroism is weak, with amber yellow tones. Electron microprobe analyses (average of 13 grains) with H2O and FeO/Fe2O3 calculated on structural grounds and normalized to 100%, gave Fe2O3 17.2, FeO 4.8, MnO 5.4, MgO 2.2, ZnO 0.5, Al2O3 9.8, P2O5 27.6, H2O 32.5, total 100 wt.%. The empirical formula, based on 3 metal apfu is (Mn2+0.37Mg0.27Zn0.03Fe2+0.33)Σ1.00(Fe3+1.06Al0. 94)Σ2.00PO4)1.91(OH)2.27(H2O)7.73. Kummerite is triclinic, P1̄, with the unit-cell parameters of a = 5.316(1) Å, b =10.620(3) Å , c = 7.118(1) Å, α = 107.33(3)°, β= 111.22(3)°, γ = 72.22(2)° and V= 348.4(2) Å3. The strongest lines in the powder X-ray diffraction pattern are [dobs in Å(I) (hkl)] 9.885 (100) (010); 6.476 (20) (001); 4.942 (30) (020); 3.988 (9) (̄110); 3.116 (18) (1̄20); 2.873 (11) (1̄21). Kummerite is isostructural with laueite, but differs in having Al and Fe3+ ordered into alternate octahedral sites in the 7.1 Å trans-connected octahedral chains.
The collective response of electrons in an ultrathin foil target irradiated by an ultraintense (
) laser pulse is investigated experimentally and via 3D particle-in-cell simulations. It is shown that if the target is sufficiently thin that the laser induces significant radiation pressure, but not thin enough to become relativistically transparent to the laser light, the resulting relativistic electron beam is elliptical, with the major axis of the ellipse directed along the laser polarization axis. When the target thickness is decreased such that it becomes relativistically transparent early in the interaction with the laser pulse, diffraction of the transmitted laser light occurs through a so called ‘relativistic plasma aperture’, inducing structure in the spatial-intensity profile of the beam of energetic electrons. It is shown that the electron beam profile can be modified by variation of the target thickness and degree of ellipticity in the laser polarization.
We calibrated portions of the radiocarbon time scale with combined 230Th, 231Pa, 14C measurements of corals collected from Espiritu Santo, Vanuatu and the Huon Peninsula, Papua New Guinea. The new data map 14C variations ranging from the current limit of the tree-ring calibration [11,900 calendar years before present (cal BP), Kromer and Spurk 1998, now updated to 12,400 cal B P, see Kromer et al., this issue], to the 14C-dating limit of 50,000 cal BP, with detailed structure between 14 to 16 cal kyr BP and 19 to 24 cal kyr BP. Samples older than 25,000 cal BP were analyzed with high-precision 231Pa dating methods (Pickett et al. 1994; Edwards et al. 1997) as a rigorous second check on the accuracy of the 230Th ages. These are the first coral calibration data to receive this additional check, adding confidence to the age data forming the older portion of the calibration. Our results, in general, show that the offset between calibrated and 14C ages generally increases with age until about 28,000 cal BP, when the recorded 14C age is nearly 6800 yr too young. The gap between ages before this time is less; at 50,000 cal BP, the recorded 14C age is 4600 yr too young. Two major 14C-age plateaus result from a 130 drop in Δ14C between 14–15 cal kyr BP and a 700 drop in Δ14C between 22–25 cal kyr BP. In addition, a large atmospheric Δ14C excursion to values over 1000 occurs at 28 cal kyr BP. Between 20 and 10 cal kyr BP, a component of atmospheric Δ14C anti-correlates with Greenland ice δ18O, indicating that some portion of the variability in atmospheric Δ14C is related to climate change, most likely through climate-related changes in the carbon cycle. Furthermore, the 28-kyr excursion occurs at about the time of significant climate shifts. Taken as a whole, our data indicate that in addition to a terrestrial magnetic field, factors related to climate change have affected the history of atmospheric 14C.
Introduction: Resuscitation of a trauma patient requires a multidisciplinary team to perform in a dynamic, high-stakes environment. Error is ubiquitous in trauma care, often related to latent safety threats (LSTs) - previously unrecognized threats that can materialize at any time. In-situ simulation (ISS) allows a team to practice in their authentic environment while providing an opportunistic milieu to explore critical events and uncover LSTs that impact patient safety. Methods: At a Canadian Level 1 trauma centre, regular, unannounced trauma ISSs were conducted and video-recorded. A retrospective chart review of adverse events or unexpected deaths informed ISS scenario design. Each session began with a trauma team activation. The on-duty trauma team arrived in the trauma bay and provided care as they would for a real patient. Semi-structured debriefing with participant-driven LST identification and ethnographic observation occurred in real time. A framework analysis using video review was conducted by human factors experts to identify and evaluate LSTs. Feasibility was measured by the impact on ED workflow, interruptions of clinical care and participant feedback. Results: Six multidisciplinary, high-fidelity, ISS sessions were conducted and 70 multidisciplinary staff and trainees participated in at least one session. Using a framework analysis, LSTs were identified and categorized into seven themes that relate to clinical tasks, equipment, team communication, and participant workflow. LSTs were quantified and prioritized using a hazard scoring matrix. ISS was effectively implemented during both low and high patient volume situations. No critical interruptions in patient care were identified during ISS sessions and overall participant feedback was positive. Conclusion: This novel, multidisciplinary ISS trauma training program integrated risk-informed simulation cases with human factors analysis to identify LSTs. ISS offers an opportunity for an iterative review process of high-risk situations beyond the traditional morbidity and mortality rounds; rather than waiting for an actual case to generate discussion and review, we prophylactically examined critical situations and processes. Findings form a framework for recommendations about improvements in equipment, environment layout, workflow, system processes, effective team training, and ultimately patient safety.
During 1990 we surveyed the southern sky using a multi-beam receiver at frequencies of 4850 and 843 MHz. The half-power beamwidths were 4 and 25 arcmin respectively. The finished surveys cover the declination range between +10 and −90 degrees declination, essentially complete in right ascension, an area of 7.30 steradians. Preliminary analysis of the 4850 MHz data indicates that we will achieve a five sigma flux density limit of about 30 mJy. We estimate that we will find between 80 000 and 90 000 new sources above this limit. This is a revised version of the paper presented at the Regional Meeting by the first four authors; the surveys now have been completed.
The subsurface exploration of other planetary bodies can be used to unravel their geological history and assess their habitability. On Mars in particular, present-day habitable conditions may be restricted to the subsurface. Using a deep subsurface mine, we carried out a program of extraterrestrial analog research – MINe Analog Research (MINAR). MINAR aims to carry out the scientific study of the deep subsurface and test instrumentation designed for planetary surface exploration by investigating deep subsurface geology, whilst establishing the potential this technology has to be transferred into the mining industry. An integrated multi-instrument suite was used to investigate samples of representative evaporite minerals from a subsurface Permian evaporite sequence, in particular to assess mineral and elemental variations which provide small-scale regions of enhanced habitability. The instruments used were the Panoramic Camera emulator, Close-Up Imager, Raman spectrometer, Small Planetary Linear Impulse Tool, Ultrasonic drill and handheld X-ray diffraction (XRD). We present science results from the analog research and show that these instruments can be used to investigate in situ the geological context and mineralogical variations of a deep subsurface environment, and thus habitability, from millimetre to metre scales. We also show that these instruments are complementary. For example, the identification of primary evaporite minerals such as NaCl and KCl, which are difficult to detect by portable Raman spectrometers, can be accomplished with XRD. By contrast, Raman is highly effective at locating and detecting mineral inclusions in primary evaporite minerals. MINAR demonstrates the effective use of a deep subsurface environment for planetary instrument development, understanding the habitability of extreme deep subsurface environments on Earth and other planetary bodies, and advancing the use of space technology in economic mining.