diff --git a/source/taup/index.rst b/source/taup/index.rst
index 804d86a5..b42c623e 100644
--- a/source/taup/index.rst
+++ b/source/taup/index.rst
@@ -26,3 +26,4 @@ TauP
    taup_curve
    taup_path
    taup_setsac
+   taup_velplot
diff --git a/source/taup/introduction.rst b/source/taup/introduction.rst
index 97bc026a..af31c40b 100644
--- a/source/taup/introduction.rst
+++ b/source/taup/introduction.rst
@@ -47,7 +47,7 @@ TauP 既提供了图形界面，也提供了命令行工具。
   - ``taup create``\ ：预计算速度模型
   - ``taup phase``\ ：通过语言描述震相的射线路径
   - ``taup slowplot``\ ：绘制速度模型的慢度-深度剖面
-  - ``taup velplot``\ ：绘制速度模型的速度-深度剖面
+  - :doc:`taup_velplot`\ ：绘制速度模型的速度-深度剖面
   - ``taup velmerge``\ ：将部分速度模型与一个标准模型合并以构成一个新的速度模型
 
 参考文献
diff --git a/source/taup/taup_velplot.jpg b/source/taup/taup_velplot.jpg
new file mode 100644
index 00000000..4f34177c
Binary files /dev/null and b/source/taup/taup_velplot.jpg differ
diff --git a/source/taup/taup_velplot.rst b/source/taup/taup_velplot.rst
new file mode 100644
index 00000000..2e3419cc
--- /dev/null
+++ b/source/taup/taup_velplot.rst
@@ -0,0 +1,26 @@
+taup velplot
+============
+
+:本节贡献者: |田冬冬|\（作者）、
+             |姚家园|\（审稿）
+:最近更新日期: 2021-06-02
+
+----
+
+``taup velplot`` 可以很方便地生成用于绘制速度-深度剖面的 GMT 脚本。
+
+使用 ``taup velplot --help`` 可查询该命令的完整用法和选项。此处仅介绍其常用用法。
+
+生成 iasp91 参考模型的速度-深度剖面绘图脚本::
+
+    $ taup velplot -mod iasp91
+
+执行以上命令生成的脚本 :file:`iasp91_vel.gmt`\ ，生成 PDF 格式的速度-深度剖面::
+
+    $ sh iasp91_vel.gmt
+
+绘制出来的图件效果如下（蓝线为 P 波速度，红线为 S 波速度）：
+
+.. image:: taup_velplot.jpg
+   :width: 400 px
+   :align: center