S U M M A R YPure shear wave data are only very rarely acquired for offshore site investigations and exploration. Here, we present details of a novel, seabed-coupled, shear wave vibrator and field data recorded by a densely populated, multicomponent ocean-bottom cable, to improve shallow soil characterization.The prototype shear wave vibrator uses vibroseis technology adopted for marine environments through its instalment on top of a suction anchor, assuring seabed coupling in combination with self-weight penetration. The prototype is depth rated to 1500 m water depth, and can be rotated while installed in the seabed. The philosophy is to acquire fully complementary seismic data to conventional P-and P-to-S-converted waves, in particular for 2-D profiling, VSP (vertical seismic profiling) or monitoring purposes, thereby exploiting advantages of shear waves over compressional waves for determining, for example, anisotropy, small-strain shear modulus and excess pore pressures/effective stress. The source was primarily designed for reservoir depths. However, significant energy is emitted as surface waves, which provide detailed geotechnical information through mapping of shear wave velocities in potentially high resolution of the upper soil units. To fully utilize pure shear wave content, a proper analysis of surface waves is paramount, due to the proximity of surface wave propagation speed with shear wave velocities.The experiment was carried out in the northern North Sea in 364 m water depth. Cable dragging was necessary to obtain close receiver spacing (2.5 m effective spacing), with total line length of 600 m. Frequency-waveform transforms reveal both Scholte and Love waves. Up to six surface wave modes are identified, that is, fundamental mode and several higher surface wave modes. The occurrence of these two different dispersive surface wave types with well-resolved higher modes allows for a unique analysis and inversion scheme for highresolution mapping of physical properties in the shallow subsurface as well as anisotropy, which is discussed in an accompanying paper. The data presented in this paper are thus a unique (long and densely populated receiver array allows for multimodal Love and Scholte surface waves from the marine environment) but challenging (marine operations) marine data set.

The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. Response to Reviewers: My comments to the reviewers are in the attached word-document, "comments reviewers final.docx". Powered by Editorial Manager® and Preprint Manager® from Aries Systems CorporationComputational Geoscience manuscript No. The former ignores capillary pressure and can be solved analytically. The latter has traditionally been formulated as an iterative integral solution, which suffers from convergence problems as the injection saturation approaches unity. Here an alternative approach is presented that solves the MSE using a pseudospectral Chebyshev differentiation matrix. The resulting pseudospectral solution is compared to results obtained from the original integral implementation and the Buckley-Leverett limit, when the capillary pressure becomes negligible. A self-contained MATLAB code to implement the new solution is provided within the manuscript.The new approach offers a robust and accurate method for verification of numerical codes solving two-phase flow with capillary pressure.

Peer-review status of attached le: eerEreviewed Citation for published item: tyeoD eyodeji nd wthisD imon eF nd xielsenD tefn nd ilrrsD itor nd fjørnr¦ D ore sF @PHIWA 9glosedEform eqution for susidene due to)uid prodution from ylindril on(ned quiferF9D tournl of hydrologyFD SUQ F ppF WTREWTWF Further information on publisher's website: httpsXGGdoiForgGIHFIHITGjFjhydrolFPHIWFHQFHRI

Abstract. Fracture attribute scaling and connectivity datasets from analogue systems are widely used to inform sub-surface fractured reservoir models in a range of geological settings. However, significant uncertainties are associated with the determination of reliable scaling parameters in surface outcrops. This has limited our ability to upscale key parameters that control fluid flow at reservoir to basin scales. In this study, we present nine 1D-transect (scanline) fault and fracture attribute datasets from Middle Devonian sandstones in Caithness (Scotland) that are used as an onshore analogue for nearby sub-surface reservoirs such as the Clair field, west of Shetland. By taking account of truncation and censoring effects in individual datasets, our multiscale analysis shows a preference for power-law scaling of fracture length over 8 orders of magnitude (10−4 to 104 m) and kinematic aperture over 4 orders of magnitude (10−6 to 10−2 m). Our assessment of the spatial organization (clustering and topology) provides a new basis for up-scaling fracture attributes collected in outcrop- to regional-scale analogues. We show how these relationships may inform knowledge of geologically equivalent sub-surface fractured reservoirs.