KMZ fault download from USGS,

US Faults data

https://earthquake.usgs.gov/hazards/qfaults/

cp *.kmz  *.zip
unzip *.zip #get the kml files.
gmt kml2gmt *.kml > data # convert to gmt xy file.

but there is something wrong with this conversion.

another way make it success

download GIS Shapefiles from https://earthquake.usgs.gov/hazards/qfaults/

In GDAL 2.2.1 the call is:

ogr2ogr -f "OGR_GMT" filename.gmt filename.shp
gmt pscoast -R-170/170/-80/80 -JM6i -B1WS -P -Ggray -K >test.ps
gmt psxy -R -J -O qfaults.gmt -W0.5p,blue  >>test.ps #plot fault geometry
#gmt psxy data -R90/130/30/60 -JM15c -B5 -Sqn1:+Lh > map.ps， to write the name of the fault data has form of
#> -L"AAA Fault"
#95.00 35.00
#-Sqn1:+Lh+f15p,blue+n0c/0.3c change the color and location of name.
#-Sqn1:+Lh+f15p,blue+n0c/0.3c+Tfaultname.txt , return faultname information on the graph.
   X   Y   angle   label# in faultname file,change angle and X Y,then use pstext adding it to the picture.