diff --git a/NEWS.md b/NEWS.md
index d09af6f..c7537eb 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,5 +1,11 @@
 # distributional (development version)
 
+## Breaking changes
+
+* The `quantile()` method for `dist_multivariate_normal()` now defaults to 
+  equicoordinate quantiles. To obtain marginal quantiles, you should use 
+  `quantile(dist, p, type = "marginal")`.
+
 # distributional 0.5.0
 
 ## New features
diff --git a/R/dist_multivariate_normal.R b/R/dist_multivariate_normal.R
index ed8d273..49a1540 100644
--- a/R/dist_multivariate_normal.R
+++ b/R/dist_multivariate_normal.R
@@ -56,7 +56,7 @@ log_density.dist_mvnorm <- function(x, at, ..., na.rm = FALSE){
 }
 
 #' @export
-quantile.dist_mvnorm <- function(x, p, type = c("marginal", "equicoordinate"),
+quantile.dist_mvnorm <- function(x, p, type = c("equicoordinate", "marginal"),
                                  ..., na.rm = FALSE){
   type <- match.arg(type)
   q <- if (type == "marginal") {