Many complex networks erase parts of their geometry as they develop, so that their evolution is difficult to quantify and trace. Here we introduce entropy measures for quantifying the complexity of street orientations and length variations within planar networks and apply them to the street networks of 41 British cities, whose geometric evolution over centuries can be explored. The results show that the street networks of the old central parts of the cities have lower orientation/length entropies - the streets are more tightly ordered and form denser networks - than the outer and more recent parts. Entropy and street length increase, because of spreading, with distance from the network centre. Tracing the 400-year evolution of one network indicates growth through densification (streets are added within the existing network) and expansion (streets are added at the margin of the network) and a gradual increase in entropy over time.

Fractures generated by internal fluid pressure, for example, dykes, mineral veins, many joints and man-made hydraulic fractures, are referred to as hydrofractures. Together with shear fractures, they contribute significantly to the permeability of fluid reservoirs such as those of petroleum, geothermal water, and groundwater. Analytical and numerical models show that – in homogeneous host rocks – any significant overpressure in hydrofractures theoretically generates very high crack tip tensile stresses. Consequently, overpressured hydrofractures should propagate and help to form interconnected fracture systems that would then contribute to the permeability of fluid reservoirs. Field observations, however, show that in heterogeneous and anisotropic, e.g., layered, rocks many hydrofractures become arrested or offset at layer contacts and do not form vertically interconnected networks. The most important factors that contribute to hydrofracture arrest are discontinuities (including contacts), stiffness changes between layers, and stress barriers, where the local stress field is unfavourable to hydrofracture propagation. A necessary condition for a hydrofracture to propagate to the surface is that the stress field along its potential path is everywhere favourable to extension-fracture formation so that the probability of hydrofracture arrest is minimised. Mechanical layering and the resulting heterogeneous stress field largely control whether evolving hydrofractures become confined to single layers (strata¬bound frac¬tures) or not (non-stratabound fractures) and, there¬fore, if a vertically intercon¬nec¬ted fracture system forms. Non-stratabound hydrofractures may propagate through many layers and generate interconnected fracture systems. Such systems commonly reach the percolation threshold and largely control the overall permeability of the fluid reservoirs within which they develop.

How much magma needs to be added to a shallow magma chamber to cause rupture, dyke injection and a potential eruption? Models that yield reliable answers to this question are needed in order to facilitate eruption forecasting. Development of a long-lived shallow magma chamber requires periodic influx of magmas from a parental body at depth. This redistribution process does not necessarily cause an eruption but produces a net volume change that can be measured geodetically by inversion techniques. Using continuum-mechanics and fracture-mechanics principles, we calculate the amount of magma contained at shallow depth beneath Santorini volcano, Greece. We demonstrate through structural analysis of dykes exposed within the Santorini caldera, previously published data on the volume of recent eruptions and geodetic measurements of the 2011–2012 unrest period, that the measured 0.02% increase in volume of Santorini’s shallow magma chamber was associated with magmatic excess pressure increase of around 1.1 MPa. This excess pressure was high enough to bring the chamber roof close to rupture and dyke injection. For volcanoes with known typical extrusion and intrusion (dyke) volumes, the new methodology presented here makes it possible to forecast the conditions for magma-chamber failure and dyke injection at any geodetically well-monitored volcano.

Natural hazards associated with volcanic edifices depend partly on how fracture resistant the edifices are, i.e. on their strengths. Observations worldwide indicate that large fluid-driven extension fractures (dikes, inclined sheets), shear fractures (landslides), and mixed-mode fractures (ring dikes and ring faults) normally propagate more easily in a basaltic edifice (shield volcano) than in a stratovolcano. For example, dike-fed eruptions occur once every few years in many basaltic edifices but once every 10^2-3 yr in many stratovolcanoes. Large landslides and caldera collapses also appear to be more common in a typical basaltic edifice/shield volcano than in a typical stratovolcano. In contrast to a basaltic edifice, a stratovolcano is composed of mechanically dissimilar rock layers, i.e. layers with mismatching elastic properties (primarily Young's modulus). Elastic mismatch encourages fracture deflection and arrest at contacts and increases the amount of energy needed for a large-scale edifice failure. Fracture-related hazards depend on the potential energy available to propagate the fractures which, in turn, depends on the boundary conditions during fracture propagation. Here there are two possible scenarios: one in which the outer boundary of the volcanic edifice or rift zone does not move during the fracture propagation (constant displacement); the other in which the boundary moves (constant load). In the former, the total potential energy is the strain energy stored in the volcano before fracture formation; in the latter, the total potential energy is the strain energy plus the work done by the forces moving the boundary. Constant-displacement boundary conditions favor small eruptions, landslides, and caldera collapses, whereas constant-load conditions favor comparatively large eruptions, landslides, and collapses. For a typical magma chamber (sill-like with a diameter of 8 km), the strain energy change due to magma-chamber inflation is estimated at the order of 10¹⁴ J (0.1 PJ). For comparison, the surface energy needed to form a typical feeder dike is of the same order of magnitude, or 10¹⁴ J. There are several processes besides magma-chamber inflation that may increase the strain energy in a volcano before eruption. Thus, during a typical unrest period with magma-chamber inflation, the added strain energy in the volcano is large enough for a typical feeder dike to form. An injected dike, however, only reaches the surface and becomes a feeder if it is able to propagate through the numerous layers and contacts that tend to deflect or arrest dikes. The strong elastic mismatch between layers that constitute stratovolcanoes not only encourages fracture arrest, but also the storage of more strain energy (than in a typical basaltic edifice/shield volcano) before fracture formation and failure. It is thus through producing materials of widely different mechanical properties that stratovolcanoes become strong and resilient.

Fractures that form when fluid pressure ruptures the rock are referred to as fluid-driven fractures or hydrofractures. These include most dykes, inclined sheets and sills, but also many mineral veins and joints, as well as human-made hydraulic fractures. While considerable field and theoretical work has focused on the geometry and arrest of hydrofractures, how they select their propagation paths, particularly in layered and faulted rocks, has received less attention. Here I propose that of all the possible paths that a given hydrofracture may follow, it selects the path of least (minimum) action as determined by Hamilton’s principle. This means that the selected path is that along which the energy transformed (released) multiplied by the time taken for the propagation is a minimum. Hydrofractures advance their tips/fronts in steps, with a time lag between the fracture front and the fluid front. In the present framework, each step is then controlled by Hamilton’s principle. The results suggest that when the hosting rock body is regarded as homogeneous, isotropic and non-fractured, hydrofracture paths are everywhere perpendicular to the trajectories of the minimum compressive (maximum tensile) principal stress σ₃ and follow the trajectories of the maximum principal compressive stress σ₁. When applied to layered and faulted rock body, the results indicate that hydrofracture paths may follow existing faults for a while, depending primarily on (1) the dip of the fault (steep faults are the most likely to be used by vertically propagating hydrofractures), and (2) the tensile strength across the fault compared with the tensile strength of the host rock along a path following the direction of σ₁. The results suggest that hydrofractures may use faults as parts of their paths primarily if the fault is steeply dipping and with close to zero tensile strength.