3D Map of Earthquakes in 2025

Plot the location and depth of worldwide earthquakes in 2025

The majority of earthquakes around the world can be found along continental plate boundaries, either due to divergent plate boundaries where new oceanic crust is generated, or at convergent plate boundaries where collisions and subduction can occur.

Query for all earthquakes in 2025 with magnitude between 2 and 10:

In[1]:=

data=EarthquakeData[All,{2,10},{{2025,1,1},{2025,12,31}}];

Count the number of earthquakes retrieved:

In[2]:=

Length[data]

Out[2]=

26953

Extract the magnitude, position and depth from the full data:

In[3]:=

props=Lookup[#,{"Magnitude","Position","Depth"}]&/@data;

Convert the position to Cartesian coordinates and the depth, both in meters:

In[4]:=

triples0=ReverseSort[{#〚1〛,GeoPositionXYZ[#〚2〛]〚1〛,QuantityMagnitude[#〚3〛,"Meters"]}&/@Normal[props]〚All,2〛];

Bin the triples into groups of 0.1 magnitude and subtract the depth, exaggerated by 2.5 times, from the surface location:

In[5]:=

triples={Round[#〚1〛,0.1],#〚2〛(Norm[#〚2〛]-2.5#〚3〛)/Norm[#〚2〛]}&/@triples0;

Calculate the range of earthquake magnitudes found in the retrieved data:

In[6]:=

range=MinMax[triples〚All,1〛]

Out[6]=

{2.,8.8}

Query for the bounding polygons for all countries, including their outlying areas:

In[7]:=

polys=EntityValue["Country",GeoVariant["Polygon","AllAreas"]];

Convert the polygons to 3D lines and apply normals that radiate from the origin:

In[8]:=

lines=Line[#〚1〛,VertexNormals#〚1〛]&/@(Flatten[GeoBoundary[polys]/.JoinedCurveIdentity]/.gp_GeoPositionGeoPositionXYZ[gp]["Data"]);

Group the earthquakes by magnitude, apply size and color and apply normals that radiate from the origin:

In[9]:=

points=KeyValueMap[{Blend[{StandardRed,StandardYellow,White},Rescale[#,range]],PointSize[#^3/20000],Point[#2,VertexNormals#2]}&,GroupBy[triples,FirstLast]];

Grouping the points optimizes rendering time by reducing the number of

Point

primitives needed from 26,953 to only 61:

In[10]:=

Length[points]

Out[10]=

Assemble a 3D graphics scene:

In[11]:=

gr=With[{t=3Pi/4},Graphics3D[{Inset[Graphics[{Gray,Circle[]}],{0,0,0},Automatic,ImageScaled[.952]],points,GrayLevel[.5],lines},BackgroundBlack,BoxedFalse,SphericalRegionTrue,ViewAnglePi/12,Lighting{AmbientLight[GrayLevel[.1]],DirectionalLight[White,ImageScaled[{0,0,1}]]},ViewPoint4{Cos[t],Sin[t],0},ImageSize450]]

Out[11]=

A video can be generated if you are patient for each frame to be rendered:

In[12]:=

AnimationVideo[Show[gr,ViewPoint4{Cos[t],Sin[t],0},ImageSize450],{t,0,2Pi},RasterSize450,DefaultDuration15]

Out[12]=

Publisher Information

Contributed by: Jeff Bryant

Wolfram Language Example Repository

3D Map of Earthquakes in 2025

See Also

Publisher Information

3D Map of Earthquakes in 2025

See Also

Related Symbols

Publisher Information