diff --git a/src/empirical/transiogram.jl b/src/empirical/transiogram.jl
index b344372..a089413 100644
--- a/src/empirical/transiogram.jl
+++ b/src/empirical/transiogram.jl
@@ -100,8 +100,7 @@ geospatial `data` along a given `direction` with band tolerance `dtol` in length
 Optionally, forward `parameters` for the underlying [`EmpiricalTransiogram`](@ref).
 """
 function DirectionalTransiogram(dir, data::AbstractGeoTable, var; dtol=1e-6u"m", kwargs...)
-  rng = MersenneTwister(123)
-  Π = partition(rng, data, DirectionPartition(dir; tol=dtol))
+  Π = partition(MersenneTwister(123), data, DirectionPartition(dir; tol=dtol))
   EmpiricalTransiogram(Π, var; kwargs...)
 end
 
@@ -115,8 +114,7 @@ tolerance `ntol` in length units.
 Optionally, forward `parameters` for the underlying [`EmpiricalTransiogram`](@ref).
 """
 function PlanarTransiogram(normal, data::AbstractGeoTable, var; ntol=1e-6u"m", kwargs...)
-  rng = MersenneTwister(123)
-  Π = partition(rng, data, PlanePartition(normal; tol=ntol))
+  Π = partition(MersenneTwister(123), data, PlanePartition(normal; tol=ntol))
   EmpiricalTransiogram(Π, var; kwargs...)
 end
 
diff --git a/src/empirical/variogram.jl b/src/empirical/variogram.jl
index f1bdf92..f492970 100644
--- a/src/empirical/variogram.jl
+++ b/src/empirical/variogram.jl
@@ -110,8 +110,7 @@ geospatial `data` along a given `direction` with band tolerance `dtol` in length
 Optionally, forward `parameters` for the underlying [`EmpiricalVariogram`](@ref).
 """
 function DirectionalVariogram(dir, data::AbstractGeoTable, var₁, var₂=var₁; dtol=1e-6u"m", kwargs...)
-  rng = MersenneTwister(123)
-  Π = partition(rng, data, DirectionPartition(dir; tol=dtol))
+  Π = partition(MersenneTwister(123), data, DirectionPartition(dir; tol=dtol))
   EmpiricalVariogram(Π, var₁, var₂; kwargs...)
 end
 
@@ -125,8 +124,7 @@ tolerance `ntol` in length units.
 Optionally, forward `parameters` for the underlying [`EmpiricalVariogram`](@ref).
 """
 function PlanarVariogram(normal, data::AbstractGeoTable, var₁, var₂=var₁; ntol=1e-6u"m", kwargs...)
-  rng = MersenneTwister(123)
-  Π = partition(rng, data, PlanePartition(normal; tol=ntol))
+  Π = partition(MersenneTwister(123), data, PlanePartition(normal; tol=ntol))
   EmpiricalVariogram(Π, var₁, var₂; kwargs...)
 end