Shear wave splitting, also called seismic birefringence, is the phenomenon that occurs when a polarized shear wave enters an anisotropic medium (Fig. 1). The incident shear wave splits into two polarized shear waves (Fig. 2). Shear wave splitting is typically used as a tool for testing the anisotropy of an area of interest. These measurements reflect the degree of anisotropy and lead to a better understanding of the area’s crack density and orientation or crystal alignment.^{[1]} We can think of the anisotropy of a particular area as a black box and the shear wave splitting measurements as a way of looking at what is in the box.
Figure 1. (a) isotropic media, (b) anisotropic media with preferentially oriented cracks.
Figure 2. Animation of shear wave splitting upon entering an anisotropic medium. Courtesy of Ed Garnero.
Figure 3. Polarization diagram of shear wave arrivals. Abrupt changes in particle motion can be explained by the arrivals of two polarized shear waves.
Introduction
An incident shear wave may enter an anisotropic medium from an isotropic media by encountering a change in the preferred orientation or character of the medium. When a polarized shear wave enters a new, anisotropic medium, it splits into two shear waves (Fig.2). One of these shear waves will be faster than the other and oriented parallel to the cracks or crystals in the medium. The second wave will be slower than the first and sometimes orthogonal to both the first shear wave and the cracks or crystals in the media. The time delays observed between the slow and fast shear waves give information about the density of cracks in the medium. The orientation of the fast shear wave records the direction of the cracks in the medium.
When plotted using polarization diagrams, the arrival of split shear waves can be identified by the abrupt changes in direction of the particle motion (Fig.3).
In a homogeneous material that is weakly anisotropic, the incident shear wave will split into two quasishear waves with approximately orthogonal polarizations that reach the receiver at approximately the same time. In the deeper crust and upper mantle, the high frequency shear waves split completely into two separate shear waves with different polarizations and a time delay between them that may be up to a few seconds.^{[2]}
History
Hess^{[3]} (1964) made the first measurements of Pwave azimuthal velocity variations in oceanic basins. This area was chosen for this study because oceanic basins are made of large, relatively uniform homogeneous rocks. Hess observed, from previous seismic velocity experiments with olivine crystals, that if the crystals had even a slight statistical orientation this would be extremely evident in the seismic velocities recorded using seismic refraction. This concept was tested using seismic refraction profiles from the Mendocino Fracture Zone. Hess found that the slow shear waves propagated perpendicular to the plane of slip and the higher velocity component was parallel to it. He inferred that the structure of oceanic basins could be recorded quickly and understood better if these techniques were used.
Ando^{[4]} (1980) focused on identifying shearwave anisotropy in the upper mantle. This study focused on shear wave splitting recorded near the Chubu Volcanic Area in Japan. Using newly implemented telemetric seismographic stations, they were able to record both Pwave and Swave arrivals from earthquakes up to 260 km beneath the volcanic area. The depths of these earthquakes make this area ideal for studying the structure of the upper mantle. They noted the arrivals of two distinct shear waves with different polarizations (NS, fast and EW, slow) approximately 0.7 seconds apart. It was concluded that the splitting was not caused by the earthquake source but by the travel path of the waves on the way to the seismometers. Data from other nearby stations were used to constrain the source of the seismic anisotropy. He found the anisotropy to be consistent with the area directly below the volcanic area and was hypothesized to occur due to oriented crystals in a deep rooted magma chamber. If the magma chamber contained elliptical inclusions oriented approximately NS, then the maximum velocity direction would also be NS, accounting for the presence of seismic birefringence.
Crampin^{[5]} (1980) proposed the theory of earthquake prediction using shear wave splitting measurements. This theory is based on the fact that microcracks between the grains or crystals in rocks will open wider than normal at high stress levels. After the stress subsides, the microcracks will return to their original positions. This phenomenon of cracks opening and closing in response to changing stress conditions is called dilatancy. Because shear wave splitting signatures are dependent on both the orientation of the microcracks (perpendicular to the dominant stress direction) and the abundance of cracks, the signature will change over time to reflect the stress changes in the area. Once the signatures for an area are recognized, they may then be applied to predict nearby earthquakes with the same signatures.
Crampin^{[6]} (1981) first acknowledged the phenomenon of azimuthallyaligned shear wave splitting in the crust. He reviewed the current theory, updated equations to better understand shearwave splitting, and presented a few new concepts. Crampin established that the solution to most anisotropic problems can be developed. If a corresponding solution for an isotropic case can be formulated, then the anisotropic case can be arrived at with more calculations. The correct identification of body and surface wave polarizations is the key to determining the degree of anisotropy. The modeling of many twophase materials can be simplified by the use of anisotropic elasticconstants. These constants can be found by looking at recorded data. This has been observed in several areas world wide.^{[7]}
Physical Mechanism
Figure 4. Schematic diagram of two orthogonal polarized shear waves traveling through an anisotropic medium.
The difference in the travel velocities of the two shear waves can be explained by comparing their polarizations with the dominant direction of anisotropy in the area. The interactions between the tiny particles that make up solids and liquids can be used as an analogue for the way a wave travels through a medium. Solids have very tightly bound particles that transmit energy very quickly and efficiently. In a liquid, the particles are much less tightly bound and it generally takes a longer time for the energy to be transmitted. This is because the particles have further to travel to transfer the energy from one to another. If a shear wave is polarized parallel to the cracks in this anisotropic medium, then it may look like the dark blue wave in Figure 4. This wave is acting on the particles like energy being transferred through a solid. It will have a high velocity because of the proximity of the grains to each other. If there is a shear wave that is polarized perpendicular to the liquidfilled cracks or elongated olivine crystals present in the medium, then it would act upon these particles like those that make up a liquid or gas. The energy would be transferred more slowly through the medium and the velocity would be slower than the first shear wave. The time delay between the shear wave arrivals depends on several factors including the degree of anisotropy and the distance the waves travel to the recording station. Media with wider, larger cracks will have a longer time delay than a media with small or even closed cracks. Shear wave splitting will continue to occur until the shearwave velocity anisotropy reaches about 5.5%.^{[7]}
Mathematical Explanation
Mathematical Explanation(Ray theory)^{[8]}
The equation of motion in rectangular Cartesian coordinates can be written as

\frac{\partial}{\partial x_i} \left[c_{ijkl}\frac{\partial U_k}{\partial x_l}\right]=\rho\frac{\partial^2 U_j}{\partial t^2}


(1)

where t is the time, \rho is the density, U_j is the component of the displacement vector U, and c_{ijkl} represents the elastic tensor.
A wave front can be described by the equation

t=\tau\left(x_i\right)


(2)

The solution to (1) can be expressed as a ray series

U_k\left(x_i,t\right)=\sum_{n=0}^\infty U_k^\left(n\right)\left(x_i\right)f_n\left(t\tau\left(x_i\right)\right)


(3)

where the function f_n\left(\vartheta\right) satisfies the relation

df_{n+1}\left(\vartheta\right)/d\vartheta=f_n\left(\vartheta\right)


(4)

Substitute (3) into (1),

N\left(U^\left(n\right)\right)M\left(U^\left({n1}\right)\right)+L\left(U^\left({n2}\right)\right)=0


(5)

where the vector operators N,M,L are given by the formula:

\begin{cases} N_j\left(U^\left(n\right)\right)=\Gamma_{jk}U_k^\left(n\right)U_j^\left(n\right) \\ M_j\left(U^\left(n\right)\right)=p_i~a_{ijkl}\frac{\partial U_k^\left(n\right)}{\partial x_l}+ \rho^{1}\frac{\partial }{\partial x_i}\left(\rho~a_{ijkl}~p_l U_k^\left(n\right)\right) \\ L_j\left(U^\left(n\right)\right)=\rho^{1}\frac{\partial }{\partial x_i}\left(\rho~a_{ijkl} \frac{\partial U_k^\left(n\right)} {\partial x_l} \right) \end{cases}


(6)

where

\Gamma_{jk}=p_i~p_l~a_{ijkl} ,\quad a_{ijkl}=c_{ijkl}/\rho ,\quad p_i=\frac{\partial \tau}{\partial x_i}


(7)

For the first order n=0, so U^\left(1\right)=U^\left(2\right)=0, and only the first component of the equation (5) is left.
Thus,

N_j\left(U^\left(0\right)\right) =\Gamma_{jk}U_k^\left(0\right)U_j^\left(0\right) ==\left(\Gamma_{jk}\delta_{jk}\right)U_k^\left(0\right)==0


(8)

To obtain the solution of (8), the eigenvalues and eigenvectors of matrix \Gamma_{jk} are needed,

Det\left(\Gamma_{jk}G\delta_{jk}\right)=0


(9)

which can be rewritten as

where the values P,Q and R are the invariants of the symmetric matrix \Gamma_{jk}.
The matrix \Gamma_{jk} has three eigenvectors: g_1,~g_2,~g_3, which correspond to three eigenvalues of G_1,~G_2, and ~G_3.

For isotropic media, G_1=\alpha^2p_ip_i corresponds to the compressional wave and G_2=G_3=\beta^2p_ip_i corresponds to the two shear waves traveling together.

For anisotropic media,G_2 \ne G_3, indicates that the two shear waves have split.
Measurement of shear wave splitting parameters
Modeling^{[9]}
In an isotropic homogeneous medium, the shear wave function can be written as

u\left(\omega\right)=A~w\left(\omega\right)exp\left[i\omega T_0\right]\cdot\hat p


(10)

where A is the complex amplitude, w\left(\omega\right) is the wavelet function (the result of the Fourier transformed source time function), and \hat p is a real unit vector pointing in the displacement direction and contained in the plane orthogonal to the propagation direction.
The process of shear wave splitting can be represented as the application of the splitting operator to the shear wave function.

\Gamma=exp\left[i\omega \delta t/2\right]\hat f \hat f+ exp\left[i\omega \delta t/2\right]\hat s \hat s


(11)

where \hat f and \hat s are eigenvectors of the polarization matrix with eigenvalues corresponding to the two shear wave velocities.
The resulting split waveform is

u_s\left(\omega\right)=A~w\left(\omega\right)exp\left[i\omega T_0\right] \Gamma\left(\phi,\delta t\right)\cdot\hat p


(12)

Figure 5. Physical explanation of \delta t and \phi . Courtesy of Ed_Garnero.
Where \delta t is the time delay between the slow and fast shear waves and \phi is the angle between the polarizations of the slow and fast shear waves. These two parameters can be individually estimated from multiple component seismic recordings (Fig. 5).
Schematic Model
Figure 6 is a schematic animation showing the process of shear wave splitting and the seismic signature generated by the arrivals of two polarized shear waves at the surface recording station. There is one incident shear wave (blue) traveling vertically along the center grey axis through an isotropic medium (green). This single incident shear wave splits into two shear waves (orange and purple) upon entering the anisotropic media (red). The faster shear wave is oriented parallel to the cracks or crystals in the medium. The arrivals of the shear waves are shown on the right, as they appear at the recording station. The northsouth polarized shear wave arrives first (purple) and the eastwest polarized shear wave (orange) arrives about a second later.^{[5]}
Figure 6. Schematic animation of shear wave splitting seismic energy recorded by the
geophone recording station.
Applications/Justification/Usefulness
Shear wave splitting measurements have been used to successfully predict earthquakes, and to map fracture networks created by high pressure fracturing of reservoirs.
According to Crampin^{[5]} shear wave splitting measurements can be used to monitor stress levels in the earth. It is well known that rocks near an earthquakeprone zone will exhibit dilatancy. Shear wave splitting is produced by seismic waves traveling through a medium with oriented cracks or crystals. The changes in shear wave splitting measurements over the time leading up to an impending earthquake can be studied to give insight to the timing and location of the earthquake. These phenomena may be observed many hundreds of kilometers from the epicenter.
The petroleum industry uses shearwave splitting measurements to map the fractures throughout a hydrocarbon reservoir. To date, this is the best method to gain in situ information about the fracture network present in a hydrocarbon reservoir.^{[10]} The best production in a field is associated with an area where there are multiple small fractures that are open, allowing for constant flow of the hydrocarbons. Shearwave splitting measurements are recorded and analyzed to obtain the degree of anisotropy throughout the reservoir. The area with the largest degree of anisotropy will generally be the best place to drill because it will contain the largest number of open fractures.^{[11]}
Case Examples
A successfully stressforecast earthquake in Iceland^{[12]}
On October 27, 1998, during a four year study of shear wave splitting in Iceland, Crampin and his coworkers recognized that time delays between split shearwaves were increasing at two seismic recording stations, BJA and SAU, in southwest Iceland. The following factors lead the group to recognize this as a possible precursor to an earthquake:

The increase persisted for nearly 4 months.

It had approximately the same duration and slope as a previously recorded magnitude 5.1 earthquake in Iceland.

The time delay increase at station BJA started at about 4ms/km and escalated to approximately 10ms/km.

10ms/km was the inferred level of fracture for the previous earthquake.
These features suggested that the crust was approaching fracture criticality and that an earthquake was likely to occur in the near future. Based on this information, an alert was sent to the Iceland Meteorological Office (IMO) on October 27 and 29, warning of an approaching earthquake. On November 10, they sent another email specifying that an earthquake was likely to occur within the next 5 months. Three days later, on November 13, IMO reported a magnitude 5 earthquake near the BJA station. Crampin et al. suggests that this is the first scientifically, as opposed to precursory or statistically, predicted earthquake. They proved that variations of shearwave splitting can be used to forecast earthquakes.
This technique was not successful again until 2008 due to the lack of appropriate sourcegeophoneearthquake geometry needed to evaluate changes in shear wave splitting signatures and time delays.^{[7]}
Volti and Crampin observed temporal increases in Band1 timedelays for 5 months at approximately 240 kilometer depth in directions N,SW and W,SW before the 1996 Gjalp Eruption in Vatnajökull Icefield. This was the largest eruption in Iceland in several decades.
The pattern of increasing shear wave splitting timedelays is typical of the increase now seen before many earthquakes in Iceland and elsewhere. The time delays just before earthquakes characteristically decrease immediately following the eruption because the majority of the stress is released at that one time. The increase in normalized timedelays in volcanic eruptions does not decrease at the time of the eruption but gradually declines at about 2ms/km/year over several. This decrease is approximately linear and there appeared to be no other significant magmatic disturbances during the period following the eruption.
More observations are needed to confirm whether the increase and decrease time delay pattern is universal for all volcanic eruptions or if each area is different. It is possible that different types of eruptions show different shear wave splitting behaviors.^{[7]}
Fluidinjection in Petroleum Engineering^{[14]}
Bokelmann and Harjes reported the effects on the shear waves of fluid injection at about 9 kilometer depth in the German Continental Deep Drilling Program (KTB) deep drilling site in southeast Germany. They observed shearwave splitting from injectioninduced events at a pilot well offset 190 meters form the KTB well. A borehole recorder at a depth of 4,000 meters was used to record the splitting measurements.
They found:

Temporal variations in shearwave splitting as a direct result of injectioninduced events.

That the initial ~1% shear wave splitting decreases by 2.5% in the next 12 hours following the injection.

The largest decrease occurred within two hours after the injection.

The splitting time to be very stable after the injection ceased.
No direct interpretation of the decrease is proposed but it is suggested that the decrease is associated with stress release by the induced events.
Limitations
Shearwave splitting measurements can provide the most accurate and in depth information about a particular region. However, there are limits that need to be accounted for when recording or analyzing shear wave splitting measurements. These include the sensitive nature of shear waves, that shear wave splitting varies with incidence and azimuth, and that shear waves may split multiple times throughout an anisotropic medium, possibly every time the orientation changes.^{[15]}
Shear wave splitting is very sensitive to fine changes in the pore pressure in the Earth's crust. In order to successfully detect the degree of anisotropy in a region there must be more several arrivals that are well distributed in time. Too few events cannot detect the change even if they are from similar waveforms.^{[7]} The Shear wave splitting varies with both incidence angle and propagation azimuth. Unless this data is viewed in polar projection, the 3D nature is not reflected and may be misleading.^{[7]} Shear wave splitting may be caused by more than just one layer that is anisotropic and located anywhere between the source and the receiver station. The shear wave splitting measurements have extensive lateral resolution but very poor vertical resolution.^{[16]} The polarizations of shear waves vary throughout the rock mass. Therefore the observed polarizations may be those of the near surface structure and are not necessarily representative of the structure of interest.^{[17]}
Common Misunderstandings
Due to the nature of split shear waves, when they are recorded in typical threecomponent seismograms, they write very complicated signatures. Polarizations and time delays are heavily scattered and vary greatly both in time and space. Because of the variation in signature, it is easy to misinterpret the arrivals and polarization of incoming shear waves.^{[18]} Below is an explanation of a few of the common misunderstandings associated with shear waves, further information can be found in Crampin and Peacock (2008).^{[7]}
Shear waves that propagate along the ray path at a group velocity have polarizations that are only orthogonal in a few specific directions. Polarizations of body waves are orthogonal in all phase velocity directions, however this type of propagation is generally very difficult to observe or record.

Polarizations of split shearwaves are fixed, parallel to cracks, or normal to spreading centers.^{[7]}^{[18]}
Even when propagating through parallel cracks or perpendicular to spreading centers or parallel to cracks, the polarizations of shear waves will always vary in three dimensions with incidence and azimuth within the shear wave window.
This statement only holds true if the fluid in the cracks is somehow removed. This may be accomplished via chemical absorption, drainage, or flow to the surface. However, these occur in relatively rare instances and there is evidence that supports the presence of fluids at depth. This includes data from the Kola deep well and the presence of high conductivity in the lower crust.

Signaltonoise ratios of shearwave splitting above small earthquakes can be improved by stacking.^{[7]}^{[18]}
Stacking seismic data from a reflection survey is useful because it was collected with a predictable, controlled source. When the source is uncontrolled and unpredictable, stacking the data only degrades the signal. Because recorded shear wave time delays and polarizations vary in their incidence angle and azimuth of radio propagation, stacking these arrivals will degrade the signal and decrease the signal to noise ratio, resulting in a plot that is noisy and hard to interpret at best.^{[7]}
Future Trends
Our understanding of shear wave splitting and how to best use the measurements is constantly improving. As our knowledge improves in this area, there will invariably be better ways of recording and interpreting these measurements and more opportunities to use the data. Currently, it is being developed for use in the petroleum industry and for predicting earthquakes and volcanic eruptions.
Shear wave splitting measurements have been used successfully to predict several earthquakes. With better equipment and more densely spaced recording stations, we have been able to study the signature variations of shear wave splitting over earthquakes in different regions. These signatures change over time to reflect the amount of stress present in an area. After several earthquakes have been recorded and studied, the signatures of shear wave splitting just before an earthquake occurs become well known and this can be used to predict future events. This same phenomenon can be seen before a volcanic eruption and it is inferred that they may be predicted in the same manner.
The petroleum industry has been using shear wave splitting measurements recorded above hydrocarbon reservoirs to gain invaluable information about the reservoir for years. Equipment is constantly being updated to reveal new images and more information.^{[7]}
See also
References

^ Aki, K; Richards, P.G. (2002). "Quantitative Seismology" (Second ed.). University Science Books, Sausalito, CA.

^ Vecsey, L., J; Babuska, V. (2008). "Shearwave splitting measurementsProblems and solutions". Tectonophysics 462: 178–196.

^ H. H. Hess (1964). "Seismic Anisotropy of the Uppermost Mantle under Oceans".

^ M. Ando; Y.Ishikawa, H. Wada (1980). "Swave anisotropy in the upper mantle under a volcanic area in Japan".

^ ^{}a ^{b} ^{c} S. Crampin; R.Evans , B. Ucer, M.Doyle, J. Davis, G. Yegorkina and A. Miller (1980). "Observations of dilatancyinduced polarization anomalies and earthquake prediction".

^ S. Crampin (1981). "A review of wave motion in anisotropic and cracked elasticmedia".

^ ^{}a ^{b} ^{c} ^{d} ^{e} ^{f} ^{g} ^{h} ^{i} ^{j} ^{k} ^{l} ^{m} S. Crampin; S. Peacock (2008). "A review of the current understanding of seismic shearwave splitting in the Earth's crust and common fallacies in interpretation".

^ V. Cerveny (1972). "Seismic rays and ray intensities in inhomogeneous anisotropic media". Geophys. JR astr. Soc.

^ P. Silver; W. Chan (1991). "Shear Wave Splitting and Subcontinental Mantle Deformation". Journal of Geophysical Research 96: 16,429–16,454.

^ R. Bale; J. Li, B. Mattocks and S. Ronen (2006). "LeastSquares Measurement of ShearWave Splitting". CSPG / CSEG / CWLS Joint Conference. Retrieved 05/12/2010.

^ E. LaBarre; T. Davis; R. Benson (March 19, 2008). "Finding the sweet spot". E&P. Retrieved 06/05/2012.

^ S. Crampin; T.Volti,R. Stefánsson (1999). "A successfully stressforecast earthquake".

^ T. Volti; S. Crampin (2003). "A fouryear study of shearwave splitting in Iceland: 2. Temporal changes before earthquakes and volcanic eruptions". Geological Society, London, Special Publications 212 (1): 135–149.

^ G. Bokelmann; H. Harjes (2000). "Evidence for temporal variation of seismic velocity within the upper continental crust". Journal of geophysical research 105: 23879–23894.

^ R. Hoar; K. Stokoe (1978). "Generation and Measurement of Shear Waves In Situ". Dynamic Geochemical Testing: 3–29.

^ M. K. Savage (February 1999). "Seismic anisotropy and mantle deformation: What have we learned from shear wave splitting?". Reviews of Geophysics 37: 65–106.

^ S. Crampin; Lovell, John H. (1991). "A decade of shearwave splitting in the Earth's crust: what does it mean? what use can we make of it? and what should we do next?". Geophysics Journal International 107 (3): 387–407.

^ ^{}a ^{b} ^{c} ^{d} ^{e} S. Crampin; Y. Gao (2006). "A review of techniques for measuring shearwave splitting above small earthquakes". Physics of the Earth and Planetary Interiors 159: 1–14.
Further reading

Crampin, S., and J.H. Lovell, 1991, A decade of shearwave splitting in the Earth’s crust: what does it mean? what use can we make of it? and what should we do next? : Geophysics Journal International, v. 107, p. 387407.

Crampin, S., and S. Peacock, 2005, A review of shearwave splitting in the compliant crackcritical anisotropic Earth: Wave Motion, v. 41, p. 5977.

Long, M. D., Hoop, M. V., and R.D. van der Hilst, 2008, Waveequation shear wave splitting tomography: Geophysics Journal International, v. 172, p. 311330.

M. Pastori, D. Piccinini, L. Margheriti, L. Improta, L. Valoroso, L. Chiaraluce, C. Chiarabba Stress aligned cracks in the upper crust of the Val d'Agri region as revealed by shear wave splitting, Geophysical Journal International, Volume 179, Issue 1, pages 601–614, October 2009

Piccinini, D, Pastori M., Margheriti L., ANISOMAT+: An automatic tool to retrieve seismic anisotropy from local earthquakes, Computers & Geosciences, Volume 56, July 2013, Pages 6268, ISSN 00983004, http://dx.doi.org/10.1016/j.cageo.2013.01.012.

Savage, M. K., February 1999, Seismic anisotropy and mantle deformation: What have we learned from shear wave splitting?: Reviews of Geophysics, v. 37, p. 65106.

Ucisik, N., l. Gudmundsson, W. Hanka, T. DahlJensen, K. Mosegaard, and K. Priestley, 2008, Variations of shearwave splitting in Greenland: Mantle anisotropy and possible impact of the Iceland plume: Tectonophysics, v. 462, p. 137148.
External links

Alfred Wegener Institute for polar and Marine Research(AWI)(Germany)

Shearwave splitting in Matlab(France)

A lot of interesting seismic images(ASU)

Information on Solids, Liquids, and Gasses
MATLAB Code for demonstration
You can download a MATLAB code and create a demonstration movie by yourself here on MathWorks website.
Figure 7 is a screen shot of the Matlab Demo output.
Figure 7. Screen shot of Matlab demo of shear wave splitting
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.