Vol49,4_5,2006 Key words gradiometric gravimeter – accelerome- ter – Whispering gallery cavity 1. Introduction The gravimetric methods are widely used in geological structure investigations, in particular in oil and mine exploration industry. This kind of research requires a very high level of instrumen- tal sensitivity and an accuracy of the order of one microGal (1 Gal = 1 cm/s2), that is 10−9 g, in or- der to achieve significant results. At this sensitiv- ity level spurious effects, not connected with un- derground mass distributions, but that are due to surface topographic irregularities or to anthropic activities or to tidal fluctuations become impor- tant. The instrumental drift, unavoidable in me- chanical sensors because of creep of the elastic elements, thermal gradients etc., is another source of errors, which are more and more im- portant as the data acquisition time is longer. Much more severe constraints rise in implement- ing a gravimeter on a ship or on a plane, because the instrument cannot distinguish between gravi- ty and kinematic acceleration vectors. Gravity gradiometric methods can avoid many of these problems. This technique is based on the measurement of one or more of the five independent components of the tensor A sapphire monolithic differential accelerometer as core sensor for gravity gradiometric geophysical instrumentation Alessandro Bertolini (1)(3)(∗), Nicolò Beverini (1) (2), Andrea De Michele (1), Francesco Fidecaro (1) and Francesco Mango (1) (1) Dipartimento di Fisica «E. Fermi», Università degli Studi di Pisa, and Istituto Nazionale di Fisica della Materia (INFM), CNR, Pisa, Italy (2) Consorzio Universitario della Spezia, La Spezia, Italy (3) Marwan Technology Srl Pisa, Italy Abstract Gradiometric gravimetry is a survey technique widely used in geological structure investigation. This work demon- strates the feasibility of a new class of low frequency accelerometers for geodynamics studies and space applica- tions. We present the design features of a new low noise single-axis differential accelerometer; the sensor is suit- able to be used in a Gravity Gradiometer (GG) system for land geophysical survey and gravity gradient measure- ments. A resolution of 1 Eötvös (1 Eö = 10−9s−2) at one sample per second is achievable in a compact, lightweight (less than 2 kg) portable instrument, operating at room temperature. The basic components of the sensor are two identical rigidly connected accelerometers separated by a 15-cm baseline vector and the useful signal is extracted as the subtraction of the two outputs, by means of an interferometric microwave readout system. The structure will be engraved in a monocrystal of sapphire by means of Computer-Numerically-Controlled (CNC) ultrasonic ma- chining: the material was chosen because of its unique mix of outstanding mechanical and dielectric properties. Mailing address: Prof. Nicolò Beverini, Dipartimento di Fisica «E. Fermi», Università degli Studi di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy; e-mail: beverini@df.unipi.it (∗) Present address: DESY, Notkestrasse 85, D 22607 Hamburg, Germany. 1095 ANNALS OF GEOPHYSICS, VOL. 49, N. 4/5, August/October 2006 1096 Alessandro Bertolini, Nicolò Beverini, Andrea De Michele, Francesco Fidecaro and Francesco Mango Γij = ∂ gj /∂ xi (that is the gradient of the gravity vector g), as the difference between the values of the gj components across a fixed distance ba- sis. This technique has demonstrated his signif- icance in the oil deposit signature and in the re- search and localization of underground cavities. The on-field implementation of the gradio- metric principle can be made by reading, with a high sensitivity gravimetric instrument, the verti- cal gravity values gz in two different point, sepa- rated by a known vectorial basis ri, and evaluat- ing the Γiz component of the tensor. However, the full exploitation of the peculiarity of the gradio- metric surveys could be achieved by using a sin- gle instrument with built-in differential readout capability, in order to efficiently reject common mode noise and spurious signal, like kinematic effects due to linear accelerations of the meas- urement platform. The basic practical implemen- tation consists of two identical accelerometers fixed rigidly together with a common sensitive axis (see fig. 1). Below the mechanical resonant frequency ω0 of the accelerometers, the difference in the vertical position of the proof masses is linearly related to the gravity gradient component Γiz ac- cording to The second term is related to static mechanical mismatch δω between the two accelerometers. Being almost impossible to build two perfectly matched accelerometers, a residual sensitivity to kinematic effects arises, polluting the signal channel by an amount proportional to the plat- form motion acceleration az. The factor ω0 / 2δω is defined as the Common-Mode-Rejection-Ra- tio (CMRR). Detecting 1 Eö gravity gradient in presence of 10−3 m/s2 linear acceleration, with an instrument with 15-cm baseline, requires a 130 dB of CMRR. Non-linearities of the accelerom- eter spring system also affect the CMRR for large environmental accelerations. This imposes suitable vibration isolation platforms, especially in airborne instruments. At present only a few gravity gradiometers have been developed. Spaceborne gravity gra- diometer devices are capable of outstanding noise levels (10−4 − 10−3 Eö) but rely on the very low environmental noise of a low acceleration satellite (at least one million times more silent than a typical light aircraft, for instance). Two ex- amples are the 6-axis superconducting differen- tial accelerometer developed by Ho Jung Paik group at University of Maryland (Chain and Paik, 1988a,b), and now being built by Oxford Instru- ments, and the GOCE Mission accelerometer (Rebhan et al., 2000). Laboratory use only appa- ratus have also been devised by doubling the core sensor of absolute gravimeters, like the supercon- ducting levitation device from GWR Inc. (Richter and Warburton, 1998), the free-fall Michelson interferometer from Micro-g Solutions Inc. (http://www.microgsolutions.com/gradiometer.pr ess.htm) and the atom interferometer developed by Kasevich group at Yale University (Snadden et al., 1998). All these instruments claim a reso- lution ranging from 0.1 to 1 Eö at one sample per second, but they are still far away from the porta- bility and the on-field applications. The most popular portable instrumentation is the Gravity Gradiometer Instrument (GGI), de- .z z r g a 2 iz i z z 1 2 0 2 0 0 2. ω ω δω ω Γ - + + b ^ l h Fig. 1. Example basic configurations of a single- component gravity gradiometer. The geometric ar- rangement of the two linear accelerometers is shown. A sapphire monolithic differential accelerometer as core sensor for gravity gradiometric geophysical instrumentation veloped primarily by Lockheed Martin for mili- tary application, but recently declassified and pre- sented now on the civil market (Goldstein and Brett, 1998). They are instruments specifically designed to be mounted on mobile platforms (air- crafts or ships). In the simplest form it consists of two matched pairs of force-balance accelerome- ters mounted orthogonally on a rotating platform (see fig. 2). Linear accelerations perpendicular to the rotation axis are modulated at the spinning frequency, while the gravity gradient effect is modulated at twice spinning frequency. In this way the useful signal is separated in the frequen- cy domain from the instrument bias and vehicle accelerations (e.g., reference ellipsoid, latitude, earth tide and isostatic effects). The signal is ex- tracted by a suitable combination of the differ- ence in readings between opposing pairs of ac- celerometers. A high CMRR is achieved by a feedback system acting on the accelerometers scale-factor, and by the effect of the rotating plat- form. The sensitivity declared is of the order of 5 Eö. The GGI is the basic element of the most complete instrument, the Bell Geospace Full Tensor Gradiometer (http://www.bellgeo.com/ tech/system.html). This instrument measures at the same time all the components of the gradient tensor, allowing a complete reconstruction of the gravity map. The performances of this instrument are very high; it is however a high cost cumber- some instrument, of large volume and weight that must be carried on a ship or a plane. More recent- ly the Bell technology has been shared for ad- vanced developments, the Falcon Gravity Gra- diometer project, with the Australian company BHP Billiton Discovery Technologies (Lee, 2001). The Falcon system has proven to be effec- tive in detecting sources with Γzz > 10 Eö rms and a full width at half maximum of 500 m of spatial resolution. At this time, the most promising portable instrument seems to be the mobile super- conducting beam-balance gravity gradiometer developed by F. Van Kann (Van Kann et al., 1994) at University of Western Australia and now being launched on the market by Gedex Inc. (for details see Matthews, 2002). The instrument is specifically designed for airborne surveys and is capable of unprecedented CMRR (exceeding 190 dB) and of a resolution better than 1 Eö at one sample per second, even in the presence of 1 m/s2 linear accelerations. However, even in this case, the size of the accessory inertial platform and of the liquid helium cooling system still limit the portability of the instrument. No available device currently fills the gap between laboratory-only apparatus and large size, costly GGs for airborne and shipborne ap- plications. In this work we present the design of a new core sensor for an advanced lightweight portable gravity gradiometer. The sensor is a differential vertical accel- erometer using state-of-art dielectric crystals machining techniques and microwave position sensing technologies to achieve 1 Eö level res- olutions operating even at room temperature. 2. The mechanical apparatus The instrument consists of three main parts: the accelerometer plate, hosting two identical ac- celerometers sensitive along the z-axis, the refer- 1097 Fig. 2. Basic scheme of a GGI (Gravity Gradiometer Instrument): two opposite matched pairs of linear ac- celerometers sit on a spinning platform. The gravity gradient signal is given by a linear combination of the accelerometer outputs at twice the spinning frequency. 1098 Alessandro Bertolini, Nicolò Beverini, Andrea De Michele, Francesco Fidecaro and Francesco Mango ence plate, hosting the motion sensing mi- crowave elements and a conductive shield. Both the accelerometer and the reference plates are completely carved out of a high quality single- crystal sapphire plate 1-cm thick by means of Computer Numerically Controlled (CNC) chisel ultrasonic machining (Bertolini et al., 2003). Ul- trasonic machining yields complex shapes in crystals, ceramics and composite materials with- out any release of thermal and mechanical stress into the work piece. Because of its geometry, in the present apparatus design, the instrument is sensitive to the off-line Γyz component of the gravity gradient tensor. Each accelerometer is a small torsion pendulum with a 0.04 kg proof mass suspended by two torsion beams. The base- line, defined as the relative position of the cen- ters of mass of the moving bodies (fig. 3), is 15- cm. A force Fz applied to the proof mass pro- duces a twist θ of the torsion pendulum given by where l, b, t are respectively the length, height and thickness of the torsion beams; G is the rigidity modulus of the material and d is the dis- bt G F ld3 z 3.ϑ tance between the center-of-mass and the rota- tion axis. By using sapphire beams 1×1 mm2 in section and 10 mm long, a resonant frequency of 20 Hz can be obtained. In this configuration, a 1 Eö gravity gradient is equivalent to a differential displacement of the GG test masses of 10−14m. Sapphire has been chosen because its excellent mechanical and dielectric properties. It com- bines high strength, low hysteresis and very low electrical losses; HEMEX™ grade sapphire sin- gle-crystal exhibit mechanical and electrical quality factors greater than 105 even at room temperature (http://www.crystalsystems.com). The monolithic structure obtained by CNC ma- chining allows the accuracy in the relative posi- tioning of the different mechanical parts to be kept within a few micrometers. The monolithic design also allows very low hysteresis, elimi- nating shear effects at the contact surfaces be- tween separate mechanical parts (stick-and- slip), while the use of a single material limits the instrument’s thermal sensitivity. Co-ma- chining of both the torsion springs effectively strongly reduces directional cross-talks due to misalignments between the components in the assembly, typical of a composite structure. In this way, a 60 dB mechanical common mode re- jection factor can be reasonably achieved just by construction. A fine trimming of the CMRR will be achieved by acting on the «flap» elec- trodes. By applying a high-voltage on both sides of a flap, a tunable negative stiffness ef- fect is realized. The trimming procedure will make use of a suitable tilting jig which eventu- ally should even be integrated in the GG system to allow for an automated balancing routine. Tilting is the easiest way to inject a 1-g level common-mode acceleration into the system and to allow in this way for a high signal-to-noise ratio in the balancing operation. We expect to gain an additional 60 dB in the instrument’s CMRR. The residual sensitivity to environmen- tal linear accelerations will be cancelled out by using a suitable vertical mechanical active iso- lation system, driven, with a proper feedback gain, by the common-mode acceleration signal itself. Both the elastic constant of the single oscil- lator and its inertia moment are affected by temperature. Such thermal effects are not a Fig. 3. Top view of the differential accelerometer mechanics. Two identical torsion pendula are carved out of a single-crystal sapphire plate by means of CNC ultrasonic impact grinding. The two disks are used as the tuning elements of two sapphire-loaded mi- crowave resonators. Two flaps with evaporated gold electrodes are used for electrostatic stiffness fine bal- ancing and cold damping of the accelerometers. 1099 A sapphire monolithic differential accelerometer as core sensor for gravity gradiometric geophysical instrumentation problem, as long as both sensors are at the same temperature. A difference of 1°C between the two sensors can however produce a sensitivity difference of the order of 10−5. In order to achieve 1 Eö sensitivity the temperature differ- ences must be controlled better than 10−5 K. It will be necessary an efficient passive isolation of the instrument container. The monolithic structure of the apparatus and the vacuum envi- ronment, already necessary in order to avoid viscous damping of the mechanical resonance, can warrant a very high thermal inertia. Semi- conductor temperature sensors will monitor the eventual thermal gradient. 3. The microwave transducer and readout system The motion of the proof-mass of each ac- cele-rometer is detected by means of a micro- wave Sapphire Dielectric Resonator Transducer (SDRT) (Blair et al., 1992). In the simplest form, a SDRT consists of two coaxial disks mounted inside a conducting shell which prevents radia- tion losses and unwanted coupling with the en- vironment. Each disk resonates with a Whisper- ing Gallery (WG) mode of the electromagnetic field. These modes are generated by total reflec- tion of the electromagnetic field on the external surface of the dielectric cylinder. By coupling two of these resonators we obtain a new res- onator whose resonant frequency ƒ is a function of the gap width between the two disks (fig. 4). With a suitable choice of the resonator geometry is possible to optimize the merit factor of the transducer Q⋅dƒ/dz, the product between the quality factor of the WG coupled mode and the tuning coefficient, which determines the dis- placement sensitivity of the transducer itself. There are two types of WG modes: E modes (Ez, Hρ, Hφ) and H modes (Hz, Eρ, Eφ), where z denotes the component along the disk rotation axis of the field vector, ρ its radial component, and φ its azimuthal component. Each mode is labeled by three numbers m, n and k, describing the number of field maximums in the azimuth, radial and longitudinal direction respectively. The principal modes (n = k = 1) allow Q factors that can be equal or even greater than the in- verse of the dielectric loss tangent of the dielec- tric material. The most sensitive resonances to changes in the gap spacing has been demon- strated to be the WG E modes, with high order azimuthal number (Blair et al., 1998). The rea- son is the E modes have a higher field density in the gap between the disks. Because of the continuity condition of the electromagnetic vector at the boundary of a dielectric, the Ez component of E mode increases from the value Ez (vacuum) in the gap spacing to εzEz(vacuum) in the dielectric, whereas the main field components of H mode remain unchanged. The tuning coefficient is also a function of the height of the disks. When the disk thickness decreases, the fraction of the electromagnetic power in the gap increases. As a consequence, the value of dƒ/dz increases. In our design the proof mass of each ac- celerometer is used as the tuning disk of the SDRT formed with the corresponding reference puck dug into the reference plate (see fig. 5). The dimensions of the disks was chosen in the way that the high order azimuthal whispering gallery E modes (m=7−12) resonate in the X- band (between 9 and 12 GHz). We have performed some preliminary tests on the SDRT using two sapphire disks with 42 mm Fig. 4. Sketch of a double-disk Whispering Gallery (WG) microwave resonator. 1100 Alessandro Bertolini, Nicolò Beverini, Andrea De Michele, Francesco Fidecaro and Francesco Mango in diameter and 5 mm in height. The sapphire crystals were of optical grade quality and without surface polish. The distance between the two disks was changed moving vertically one of them mounted on a micrometer screw. A microwave sweeper was used as tuneable microwave source and a diode as square law detector. Table I reports the experimental and theoret- ical value of the resonance frequencies of some modes when the two disks are in contact (zero distance), and the Q factor of the resonances. The theoretical frequencies are calculated using the method described in Blair et al. (1998). Fig- ure 6 plots experimental and theoretical values of resonance frequency and of the tuning coef- ficient dƒ/dz, as a function of the distance be- tween the two disks for the E11,1,1 mode. These results demonstrate that a gap spacing of 100 µm ensures for the E11,1,1 mode in our geometry a high tuning coefficient, around 6 MHz/µm, and a quality factor of the order of Q ≈ 50 000- 60 000 even at room temperature. Higher Q val- ues can be expected with sapphire crystals of better quality. Radiation loss shielding is achieved by gold plating the inner surface of the sapphire plates. The readout system will be realized following the fully differential interferometric scheme in- troduced by Blair et al. (1994), which has been demonstrated to be capable of displacement resolutions better than 10−14 m. The idea is to use one of the two resonators to build a loop os- cillator and excite the second one. The differen- tial acceleration is detected by measuring the frequency mismatch between the two SDRT with a balanced detector. A practical implemen- tation is sketched in fig. 7. A stable microwave signal is generated by a loop oscillator includ- ing one of the gradiometer resonators as the ref- erence cavity, a phase shifter and a low noise high-gain microwave amplifier. A power di- vider splits the microwave power between the two resonators used in trasmission mode. Cou- pling is achieved by means of two circulators. The power transmitted by the SDRTs is com- pared by means of a 180° power combiner, whose output is demodulated in-phase with re- spect to the loop oscillator in a low-noise Dou- ble-Balanced Mixer (DBM). The low-pass fil- Table I. Experimental and theoretical resonance frequencies of a Sapphire Dielectric Resonator made of two disks of 42 mm diameter and 5 mm height in contact. Mode Exp (GHz) Theo (GHz) Q E7,1,1 8.235 8.262 14000 E8,1,1 8.889 8.918 18000 E9,1,1 9.551 9.581 24000 E10,1,1 10.221 10.250 26000 E11,1,1 10.897 10.926 55000 E12,1,1 11.578 11.606 58000 Fig. 5. Section view of the differential accelerometer. The accelerometer and the reference plates are surround- ed by a metallic shielding enclosure. The two pairs of control electrodes allow for fine trimming of the CMRR and cold damping of the accelerometer motion. A sapphire monolithic differential accelerometer as core sensor for gravity gradiometric geophysical instrumentation tered output of the DBM is a voltage propor- tional to the sensed differential acceleration, in the frequency band DC-20 Hz. The balanced in- terferometric scheme allows for carrier sup- pression, enhancing the resolution, especially at low frequencies. Otherwise the sensitivity would be limited by the loop oscillator phase noise and by the frequency discriminator noise. By considering a reasonable dynamic range of 130 dB for the mixer, the frequency mismatch between the two microwave resonators must be limited within 100 kHz to allow the gradiome- ter to operate. A trimming procedure will be then necessary, after the machining, to suppress the bias effect due to shape mismatch in the di- electric resonators, producing a static frequency offset. The high mechanical Q level calls for a mechanism that limits the dynamics of the os- cillators. This will be done by an active damp- ing feedback control as used in sensitive detec- tors like the VIRGO gravitational antenna (Ac- ernese et al., 2004). Each sensor will be fitted with an actuator consisting in a small magnet fixed on the upper face of the moving arm, driv- en by a coil fixed to the frame. A laser beam is reflected by the moving arm and sent to a posi- tion sensitive detector. A current signal propor- 1101 Fig. 7. Schematic of the interferometric microwave readout system. One resonator is used as the refer- ence cavity for a X-band loop-oscillator. Its output is used to interrogate the other dielectric resonator and the frequency difference is read by a double-bal- anced mixer. Fig. 6. Resonance frequency and tuning coefficient of the E11,1,1 mode in function of the distance between the two disk of 42 mm diameter and 5 mm height (points-experimental values; line-theoretical values). 1102 Alessandro Bertolini, Nicolò Beverini, Andrea De Michele, Francesco Fidecaro and Francesco Mango tional to the moving arm velocity is injected. It is possible in this way to reduce the dynamic of the oscillation by orders of magnitude. The minimum detectable gravity gradient per bandwidth unit can be evaluated by comparing the expected signal with the noise. By consider- ing thermal noise due to the Brownian motion of the accelerometer test masses and the phase noise of the loop oscillator, we find the equation where ω0 is the mechanical resonant frequency, l is the accelerometer baseline, m is the test mass, Q is the mechanical quality factor, d f /dx is the tuning coefficient, Sφ is the loop oscillator phase noise spectral density and ε is the balancing fac- tor of the interferometric system. At room tempe- rature typical expected parameters (Tobar et al., 1993) are Q ≈104, Sφ≈10 −3/f 2 rad2/Hz and ε ≈ ≈ 0.01, yielding Γmin ≈ 1 Eö for 1-s of integra- tion time. 4. Conclusions We have presented the design features of the core differential accelerometer for new low cost lightweight gravity gradiometer instruments suitable for land geophysical surveys. The com- bination of the monolithic CNC machined me- chanics and of a fine electrostatic stiffness bal- ancing would produce a very high mechanical CMRR, of the order of 120 dB, allowing 1-Eö level resolutions to be achieved. At present the main components of the first prototype are being developed and fixed. The core mechanics is being built in the ultrasonic machining facility of University of Pisa, using a low cost, lower purity grade HEMLITE™ sap- phire. The microwave readout system design is being done by investigating on the WG modes of the single transducer. After the integration of the core sensor with the necessary sub-systems like the temperature control, the gimbal level- ing frame, the balancing jig etc., a first labora- ( ) ( ) d d l x f S mQ Q k T4 B 0 2 2 2 0 2 2 0 2 0 2 $ $ ω ω π ω ε ω ω ω ωω ω Γ = + - + φ f ^ _p h i7 A tory prototype of Γyz gradiometer will be tested. Moreover new Γzz sensors, less sensitive to the alignment with respect to the vertical direction, will be designed. Acknowledgements This work has been partially supported by «Fondazione della Cassa di Risparmio di Pisa». REFERENCES ACERNESE, F. et al. (2004): A local control system for the test masses of the VIRGO gravitational wave detector, Astroparticle Phys., 20, 617-628. BERTOLINI, A., N. BEVERINI, G. DAL LAGO, R. DESALVO, F. FIDECARO, F. FRANCESCONI, M. FRANCESCONI and G. OMICCIOLI (2003): A 5-axis CNC ultrasonic cutting machine: design and preliminary test, LIGO Technical Note T020198-00. BLAIR, D.G., E.N. IVANOV and H. PENG (1992): Sapphire di- electric resonator transducers, J. Phys. D, 25, 1110. BLAIR, D.G., E.N. IVANOV and H. PENG (1994): An ultra- high sensitivity sapphire transducer for vibration meas- urements, J. Phys. D, 27, 1150. BLAIR, D.G., E.N. IVANOV, M.E. TOBAR and B.D. CUTHBERT- SON (1998): Sensitivity and optimization of a high-Q sap- phire dielectric motion-sensing transducer, IEEE Trans. Ultrason. Ferroelect. Freq. Contr., 45, 1303-1313. CHAIN, H.A. and H.J. PAIK (1988a): Superconductive grav- ity gradiometer for sensitive gravity measurements, I. Theory, Phys. Rev. D, 35, 3551. CHAIN, H.A. and H.J. PAIK (1988b): Superconductive grav- ity gradiometer for sensitive gravity measurements, II. Experiment, Phys. Rev. D, 35, 3572. GOLDSTEIN, M.S. and J.J. BRETT (1998): Precision gravity gradiometer/AUV system, in Proceedings of the 1998 Workshop on Autonomous Underwater Vehicles, AUV’98, 167-174. LEE, J.B. (2001): FALCON gravity gradiometer technology, Exploration Geophys., 32, 247-251 (available on line at http://falcon.bhpbilliton.com/docs/falcon_gravity_gra- diometer_technology.pdf). MATTHEWS, R. (2002): Mobile Gravity Gradiometry, PhD Thesis (University of Western Australia), (unpublished). REBHAN, H., M. AGUIRRE and J. JOHANNESSEN (2000): The gravity field and steady-state ocean explorer mission - GOCE, ESA Earth Observation Q., 66, 6-11. RICHTER, B. and R.J. WARBURTON (1998): A new generation of superconducting gravimeters, in Proceedings of the 13th International Symposium on Earth Tides, Série Géophysique, Royal Observatory of Belgium, Brussel, 545-556. SNADDEN, M.J., J.M. MCGUIRK, P. BOUYER, K.G. HARITOS and M.A. KASEVICH (1998): Measurement of the Earth’s gravity gradient with an atom interferometer-based grav- ity gradiometer, Phys. Rev. Lett., 81, 971-974. TOBAR, M.E., A.J. GILES, S. EDWARDS and J. SEARLS (1993): 1103 A sapphire monolithic differential accelerometer as core sensor for gravity gradiometric geophysical instrumentation High Q TE stabilised sapphire microwave resonators for low noise applications, in Proc. IEEE Frequency Control Symposium, U.S.A. VAN KANN, F.J., M.J. BUCKINGHAM, C. EDWARDS and R. MATTHEWS (1994): Performances of a superconducting gravity gradiometer, Physica B, 194 and 196, 61-62. (received January 28, 2006; accepted July 20, 2006)